Decadal Urban Land Use/Land Cover Changes and Its Impact on Surface Runoff Potential for the Dhaka City and Surroundings Using Remote Sensing

: Rapid urban growth processes give rise to impervious surfaces and are regarded as the primary cause of urban ﬂooding or waterlogging in urban areas. The high rate of urbanization has caused waterlogging and urban ﬂooding in many parts of Dhaka city. Therefore, the study is undertaken to quantify the changes in land use/land cover (LULC) and urban runoff extent based on the Natural Resources Conservation Service (NRCS) Curve Number (CN) during 1978–2018. The ﬁve-decadal LULC has been analyzed using three-generation Landsat time-series data considering six different classes, namely agriculture, built-up, wetland, open land, green spaces, and water bodies for the years 1978, 1988, 1998, 2007, and 2018. Signiﬁcant changes in LULC for the study area from 1978–2018 are observed as 13.1%, 4.8%, and 7.8% reduction in agricultural land, green spaces, and water bodies, respectively, and a 22.1% increase in the built-up area is estimated. Within Dhaka city, 14.6%, 16.0%, and 12.3% reduction in agricultural land, green spaces, and water bodies, respectively, and a radical increase of 41.9% in built-up area are reckoned. The decadal runoff assessment has been carried out using the NRCS-CN method, considering an extreme rainfall event of 341 mm/day (13 September 2004). The catchment area under very high runoff category is observed as 159.5 km 2 (1978) and 318.3 km 2 (2018), whereas, for Dhaka city, the setting is dynamic as the area under the very high runoff category has increased from 74.24 km 2 (24.44%) to 174.23 km 2 (57.36%) in years 1978 and 2018, respectively, and, mostly, the very high runoff potential areas correspond to the dense built-up surfaces.


Study Area
Dhaka, the capital of the People's Republic of Bangladesh, is the center of all development activities of the country, the home to numerous state and diplomatic institutions. The watershed area considered in this study encompasses Dhaka city and surroundings and lies between 23°36′5″N to 23°55′49″N latitude and 90°8′36″E to 90°34′46″E longitude (Figure 1a). The study area is spread over 846 km 2 , which is 0.57% of the entire geographical area of Bangladesh. The greater Dhaka city is surrounded by the Buriganga river to the southwest; the Turag river to the northwest; the Balu river to the northeast, and the Shitalakshya river to the south. The Dhaka district is bounded by Gazipur, Tangail, Munshiganj, Rajbari, Narayanganj, and Manikganj districts [77]. Tropical vegetation and moist soils characterize the land, which is flat and close to sea level. Topographically, the catchment area is flat with surface elevation ranging from 1 m to 22 m, and, for Dhaka city, it varies between 1 m to 14 m [78]. This leaves Dhaka city susceptible to flooding during the monsoon seasons owing to heavy rainfall and cyclones. Dhaka is one of the most densely populated megacities in the world, with a density of 68,561 people per km 2 [79]. The 400 years of urban growth scenario for Dhaka have been shown by Kabir [80].
The temporal Landsat images at a decadal time interval (1978, 1988, 1998, 2007, and 2018) covering the entire watershed were used along with the Shuttle Radar Topography Mission (SRTM) 30 m digital elevation model (DEM) in the present study. The specifications of the satellite data used The temporal Landsat images at a decadal time interval (1978, 1988, 1998, 2007, and 2018) covering the entire watershed were used along with the Shuttle Radar Topography Mission (SRTM) 30 m digital elevation model (DEM) in the present study. The specifications of the satellite data used are given in Table 1. In addition, a variety of other datasets are used for demographic (population), rainfall analysis and soil map for hydrologic soil group preparation ( Table 2).  Table 2. Specifications of other data used.

Type of Data Source
Demographic data World Population Review Rainfall data Bangladesh Meteorological Department Soil map Soil Resources Development Institute, Bangladesh SRTM 30 m spatial resolution DEM is used for the watershed delineation for the study area using hydro-processing techniques with an area threshold of 3 km 2 . The methodology flow chart (Figure 2) describes the entire workflow sequentially to the quantification of the final runoff potential produced from the catchment.
Remote Sens. 2020, 12, x FOR PEER REVIEW 4 of 23 are given in Table 1. In addition, a variety of other datasets are used for demographic (population), rainfall analysis and soil map for hydrologic soil group preparation (Table 2). SRTM 30 m spatial resolution DEM is used for the watershed delineation for the study area using hydro-processing techniques with an area threshold of 3 km 2 . The methodology flow chart ( Figure 2) describes the entire workflow sequentially to the quantification of the final runoff potential produced from the catchment.

Land Use/Land Cover (LULC)
LULC plays a dominating role in direct runoff estimation and affects infiltration mechanisms, hence considered as the fundamental input for the NRCS-CN method [51,60,63,81,82]. Landsat cloudfree time-series imagery were selected for preparing five decadal LULC layers for the years 1978, 1988, 1998, 2007, and 2018 due to its long-term availability and cost-effectiveness. All the images were obtained from the United States Geological Survey website [83]. The acquisition dates of the collected data vary in-between late January to late February because of seasonality, cloud, and fog-free data

Land Use/Land Cover (LULC)
LULC plays a dominating role in direct runoff estimation and affects infiltration mechanisms, hence considered as the fundamental input for the NRCS-CN method [51,60,63,81,82]. Landsat cloud-free time-series imagery were selected for preparing five decadal LULC layers for the years 1978, 1988, 1998, 2007, and 2018 due to its long-term availability and cost-effectiveness. All the images were obtained from the United States Geological Survey website [83]. The acquisition dates of the collected data vary in-between late January to late February because of seasonality, cloud, and fog-free data and finally to ensure data for maximum closest dates. The supervised classification technique was used for LULC classification using the maximum likelihood classification algorithm. Contextual refinement was carried out manually for each classified map. LULC classes have been categorized into six different classes (Table 3) for the study area as follows: agriculture, built-up, wetland, open land, green spaces, and water bodies [7,8,72,84,85]. The accuracy assessment was performed on a GIS platform. Firstly, a total of 600 random points was generated for all six classes. Afterwards, with the help of training samples, very high resolution Google images (using Google Earth Pro platform), and various secondary sources, the ground truth points were finalized for each year from 1978-2018. Finally, the classification accuracies were assessed as 87

Hydrologic Soil Group
Soil Resources Development Institute has recognized a total of 500 soil series in Bangladesh and further reclassified into 21 general soil types. The different soil textures present in the study area are peaty, loamy fine sand, sandy loam, clay loam, loam, mucky clay, silt, and silty clay loam. However, most of the area are surrounded by clay and loamy soil. The hydrologic soil group (HSG) is an essential input for the NRCS-CN model. According to United States Department of Agriculture (2009) [88] and Musgrave (1955) [86], soil reclassification has been done to produce the HSG layer considering soil characteristics and their potential towards surface runoff (Table 5). It is seen that Dhaka city contains all the four types of HSG classes, i.e., A, B, C, and D, where the soil type 'D' denotes high runoff potential, type 'C' as high to moderate runoff potential soil group, type 'B' as moderate to low runoff potential and type 'A' as low runoff potential ( Figure 3). The study area is dominantly under HSG class 'C', followed by class 'B' and 'D'. It is observed that most of the urbanization has happened along the HSG class 'D' (Figure 3), which also signifies high runoff potential. Cronshey (1986) has suggested a list of curve numbers that can be assigned to different types of land uses with varying HSG classes [87]. Table 5. Hydrological Soil Group (HSG) characteristics [87,89].

HSG Characteristics
A Low runoff potential. Soils having high infiltration rates even when thoroughly wetted and consisting chiefly of deep, well to excessively drained sands or gravels. B Soils having moderate infiltration rates. e.g., shallow loess, sandy loam. C Soils having slow infiltration rates. e.g., clay loams, shallow sandy loam. D High runoff potential. Soils having slow infiltration rates, soils with a permanent high water table, soils with a clay-pan or clay layer at or near the surface, and shallow soils over nearly impervious material.
Remote Sens. 2020, 12, x FOR PEER REVIEW 6 of 23 considering soil characteristics and their potential towards surface runoff (Table 5). It is seen that Dhaka city contains all the four types of HSG classes, i.e., A, B, C, and D, where the soil type 'D' denotes high runoff potential, type 'C' as high to moderate runoff potential soil group, type 'B' as moderate to low runoff potential and type 'A' as low runoff potential ( Figure 3). The study area is dominantly under HSG class 'C', followed by class 'B' and 'D'. It is observed that most of the urbanization has happened along the HSG class 'D' (Figure 3), which also signifies high runoff potential. Cronshey (1986) has suggested a list of curve numbers that can be assigned to different types of land uses with varying HSG classes [87].  [87,89].

HSG Characteristics A
Low runoff potential. Soils having high infiltration rates even when thoroughly wetted and consisting chiefly of deep, well to excessively drained sands or gravels. B Soils having moderate infiltration rates. e.g., shallow loess, sandy loam. C Soils having slow infiltration rates. e.g., clay loams, shallow sandy loam. D High runoff potential. Soils having slow infiltration rates, soils with a permanent high water table, soils with a clay-pan or clay layer at or near the surface, and shallow soils over nearly impervious material. These LULC and HSG maps along with the look up table of CN-II (Table 3) [61,87-89] values were used to derive the CN-II map ( Figure 4). However, the selection of CN-II values is based on the standard table; they can be subjective and usually calibrated when we use them in an event based rainfall-runoff or a daily time scale hydrological model. This study has not performed full scale  table of CN-II (Table 3) [61,[87][88][89] values were used to derive the CN-II map ( Figure 4). However, the selection of CN-II values is based on the standard table; they can be subjective and usually calibrated when we use them in an event based rainfall-runoff or a daily time scale hydrological model. This study has not performed full scale hydrological or flood modeling, but has attempted to estimate changes in runoff potential using CN-II based runoff values. The detailed calibration of surface runoff can be taken as a separate study, where observed river flow data can be used for calibration of CN values.

Antecedent Moisture Condition (AMC)
Antecedent moisture condition (AMC) refers to the moisture content present in the soil at the beginning of the rainfall-runoff event under consideration [89]. It is well known that AMC governs initial absorption and infiltration. The AMC of a particular watershed is taken under consideration, so that the already present moisture content in the soil, which plays a significant role in generating runoff can be estimated. There are three levels of AMC, as recognized by NRCS, as shown in Table 6 [89]. In the present study, the AMC-II condition has been applied to estimate runoff extent.

Antecedent Moisture Condition (AMC)
Antecedent moisture condition (AMC) refers to the moisture content present in the soil at the beginning of the rainfall-runoff event under consideration [89]. It is well known that AMC governs initial absorption and infiltration. The AMC of a particular watershed is taken under consideration, so that the already present moisture content in the soil, which plays a significant role in generating runoff can be estimated. There are three levels of AMC, as recognized by NRCS, as shown in Table 6 [89]. In the present study, the AMC-II condition has been applied to estimate runoff extent.

AMC Type Condition
AMC-I Soils are dry but not to wilting point. Satisfactory cultivation has taken place AMC-II Average conditions AMC-III Sufficient rainfall has occurred within the immediate past five days. Saturated soil conditions prevail.

Overview of NRCS-CN
The NRCS-CN (Natural Resource Conservation Service-Curve Number) is a very popular and widely used method of estimating rainfall excess from rainfall [90] and relies on only one parameter, i.e., curve number (CN) [89]. The NRCS-CN model has been used in this study to estimate the urban runoff potential using Landsat time-series (1978-2018) data. The model presents a simplified procedure for estimating runoff potential in small watersheds, which have been presented by using the CN method (Table 3), where the CN values are assigned to various LULC classes. Thus, it has become very useful for urban areas as the runoff potential is very high even in very small watersheds [63,[91][92][93][94][95]. Therefore, after integrating all parameters, i.e., LULC, CN II map, and rainfall, the runoff is calculated for the entire watershed area as given in Equation (1): Where Q is the runoff, P is rainfall, I a is the initial abstraction, and S = maximum retention capacity as given in Equation (2): Here, the rainfall P was taken from the maximum daily rainfall measured by the Bangladesh Meteorological Department. Two different maximum rainfall events have been used to quantify the surface runoff based on the observed rainfall data starting from year 1953. Therefore, until 1978, 326 mm/day was observed as a maximum daily rainfall event. In addition, for the next two decades, it remained as the maximum rainfall; therefore, 326 mm/day rainfall was used as input also for the years 1988 and 1998. The next decade onwards, 341 mm/day rainfall was observed in the year 2004, which was used as input for the next two years, i.e., 2007 and 2018.

Rainfall Data Analysis
In the present study, daily rainfall data for 65 years (1953-2018) from the Bangladesh Meteorological Department for Dhaka station have been used to analyze the rainfall trend. It is found that maximum daily rainfall of 341 mm/day has occurred on 13 September 2004. The annual maximum rainfall variation in a day (1953-2018) is shown in Figure 5, and it should also be noted that, for the year 1974, the data were not available. The rainfall analysis has been performed based on maximum daily rainfall of the past 65 year's rainfall patterns to know the return period of such high-intensity rain events so that the corresponding flood hazard can be estimated. The climate change impact on flood frequency, magnitude, and runoff variability has been assessed by a number of researchers. Mohammed et al. observed that floods (extreme discharges) and hydrological droughts (low flows) will become more frequent and severe in terms of incidence and extent in the coming days due

Decadal LULC Change Analysis
LULC changes have been investigated for the years 1978, 1988, 1998, 2007, and 2018, using threegeneration Landsat time-series data considering six different classes, namely agriculture, built-up, wetland, open land, green spaces, and water bodies. The classified LULC maps for corresponding years are shown in Figure 6. The change matrix among inter-class and intra-class following the considered time-series has been performed.

Decadal LULC Change Analysis
LULC changes have been investigated for the years 1978, 1988, 1998, 2007, and 2018, using three-generation Landsat time-series data considering six different classes, namely agriculture, built-up, wetland, open land, green spaces, and water bodies. The classified LULC maps for corresponding years are shown in Figure 6. The change matrix among inter-class and intra-class following the considered time-series has been performed.

Decadal LULC Change Analysis
LULC changes have been investigated for the years 1978, 1988, 1998, 2007, and 2018, using threegeneration Landsat time-series data considering six different classes, namely agriculture, built-up, wetland, open land, green spaces, and water bodies. The classified LULC maps for corresponding years are shown in Figure 6. The change matrix among inter-class and intra-class following the considered time-series has been performed.    Figure 7a illustrate the decadal LULC change analysis for the entire study area for the chosen years. In addition, the LULC changes over the decades, i.e., 1978-2018, 1978-1988, 1988-1998, 1998-2007, and 2007-2018 have been performed. It has been found that agricultural land has the largest share among all LULC classes, i.e., 41961.6 ha (49.57%) and 30890.8 ha (36.51%) during the years 1978 and 2018, respectively. It has declined at the rate of 13.1%, 6.6%, 4.3%, and 6.6% in the years 1978-2018, 1978-1988, 1998-2007, and 2007-2018 (1978-1988, 1988-1998, and 1998-2007), the decline is at the rate of 2.1%, 7.7%, and 1.6%, respectively. However, during the decade 2007-2018, it has increased at the rate of 3.6% due to various initiatives undertaken by the Government to revitalize rivers.   In the year 1978, both the classes shared 32.3% and 18.3%, which was the second and third largest share, respectively, of the city area.

Decadal Runoff Assessment
The entire catchment is much larger as compared to the Dhaka city, which contains natural flow from the upstream region and encompasses the city area. Thus, runoff potential analysis has been performed for the entire catchment to identify the vulnerable section of urban areas based on rainfall-runoff analysis and hydrological soil cover variability. The NRCS Curve Number method has been used for estimating the runoff potential and carried out for daily maximum rainfall observed as 341 mm/day on 13 September 2004, for the hydrologic-soil-cover complex scenario of years 2007 and 2018; and 326 mm/day on 14 July 1956, for the years 1978, 1988, and 1998. Figure 8 shows that the decadal storm water runoff-based inundation is relatively high in most urbanized parts of the entire catchment. The runoff potential varies from 164 mm/day to 341 mm/day and classified into five classes as follows: very high (315-341 mm/day), high (275-315 mm/day), moderate (240-275 mm/day), low (210-240 mm/day) and very low (164-210 mm/day). Table 9 is showing the areal extent in square km and the percentage of five distinct runoff classes for the entire catchment. The inundation area under the very low runoff category is found to be varying between 0.2-1% of the total area during 1978-2018, whereas the area under low category ranges between 0.4-4.6% with a sudden increase in area by 40 km 2 and 37.6 km 2 in 1988 and 1998, respectively. Contrarily, the very high runoff potential area has increased from 159.5 km 2 (18.3%) to 318.3 km 2 (36.5%) for the year 1978 and 2018, respectively. It is also observed that the inundated area has mostly occupied the built-up surfaces. The area under the high runoff class has also increased throughout the period except a sudden drop noticed in the year 1998. It is observed that the high runoff potential area (km 2 ) and mean runoff of the study area is increasing with respect to the built-up area (%) (Figure 9). In addition, the mean runoff has increased from 277 to 301 mm/day during the year 1978 to 2018, and the built-up area has increased by 22.1% of the total area during this period.

Decadal Runoff Assessment
The entire catchment is much larger as compared to the Dhaka city, which contains natural flow from the upstream region and encompasses the city area. Thus, runoff potential analysis has been performed for the entire catchment to identify the vulnerable section of urban areas based on rainfallrunoff analysis and hydrological soil cover variability. The NRCS Curve Number method has been used for estimating the runoff potential and carried out for daily maximum rainfall observed as 341 mm/day on 13 September 2004, for the hydrologic-soil-cover complex scenario of years 2007 and 2018; and 326 mm/day on 14 July 1956, for the years 1978, 1988, and 1998. Figure 8 shows that the decadal storm water runoff-based inundation is relatively high in most urbanized parts of the entire catchment. The runoff potential varies from 164 mm/day to 341 mm/day and classified into five classes as follows: very high (315-341 mm/day), high (275-315 mm/day), moderate (240-275 mm/day), low (210-240 mm/day) and very low (164-210 mm/day). Table 9 is showing the areal extent in square km and the percentage of five distinct runoff classes for the entire catchment. The inundation area under the very low runoff category is found to be varying between 0.2-1% of the total area during 1978-2018, whereas the area under low category ranges between 0.4-4.6% with a sudden increase in area by 40 km 2 and 37.6 km 2 in 1988 and 1998, respectively. Contrarily, the very high runoff potential area has increased from 159.5 km 2 (18.3%) to 318.3 km 2 (36.5%) for the year 1978 and 2018, respectively. It is also observed that the inundated area has mostly occupied the built-up surfaces. The area under the high runoff class has also increased throughout the period except a sudden drop noticed in the year 1998. It is observed that the high runoff potential area (km 2 ) and mean runoff of the study area is increasing with respect to the built-up area (%) (Figure 9). In addition, the mean runoff has increased from 277 to 301 mm/day during the year 1978 to 2018, and the built-up area has increased by 22.1% of the total area during this period.     Decadal Runoff Scenario of Dhaka City Table 10 shows that Dhaka city has experienced very high urban runoff over the decades, mostly inundating the built-up areas as shown in Figure 10   Decadal Runoff Scenario of Dhaka City Table 10 shows that Dhaka city has experienced very high urban runoff over the decades, mostly inundating the built-up areas as shown in Figure 10

Validation
The study has been undertaken for the quantification of urban runoff over five decades. Dewan and Yamaguchi [78] have reported that flood inundation of the greater Dhaka for the years 1988, 1998, and 2004 were 47.1%, 53%, and 43%, respectively. It is also revealed that rapid LULC change plays a vital role in escalating the flood processes, which also signifies the findings of the present study. Barua and Ast have also shown the monsoon flood inundation for the year 1998 [98]. Similarly, Sayed and Haruyama 2016 have reported that 70% of the Greater Dhaka area is under moderate to very high flood hazard zone [99]. Several national and international electronic media have also shown the on-street historic flood scenario due to sudden and extreme rainfall events every year ( Figure 11). Furthermore, to understand the efficacy of the storm water drainage network for Dhaka city, a study

Validation
The study has been undertaken for the quantification of urban runoff over five decades. Dewan and Yamaguchi [78] have reported that flood inundation of the greater Dhaka for the years 1988, 1998, and 2004 were 47.1%, 53%, and 43%, respectively. It is also revealed that rapid LULC change plays a vital role in escalating the flood processes, which also signifies the findings of the present study. Barua and Ast have also shown the monsoon flood inundation for the year 1998 [98]. Similarly, Sayed and Haruyama 2016 have reported that 70% of the Greater Dhaka area is under moderate to very high flood hazard zone [99]. Several national and international electronic media have also shown the on-street historic flood scenario due to sudden and extreme rainfall events every year ( Figure 11). Furthermore, to understand the efficacy of the storm water drainage network for Dhaka city, a study using a Storm Water Management Model (SWMM) can be undertaken.

Conclusions
The present study aimed to analyze the spatio-temporal variability of urban runoff characteristics in parts of Dhaka city corresponding to extreme rainfall events in comparison to decadal LULC changes from 1978-2018. It is observed that within catchment area during 1978-2018, 11,070.8 ha of agricultural land has diminished (13.1%); 18684.1 ha of built-up area has increased (22.1%); 2406.3 ha of the wetland has increased (2.8%); 671.9 ha of open land has increased (0.8%); 4076.0 ha of vegetation/ green spaces has diminished (4.8%), and 6628.1 ha of water bodies has diminished (7.8%). On the other hand, the LULC change scenario observed for the Dhaka city over the period 1978-2018 reveals that 4300.5 ha of agricultural land has diminished (14.6%); 12,345.5 ha of built-up area has drastically increased (41.9%); 406.4 ha of the wetland has increased (1.4%); 95.7 ha of open land has decreased (0.3%); 4723.7 ha of vegetation/ green spaces has diminished (16.0%), and 3635.7 ha of water bodies has diminished (12.3%), which is more alarming for the city to tackle current urban storm water runoff.
As it is evident from the current analysis of rainfall, urban growth, and LULC change of the last 40 years, the number of very high rainfall events have increased (one event of daily rainfall >300 mm

Conclusions
The present study aimed to analyze the spatio-temporal variability of urban runoff characteristics in parts of Dhaka city corresponding to extreme rainfall events in comparison to decadal LULC changes from 1978-2018. It is observed that within catchment area during 1978-2018, 11,070.8 ha of agricultural land has diminished (13.1%); 18684.1 ha of built-up area has increased (22.1%); 2406.3 ha of the wetland has increased (2.8%); 671.9 ha of open land has increased (0.8%); 4076.0 ha of vegetation/ green spaces has diminished (4.8%), and 6628.1 ha of water bodies has diminished (7.8%). On the other hand, the LULC change scenario observed for the Dhaka city over the period 1978-2018 reveals that 4300.5 ha of agricultural land has diminished (14.6%); 12,345.5 ha of built-up area has drastically increased (41.9%); 406.4 ha of the wetland has increased (1.4%); 95.7 ha of open land has decreased (0.3%); 4723.7 ha of vegetation/ green spaces has diminished (16.0%), and 3635.7 ha of water bodies has diminished (12.3%), which is more alarming for the city to tackle current urban storm water runoff.
As it is evident from the current analysis of rainfall, urban growth, and LULC change of the last 40 years, the number of very high rainfall events have increased (one event of daily rainfall >300 mm during 1953-2001; two events >300 mm rainfall from 2002-2018) ( Figure 5) as well as exponential urban and population growth. In the future, due to increasing build up land use and the chance of occurring more high rainfall events, the urban flooding is bound to increase.
The runoff assessment using NRCS-CN reveals the significant increase during the above period and mostly in the urbanized areas. The very high and high runoff areas have increased within catchment from 159.5 km 2 (18.3%) and 144.2 km 2 (16.6%) in the year 1978, to 318.3 km 2 (36.5%) and 439.1 km 2 (50.4%) in the year 2018, respectively, which is spread throughout the urban and peri-urban areas; and very low and low runoff area reveal insignificant trend along with the significant contribution of moderate runoff area i.e., 527 km 2 (60.6%), 510.1 km 2 (58.5%) and 564.3 km 2 (64.7%) for the year 1978, 1988, and 1998, respectively. In Dhaka city, the scenario is found to be very alarming as very high runoff has increased from 74.24 km 2 to 174.23 km 2 sharing 24.44% and 57.36% of the city area for the year 1978 and 2018, respectively, which mainly coincide with dense built-up surfaces; on the contrary, low and very low runoff area reveals an insignificant trend sharing a very small area, similar to catchment area; the high runoff area is found to be almost the same in years 1978 and 1988, whereas, from 1998 to 2018, the inundation area has significantly increased from 91.93 km 2 to 126.35 km 2 , sharing 30.27%, 41.59%, respectively; and moderate runoff area exhibited an increasing trend in the year 1978-1998 and an insignificant trend in the year 2007-2018.
Therefore, the study demonstrates the urban storm water runoff scenario in parts of Dhaka city and its surrounding area based on the NRCS-CN model considering LULC change analysis over the decades. One of the limitations of the present work is considering a single urban built-up LULC class, which has resulted in very high CN value. This class can be improved with level 3 or 4 urban LULC classes using high-resolution remote sensing data, giving better control over CN values in this area. Therefore, it is recommended that further studies can be undertaken with improved datasets such as high-resolution remote sensing data and DEM, high temporal rainfall data (15 minutes' interval rainfall data is recommended). Additionally, the storm water modeling software, such as SWMM or MIKE URBAN, and other modeling techniques can be used for comparison.