Computerized constructed wetlands allocation model (based on ecosystem services demand)

Egypt’s urban communities face many threats, including, pluvial floods, heat waves, and lack of publicly accessible urban green spaces. Nature‑based solutions such as constructed wetlands (CWs) present a promising solution that could offer a wide range of ecosystem services (ES). However, the adoption of CWs is challenged by the lack of local planning guides and uncertainty about potential benefits. There are various models and tools available for quantifying and valuing ES, however, many of them are either highly complex or require extensive data and expertise. The aim of this paper is to develop a GIS‑based multi‑criteria decision model to select CW sites based on the supply and demand of ES. The model is to operate on three main stages: (i) demand: based on the need for risk reduction or benefit provisioning, (ii) potential sites (PSs): identify and score suitable sites for establishing a CW, and (iii) supply: define the service benefiting areas (SBA). An experimental approach is used, where the model is tested on New Damietta, an Egyptian Mediterranean city, proving the model is a reli‑ able decision‑making tool during preliminary urban planning stages due to its practi‑ cality, flexibility, and reasonable data requirements.


Introduction
People living in urban areas exceed 50% of the world's population and continue to grow [1].However, urban expansion exceeds urban population growth, particularly in lowelevation coastal zones [2].Between climate change, and anthropogenic activities altering the natural ecosystems, urban coastal communities face growing risks including sea level rise, erosion, degradation of natural ecosystems (e.g., wetlands), stormwater runoff, air and water pollution, and psychological stresses [3,4].
To overcome these risks and establish more sustainable and resilient communities, previous literature highlighted that urban designers and planners must shift away from conventional planning and urban management approaches.Consequently, various new planning and management-related approaches appeared such as blue-green infrastructure, water-sensitive urban design, and nature-based solutions [5].These approaches call for smarter and more resilient management of urban systems, particularly the urban water cycle [6].
Moreover, these approaches pursue multifunctionality in the form of Ecosystem Services (ES) Provision [7].ES could be defined as the direct and indirect benefits to humans provided by planned and/or natural ecosystems [8].To supply ES, nature-based solutions veer away from heavily engineered solutions to softer, greener, and more natural ones that reintegrate nature into the urban realm.Such solutions include rain gardens, green roofs, bio-swales, retention ponds, and constructed wetlands [9], which is the focus of this study.
This paper adopts a comprehensive research methodology aiming to achieve the paper objective comprising three main stages.The first stage is a deductive research based on a review of the relevant literature discussing the concept of CWs, their ES, and the factors impacting the demand for and delivery of these services.Based on the results of the literature review, criteria were deduced for the modelling of ES demand and supply.A number of criteria were defined based on data availability, applicability and reliability of results.
In the second stage, the research devised a tool that could help decision-makers prioritize sites for CWs to ensure optimal ES delivery.A correlation analysis was conducted, integrating the previous criteria into a multi-criteria model with specific ranges and classes set for each criterion based on data from the literature.The criteria used for modeling ES demand were based on the need for risk reduction or benefit provisioning.
While service supply extent of potential CWs sites was used as the criterion for modelling supply of each ES.
The third stage is the empirical phase, where an experimental research design was adopted for investigating the model applicability.It was verified through a scenario analysis applied to a selected case study area, the city of New Damietta, Egypt.The results of the spatial modelling were analyzed through descriptive statistics.The outcome of this study is an innovative spatial multi-criteria decision model for allocating CWs within the urban fabric to optimize ES delivery.

Constructed wetlands
Wetlands are transitional zones between aquatic and terrestrial ecosystems, which are covered by water or have waterlogged soils for a significant part of the vegetation-growing period [10].Constructed wetlands (CWs) are engineered landscape elements that exploit natural processes (natural vegetation, soil, and micro-organisms) to achieve the required functions [11].CWs have the potential to supply multiple ES such as runoff and wastewater treatment, temperature regulation, and recreational services [12].Moreover, CWs are resilient and economically viable during the construction, operation, and maintenance phases [13].
Due to the high feasibility and extensive potential ES supplied by CWs, urban designers and landscape architects have shown growing interest in their integration into the urban landscape [10].The provision of ES is strongly correlated with CW design characteristics including hydrology, size, bed media, and macrophytes [14].Contextual properties such as climate, geology, topography, hydrology, and adjacent land use not only impact CW performance but also determine the suitability of the system within its context [11,15].Thus, CWs design is a complex endeavor that requires interdisciplinary efforts to ensure reliable results.This study is delivered from an urban planning perspective; therefore, it investigates the spatial flow of ES from service provisioning areas (potential CW sites in this case) to service benefiting areas (SBA), which is the service supply extent [16].
While CW characteristics and contextual properties influence the "supply" of ES, it is vital to recognize where "demand" for these ES is highest to ensure optimum impact [17].Consequently, this study aims to develop a spatial multi-criteria decision model to guide CWs site selection to facilitate optimum ES delivery, with a focus on three pressing risks facing urban coastal communities in Egypt: (a) runoff control, (b) climate regulation, and (c) recreational services.

The constructed wetlands allocation model
The suggested spatial model utilizes geographic information system (GIS) and remote sensing, which are widely used tools in urban-related studies used for entering, processing, integrating, displaying, and analyzing data from various sources [18].The suggested GIS-based planning model applies a multicriteria decision approach which will help decision-makers, urban designers, and planners prioritize sites where CWs would be most effective and facilitate on-ground implementation.The developed model operates through three stages (Fig. 1): i. Demand: map the levels of demand for each ES.ii.Potential sites: identify potential sites (PSs) within an urban context and score their suitability for establishing a CW.iii.Supply: define service benefiting area (SBA) for each ES provided by the PSs.

Model first stage: demand
The first stage of the model is to map demand for ES based on the need for risk reduction or benefit provisioning.The following sections describe the criteria used for mapping demand for each of the three ES currently under investigation.All criteria are classified into five classes and scored as follows: very low (0.2), low (0.4), medium (0.6), high (0.8), and very high (1).

Runoff control demand
Runoff generation is mainly driven by lower infiltration rates due to extensive use of impervious surfaces [14].Despite Egypt's relatively low precipitation (1 to 250 mm), its cities are vulnerable to runoff and flood risks [19].Demand for runoff control could be determined by the anticipated level of runoff risk which is a function of runoff hazard and vulnerability [20].For the purposes of this study hazard will be mapped based on runoff depths calculated through the NRCS-CN method, and vulnerability will be mapped based on four criteria identified from the literature: i. Elevation, lower altitude lands are naturally more vulnerable to runoff [21].ii.Slope strongly impacts surface runoff, especially runoff velocity, therefore, steeper slopes mean higher runoff vulnerability [22,23].iii.Distance to streams, where the closer a location is to natural streams the more vulnerable it is to runoff risk [24].iv.Land use depicts activities taking place within a particular spatial context and has been used as a runoff vulnerability parameter in various studies [25].
Runoff depth (hazard): the NRCS-CN method is utilized to analyze the rainfall-runoff processes in agricultural and urbanized watersheds [26,27] and estimate the proportion of precipitation that becomes surface runoff for a specific rainfall event [22,23].It is among the most used surface rainfall-runoff models [28], due to the simplicity of its structure and its dependency on a few basic parameters [28,29].
The NRCS-CN method relies on three data inputs: (i) precipitation data for a selected storm or event, (ii) land cover data (LC), and hydrologic soil group data (HSG) [26].Surface runoff depth is calculated through Eq. ( 1): Where Q is the runoff (mm), P is the precipitation depths in (mm), S is the potential retention capacity (mm), λ is the dimensionless initial loss coefficient with an average value of 0.2 [28].However, when λ is small, it is assumed that there is no surface runoff.The value of S could vary significantly between 0 and 25,146) and is estimated using the average CN for each subbasin as per Eq. ( 2): where curve number (CN) is a dimensionless parameter [28] that denotes the surface runoff potential of a certain land cover (LC) and hydrologic soil group (HSG) maps.A weighted CN, which is a weighted average of the CNs for each basin of the study area, is then calculated using Eq. ( 3): where CN w is weighted CN; CN i and A i are the curve number and area (respectively) for each LC-HSG polygon.

Climate regulation demand
Urban areas commonly exhibit higher temperatures than adjacent unurbanized areas and urban coastal areas are no exception [30,31].This is triggered by increased artificial surfaces, low evapotranspiration, fossil fuel combustion, and low ventilation due to high building density [23].Higher temperatures negatively impact human health, inhibiting outdoor recreation and activity, while also increasing energy and water consumption [3, (1) Ai 32].Based on the literature, three proxy criteria were identified to map climate regulation demand using remotely sensed data: i. Land surface temperature (LST) assesses thermal emissions radiated by the retention of solar energy by different surfaces [33].Various studies have found a strong correlation between LST and air temperature, and thus have been used to map cooling effects of blue-green infrastructure [15,34,35].ii.Normalized difference vegetation index (NDVI) is a common indicator of vegetation density.NDVI values range from + 1.0 to − 1.0.where values close to zero (− 0.1 to 0.1) generally correspond to barren land, while high values indicate dense vegetation.Vegetation cools urban areas through shading and evaporative cooling, therefore, areas with less vegetation are assumed to have a higher demand for cooling [23,30].iii.Normalized difference moisture index (NDMI) is a reliable indicator of water stress in vegetation.NDMI values range from + 1.0 to − 1.0, where negative values could mean water stress or bare soil, while high positive values may indicate waterlogging or dense canopy cover.Water stress in plants impacts their ability to cool through transpiration.Therefore, water-stressed areas will have a higher demand for cooling and climate regulation.

Recreational services
The demand for recreational services provided by CWs is not associated with a need to mitigate a certain risk but rather the potential to provide a unique experience to users.Therefore, the demand for this service is assumed to be constant throughout the study area where the goal is to provide residents with access to this service.

Model second stage: potential sites and site suitability score
After mapping demand for ES, the next step is to identify potential or suitable sites for CWs.Potential sites (PSs) are land parcels that are suitable for or fulfill the requirements for the implementation of CWs [11].Identifying PSs within an urban context could be challenging, especially in high-density areas.Based on the literature, five criteria were designated to assess the suitability of PSs for constructing a wetland: i. Area: there is no such thing as a typical CW [36].Russell et al. [11] suggest that the average size of a CW in an urban area is likely to be around 1000-3000 m 2 , while the overall needed space will be 2-3 times larger than this.ii.Soil salinity must be considered when selecting CWs sites as soil salinity impacts vegetation selection and could impact water treatment efficiency.iii.Slope, where sites with low or gradual slopes are preferred since they simplify the design and construction processes and minimize costs [36].Secondly, they provide better natural storage capacities, easily maintain permanent pool volumes, and offer longer infiltration times [37].
iv. Site accessibility is essential for construction, maintenance, and monitoring as well as potential future visitors.v. Soil permeability must be considered when selecting CWs sites to determine infiltration rate and storage capacity [36,38].Highly permeable soils will make it difficult to maintain the suitable soil saturation needed for efficient hydrological performance [36].

Runoff control supply
Runoff control service benefiting areas (SBAs) are mapped by delineating the sub-basin for each PS.This means using topographic data to define the area (sub-basin) where all the water (including runoff ) eventually flows through a particular outlet (in this case PSs).Risk reduction is then determined by extrapolating the number of residents within the SBA of each site.

Climate regulation supply
The climate regulation SBA is mapped according to the cooling range, which is the maximum distance a CW can cool [35].Cooling range depends on multiple parameters including wetland size, vegetation, urban form as well as environmental conditions [15,34].Various studies have arrived at different threshold values for cooling ranges of CWs varying from 60 to 500 m [15,34,35].This study will assume maximum potential and apply a uniform cooling range of 500 m.

Recreational services supply
Methods used to quantify recreational services usually use proxies such as number of people within a given distance from green spaces [39].The Egyptian National Organization for Urban Harmony [40] indicates that residents should have access to a 12,000-30,000 m 2 open area within a 1000 m distance at the district level, and 4000 m 2 area within a 400-m distance at the neighborhood level.Therefore, 400 and 1000 m were used to map recreational SBAs of CWs.

Methods
As previously stated, the suggested multi-criteria model aims to simulate the flow of ES by connecting demand and supply.The demand for each ES is simulated using multiple criteria retrieved from satellite and remotely sensed data.While supply is mapped as the service benefits area (SBA) modeled based on evidence found in the literature.Table 1 summarizes the criteria used in the suggested multi-criteria model, as well as the range, classification, and references for each criterion.

Case study analysis
This study implemented an experimental approach through a case study analysis of New Damietta, which is an Egyptian city (31°26′46′′N, E 31°40′6′′E) bordered from the north by the Mediterranean and by agricultural lands from the east-west and south.The city is close to two major Egyptian Coastal Wetlands; Manzala and Burullus (30 and 80 km respectively).The city has a Mediterranean climate with short rainy winters (November to April) and extended dry summers [41] (El-Bokl et al. 2015).It is divided into 6 residential districts, a coastal tourism and recreational zone, and an industrial zone [42] (Fig. 2c).New Damietta was selected as a case study mainly due to the uniqueness of its predevelopment condition which still affects the city to this day.Located at the fringe between the Nile Delta and the Mediterranean, the area used to be predominantly Sabkha, which is gradational between land and the intertidal zone where evaporitesaline minerals accumulate due to arid climate.Most of the Sabkha areas were dried to enable the construction and development of the city.However, lately, there have been plans to rehabilitate and fortify the coastal area of the city as part of a national Integrated Coastal Zone Management (ICZM) plan [43].This includes building lowcost dike systems across the coast to protect the city from flooding and erosion, as well as reclaiming some of the areas' natural land cover and providing recreational services [43].Therefore, this study aims to build on these plans and investigate the potential for integrating nature-based solutions, particularly CWs, further into the heart of the city not just the coastal zone.The next sections use satellite images and  remotely sensed data to map the demand for ES, followed by supply options.The analysis was conducted in ArcGIS© Desktop 10.5.

Runoff control demand
Runoff control demand was modeled by mapping runoff vulnerability and hazard criteria through multiple steps.First, digital elevation model (DEM) data from the Shuttle Radar Topography Mission (SRTM) was accessed from the United States Geological Survey (USGS) website for the date (11-8-2022) and processed by filling the sinks to eliminate minor imperfections in the data.The resulting data layer was the mapped Elevation criterion (Fig. 3a).Secondly, the slope was mapped from the previously processed DEM data using the "surface-slope" tool (Fig. 3b).Thirdly, the distance to streams was modeled by using the processed DEM data to generate "flow direction", "flow accumulation", and finally "stream order", and finally using the "Euclidean distance" tool to generate multiple buffer zones around the streams which were classified into the five ranges previously defined (Fig. 3c).Fourthly, land use data were obtained in shapefile format and reclassified into five classes as previously defined (Fig. 3d).All four runoff vulnerability criteria were classified into five classes and scored between 0.2 (very low) and 1 (very high).Then, total runoff vulnerability was computed using the "raster calculator" tool as the sum of the four vulnerability criteria and again reclassified into five classes (Fig. 3e).To model runoff depth (hazard), first, land cover and soil data (based on Ali and Shalaby, 2015) were combined and CN values were assigned for each pixel (Fig. 3f ).Second, the "basin" tool in ArcGIS was used to delineate major basins, and weighted CN was calculated for each (Fig. 3g).Third, daily rainfall data from 1981 to 2023 was retrieved from NASA's Power Data Access Viewer and used to calculate rainfall events with 2-and 10-year return period (14.33 and 33.54 mm respectively).These scenarios were used in the NRCS-CN equation to calculate runoff depth.Next, the runoff depth layer was reclassified into five classes and scored between 0.2 and 1.Finally, demand for runoff control was calculated by multiplying total vulnerability and hazard scores and then normalizing the resulting data (Fig. 3h).

Climate regulation demand
Landsat 8 data retrieved for (8-11-2022) from the USGS website was used to model LST, NDVI, and NDMI.The three data layers were classified into five classes according to the previously defined ranges.Demand for climate regulation was calculated using a raster calculator as the sum of the three demand criteria and then the resulting layer was normalized (Fig. 4).

Second stage: identifying potential sites and site suitability score
For the current study, a preliminary screening was conducted to identify available land parcels.Firstly, land use data was used to isolate parcels assigned: green areas, recreational, sports, vacant, or unassigned.Secondly, only parcels of 9000 m 2 or more were selected.Consequently, 20 PSs were identified, and the suitability score criteria were calculated for each (Fig. 5).Landsat 8 data was used to calculate Normalized Difference Salinity Index Soil as an indicator of soil salinity (Fig. 5a).The previously mapped slope layer was used but reclassified to match the ranges defined for site suitability (Fig. 5b).
Accessibility was mapped based on number of sides with access to roads and the road grade (Fig. 5c).Soil permeability was mapped using soil data from Ali and Shalaby (2015) (Fig. 5d).Then "zonal statistics" tool was used to extract mean values for each PS.Each site suitability criterion was divided into five classes and assigned scores from 0.2 (very low suitability) to 1 (very high suitability).Finally, the total site suitability score was calculated as the sum of the five site suitability criteria.

Third stage: mapping supply
Supply for each of the three ES under investigation was modeled by defining the SBAs of each PS.Firstly, for modeling SBAs of runoff control, the "Watershed" tool in Arc-GIS was used and the PSs layer was used as input in the "pour-points" field.The output was the delineated sub-basin for each site (Fig. 6a), which is the area where each site is expected to collect runoff based on natural surface flow.Secondly, for modeling SBAs of climate regulation, the "buffer" tool was used to draw the presumed 500 m cooling range of each PS (Fig. 6b).Finally, for modeling SBAs of recreational services, the road shapefile for New Damietta was used to create network analysis layer, then the "Service Area" tool from ArcGIS "Network Analyst extension" was used to define areas within 400 and 1000 m distance from each PS (Fig. 6c).

Runoff control: connecting demand and supply \* MERGEFORMAT
Table 2 shows the runoff depth (mm) for the 2-and 10-year return period rainfall events calculated using the NRCS-CN method for each sub-basin (SB) in the study area is displayed in Table 2.Moreover, the runoff volume (m 3 ) generated in each SB was calculated by multiplying the runoff depth with the area of each sub-basin.It appears that SBs 1, 2, 6, and 10 have the highest runoff depth, while SBs 5, 7, and 8 have the lowest.This corresponds to these SBs having the highest and lowest CN w respectively, where CN strongly impacts the runoff generation potential of each land cover type.The analysis also suggests that a day rainfall event of 14.33 mm (2-year return period) would generate 93,533 m 3 of runoff across the entire city.While an event of 33.54 mm of rainfall (10year return period) could generate 481,165 m 3 of runoff.By combining the runoff control demand map (Fig. 3h) with the population density map, it was possible to determine the number and percentage of the population under each demand class.\* MERGEFORMAT Fig. 7 demonstrates the percentage of the population and area of the city falling under each of the five runoff demand classes.The results suggest that the majority of the population lives in high (scoring 0.8 for runoff control demand level) and medium (scoring 0.6 for runoff control demand level) demand  Similarly, approximately 65% of the city's area is labeled as high and medium-demand zones.However, it is worth noting that most of the areas under high and very high runoff demand are concentrated in the industrial zone in the southern part of the city.While the lowdemand zones are in the areas under development where there is still no or low population density (Fig. 3h).
As for runoff control supply, by delineating the sub-basin for each potential CW site, it was possible to determine the potential runoff collected by each PS.\* MERGEFOR-MAT Fig. 8 shows the percentage runoff volume collected by each site for the 2-and 10-year events.The analysis suggests that the PSs could collect over 50% of the total runoff generated over the entire area of the city by the 2 and 10-year rainfall events (49,556 and 253,529 m 3 respectively).Sites 11 and 19 offer the highest runoff control, due to their location along high-level streams.While sites 3, 4, 6, 7, and 15 offer the least runoff control.However, it is worth noting that more detailed hydraulic modeling could impact the results.Also, using other runoff collection and transfer measures (such as bioswales) in conjunction with CWs could yield even better results.
Moreover, \* MERGEFORMAT Fig. 8 reveals the percentage reduction in population within each runoff control demand class by each PS.The analysis suggests that the combined impact of the 20 PSs could protect 36% of the population residing in very low runoff demand areas, 61.6% of those living in low-demand areas, 81% of those in medium and high-demand areas, and 33% within very high-demand zones.A breakdown of the percentage runoff control demand class can be found in Supplementary data.

Climate regulation: connecting demand and supply
Based on the results of the climate regulation demand analysis, \* MERGEFORMAT Fig. 9 demonstrates the percentage of the population and city area falling under each of the five demand classes.Where 2.6% of the population live in very high-demand areas, 41.7% in high-demand zones, 38.5% live in medium-demand zones, and the remaining live in low and very low-demand areas (4.4% and 12.8% respectively).
As for climate regulation supply, the SBA map was overlapped with the demand and population density maps to determine the percentage reduction in population within each climate regulation demand class ( \* MERGEFORMAT Fig. 10).The total combined impact of the 20 PSs was calculated by merging all the SBAs to avoid double-counting of population in zones where SBAs overlap.The results suggest that the combined effect of the PSs could cause a 76% cooling demand reduction in low-demand areas, a 71% reduction in low and medium-demand areas, and over 60% reduction in high and very highdemand areas.

Recreational services: connecting demand and supply
The network analysis area was combined with the population density area to demonstrate the percentage of the population within a 400-and 1000-m walking distance from each PSs ( \* MERGEFORMAT Fig. 11).The results suggest that the PSs could offer recreational opportunities for over 188 thousand residents within a 400-m walking distance (accounting for 37.3% of the population).In addition, almost 450 thousand residents are within 1000 m distance from the PSs (accounting for 88.7% of the population).

Site suitability score
Site suitability score was calculated using the sum of the five site selection parameters ( \* MERGEFORMAT Fig. 12).The analysis shows that the normalized total suitability score of the 20 PSs range between 0.4 and 0.8.Site ( 16) is the only one to score 0.8, while seven other sites scored 0.7.Only two sites had a low score of 0.4 (sites 1 and 7), these sites had small areas, high salinity, and high slope, which impacted their score.Four sites (9:12) had a medium suitability score of 0.5, where all four had small areas but high to medium salinity with an average slope.The remaining six sites scored 0.6.High-scoring sites should be prioritized when making decisions for establishing CWs within the city, coupled with the potential ES provision of each site.

Conclusions
This study discussed modeling the demand and supply of ecosystem services and their challenges, especially with the lack of field data and accurate information.Moreover, the study highlighted that despite the potential benefits and services offered by these green and blue systems it is difficult to accurately predict the performance of Naturebased solutions such as CWs, specifically in areas with no predecessors or prototypes to guide the design and planning processes.Therefore, the aim of this study was to introduce a simple and resilient model as a tool to assist planners and decision-makers in evaluating CWs to offer maximum benefit as urban planners and designers are the first link towards encouraging the uptake of these systems.
The suggested computerized constructed wetlands allocation model (CCWA) provided a GIS-based multi-criteria decision model that utilizes widely available remotely sensed and easily accessible satellite data to map and correlate demand and supply of ES.The model was tested on an Egyptian coastal city; New Damietta, to map three ES: runoff control, climate regulation, and recreational services.The case study findings showed that the majority of the population live in medium-to-very high demand zones for runoff control (83.5%) and that the 20 PSs could cut medium and high demand by 81% each, and very-high demand by 33%.Similarly, most of the population lives in medium-to-very high demand zones for climate regulation (82.8%), and the 20 PSs could reduce medium, high, and very high demand by 71, 61.4, and 68.8% respectively.Finally, 37.3 and 88.7% of the population are within 400 and 1000 m distance from the PS respectively.Moreover, the findings suggest that planning a network of CWs that covers the whole city is viable for optimum benefits.
In conclusion, the model was found to be practical where it can cover more ES and allow for more data layers to be added for more accuracy.The findings suggested that this model could be used and generalized to prioritize CWs or other nature-based solutions.However, further studies are recommended to validate the model and expand the list of ES covered by it.

2 ) 2 ( 1 Fig. 2 a
Fig. 2 a New Damietta land use plan, b city location, and c city main districts.Source: Prepared by authors using ArcMap 10.5

Fig. 3
Fig. 3 Runoff vulnerability layers for New Damietta: a elevation, b slope, c distance to streams, d land use.And total runoff vulnerability (e), curve number grid (f), runoff hazard per basin (g), and demand for runoff control (h).Source: Prepared by authors using ArcMap 10.5

Fig. 4
Fig. 4 Climate regulation demand criteria: a LST, b NDVI, c NDMI, d climate regulation demand.Source: Prepared by authors using ArcMap 10.5

Fig. 5
Fig. 5 CWs site suitability criteria: a soil salinity, b slope, c accessibility, d soil permeability, and e preliminary and final PSs for CWs.Source: Prepared by authors using ArcMap 10.5

Fig. 6 a
Fig. 6 a Runoff control SBAs, b climate regulation SBAs, and c recreational SBAs (400-1000 m walking distance) for each PS.Source: Prepared by authors using ArcMap 10.5

Fig. 9
Fig. 9 Population and area (%) under the five climate regulation demand classes.Source: Prepared by authors

Table 2
Runoff depth (Q-mm) and volume (Vol-m 3 ) for each sub-basins (SB) for the 2-and 10-year return period rainfall events.Source: Prepared by authors based on data extracted from ArcMap 10.5 analysis Population and area percentage per runoff control demand class.Source: Prepared by authors zones.Where 7.5% of the population lives in very high runoff demand areas, 40.5% in high runoff demand zones, 35.5% live in medium demand zones, and the remaining live in low and very low runoff demand areas (10.2% and 11.9% respectively).