Distributive Justice and Urban Form Adaptation to Flooding Risks: Spatial Analysis to Identify Toronto's Priority Neighborhoods

Empirical evidence points out that urban form adaptation to climate-induced flooding events—through interventions in land uses and town plans (i. e., street networks, building footprints, and urban blocks)—might exacerbate vulnerabilities and exposures, engendering risk inequalities and climate injustice. We develop a multicriteria model that draws on distributive justice's interconnections with the risk drivers of social vulnerabilities, flood hazard exposures, and the adaptive capacity of urban form (through land uses and town plans). The model assesses “who” is unequally at-risk to flooding events, hence, should be prioritized in adaptation responses; “where” are the high-risk priority areas located; and “how” can urban form adaptive interventions advance climate justice in the priority areas. We test the model in Toronto, Ontario, Canada, where there are indications of increased rainfall events and disparities in social vulnerabilities. Our methodology started with surveying Toronto-based flooding experts who assigned weights to the risk drivers based on their importance. Using ArcGIS, we then mapped and overlayed the risk drivers' values in all the neighborhoods across the city based on the experts' assigned weights. Accordingly, we identified four high-risk tower communities with old infrastructure and vulnerable populations as the priority neighborhoods for adaptation interventions within the urban form. These four neighborhoods are typical of inner-city tower blocks built in the 20th century across North America, Europe, and Asia based on modern architectural ideas. Considering the lifespan of these blocks, this study calls for future studies to investigate how these types of neighborhoods can be adapted to climate change to advance climate justice.

marginalized groups unable to cope with flood hazards Michael et al., 2019). For instance, there is evidence that low-income neighborhoods contain a higher percentage of impervious surfaces than affluent neighborhoods due to a lack of green spaces (Bautista et al., 2015;Garcia-Lamarca et al., 2021), leading to their inadequate adaptive capacity.
These risk inequities are rooted in the uneven patterns of urban development based on economic rationales that have long prioritized infrastructure investments in high-value real estate, leading to decades of disinvestments in hazard-exposed and impoverished yet vulnerable neighborhoods (Herreros-Cantis et al., 2020). The prevalence of climate change further extended the rationales underlying inequities, hence, exacerbated vulnerabilities and exposures through land use planning , and we argue the town plans' design, where the town plan is defined as the streets and their networks and the arrangements of the building footprints and urban blocks (Conzen, 1960). Henceforth, urban form refers to land uses and the town plan -two of the three Conzen's (1960) urban morphology components 1 . For example, when retreat is adopted as a land use adaptation measure for flood-prone areas, it often entails the forced relocation of marginalized communities to sites far away from their social networks and livelihoods, hence worsening their vulnerabilities (Henrique and Tschakert, 2019).
Despite such unequal outcomes, there is a deficit of empirical studies that propose methodologies to measure how the adaptive capacity (or adaptation) of urban form is connected to the differential vulnerabilities (i.e., different sensitivities to risks), exposures, and risk inequities (Mohtat and Khirfan, 2021). This deficit is attributed to the nascence of theoretical studies that connect urban form with adaptation and adaptive capacity in general (Dhar and Khirfan, 2017;Sharifi, 2019c), and with climate justice in particular (Mohtat and Khirfan, 2021).
To identify how adaptation interventions can be distributed to avoid flood risk inequities, hence advance climate justice, this study draws on Rawls's (1971) distributive justice, referring to the just spatial distribution of resources to maximize benefits to the disadvantaged. We operationalize Dhar and Khirfan's (2017) framework for measuring urban form's adaptive capacity to investigate the spatial distribution of adaptation interventions, hence urban form's adaptive capacity, and explore this adaptive capacity's connections to differential vulnerabilities and hazard exposures. Accordingly, we develop a multicriteria model that includes indicators and variables to identify the spatial distribution patterns of risk drivers: social vulnerabilities, flood hazard exposures, and areas with a low adaptive capacity of urban form. Our model assesses specifically "who" are unequally at-risk to flooding events, hence should be prioritized for adaptation interventions; "where" are the high-risk priority areas located; and "how" urban form adaptive interventions may advance climate justice in these priority areas.
We test this model, which can be applied in any city within Canada and beyond, in Toronto, Ontario, Canada, where there are indications of increased frequency and intensity of flood 1 The third component is the three-dimensional built form. events combined with the disparities in social vulnerabilities (Feltmate and Thistlethwaite, 2012;Rincón et al., 2018). We aim to identify how social vulnerabilities, flood exposures, and adaptation interventions within the urban form are distributed in Toronto? Based on this, which neighborhoods are experiencing the highest risks of floods and need to be prioritized in adaptation? And how can we identify these priority neighborhoods?
To answer these queries, we developed a survey that asked Toronto-based flooding experts to weigh the importance of risk drivers of our multicriteria model and their associated indicators in triggering flood risks in Toronto. We then overlaid the values of the risk drivers and their indicators in ArcGIS, using the experts' assigned weights. The results reveal that flood risks are disproportionately distributed in four tower neighborhoods with old infrastructure, where low-income, racialized, and migrant populations concentrate, namely: Thorncliffe Park, Flemingdon Park, North St. James Town, and Black Creek.

FROM DIFFERENTIAL VULNERABILITIES TO CLIMATE JUSTICE IN URBAN FORM ADAPTATION TO FLOODING RISKS
Vulnerability, or people's susceptibility to being adversely affected by shocks, stresses, and hazards (Adger, 2006b;Gallopín, 2006), is not equal but differential. Differential vulnerability entails that some social groups undergo greater human, livelihood, and financial losses due to their exposure to stresses and lack of coping capacity (Suarez, 2002;. Evidence on differential vulnerabilities abounds globally: from the proximity of racial neighborhoods to contaminated sites and the ensuing negative impacts on the health of their residents in the USA, to the lack of low-income communities' access to potable water and sanitary services, hence, their sensitivity to droughts in the Philippines (Bautista et al., 2015;Porio et al., 2019). Differential vulnerabilities are rooted in the historical capitalist processes of urban development and their embedded domination and oppression patterns that shape inequity in the spatial distribution of urban assets (e.g., housing, land, green space, and infrastructure) and entitlements among socially different groups, whether across income, race, gender, or ethnicity, among others (Sen, 1982;Ribot, 2014;Michael et al., 2019). With the emergence of climate change as an urban crisis, the historical disinvestments in disenfranchised neighborhoods and the systematic exclusions of the disadvantaged from power structures place vulnerable groups in unsafe living conditions, exacerbating their vulnerabilities and exposures to different hazards, including flooding events (Blaikie et al., 2005;Michael et al., 2019). Additionally, efforts to mitigate climatic hazards, such as through adaptation, align with the uneven historical mechanisms of urban development, prioritizing the protection of urban economies over climate justice through selective investment in vital urban infrastructure and wealth reproduction systems Rice, 2019, 2020).

Flood Risks and Climate Justice Challenges
Changing precipitation rates combined with the increase in the density of urban impervious surfaces, old and overburdened drainage systems, and urban population, particularly in lowlying areas, intensify the risk of loss of lives and livelihoods and damage to properties and infrastructure from rainfall run-off and river flooding events (Faccini et al., 2018;O'donnell and Thorne, 2020;Sohn et al., 2020). Yet, individuals experience these flood risks differentially, depending on three contextspecific risk drivers: social vulnerabilities, low adaptive capacity, and exposure to flooding hazards. In fact, empirical evidence indicates that social vulnerabilities are associated with inequities in flood hazard exposures and access to adaptive capacity, triggering inequities in the spatial distribution of risks across the lines of race, income, and ethnicity, among others (Suarez, 2002;Islam and Winkel, 2017;Herreros-Cantis et al., 2020).
The uneven processes of urban development have forced marginalized groups with economically precarious and socially unstable conditions to live in deteriorating settlements, prone to power outages and infrastructure failures in the face of hazards (Walker and Burningham, 2011;Graham et al., 2016). Many of these settlements are located in low real estate value and precarious sites, like low-lying areas, floodplains, and industrial zones with impervious surfaces, which increase their exposure to flooding events. The lack of land tenure rights and informality in the Global South and discriminatory policies and zoning laws based on market rules in the Global North have led, over time, to the systematic disinvestment in these vulnerable and flood-prone neighborhoods (Chakraborty et al., 2014;Borie et al., 2019;Michael et al., 2019). Among the residents of these neighborhoods are new immigrants with language and employment barriers who lack community connections and citizenship entitlements, including election rights, to influence the formal urban governance structures and local decisionmakers; hence, they are often excluded from flood awareness, warning, and management programs (Donner and Rodríguez, 2008;Dodman et al., 2019;Turhan and Armiero, 2019). Additionally, the employment of these vulnerable groups in lowpaying service jobs, their everyday struggles for basic needs like food, and their lack of housing ownership render them financially unable to adopt flood protective behavior, such as buying insurance and retrofitting their flimsy settlements Herreros-Cantis et al., 2020;Ziervogel, 2020).
With their lack of preparedness, disenfranchised and marginalized vulnerable groups are more at risk of losing life, assets, and income due to flood hazards than the affluent groups in society (Collins et al., 2018;Kim et al., 2018). In addition, they have fewer opportunities for recovery, reconstruction, and relief due to their lack of access to personal wealth and timely and adequate assistance programs such as loans and emergency services (Rufat et al., 2015;Graham et al., 2016;Thomas and Warner, 2019). Hence, their frequent experience of risks worsens their existing vulnerabilities, reproduces new ones, and reduces their capacity to cope with future hazards.

Climate Justice Challenges in Urban Form Adaptation
Climate change adaptation refers to "the process of adjustment to actual or expected climate and its effects . . . to moderate or avoid harm or exploit beneficial opportunities" while adaptive capacity is the ability of humans, institutions, and systems to adapt to climatic effects (IPCC, 2014, p. 5). Urban form adaptation entails physical interventions in the built environment and functions to minimize risks by improving the adaptive capacity of urban form to reduce vulnerabilities and exposures, thereby coping with, surviving, and recovering from hazards (Dhar and Khirfan, 2017). Specifically, improving the adaptive capacity of town plans and land uses can enhance urban form's flexibility to absorb unknown climatic events with uncertain patterns, such as flooding ensuing from extreme precipitations. This improved adaptive capacity can ensure that the urban form maintains its functions and structure, contributing to urban form resilience (Dhar and Khirfan, 2017;Khirfan and El-Shayeb, 2020).
Khirfan and El-Shayeb (2020) connect urban form adaptation and resilience by drawing on Meerow et al. (2016, p. 39) definition of resilience: "the ability of an urban system-and all its constituent socio-ecological and socio-technical networks across temporal and spatial scales to maintain or rapidly return to desired functions in the face of a disturbance, to adapt to change, and to quickly transform systems that limit current or future adaptive capacity." Accordingly, adaptation (and adaptive capacity) is among the three pathways to resilience, along with persistence and transformation. Framing urban form adaptation under the umbrella of resilient planning has rendered resilience central to flood adaptation policies and projects to enhance the flexibility and adaptability of urban forms to increased rainfall events (Lennon, 2015;Graham et al., 2016;Shi, 2020;Shokry et al., 2020). This is evident in the shift in land use policies to integrate large-scale green projects that absorb and dissipate rainwater run-off (Anguelovski et al., 2019;Shi, 2020) and urban design interventions that incorporate resilient water-sensitive infrastructure in town plans to infiltrate, harvest, and convey rainwater (Watson and Adams, 2010;Matos Silva and Costa, 2016).
Despite its benefits for urban form adaptation to climate change-induced floods, resilient planning risks ignoring the underlying causes behind risk inequities and differences in adaptive capacity, perpetuating the historic uneven processes of urban development (Meerow et al., 2019). In particular, the application of resilience planning as a development agenda in the last decades has capitalized on branding cities as climateand flood-adaptive sites to encourage investments by the tourism industry, real-estate developers, and the new sustainability class (Connolly, 2019;Anguelovski et al., 2020;Garcia-Lamarca et al., 2021). The exclusionary controls over the types and locations of investments can exacerbate vulnerabilities, whether through increased land values/rents and ensuing climate gentrification (Chu et al., 2017;Shi, 2020) or through the forced relocation of marginalized groups to clear space for large-scale projects (Henrique and Tschakert, 2019). Accordingly, enhancing urban form resilience through flood-adaptive land use planning and town plan design risks excluding vulnerable neighborhoods that already lack sufficient adaptive capacity .

THEORETICAL FRAMING: HOW TO COMBINE URBAN FORM ADAPTATION WITH DISTRIBUTIVE JUSTICE
To investigate how the adaptive capacity of urban form is connected to differential vulnerabilities and exposures to floods, and how urban form adaptation responses should be distributed to advance climate justice, the theoretical framework of this study combines Dhar and Khirfan's (2017) urban design resilient index (UDRI) the distributive justice notion.

The UDRI Framework
We draw on Dhar and Khirfan's (2017) UDRI framework to assess and compare the adaptive capacity of urban form in different urban neighborhoods to identify the disadvantaged ones whose adaptive capacity is also low. We focus on this framework because it is clear, comprehensive, and generalizable; it is also applicable to Conzen's (1960) urban form components, particularly land uses and town plans. The framework includes seven concepts that impact the resilience, hence the adaptive capacity of urban form across functional, spatial, physical, and temporal dimensions; they are harmony with nature, polyvalency, heterogeneity, connectivity, indeterminacy, latency, and modularity (refer to Table 1 for definitions).
While Dhar and Khirfan (2017) developed their UDRI framework for measuring the resilience of urban form at the neighborhood scale, this study applies it at the urban scalethat is, for the entire city. To facilitate this, we draw on only four of the seven concepts in the UDRI, namely: harmony with nature, polyvalency, heterogeneity, and connectivity, for which we found empirical evidence of their application at the city scale (see Table 1). Furthermore, we add a fifth concept, flexibility, due to the numerous theoretical and empirical debates regarding its application in assessing the general resilience of urban form at the city scale (Roggema, 2014;Sharifi, 2019a;Freire and Monteiro, 2020), particularly with regards to flooding events (Sharifi, 2019c).
While Dhar and Khirfan (2017) applied their resilient concepts to all three of Conzen's (1960) urban form components (i.e., land uses, town plans, and the threedimensional (3D) built form), we apply the five concepts only to land uses and town plans. Our reason for this is the lack of data that facilitate measuring the adaptive capacity of the 3D urban form elements to flood risks at the city scale.
Accordingly, beginning with land uses, we consider that their adaptive capacity can be enhanced through the configurational characteristics of harmony with nature, heterogeneity, and polyvalency. Land uses in harmony with nature have a minimal impact on the natural environment and can mitigate climatic hazards by strengthening ecosystem functions. One of the prevalent ways to enhance harmony with nature through land uses for adaptation to flooding is by integrating green and blue infrastructure (GBI). GBI refers to an interconnected network of natural (e.g., lakes, streams, and parks) and seminatural ecosystems (e.g., community gardens and green roofs) that benefit humans through providing ecosystem services (Bolund and Hunhammar, 1999;Mohtat and Khirfan, 2021). GBI can mimic natural hydrological processes such as infiltration, evapotranspiration, retention, detention, and slow flow (Liu et al., 2019) that collectively promote nature-based solutions for adaptation (IPCC, 2022). Therefore, several studies have introduced GBI as a decentralized approach for managing the excess rainwater and regulating flooding, which can supplement the centralized urban drainage gray infrastructure (Abebe et al., 2018;Li et al., 2020). Heterogeneous land uses, through the variation of their types over a spatial unit, facilitate the spread and dissipation of hazards across space. For instance, urban forms that include a rich combination of land use types with different porosity (e.g., open spaces, industrial uses, green spaces, and residential uses) are better able to dissipate rainwater run-off (Cadenasso et al., 2013;Dhar and Khirfan, 2017;Zhou et al., 2017). Polyvalent land uses allow a change in functions without significant physical changes to accommodate hazards (Dhar and Khirfan, 2017). For example, recreational spaces adjacent to rivers can become spaces that temporarily accommodate floods (Macintosh, 2013).
As for town plans, we consider that their adaptive capacity increases when they are flexible and connected. Flexibility refers to the urban form's ability to integrate future changes and interventions for adaptation; hence, it bears some similarities with Dhar and Khirfan's (2017) latency and indeterminacy concepts (Table 1). However, the difference is that Dhar and Khirfan (2017) used latency and indeterminacy at the micro scale in relation to small spaces adjacent to streets and intersections while Sharifi (2019a,c) used flexibility to elaborate on the characteristic of adaptable urban form at the larger scale, such as urban blocks. Therefore, we replace latency and indeterminacy with flexibility. Flexibility and connectivity often go hand in hand. Flexible town plans facilitate accommodating adaptive interventions and incorporating land modification regulations. For example, integrating green spaces in fined-grained urban blocks is easier and more cost-effective than large-grained ones with little connectivity (Salat, 2017;Sharifi, 2019a,c). Connectivity enhances the town plans' permeability by increasing the contact between blocks with streets. It, therefore, accelerates access to buildings and emergency management in the advent of intense rainfall events leading to run-off flooding (Sharifi and Yamagata, 2014;Sharifi, 2019b).

Distributive Justice
Distributive justice refers to the just spatial/temporal distribution of resources to maximize benefits to the most vulnerable (Rawls, 1971;Adger, 2006a;Shi et al., 2016). Building on Rawls's (1971) liberty and maximization rules, distributive justice gives those with the greatest need the right to equal access to resources and the priority in their spatial allocation (Sen, 1992;Adger, 2006a). We draw on the distributive justice notion to identify how different the neighborhoods' urban forms are shaped in  Dhar and Khirfan (2017, p. 83-84) and the concepts from this framework that this study uses.
The UDRI framework (Dhar and Khirfan, 2017) Evidence on how to apply concepts to: The concepts used in our theoretical framework

Theoretical Framework
Our theoretical framework connects these disparate notions whereby the identification of high-risk neighborhoods ensues from the simultaneous presence of four flood risk drivers: (1) exposures to flooding hazards; (2) social vulnerabilities; (3) low adaptive capacity of land uses; and (4) low adaptive capacity of town plans (Figure 1). Drawing on our interpretation from the UDRI framework, we assess the adaptive capacity of land uses based on their degree of harmony with nature, heterogeneity, and polyvalency while evaluating the adaptive capacity of town plans in terms of their connectivity and flexibility.

METHODOLOGY
To operationalize our theoretical framework, our methodology started with developing a conceptual framework that includes indicators and variables for measuring the four flood risk drivers (hereafter, we dub this conceptual framework "the multicriteria model"). We then conducted overlay analysis in ArcGIS using the experts' assigned weights.

Conceptual Framework: Our Proposed Multicriteria Model
Several studies propose multicriteria models (MM) to identify the spatial distribution of flood risks and their drivers, hence the priority areas for adaptation responses. In most existing studies, MMs include physical factors that cause flood hazards and exposures, such as slope, elevation, rainfall, and soil types (Lin et al., 2019;Ogato et al., 2020). However, less attention is paid to the unequal spatial distribution of adaptation interventions, hence differences in the adaptive capacity of land uses and town plans across neighborhoods and their connections to differential vulnerabilities and exposures. The few empirical studies that connect urban form adaptation with differential vulnerabilities and flood exposures also consider the unequal access of vulnerable groups to GBI to identify priority areas for just adaptative interventions-see Meerow and Newell (2017) and Li et al. (2020). However, they overlook frameworks like the UDRI that take into account the configurational characteristics of resilient urban form. Thus, we operationalized our theoretical framework to address this deficit by proposing a MM whose indicators and variables tackle the spatial distribution of four co-existing risk drivers: (1) flood hazard exposures; (2) social vulnerabilities; (3) low adaptive capacity of land uses (due to a lack of harmony with nature, heterogeneity, and polyvalency); (4) low adaptive capacity of town plans (due to a lack of flexibility and connectivity) -(see Figure 1 and Table 2). Our MM, in total, includes 38 variables, which measure 15 indicators per neighborhood as the unit of analysis, whereby the City of Toronto has defined the neighborhood's boundaries since the 1990's to facilitate collecting data, planning, and analysis 2 (City of Toronto, 2019). We mapped all the variables in ArcGIS. We normalized the variables' values from zero to ten (using linear scale transformation) to make them comparable to and combinable with each other (Meerow and Newell, 2017;Lin et al., 2019;Li et al., 2020). We calculated the average of variable values to map each indicator. The following sections explain in detail each risk driver's indicators and variables and the data sources.

Exposure to Flood Hazards
To identify the exposure of each neighborhood to flood hazards, our MM proposes two indicators: "proximity to flood plains" (Chakraborty et al., 2014;Lyu et al., 2016) and "run-off coefficients" (Thompson, 2006;Meerow and Newell, 2017;Li et al., 2020). We measured the proximity to floodplains by calculating the percentage of land covered by floodplains in each neighborhood using the Floodplain Mapping Index data (TRCA, 2020b) and the Intersect Analysis tool in ArcGIS. We estimated the average Run-off coefficients for each neighborhood, using Thompson's (2006) rational method (also see Li et al., 2020). We first estimated the average area of land covered by land use categories in the rational approach, using the land use data (City of Toronto, 2020e) and the Intersect Analysis tool in ArcGIS ( Table 3). We then multiplied the percentage values with their relevant coefficient amount to calculate the average estimated amount for each neighborhood.

Social Vulnerabilities
We adopted Chakraborty et al.'s (2014) social vulnerability index to measure the vulnerability of Canadians to floods. Chakraborty et al. (2020) developed this index's indicators and variables based on theoretical debates, policy documents, and Canadians' demographic characteristics across census tracts. They used several statistical approaches to assure the index's generalizability, validity, and replicability. Hence, it is reliable enough to represent Canadians' socioeconomic characteristics, making it an appropriate tool for measuring social vulnerabilities across Toronto in this study. Accordingly, we considered "age, " "gender, " "wealth, " "ethnicity, race, and immigration status, " "employment status, " "family structure, " "education, " and "builtenvironment conditions" as social vulnerability indicators (see the full list of indicators and variables in Table 2). We extracted all the variable values from the Neighborhood Profiles, which the City of Toronto has built based on the 2016 census data (City of Toronto, 2019).

The Adaptive Capacity of Land Uses
As our theoretical framework indicates (Figure 1), when it comes to assessing the adaptive land uses, this study draws on the three indicators of "harmony with nature, " "polyvalency, " and "heterogeneity" ( Table 2). FIGURE 1 | The theoretical framework. *Flexibility is not included in the UDRI framework. Dhar and Khirfan (2017) have proposed that the larger the amounts of land covered by natural porous surfaces, such as GBI, the higher harmony with nature of land uses. Thus, we considered the percentage of land covered by green and blue spaces and the density of street trees as variables for measuring the harmony with nature and the adaptive capacity of land uses.
To measure heterogeneity, or the spatial differentiation of land uses, we calculated the values of variables proposed by Cadenasso et al. (2013) per neighborhood. These variables include: (1) the number of land use patches; (2) patch richness, in reference to the number of different land use patches such as commercial, residential, and institutional; (3) the frequency of different patch types, referring to the number of times each land use patch appears in the urban landscape ( Table 2). Note that the more the variables' values, the higher the urban form's ability to spread and mitigate climatic hazards like floods.
Several studies have referred to open spaces and mixeduse developments as polyvalent (or multifunctional) land uses The percentage of 65 years old and above population +

Gender
The percentage of females (15 years old and above) who participate in the labor force a The percentage of people (25-64 years old) whose highest degree is a secondary (high) school diploma or equivalency certificate.

Built-environment conditions
The percentage of households living in homes with need for major repair + The percentage of households with more than one person per room + The percentage of Labor Force (above 15) whose main mode of commute to work is public transportation g +  The number of Utility and Transportation patches -*Flexibility is the additional fifth concept of our theoretical framework. a Structural gender inequality causes female workers to suffer more than their male counterparts from unstable working conditions and low income (Kalev and Deutsch, 2018), reducing their access to assets to cope with risks. b Visible minority refers to "persons, other than Aboriginal peoples, who are non-Caucasian in race or non-white in color" (Chakraborty et al., 2020, p. 4). c First generation refers to "persons born outside Canada. For the most part, these are now, or once were, immigrants to Canada" (Chakraborty et al., 2020, p. 4). d Aboriginal identity relates to "persons who are First Nations (North American Indian), Métis or Inuk (Inuit) or those who are Registered or Treaty Indians or those who have membership in a First Nation or Indian band" (Chakraborty et al., 2020, p. 4). e Male not in the labor force refers to male persons "who are unwilling or unable to offer or supply labor services under conditions existing in their labor markets (including persons who were full-time students currently attending school)" (Statistics Canada, 2008). Many cultures consider males as the main persons who financially support families. Therefore, this variable can indicate the unfavorable financial condition of households, which reduces their access to adaptive resources. f Unemployed persons are those "without work, are available for work and are actively seeking work" (Statistics Canada, 2008). g Flooding events can damage public transportation infrastructure (such as subways), leading to the closure of public transit systems and delays (Nirupama et al., 2014) and adversely affecting those who depend on them. h The Canadian building codes before the 1980's were not strict enough to include emergency conditions (Archer, 2003).
that can accommodate floods and provide space for erecting emergency shelters (Macintosh, 2013;Roggema, 2014;Sharifi, 2019a). Therefore, we considered the percentages of areas covered by these land use types per neighborhood as variables to measure polyvalency. For all the indicators, we used the Zoning By-Law data provided by the City of Toronto (2020e). We used the Intersect Analysis and Summary Statistics tools in ArcGIS to map all the indicators. In addition, the Dissolve and Merge tools in ArcGIS were used for analyzing the third indicator.

The Adaptive Capacity of Town Plans
Building on our theoretical framework, we draw on two indicators of "flexibility" and "connectivity" ( Table 2) to measure the adaptive capacity of town plans.  (Thompson, 2006;Li et al., 2020). Salat (2017) and Sharifi (2019c) proposed that fine-grained blocks and building footprints are more flexible than largegrained blocks to accommodate changes, such as through smallscale adaptive interventions for incremental adaptation at a lower cost. Furthermore, they can accelerate emergency responses in the advent of flooding disasters by providing opportunities for multi-use developments and enhancing access points at street edges. Thus, we compared the flexibility of town plans in different neighborhoods by calculating the average size of their blocks and building footprints, whereby the smaller the size, the higher the flexibility. Sharifi and Yamagata (2014), Feliciotti et al. (2016), and Dhar and Khirfan (2017) argued that the connectivity of town plans promotes the accessibility of blocks and buildings through street networks, thus facilitating evacuation planning, emergency search, and rescue activities in the advent of flooding disasters. As Feliciotti et al. (2016) proposed, the higher the number of three-and four-way intersections, the higher the connectivity. Therefore, we used the average density of street networks' intersections per neighborhood as the variable for measuring the connectivity of town plans; in other words, the higher the density, the higher the connectivity.

Land use categories Coefficient
For both indicators, we used the data provided by the City of Toronto (2020b,d). Furthermore, we used ArcGIS for the Intersect Analysis and Summary Statistics tools to produce the indicators' maps ( Table 2).

Weighted Overlay Analysis Through ArcGIS
Since risks result from intersections among multiple drivers with unequal importance, the existing GIS-based multicriteria approaches on flood risk mapping often involve weighted overlay analysis. Qualitative and mixed-method research studies like this one often use experts' judgments for weightings, such as through different approaches of rating and ranking -see: Meerow and Newell (2017), Rincón et al. (2018), Li et al. (2020). This weighting approach facilitates quantifying immeasurable data and responds to the challenges of data scarcity (Wang et al., 2011;Lin et al., 2019).
For this study, we conducted an online survey (using Qualtrics) to seek the experts' opinions regarding the weights of flood risk drivers (Figure 1) and their associated indicators ( Table 2). Our survey population comprised Toronto-based planning experts who have experience in at least one of the fields of urban flood management, climate change adaptation, and/or adaptive urban form. We found these experts through a systematic search on Google, LinkedIn, and LinkedIn Premium. Our search yielded 392 relevant experts, working variously in four academic, 13 non-governmental, 27 governmental, and 44 private organizations. We shared the survey link with these experts through email and/or LinkedIn messaging from April to the end of June 2021. The survey eventually yielded 120 responses (31% response rate).
To ask the survey participants to weigh the flood risk drivers, we drew on the Analytic Hierarchy Process (AHP), which is a rational, accurate, cost-effective, and easy-to-use approach for measuring the importance of immeasurable elements through pair-wise comparisons (Lin et al., 2019). First proposed by Saaty (1990) for quantifying the weights of decision criteria, AHP became a popular approach for subjective evaluation of flood risk drivers in GIS overlay analysis -(see: Lin et al., 2019;Li et al., 2020;Ogato et al., 2020). Building on this approach, we asked the expert participants to pair-wisely compare the relative importance of the four risk drivers regarding the exacerbation of flood risks in Toronto with a scale that ranges from 1 (equal importance) to 9 (extremely more important) -see Saaty (1990) and Figure 2.
To interpret the data, we created a pair-wise comparison matrix (i.e., [C]) for each participant based on the fundamental AHP scale suggested by Saaty (1990): where m ij is the component of the normalized matrix. We eventually obtained the weight of the i th concept (1 ≤ i ≤ 4) as the average of each row in the normalized matrix: To evaluate the consistency of the survey responses, we calculated the Consistency Index as follows: where λ max is the maximum eigen value of the pair-wise comparison matrix [C] and n is the number of concepts that are compared (i.e., n = 4). According to Saaty (1990), a consistent matrix has a Consistency Index of <10% of the Consistency  Ratio (C.R.); where the value of the C.R. for a matrix with the size of four is proposed to be 0.9 by Saaty (1990). Considering these consistency criteria, we filtered the responses and calculated the weight of concepts corresponding to each participant. The final weights of concepts are the average of weights obtained for each participant.
While AHP approach is reliable for weighting the four risk drivers, it may become a lengthy task for weighting the 16 indicators due to a large number of pair-wise comparisons . Accordingly, we measured the weights of indicators through direct rating (DR) where expert participants assigned a weight (from 0 to 10) to the impact of indicators on each risk Using the weights assigned by the experts, we began our overlay analysis in two steps: overlaying the indicator maps to map their associated concepts and overlaying the risk drivers' maps to draw the final flood risk map. We used the Union Analysis tool and the weighted sum average function in ArcGIS to complete the weighted overlay analysis for both steps. Note that we normalized all the concepts and the final risk map values from 0 to 10, using linear scale transformation (see Figure 3).

THE CASE STUDY: TORONTO
Toronto is Ontario's capital and Canada's foremost economic hub. Toronto spreads over 633.5 km 2 , and its population totals 2.73 million (in 2016), 50% of which are visible minorities, which makes it the most populous city in Canada and one of the most multicultural cities in the world (Filion et al., 2015;Statistics, 2017). The city's location within the Lake Ontario Watershed and its exposure to moist air masses and high precipitation rates have caused several historical flooding events that caused a loss of lives and damages to properties and infrastructure in 1878, 1954(after Hurricane Hazel), 1976, 2005(Nirupama et al., 2014Rincón et al., 2018;TRCA, 2021).
More importantly, there is evidence of increased precipitation rates due to global climate change in this city. Feltmate and Thistlethwaite (2012) mentioned that six 50-year and two 10year precipitation events had been recorded during just 15 years -from 1996 to 2011. Over the last decade, governmental, non-governmental, academic, and private organizations at the different municipal, regional, and provincial levels have developed stormwater management plans, policies, and/or strategies to promote climate change adaptation .
In Toronto, the City's Water Division oversees developing policies and implementing strategies for stormwater management and climate change adaptation. After the approval of the Climate Change Action Plan in 2007, the City of Toronto published its first climate adaptation strategies, including actions on flood protection and emergency management, in a document titled "Ahead of the Storm: Preparing Toronto for Climate Change" (City of Toronto, 2008). Following this document, the City continued to work on its first Resilience Strategy, which includes 50 major plans, such as the Basement Flooding Protection Program and Wet Weather Flow Master Plan and Management Guidelines (City of Toronto, 2017Toronto, , 2020c. The City of Toronto collaborates with other organizations as well. Vertically, it works with regional and provincial governmental organizations such as the Toronto and Region Conservation Authority (TRCA). Being one of the 36 conservation authorities in Ontario, the TRCA receives funds from municipalities to offer them information on flood mapping, educational workshops, awareness programs, and low impact development (LID) design guidelines (TRCA and CVC, 2010;Henstra and Thistlethwaite, 2017;TRCA, 2020a,b). At the provincial level, the City receives advice from the Ministry of Environment, Conservation, and Parks (MECP) and the Ministry of Environment and Climate Change (MECC), among others (City of Toronto, 2020c;Henstra et al., 2020). Horizontally, non-governmental organizations (e.g., Toronto Environmental Alliance), private firms (e.g., Metrolinx), and academic institutions (e.g., Intact Center for Climate Change Adaptation) assist the City of Toronto in conducting feasibility assessment projects and developing strategies and standards (Mah et al., 2018;Metrolinx, 2018;City of Toronto, 2020c).

RESULTS: HOW ARE FLOOD RISKS DISTRIBUTED?
Mapping the Risk Drivers

Exposure to Flood Hazards
In terms of exposure to flood hazards, the survey results show that the average weights of proximity to flood plains and run-off coefficients are 0.47 and 0.53 ( Figure 4A). Accordingly, experts believe that the run-off coefficient has a slightly higher impact on exposure to flood hazards in Toronto than proximity to flood plains. When we overlay the indicator values in ArcGIS, using their assigned weights, the results show that Flemingdon Park neighborhood, followed by West Humber-Clairville, and Morningside, are the most exposed to flood hazards ( Figure 5A and Table 4).

Social Vulnerabilities
With regards to social vulnerabilities, the survey results show that wealth and built-environmental conditions (weighted at 0.16 each) have the greatest impact on social vulnerabilities, while gender (weighted 0.09) is the least impactful ( Figure 4B). In addition, the overlay analysis of indicators' values by using their weights shows the disproportionate spatial distribution of social vulnerabilities within the city. In this respect, Black Creek neighborhood followed by Oakridge and Thorncliffe Park have the highest social vulnerability to floods ( Figure 5B and Table 4).

The Adaptive Capacity of Land Uses
When it comes to the adaptive capacity of land uses, the survey results reveal that harmony with nature followed by polyvalency (weighted 0.36 and 0.34, respectively) have the highest impacts. In contrast, heterogeneity (weighted 0.30) maintains the minimum impact on land uses ( Figure 4C). After overlaying these indicators' values (using their assigned weights), the results show that land uses in the North St. James Town neighborhood followed by Willowdale West and South Riverdale have the lowest adaptive capacity ( Figure 5C and Table 4).

The Adaptive Capacity of Town Plans
Last, with regards to the adaptive capacity of the town plans, the survey results reveal that flexibility and connectivity (weighted 0.49 and 0.51, respectively) have relatively similar impacts on the adaptive capacity of town plans ( Figure 4D). The results of our weighted overlay analysis using ArcGIS show variation among Toronto's neighborhoods in terms of the adaptive capacity of their town plans. As shown in Figure 5D and Table 4, Thorncliffe Park followed by Humber Summit, York University Heights, and Downsview-Roding have the lowest adaptive capacity in their town plans.

Mapping the Final Flood Risk Map: Identifying the Priority Neighborhoods
To map the final flood risk map and to identify which of Toronto's neighborhoods should be prioritized for adaptation interventions, we overlaid the maps of risk drivers (Figure 5) using the weights assigned by the experts. As the experts' survey results show (Figure 6), social vulnerabilities (0.32) have the highest impact on flood risks in Toronto, while exposure to floods (0.22) has the lowest impact. In addition, the experts believe that the adaptive capacity of land uses (0.23) and town plans (0.23) have similar impacts on flood risks. The results of the weighted overlay analysis reveal that Thorncliffe Park followed by Flemingdon Park, North St. James Town, and Black Creek are four neighborhoods that are disproportionately at risk from flooding, hence, must be prioritized in urban form adaptation interventions (Figure 7 and   2008; Hassen, 2021). The unfavorable conditions of the built environment, the concentration of poverty, and the impervious surface materials with high run-off coefficients are the main reasons behind the vulnerability to increased precipitation and exposure to increased flooding (Figure 8).
More importantly, our analysis of indicator values in Figure 8 shows inadequate adaptive capacities of these neighborhoods' land uses and town plans. The high-rise developments and the separation of land uses have resulted in a lack of landuse heterogeneity as well as urban form's large-grained blocks and disconnected streets (Figure 8). Although the "towers in the park" urban form includes ample open green spaces, other factors reduce the urban form's flexibility to incorporate future changes and its ability to spread run-offs. This is due to the discriminatory policies, lack of maintenance, and, more recently, infill development that have decreased and continue to reduce the quantity and quality of open green spaces and the land uses' harmony with nature. For example, North St. James Town has the lowest area of green space per person in Toronto since new apartment complexes have replaced open green spaces between the towers over time (Nguyen, 2014;Hassen, 2021) -(see Figure 9). Additionally, in Thorncliffe Park, Metrolinx 3 plans to replace some open spaces and business buildings with the train yards of the Ontario Line (Aecom Canada Ltd., 2021). These plans have raised the concerns of grassroots environmental justice activists (Savetpark Community, 2021).
Yet, social vulnerabilities remain the most critical in triggering flood risks, particularly the lack of access to wealth when combined with unfavorable built-environmental conditions (Figure 6). Similarly, the run-off coefficient and harmony with nature are the most important indicators of flood exposure and adaptive capacity (Figure 4). Accordingly, we call for future theoretical and empirical studies to investigate how GBI interventions and nature-based solutions can address the root causes of vulnerability in tower communities in Toronto and elsewhere while advancing just adaptation to flooding. Furthermore, we propose that future research explores how low-income and disadvantaged communities and marginalized groups can participate and integrate their needs in the design and implementation of small-scale GBI in a way that curbs, if not altogether avoids, gentrification by maintaining housing affordability.
More importantly, our findings show that the need to prioritize tower communities for just adaptation to changing climate may not be specific to Toronto but applies globally. As they age and dilapidate over time, tower buildings that once were modern 20 th century housing types have become the 21 st century's affordable housing enclaves for low-income, marginalized, and vulnerable communities, particularly in Western Europe and North America. Over the last two decades, municipalities around the world have proposed strategies to advance social equity and to improve the conditions of the built-environment in similar tower buildings, whether through renovations, public realm improvements, mixed-use developments, and/or the integration of urban agriculture (E.R.A. Architects University of Toronto, 2008;Benko et al., 2018;Veschambre, 2018). Some of these improvements include climate mitigation strategies (i.e., decreasing greenhouse gas emissions through improved energy efficiency) (Aragon et al., 2018;Seebauer et al., 2019). Yet, there is a need for studies that inform both research and policy on the adaptation of tower neighborhoods to climatic events including flooding through participatory processes that are grounded in context-specific needs and the local communities' lived experiences as well as the knowledge of local experts.

CONCLUSION
This study proposed a multicriteria model whose variables and indicators assess the spatial distribution of social vulnerabilities, flood hazard exposure, and urban form's adaptive capacity to facilitate an assessment of "who" are unequally at-risk to flooding events, hence, should be prioritized in adaptation interventions; "where" are the high-risk priority areas located; and "how" urban form adaptive interventions may advance climate justice in these locations. Specifically, this model changes how risk inequalities are understood by combing sociodemographic indicators with five configurational characteristics of resilient and adaptive land uses and town plans: harmony with nature, heterogeneity, polyvalency, flexibility, and connectivity. We tested the model in Toronto, through weighted overlay analysis using ArcGIS and an online survey of 120 Torontobased flooding experts, to identify how social vulnerabilities, flood exposures, and adaptation interventions are distributed within Toronto's urban form. This information enables us to identify which neighborhoods are experiencing the highest risks of floods.
The results reveal the uneven spatial distribution of flood risks, hence, identify four neighborhoods that should be prioritized for adaptation interventions: Thorncliffe Park, Flemingdon Park, North St. James Town, and Black Creek. Indeed, these are innercity, high-density tower communities with old infrastructure and low-income, racialized, and migrant populations -typical of the 20 th century modern tower block communities dotted across North America, Europe, and Asia. This study was part of a bigger project. Building on the experts' surveys, the following steps include working with the vulnerable communities through participatory and interactive processes to develop small-scale adaptive GBI solutions grounded on place-based experiences, representing the neighborhood residents' everyday lived experiences. Surely, as more empirical studies investigate the root causes of climate related risks in tower communities beyond Toronto and Canada, we will learn more about why certain communities will need to be prioritized in adaptation interventions and how we can work with them to advance just climate solutions that are grounded in the communities' contextspecific needs.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by the University of Waterloo Research Ethics Board (ORE #42887). The patients/participants provided their written informed consent to participate in this study.

AUTHOR CONTRIBUTIONS
NM: research design, investigation, conceptualization, methodology, data curation, formal analysis, literature review, and writing-review and editing. LK: research design, conceptualization, writing-review and editing, supervision, funding acquisition, and project administration. Both authors contributed to the article and approved the submitted version.

FUNDING
The funding for this research project was provided by the Social Sciences and Humanities Research Council of Canada (SSHRC) under file number 435-2016-0243. Dr. Amrita Daniere funded the open access publication expenses of our paper. The funders had no role in our research, including manuscript preparation, data collection, and analysis.

ACKNOWLEDGMENTS
We are sincerely grateful to the experts who shared their valuable insights through surveys. We would express our deep gratitude to Dr. Alireza Najma for his assistance with analyzing the survey results. Last but not least, we would like to express our special thanks to Dr. Amrita Daniere, who kindly funded the open access publication expenses of our paper.