Analyzing Land Cover Change and Urban Growth Trajectories of the Mega-Urban Region of Dhaka Using Remotely Sensed Data and an Ensemble Classifier

Accurate information on, and human interpretation of, urban land cover using satellite-derived sensor imagery is critical given the intricate nature and niches of socioeconomic, demographic, and environmental factors occurring at multiple temporal and spatial scales. Detailed knowledge of urban land and their changing pattern over time periods associated with ecological risk is, however, required for the best use of critical land and its environmental resources. Interest in this topic has increased recently, driven by a surge in the use of open-source computing software, satellite-derived imagery, and improved classification algorithms. Using the machine learning algorithm Random Forest, combined with multi-date Landsat imagery, we classified eight periods of land cover maps with up-to-date spatial and temporal information of urban land between the period of 1972 and 2015 for the mega-urban region of greater Dhaka in Bangladesh. Random Forest—a non-parametric ensemble classifier—has shown a quantum increase in satellite-derived image classification accuracy due to its outperformance over traditional approaches, e.g., Maximum Likelihood. Employing Random Forest as an image classification approach for this study with independent cross-validation techniques, we obtained high classification accuracy, user and producer accuracy. Our overall classification accuracy ranges were between 85% and 97% with kappa values between 0.81 and 0.94. The area statistics derived from the thematic land cover map show that the built-up area in the 43-year study period expanded quickly, from 35 km2 in 1972 to 378 km2 in 2015, with a net increase rate of approximately 980% and an average annual growth rate of 6%. This growth rate, however, was higher in peripheral areas, with a 2903% increase and an annual expansion rate of 8%, compared to a 460% increase with an annual growth rate of 4% in the core city area (Dhaka City Corporation). This huge urban expansion took place in the north, northwest, and southwest regions of Dhaka, transforming areas that were previously agricultural land, vegetation cover, wetland, and water bodies. The main factors driving the city towards northern corridors include flood-free higher land, the availability of a transportation network, and the agglomeration of manufacturing-based employment centers. The resulting thematic map and spatial information produced from this study therefore serve to facilitate a detailed understanding of urban growth dynamics and land cover change patterns in the mega-urban region of Dhaka, Bangladesh.


Introduction
Although the overwhelming majority of Bangladesh's population lives in rural settings with a low level of urbanization compared to other Asian countries (currently 34%), the megacity of Sustainability 2018, 10, 10; doi:10.3390/su10010010www.mdpi.com/journal/sustainabilityDhaka has experienced massive urbanization in recent years, which sharply contrasts with the rate of urbanization in other nations.The population of Dhaka, for example, rapidly inflated from nearly 1 million people in 1972 to 17 million or more in 2015 [1][2][3], making it one of the most densely inhabited and most quickly urbanized cities in the world [4].Various estimates suggest that if the current trend of population growth continues, together with the ever-increasing economic growth and industrialization that is expected to occur over the next 10-15 years, Dhaka will emerge as the sixth-largest urban cluster on the planet [3].The Dhaka mega-urban center is projected to reach an annual population growth rate of 2.98% by 2030, with the total population in the city territory surpassing 27.37 million by this time [1,3], and exceeding the growth rate of other major urban centers such as Beijing, Shanghai, and Mexico City.This immense and swift urbanization driven by population growth in Dhaka is primarily inflated by a great influx of rural migrants.Rural dwellers immigrate to the urban center either looking for work in a booming industrial sector (i.e., such as readymade garments around the city) [5] or to avail themselves of civic facilities such as health care clinics, education, and so on.Such growth has also resulted from the increasing frequency of natural disasters like floods, erosion, increasing salinization, and concomitant loss of agricultural productivity in the surrounding countryside, which has previously caused large-scale rural-to-urban migration.
The current urbanization pattern with its fast and uncontrolled expansion, has, however, put Dhaka's critical environment, ecology, and land resources at further risk, especially with the presence of highly eco-sensitive wetland areas in close proximity to the city [6].For example, water bodies and wetlands surrounding Dhaka are facing destruction as these are being filled up to construct multi-storied buildings and other real estate developments.Coupled with an unbearable traffic gridlock, high pollution levels, and insufficient infrastructure and housing development [2], this ranks Dhaka as one of the least livable cities on Earth [7].
There have been overwhelming indications that such rapid urbanization with increasing anthropogenic activities driven by population growth has been the major driver of changing land cover and the resultant escalated impact on global climate change [8], with other long-lasting effects on many biophysical processes [9] Such unplanned urbanization and commensurate land cover change is the primary cause of habitat transformation [10,11], losses of biodiversity, and the breakdown of ecosystem services [10].Noticeable impacts on the adjacent environment of urbanized land cover include broad exploitation of productive agricultural land [11,12], wildlife habitat degradation and isolation [9], hydrological cycle disruption [13], and degradation of green cover patches [14].Furthermore, urbanization destroys soil productivity by changing the structure and function of agricultural landscapes [15,16], causing serious threats to fresh food production while simultaneously increasing the risk to long-term food security in the region and beyond [11].
The megacity of Dhaka shelters large proportions of the urban population and is responsible for a major share of the country's economic activities as well as the consumption of food and natural resources.The city is the major engine of Bangladesh's economic growth, which yields about one-third of the country's total GDP [17].Due to its massive population concentrations and rapid development dynamics, Dhaka is gaining increasing importance as an economic, political, and administrative hub for the whole nation.Despite its rapid population increase and significant urban expansion, the volume of scientific research and analysis of the urbanization trends, patterns, and processes of Dhaka in relation to their spatial and temporal context has yet to be produce conclusions.Earlier studies carried out by Dewan et al. (2009) [18] and Ahmed et al. (2015) [19] aimed to explain land cover change in Dhaka; however, their attempts were limited in terms of aerial extension, temporal frequencies, and available methodology.In addition, accurate information on the physical dimensions of urban areas is often confusing and difficult to obtain.
Given the rapid rate of urban expansion with high inter-and intra-annual variability of land use practices in developing countries like Bangladesh, a higher frequency and a longer temporal sequence of urban growth analysis is crucial, such that decadal changes are no longer enough.Rather, a much higher frequency of analysis is preferred.Such consistent monitoring enables concerned authorities to have a clear picture of the current state of various land use types, and the physical direction and the complexity of their development dynamics over time in response to evolving economic, social, and biophysical conditions.Thus, a better understanding of land and environmental resources and their changing patterns is crucial for a sustainable urban future for Dhaka.
As remote sensing represents a major source of current and historical urban information by providing spatially consistent coverage of large areas with fine spatial resolution and temporal frequencies [20], the aim of this study is, therefore, to quantify the full continuum of urban growth in Dhaka using open-source Landsat imagery.Using the ensemble classifier Random Forest with multi-date remote sensing imagery, we studied land cover change and urban expansion of the greater Dhaka for the period of 1972 to 2015.The land cover map and spatial metrics generated from this study may serve as a valuable tool for planners who need to better understand and more accurately characterize urban processes and territorial planning regarding urban growth in Dhaka's region.

Study Area
Dhaka has the unique position of being at once the oldest, historically largest, and most centrally located city in Bangladesh (23.3-24.6 • N and 90.1-90.3• E) (Figure 1).The current study area is a polycentric urbanized area, comprising a number of cities and growth centers that are physically contiguous with and relate to one another through population growth, economic activities, and physically connections (developed road networks and urban expansion).Dhaka and neighboring areas are witnessing a surge of industrial and economic activities and has become an attractive destination for direct foreign investment in Bangladesh.It is a powerful engine of economic growth, providing important economic, social, and cultural opportunities for the entire nation.The total aerial expansion of this mega-urban region as selected for this study is 1528 km 2 , which is called the Detailed Urban Planning Area of Greater Dhaka (DAP).
Historically, the city witnessed a lot of different rulers and kingdoms, which led to its development as a widely-connected urban center.In the 7th to 8th centuries, Dhaka was under the Buddhist kingdom, then controlled by Hindu rulers.From 1299 to 1608 the city was governed by Turks and Pathans, and then dominated by Mughals beginning in 1608 [21].From the beginning of the 17th century Dhaka was governed by a succession of European settlers, first by the Portuguese, then the Dutch, French, and later the English colonial powers [22].
Originally Dhaka City was a small town lying between the rivers Buriganga in the south and Dulai Khal in north, with its center near the present Bangla Bazar [23].Later the city expanded northwards from the banks of the Buriganga River due to the physical constraints imposed by vast swampy areas to the east and west of the city.This barrier still exists but is no longer as effective as in the past since much of this swampy land has already been filled in or is targeted for future urban development and infilling by land developers, co-operative societies, or individuals.As a result, while the city is still expanding northwards, it is now also experiencing growth in all directions and mostly in chaotic and unplanned ways.The 400-year-old city first came within the framework of planning in 1840 with the formation of a municipal committee [24].A more formal organization of urban development was initiated in 1959 with the adoption of a master plan [23].Dhaka was first electrified in 1878 and connected by railway with the southeast-neighboring township of Naryanganj in 1885.
Dhaka City and its surrounding areas comprise six municipalities, Kadamrasul, Gazipur, Naryanganj, Sidirganj, Savar, and Tongi, and a number of connecting small urban areas.They are physically contiguous, and have grown into one of the world's largest megacities, with a very rapid rate of urban expansion.There are five major river systems rolling across the megacity, i.e., the Buriganga-Dhaleshwari to the south, Bansi-Dhaleshwari to the west, Turag Rivers to the North-West, Tongi Khal to the North, and the Shitallakhya-Balu to the east and southeast.These rivers have been the most important means of travel and communication for a long period.Dhaka and its neighboring areas are composed of alluvial terraces of the southern part of the Pleistocene uplands known as the Madhupur tract and low-lying areas around the confluence of the Meghna and Lakhya rivers.Geomorphic classifications reveal that a relatively young floodplain constitutes the largest area, followed by the higher terraces of the Pleistocene period.Low-lying swamps and marshes located in and around the city are examples of other major geomorphic features.Geomorphic classifications reveal that a relatively young floodplain constitutes the largest area, followed by the higher terraces of the Pleistocene period.Low-lying swamps and marshes located in and around the city are examples of other major geomorphic features.
The elevation in the area ranges from 0 to 24 m above mean sea level.Except for the Pleistocene terraces on the northern edge, 17% of the land is 0-5 m above mean sea level, which accounts for 250 km 2 of the total area.In addition, around 21% of the land is 6-7 m (314 km 2 ), and the rest (53%) of the land is 9-24 m above the mean sea level.Dhaka enjoys a fairly equitable tropical monsoon climate, with an annual average temperature of 25 • C (77 • F) and monthly means varying between 18 • C (64 • F) in January and 29 • C (84 • F) in August.Annual rainfall of Dhaka City varies from 1169 to 3028 mm and the annual average rainfall of 2076 mm [25].Nearly 80% of the rainfall occurs during the monsoon season, which lasts from May to the end of September.

Data and Processing
This study uses a series of Landsat imagery, geospatial information, and multilayer census data in addition to fieldwork and survey data.Considering the local climate patterns, especially the effect of the monsoon on image spectral signatures and landscape configuration across the study area, we downloaded eight periods of cloud-free Landsat images captured in the dry season between December and January for the following years: MSS 1972, 1980, TM 1990, 1995, 2005, 2010, ETM+ 2000, and OLI-TIRS 2015.Along with Landsat data, a supplementary night light imagery, OLS (Operational Line-scan System) sensor, 1 km 2 spatial resolution, version 4, Defense Meteorological Satellite Program/Operational Linesca System (DMSP-OLS) Nighttime Lights composite images were downloaded from NOAA (https://www.ngdc.noaa.gov/).Since urban areas or developed regions reveal higher light reflectance during the night time [26], this DMSP-OLS night light signature is believed to be an important driver in our regression-based model calibration.As DMSP-OLS data are only available from 1992 to 2013, spectral signature of these imagery for 1992, 1995, 2000, 2005, 2010, and 2013 were used for model calibration for the years 1990,1995,2000,2005,2010, and 2015, respectively.In addition, a continuous surface elevation map was generated using a Shuttle Radar Topography Mission (SRTM) 30 m spatial resolution image acquired in September 2014.All downloaded images were visually inspected for geometric precision.Recent GIS data including planning boundaries, different levels of administrative boundaries, road networks, water bodies, and other physical features were collected from the Dhaka City Corporation, Urban Development Directorate, and RAJUK.Given the classifier model such as Random Forest employed in this study, only the visible, near infrared and thermal bands from Landsat imageries were stacked, mosaicked, and later subsetted using the study area boundary.Several spectral and spatial image enhancement techniques were applied to improve image signatures.Since remotely sensed data vary in spatial, geometric, radiometric, spectral, and temporal resolutions, all images were geometrically corrected, resampled, and transformed to the local coordinate system: Bangladesh Transverse Mercator Projection (BTM).As MSS images and DMSP-OLS come with 60-m and 1-km resolution, respectively, we artificially resampled MSS and DMSP-OLS imagery to 30-m spatial resolution so that all the images used in this study have identical spatial resolutions, thus enabling our mapping and analysis to continue.All preliminary image processing such as mosaicking, subsetting, spectral signature analysis, image filtering, and enhancement were carried out using ERDAS Imagine image processing software.The main image classification algorithm Random Forest was executed in R-studio.Finally, to reduce potential salt-and-pepper effects, a 3 × 3 majority filter was applied before the classified land covers were used for further analysis [27,28].To demonstrate the spatial configuration of urban expansion, eight transect zones, forming fan-shaped areas, were created with respect to the pre-determined urban center.The fan-shaped transects were drawn in GIS, referenced at the location of the RAJUK building in the center of the old Dhaka.Then using a GIS overlay, the built-up areas that fall within each zone were aggregated and the total area of each category within each sector was calculated.These transect zones and sector analysis methods are effective for characterizing the quantity of both the spatial distribution and orientation of urban expansion.Finally, these values were displayed as spider graphs allowing the spatial distribution of built-up areas to be visualized and illustrated.

Land Cover Classification Scheme
Accurate information on land use and its human interpretation using satellite-derived data is critical given the complex multi-scale interactions between socioeconomic, demographic, and environmental factors occurring at multiple temporal and spatial scales [29,30].These multifaceted interactions with their periodically changing patterns and resultant effects on anthropological and ecological wellbeing are required to be documented and utilized for planning, such as for the best uses of critical land and use of available environmental resources [31].Land cover classification using open-source Landsat imagery has been widely employed to document local to large-scale global landscape changes and their anthropogenic disturbances [31,32].Acquiring accurate land cover information for such studies, however, can still be a challenging task, particularly in urban areas, due to extremely heterogeneous landscapes and the frequent problem of mixed pixels [33,34].For example, sand or bare soil is often confused with bright impervious surfaces, e.g., concrete, due to their similar spectral reflectance [34,35].To deal with such a challenging task as land cover classification across a diverse urban landscape, a wide range of parametric and non-parametric methods for analysis of satellite-derived imagery are in practice and continue to be proposed.Since increasing accuracy in land cover classifications is essentially important for remote sensing analysis, data mining technology and machine learning algorithms such as Support Vector Machine (SVM), Artificial Neural Network (ANN), Decision Tree (DT), and ensembles classifier such as Random Forest (RF) are used in remote sensing studies and have shown a large increase in recent times [34,36].In particular, Random Forest [37] has recently received considerable attention in remote sensing image classification due to its ability to successfully handle high data dimensionality and multicollinearity with computational efficiency [38], often producing higher classification accuracy than traditional methods [39,40].The land cover classification method for this study was therefore selected to be Random Forest, due to it being such a powerful machine learning ensemble and non-parametric statistical learning technique.
The model with associated packages was loaded and executed in R open-source statistical computing and graphics software (https://www.r-project.org/).Since the aim of this study was to examine the land cover change and spatial extent of built-up areas, we adopted a simple classification system, partly derived from Anderson's [41] first-order hierarchical classification system.In this way, six representative land cover categories were generated (Table 1) using expert knowledge of the study area and the extensive field survey undertaken in July 2016.

Random Forest Classification Method
Random forest is a machine-learning algorithm, with classification and regression techniques (developed by Breiman, 2001) based on a combination of a large set of decision trees.The basic premise of this approach is that a set of classifiers performs better than an individual classifier.Simply, the model works by randomly picking a sample of observations and a sample of variables (typically, a few hundred to several thousand depending on the size and nature of the training set) with replacement to produce a large number of individual decision trees [42].These individual trees are then aggregated using bootstrap aggregation sampling techniques called bagging, from which the final class assignment is determined by applying a majority voting rule [40,42].The principal of random forest is to improve variance reduction of bagging by reducing the correlation between the trees without increasing the variance too much [37].This is achieved in the tree-growing process through combinations of many binary decision trees.Using bootstrapping techniques, each tree is trained by selecting a random set of variables coming from the training sample and choosing randomly at each node a subset of explanatory variables [36].This bootstrapping approach allows the model to perform better as it reduces the variance of the model without increasing the bias.The idea of bootstrapping is that the prediction of a single tree is highly sensitive to noise while the average of many trees is not, if they are uncorrelated.Observations in the original training data that are not used to fit the respective tree or happen in the bootstrap sample can be used for an unbiased estimate of the prediction error or the out-of-bag (OOB) error [38].Due to its ability to implicitly deal with missing values, and high-variance data, the algorithm has gained increasing popularity in recent times, with various applications in different fields of remote sensing [43,44].It is a non-parametric approach, typically involved with the identification of multiple classes and the capability of determining variables of importance [45,46].

Model Specification/Parameter Setting
There are two parameters in Random Forest that the user is required to specify; first, the number of input variables randomly chosen at a specific split (mtry) and second, the number of trees in the forest (ntree).For classification, the default value of m is √ p and the minimum node size is one.For regression the default value of m is p/3 and the minimal size of the terminal nodes of the trees was set to the default value of five.The node size allows for specification of the minimum number of observations in a node.
The selection of parameter m influences the final error rate [45].If m is increased, both the correlation between the trees and the classification accuracy of individual trees in the forest are increased [47].While the number of trees (ntree) can be as many as possible (default value 100) due to the RF model being fast and not overfitting, in practice the best values for these parameters will depend on the problem and should be treated as tuning parameters.In this study, the number of regression trees grown (ntree) was adjusted based on the out-of-bag (OOB) estimated error, and the optimum number of random variables to be tested in each tree (mtry) was estimated using the tune RF function in the Random Forest package [48].Initially, we built a RF model based on the complete set of data including all spatial, spectral, and topographic variables with varied parameter settings.Then, using the variable of importance facilities, we ranked them in order according to their importance [49].Using backward elimination processes, we eliminated the least important variables from the bottom of this list.Generally, we looked for a large break between variables in the variable of importance plot to decide how many important variables to choose [42].For each forest model iteration, the number of regression trees grown (ntree) was set at 2000 and the mtry value was optimized based on the lowest OOB error.The lowest OOB error in our final model iteration were achieved when the mtry values were tuned between 5 and 6 on 20 variables out of 23.

Collection of Reference Data and Variable Selection
Owing to its diverse physiography and extremely heterogeneous land cover characteristics resulting from significant anthropogenic influence, it is challenging to accurately classify land cover properties in the study area.As a result, sufficient training data with close investigation of image pixel behaviors for each land cover category can minimize classification errors [50].The primary reference data for this study were collected from field surveys using handheld GPS with predefined training site map locations associated with high-resolution Google Earth imagery during July 2016.Initially, training sites were manually identified on images with minimum areas of 10/10 pixels for each land cover category.As Jensen [51] recommended that the number of training pixels should at least be equal to 10 times the number of variables used in the classification model for a parametric classification approach, we collected approximately 400 samples across the study areas.Since most GPS-based ground training data collected were biased to locations across from/close to the road network, an additional 600 to 700 (varies on images) stratified random samples were also taken as reference data [28].As several studies suggest that nonparametric machine learning classifiers such as Random Forest required a large number of reference data to attain the most favorable outcome [48,52], we combined ground training data and randomly chosen samples in order to get enough information on the land cover, while maintaining the spatial distribution of our reference data throughout the study area.Later, large portions of these samples were used as input variables for the Random Forest classifier and some were retained and used for validation for the thematic map of 2015.Reference data for the outputs for 2010, 2005, and 2000 were manually collected using Google Earth imagery, topographic maps from 2010, and Landsat imagery.However, for the outputs for 1995, 1990, 1980, and 1972, on-screen digitizing techniques were adopted and nearly 1000 training samples were generated across the study area.Using these training samples, we extracted spectral properties of selected bands and indices and later used these as input variables for the RF model calibration and predictions.
Since non-parametric-based regression models performed well once they were calibrated with sufficient variables [45], we developed 23 variables including three topographic variables, i.e., elevation, slope, and aspect, and 17 spectral variables from Landsat imagery.In addition, we computed several spatial variables such as road density, distance to the main road, and distance to the river.The topographic variables were mainly extracted from digital elevation data of SRTM imagery, while spectral variables were extracted from visible bands, NIR, and thermal bands of Landsat imagery.The computed indices using spectral reflectance of these band such as NDVI, NDBI, MSI, SAVI, Tasseled Cap Greenness, Brightness, Wetness, and PCA 1-3.The details of these variables and their computed techniques are given in Table 2.If the number of variables is too large, as is often the case for multi-source studies, a RF can be applied only to those variables that have been identified as the most important and that contribute most to increased accuracy [45,53].As a result, we calibrated our model based on a complete dataset, then used the variable importance results to tune parameter settings such as increasing the number of mtry.Using backward elimination processes, we curtailed the least important variables that were ranked far down the list of the variables of importance graph (Figure 2) [48].In this way, the final RF models were developed based on a reduced number of datasets for every single period of image classification.As a large number of trees is recommended when using RF algorithms to stabilize the Mean Squared Error over many iterations [50,54], in each iteration process, the RF algorithm was run to grow 2000 trees with five to six of the 20 explanatory variables randomly selected at each node as potential variables to base the split on.

Accuracy Assessment
In Random Forest a cross-validation test is estimated internally during the bootstrapping process as 70% of training data are used for the tree growing process and the remaining 30% are used to estimate OOB error [42,53].As a result, separate validation does not get an unbiased estimate of the test error [49].However, self-testing is possible in order to assess model performance.In order to assess the accuracy of the land cover maps generated from the Landsat imagery, an independent sample was extracted using random sampling techniques in GIS.Congalton and Green (2009) [55] pointed out that an adequate number of validation data points for each class category is required to generate a statistical error matrix and accuracy assessment.In total 500 validation points were drawn using stratified random sampling to meet the demand of Congalton and Green's (2009) rule-of-thumb of a minimum of 50 samples per class.Later these samples were cross-checked and validated class by class using land cover maps, Google Earth high-resolution imagery, and true color Landsat imagery.
The overall accuracy, user's accuracy, producer's accuracy, and the non-parametric Kappa coefficient were then derived from the confusion matrices to assess the classification accuracies.The overall classification accuracy was 85-97%, with kappa values ranging between 0.81 and 0.93.The highest classification accuracies were reported for 2010 and 2015 at 97.48 and 94.05%, with overall Kappa statistics of 0.94 and 0.93, respectively (see Table 3).The lowest overall accuracies were measured for 2000, 1995, and 1980 with 85%, 86%, and 86%, respectively.
The user accuracy for some classes, such as developed land, water bodies and fallow land were higher, over 90%, while these were lower for wetland and agricultural classes as they were often misclassified for 1972-1990.These misclassifications, however, were greatly reduced with the later map due to availability of land cover information such as Google Earth imagery and a land use map of the region.The confusion matrix obtained from OOB errors for the eight land cover maps shows variance explained with the model was between 95% and 98%.The most important variables according to the values of Mean Decrease Accuracy were varied but, in general, DMSP-OLS night-light reflectance, elevation, the red band, the NIR band, the thermal band, and tasseled cap greenness were of greatest importance (Figure 2).On the other hand, slope and aspect were the least promising indicators in our regression model and so we eliminated them from the final model iteration.

Land Cover Change and Urban Expansion
The spatial pattern of urban expansion between 1972 and 2015 suggests that Dhaka and its neighboring areas have witnessed notable urban expansion in the peri-urban and suburban areas, accompanied by a significant recent increase in population growth.Figure 3 presents the land use and land cover maps of greater Dhaka for the eight time periods (e.g., 1972, 1980, 1990, 1995, 2000, 2005, 2010, and 2015).Each growth center in the following discussion is presented separately to identify areas with significant land cover transformation and facilitate a more detailed understanding of urban pattern and trajectories across the region.
Through a comparison of the land cover maps generated in this study, along with close investigation of historical topographic maps and the empirical literature, it has been observed that there are three primary growth centers in the study area, specifically Savar in the northwest, Naryangang in the east, and Dhaka in the center, located on the banks of the rivers Dhaleswari, Lakkhya, and Buriganga, respectively [23,56].These three mighty rivers fed the initial growth of these cities as they provided year-round navigable channels for trade, transportation, and commerce for the entire region.The area where the primary township developed was at a comparatively higher elevation than the surroundings (see Figure 1), and was thus safe from erosion and flooding and easily serviced by a water linkage.
Dhaka's urban growth gained momentum only after 1950 when a large influx of Muslim migrants from India [22,57] resulted in rapid population growth in the city territory, from 335,935 in 1951 to 556,712 in 1961, with a net growth rate of 65% [58].Such a huge and unexpected population explosion created enormous pressure on Dhaka's existing housing and urban infrastructure and thus demanded land for further development that is away from the old city but readily available for housing and infrastructure development.As a result, major urban development projects during that time were initiated in the immediate north of the old city such as in Banani, Gulshan, and Baridhara, with 140 hectares, 290 hectares, and 150 hectares, respectively, acquired for high-class residential developments [23].In addition, another 500 acres of agricultural land was acquired in the northwest in Dhanmondi for high-class residential development, following the recommendation of a sub-committee created by the East Pakistan planning division in 1948.Although this residential plan was only partially executed, the impacts of these housing developments have been quite significant in terms of the physical structure and spatial dimension of the city.Since then the city has expanded northward and high-class residential areas are constantly attempting to keep themselves at the northern periphery of the city (Figure 3).
The actual urban development in Dhaka, however, has increased drastically since 1971, when Dhaka became the national capital of the new sovereign state of Bangladesh.Even in the 1960s Dhaka was experiencing a very rapid population growth of over 6% annually, and this growth became much faster after 1971, with the population reaching 1,772,000 in 1974 from only 557,000 in 1961 [58].The thematic map of 1972 representing built-up areas during that time in Dhaka was approximately 2750 ha or 9% of the total DCC areas (Table 4).Its immediate neighboring cities such as Gazipur, Naryanganj, and Savar were small townships with estimated built-up areas of approximately 243 ha, 238 ha, and 30 ha, respectively.Agricultural land (60%), lowland (21%), and vegetation cover (13%) were the dominant land cover types in 1972, and accounted for 96% of the study area.While in the DCC area those land cover classes were still higher contributors, their amounts were lowered to 47% for agricultural, 28% for lowland, and 14% for vegetation cover (Figure 4).but it was still at a startling rate in the following period.For example, between 1980 and 1990, Dhaka's population doubled in size from 3,248,000 in 1980 to 6,619,000 in 1990, with an estimated population growth rate of 104% [61] and built-up growth rate of 46% during the same period.Major built-up area in the DCC area in the early 1970s still occupied the core city, but also expanded beyond the old city limits.In the period between 1972 and 1980, the built-up area in the DCC area clearly headed towards the west and north of the existing growth center (Figures 3 and 5).The total built-up area had by this time increased by 72%, while the population in the DCC areas increased to 3,248,000, with a net growth rate of 120% [59,60].Later population growth slowed down, but it was still at a startling rate in the following period.For example, between 1980 and 1990, Dhaka's population doubled in size from 3,248,000 in 1980 to 6,619,000 in 1990, with an estimated population growth rate of 104% [61] and built-up growth rate of 46% during the same period.
By this time, Dhaka experienced dense development and clear expansion towards the north, northeast, and northwest by engulfing remnants of open and undeveloped land in core city areas (Figure 5).Developed land, as seen in this period (1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990), clearly expanded to Adabar, Mirpur-1, Mirpur-11 and Airport in the north and further north to the banks of the Turag River.A visible amount of developed land, as seen in 1990 maps of Tongi bazar (Gazipur) and the surrounding industrial areas, were the main impetus for urban development (Figure 3).The expansion of the city to the east and west was limited due to the presence of a huge depression between Dhaka and Savar, which resulted in major development having to occur in other directions, such as the rapid expansion along the north-south corridor of Dhaka City.While beyond the city limits to the north and northwest crosses into depression areas, vast lands were rural and agriculture in use, actively working as hinterland for Dhaka and yet to be incorporated within the full continuum of urban development by 1990.Although built-up areas were small patches of the whole study area during the early period of this study, they rapidly increased in size in the following years, due in part to the booming business of export-oriented garment industries in Bangladesh.This type of business expansion had enormous impacts on urban growth in and around Dhaka as it contributed to industrialization and job opportunities in the cities, and spurred the mass movement of the surplus population away from the countryside [23].As elevated land in the city has already been urbanized during the rapid urban development since the 1980s, such industrial expansion coupled with residential development could only expand to the north (Tongi-Gazipur) or by reclaiming the vast depression land in the east and west of the core city.As a result, arable and wetland areas, mainly located adjacent to the existing urban edge, have become the target of urban expansion, particularly since the 1990s.As such, land filling was triggered in the outskirts of the capital, engulfing canals, paddy fields, waterways, and the vast depression areas in the east and west of the city.As observed from land cover maps of 1990, major landfilling areas are identified in several locations such as in Khilkhet (Bashundhara Residential Area) and Mirpur (Duaripara and Rupganj).As a result, landfill area increased its proportional share from 475 hectares in 1980 to 801 hectares in 1990 in the DCC area (Table 4).As population pressure with industrial development continued to increase in the city limits, landfilling battles on the city outskirts concurrently increased at a rapid rate in the following decades.For example, fallow land reached 2144 hectares in 2000 and 4067 hectares in 2015 with a net growth rate of 45% and 127%, respectively (Table 4).The major landfilling areas are still identified in the eastern part of the city, such as Khilkhet in the Bashundahra Residential Area and other new locations at Aftabnagar, Banasree, and Karial (Figure 5).Between 2000 and 2015, the study area observed landfills in all directions including the existing location, and continuous expansion of the area that had been previously filled up.However, this time landfilling extended to the southeast along the Dhaka-Naryanganj highway and to the east in Keraniganj, crossing the river of Buriganga.
The initial urban development in Savar was on the northeast bank of the Dhaleswari River and the small township was located between the present Namabazar and Savar Thana (Figure 6).The whole region to the north, east, and south was largely agricultural or marshy land.As estimated from the thematic map of 1972, approximately 14,580 hectares or 58% was agricultural and 6954 hectares or 27% of the total land area in Savar was lowland (Table 5).An analysis of the land cover map of the Savar area reveals that the built-up area had increased between 1972 and 1980 by approximately 211%, while the population growth was around 92% (estimated from the 1974 to 1981 census) [59].The trend of urban growth in Savar continued in the following years, but the rate was substantially higher than in the past.For example, in the period between 1990 and 1995 the total urban growth in Saver was approximately 325%, with an annual increase rate of 34%.Later developed land doubled in size as it increased from 628 hectares in 1995 to 1274 hectares in 2000, with a net increase rate of 103%.Although developed land continued to increase, the rate of growth was higher in the period of 2005-2010; the total built-up area reached 4660 hectares in 2010, with a net growth rate of approximately 150%.Population growth in the same period was highest at around 130% between the 2001 and 2011 census periods [62,63].
In 1974 Savar was a city with 131,429 inhabitants [58], but urban development pressure in Savar resulted from the establishment of a number of planned large-scale industrial zones, such as the Dhaka Export Processing Zone in Dhamsona in the early 1990s and heavy industrial plants in Nayarhat in the mid-1970s.As a result, the population swelled, due mainly to a great influx of rural migrants looking for work in the booming industrial sectors in Savar [23,57].In addition, a number of establishments with growing tourism importance such as Jahanginrangar University, the Savar Cantonment, the Atomic Energy Center, and the National Memorial, have further reinforced development of Savar.
The southeastern city of Naryanganj was mainly known in the past as a salt trading center and  The trend of urban growth in Savar continued in the following years, but the rate was substantially higher than in the past.For example, in the period between 1990 and 1995 the total urban growth in Saver was approximately 325%, with an annual increase rate of 34%.Later developed land doubled in size as it increased from 628 hectares in 1995 to 1274 hectares in 2000, with a net increase rate of 103%.Although developed land continued to increase, the rate of growth was higher in the period of 2005-2010; the total built-up area reached 4660 hectares in 2010, with a net growth rate of approximately 150%.Population growth in the same period was highest at around 130% between the 2001 and 2011 census periods [62,63].
In 1974 Savar was a city with 131,429 inhabitants [58], but urban development pressure in Savar resulted from the establishment of a number of planned large-scale industrial zones, such as the Dhaka Export Processing Zone in Dhamsona in the early 1990s and heavy industrial plants in Nayarhat in the mid-1970s.As a result, the population swelled, due mainly to a great influx of rural migrants looking for work in the booming industrial sectors in Savar [23,57].In addition, a number of establishments with growing tourism importance such as Jahanginrangar University, the Savar Cantonment, the Atomic Energy Center, and the National Memorial, have further reinforced development of Savar.
The southeastern city of Naryanganj was mainly known in the past as a salt trading center and for having the largest jute mill, high numbers of cotton mills, and other large-scale industrial plants.In its initial period, urban development was limited to the eastern and western banks of the Lakkhya river, as seen in the land cover map of 1972 (Figure 7): approximately 238 hectares of land were covered with built-up area (Table 6).The outer peripheral expansion of this city in the south was up to the Kashipur River and towards the north up to the Adamjee EPZ.Later, expansion of huge cotton industries in the early 1990sn along with other large-scale industrial development, meant the city underwent rapid urbanization to accommodate the growing number of workers and associated infrastructure development.As a result, between 1980 and 1990 the built-up area increased to 83%.This time the city is expanding towards the northwest, following the old axis of the Dhaka-Naryanganj highway, but, due to physical constraints such as the river Buriganga to the south and the flood-prone swampy land to the north, urban development away from the road was slow.Nevertheless, due to the aftermath of the 1988 floods, this swampy land in the northwest outskirts of Naryanganj city and southeast periphery of Dhaka City, literally called Dhaka Naryanganj Demra (DND) flood protection dam, originally meant to protect the land from flooding and improve agricultural production, became the main target of residential development [64].As a result, a large area in the DND area experienced rapid and chaotic urban development, which resulted in rapid land conversion from agricultural-lowland to residential development (Figure 7).Built-up area, therefore, increased from 612 hectares in 1990 to 1720 hectares in 2000, a net growth rate of 181%.The developed land in Narayanganj is expanding much faster in recent times due to the establishment of several planned industrial estates in Enayetnagar, Fatullaha, Adamje, Kachpur, and Panchabati, which has made the city one of the largest urban agglomerations in the region.
Beyond the Buriganga River, Keraniganj was primarily rural and dominated by various agricultural types.The built-up area estimated from the land cover map of 1972 (Figure 3) was  Nevertheless, due to the aftermath of the 1988 floods, this swampy land in the northwest outskirts of Naryanganj city and southeast periphery of Dhaka City, literally called Dhaka Naryanganj Demra (DND) flood protection dam, originally meant to protect the land from flooding and improve agricultural production, became the main target of residential development [64].As a result, a large area in the DND area experienced rapid and chaotic urban development, which resulted in rapid land conversion from agricultural-lowland to residential development (Figure 7).Built-up area, therefore, increased from 612 hectares in 1990 to 1720 hectares in 2000, a net growth rate of 181%.The developed land in Narayanganj is expanding much faster in recent times due to the establishment of several planned industrial estates in Enayetnagar, Fatullaha, Adamje, Kachpur, and Panchabati, which has made the city one of the largest urban agglomerations in the region.
Beyond the Buriganga River, Keraniganj was primarily rural and dominated by various agricultural types.The built-up area estimated from the land cover map of 1972 (Figure 3) was approximately 138 hectares, while agricultural was around 12,305 hectares or 70% of the total land in Keraniganj.However, these areas have undergone rapid urbanization in recent times due to several connecting bridges over Buriganga, which established a quick gateway to Dhaka and the rest of the south Keraniganj.The major urban development in Keraniganj is now found in the immediate south bank of Buriganga, in such locations as Zinzira, Ragunathpur, Mandail, Kalindi, Barisur, Dakpara, Kaliganj, and up towards Mirerbag.In addition, the area also houses a number of small-and medium-scale manufacturing industries, warehouses, and lower-income residences.These manufacturing industries, along with ongoing massive-scale residential development projects in Ekuria, are making this area urbanize rapidly in recent times.As shown in the land cover map for 1990 to 2015 (Figure 3), the total built-up area in Keraniganj increased by 764%, while fallow land increased by 1052% (Table 7).The area statistics derived from the land cover maps reveal that urban growth in the northern outskirts of Dhaka city in Gazipur (Figure 8) increased significantly, from 243 hectares in 1972 to 7640 hectares in 2015, with a total growth rate of 3038% (Table 7).This built-up growth rate is substantially higher than any other urban center studied here.Among the eight time periods, the highest urban growth in the region was recorded in the 1972-1980 and 2005-2010 periods at the rate of 145% and 100%, respectively.This growth can be attributed to both industrial development and expansion and human population growth in the region.Specifically, this area contains five universities, a high-tech park, a number of government, semi-government, and private offices, and the largest concentration of garment manufacturing industries in the countree.Due to the availability of flood-free, low-cost land with close proximity to the core city area, Gazipur has become the main target of both residential and industrial development in recent times.As a result, arable agricultural land and lowland cover types have been converted to developed land very quickly.As observed in the 43-year period, the total agricultural land dropped to 6784 hectares by 2015 from 25,186 hectares in 1972, a 73% decrease.The land cover statistics also suggest that vegetation cover in the region has increased over the 43-year study period; however, except for the forest cover in the north, this increase in vegetation is actually an increase in peri-urban settlement.These settlements are covered with homestead vegetation, as is typical across Bangladesh.
types have been converted to developed land very quickly.As observed in the 43-year period, the total agricultural land dropped to 6784 hectares by 2015 from 25,186 hectares in 1972, a 73% decrease.The land cover statistics also suggest that vegetation cover in the region has increased over the 43year study period; however, except for the forest cover in the north, this increase in vegetation is actually an increase in peri-urban settlement.These settlements are covered with homestead vegetation, as is typical across Bangladesh.In the smallest administrative unit, Bandar, which is located on the eastern bank of the Shitallakhya River, the dominant land cover change is an increase in developed land (Figure 3).While the dominant land use is agriculture (41%), the land cover maps for the last 43 years show that developed land increased from 98 hectares in 1972 to 1532 hectares in 2015, with a net growth rate of 910%.The greatest increase is observed between 1995 and 2000, when the total built-up growth was 137% (Table 7).Within this conversion to developed land, the largest portion is industrial and portrelated infrastructure development, which has expanded along the eastern bank of the Shitallakhya River in a north-south direction.This built-up expansion was mainly at the expense of agricultural In the smallest administrative unit, Bandar, which is located on the eastern bank of the Shitallakhya River, the dominant land cover change is an increase in developed land (Figure 3).While the dominant land use is agriculture (41%), the land cover maps for the last 43 years show that developed land increased from 98 hectares in 1972 to 1532 hectares in 2015, with a net growth rate of 910%.The greatest increase is observed between 1995 and 2000, when the total built-up growth was 137% (Table 7).Within this conversion to developed land, the largest portion is industrial and port-related infrastructure development, which has expanded along the eastern bank of the Shitallakhya River in a north-south direction.This built-up expansion was mainly at the expense of agricultural land as the total area of agricultural land declined at a rate of approximately 34%, from 3549 hectares to 2362 hectares, across the period.The landfill area during the same period also increased substantially, from only 2 hectares in 1972 to 175 hectares in 2015.Much of these landfills are comprised of brickfields, which are expanding east and south out from the center.
Among the other three connecting growth centers (i.e., Soanrgaon, Rupganj, and Kaliganj) in this study, increasing cover of developed land was documented in Soanargaon and Rupganj, both of which are located in the southeast and east of the study area (Figures 1 and 3).The development of Sonargaon was historically significant as it was the capital of Bangladesh in the 15th century.The urban development of this historic city gained momentum in recent times due to the significant increase of tourism-related services, educational institutions, and various private residential development projects.As such, developed land in Sonargaon has increased rapidly between 2000 and 2015, from 185 hectares to 851 hectares, with a net growth rate of 360%.Similarly, in Rupganj, urban area markedly increased from 197 hectares in 2000 to 1287 hectares in 2015, a net growth rate of 553% (Table 7).Additionally, the landfill area in Rupganj has increased since 1990, with 738 hectares in 2000 from only 30 hectares in 1990.This growth remained constant in the following year but at an accelerated rate, i.e., it reached 1840 hectares in 2010 and 3082 hectares in 2015 with total growth at a rate of 149% and 67%, respectively.Kaliganj in the northeast of this study area is a predominantly agricultural, rural, and peri-urban settlement (53% in total).The land use survey of 2006 estimates that 23% of the total land of Tongi, Gazipur, Kaliganj, and Rupganj was residential, including urban housing, rural homesteads, and informal housing; less than 1% is designated for industrial and commercial uses and 60% land is used for agricultural purposes [65].The recent land use survey conducted by RDP, RAJUK 2013 [66] estimated that residential land use comprises 36.5% of the current study area, while industrial and commercial uses comprise 2.6% and 1%, respectively.The total built-up area including residential, industrial, public facilities, road networks, commercial, and mixed use comprises 44% of the survey area, which closely aligns with our findings of 45% including urban and rural settlements (except for the small patches of forest in the north, vegetation classes in our study correspond with rural settlements, which account for almost 20%).

Urban Growth Trajectory Analysis
The spatial pattern of urban expansion between 1972 and 2015 suggests that Dhaka and its surrounding areas have experienced rapid urbanization, which has resulted in significant changes in land use over the 43-year period studied here.The spatial trajectory of urban expansion over the study period suggests that north, northwest, and southeast were the major directions experiencing urban development (Figure 9).By contrast, southern and western expansion are relatively low since the southern and western zones are dominated by low-lying marshy land, which clearly constrained the growth of built-up areas in these directions.The urban growth was estimated from the thematic map over the 43-year period in the direction toward the north (35% expansion), along the Dhaka-Mymensingh highway (Tongi-Gazipur) and northwest (30% expansion), the Savar-Ashulia corridor along the Dhaka-Aricha highway, and also toward the southwest (25.6%) along the Dhaka-Narayanganj and Dhaka-Chittagong highways (Figure 9).The development towards the north was significant due to comparatively flood-free land in that direction and particularly a number of development projects, such as shifting the airport from its old site at Tejgaon to its new site north of Kurmitola, which is about 10 km north of the old airport.Another major development project is the Uttara Model Town located immediately north of the new airport.North of Uttara stands the industrial town of Tongi and beyond Tongi and Joydebpur has also been partly filled up with the establishment of two institutions of higher learning, the National University, and the Islamic University of Technology.Due to population growth and the subsequent expansion of trade and commerce, the current development of the city is radiating southwest, north, northeast, west, and northwest from the banks of the Burignaga River and extending into rural areas and beyond the outskirts.Dhaka's recent growth has been linked with Narayanganj and Gazipur City, such that there are no gaps between Dhaka and those two connecting cities.
The spatial information derived from the time-series analysis land cover maps over the 43-year period shows that a continuous increase in built-up land occurred, at a rate much faster than the urban population growth in all connecting cities in this study.For example, between 1972 and 1980 and 1980 and 1990, the annual population growth rate was higher at 11.4% and 7.3%, respectively, while the urban expansion rate was only 7% and 3%, respectively (Figure 10B).Later population growth continued to fluctuate around 3.5% per annum, but the urban expansion rate maintained Another major development project is the Uttara Model Town located immediately north of the new airport.North of Uttara stands the industrial town of Tongi and beyond Tongi and Joydebpur has also been partly filled up with the establishment of two institutions of higher learning, the National University, and the Islamic University of Technology.Due to population growth and the subsequent expansion of trade and commerce, the current development of the city is radiating southwest, north, northeast, west, and northwest from the banks of the Burignaga River and extending into rural areas and beyond the outskirts.Dhaka's recent growth has been linked with Narayanganj and Gazipur City, such that there are no gaps between Dhaka and those two connecting cities.The spatial information derived from the time-series analysis land cover maps over the 43-year period shows that a continuous increase in built-up land occurred, at a rate much faster than the urban population growth in all connecting cities in this study.For example, between 1972 and 1980 and 1980 and 1990, the annual population growth rate was higher at 11.4% and 7.3%, respectively, while the urban expansion rate was only 7% and 3%, respectively (Figure 10B).Later population growth continued to fluctuate around 3.5% per annum, but the urban expansion rate maintained stable growth trends at around 6% until the end of the study period.Another important feature of this urban growth, as presented in Figure 10, is that peripheral growth was higher throughout the study period than the core city area (Figure 10A).In the 43-year period the total peripheral growth rate was 2903%, while that of the core city area was 459%.
In this study we utilized eight dates of imagery to develop our time series with different sensors and spatial resolutions.Since these images were not captured at the same time, the water level in permanent reservoirs such as ponds, rivers, canals, and seasonal wetlands was variable.As can be seen, actual quantities of water bodies and wetlands do fluctuate in different periods; however, the general trend is that wetlands in peri-urban areas were encroached on heavily by agricultural pressure, brickfield expansion, or landfilling.The areas where these processes were taking place are identified in the land cover maps.A possible explanation for this rampant urban development and encroachment in the city outskirts is the lack of space in the core area for further development, combined with the somewhat inadequate public policies concerning these wetlands and their environmental protection.The Greater Dhaka Area is a bustling megacity with more than 18 million people, making it one of the most populated cities in the world, with a density of 23,234 inhabitants per square kilometer within the core city area (DCC).Like other megacities in the developing world, such huge population growth in association with inadequate urban services and lack of authorities' reign over growth management has led to an increase in environmental degradation, exponentially growing slums, and increases in pollution, poverty, disease, crime, land grabbing, corruption, and traffic congestion.Fifty-five government agencies work to provide civic services in Dhaka, all typically following their own trajectories in decision-making, which often leads to chaos among the differing authorities, rarely reflecting citizens' wishes in the implementation of development plans.In terms of driving urban expansion, existing planning guides played little role compared to individual and corporate agents such as real-estate developers, industrial entrepreneurs, investors, and other vested groups.Real-estate developers and manufacturing investors, such as export-oriented garment manufacturers in the city outskirts, have been instrumental in shaping urban development in Dhaka.Hence, this study's findings highlight that a rapidly growing city like Dhaka should be managed through the proactive protection of natural and ecologically sensitive areas, along with the involvement of landowners, locals, land developers, and multilevel stakeholders based on a shared long-term vision concerning the most desirable and sustainable urban growth for Bangladesh.The Greater Dhaka Area is a bustling megacity with more than 18 million people, making it one of the most populated cities in the world, with a density of 23,234 inhabitants per square kilometer within the core city area (DCC).Like other megacities in the developing world, such huge population growth in association with inadequate urban services and lack of authorities' reign over growth management has led to an increase in environmental degradation, exponentially growing slums, and increases in pollution, poverty, disease, crime, land grabbing, corruption, and traffic congestion.Fifty-five government agencies work to provide civic services in Dhaka, all typically following their own trajectories in decision-making, which often leads to chaos among the differing authorities, rarely reflecting citizens' wishes in the implementation of development plans.In terms of driving urban expansion, existing planning guides played little role compared to individual and corporate agents such as real-estate developers, industrial entrepreneurs, investors, and other vested groups.Real-estate developers and manufacturing investors, such as export-oriented garment manufacturers in the city outskirts, have been instrumental in shaping urban development in Dhaka.Hence, this study's findings highlight that a rapidly growing city like Dhaka should be managed through the proactive protection of natural and ecologically sensitive areas, along with the involvement of landowners, locals, land developers, and multilevel stakeholders based on a shared long-term vision concerning the most desirable and sustainable urban growth for Bangladesh.

Conclusions
Using multi-sensor time series satellite data and an ensemble classifier "Random Forest" with socioeconomic, demographic, and biophysical information, this study quantifies spatiotemporal land cover change and captures the full continuum of urban expansion in the megacity of Dhaka, Bangladesh, over a 43-year period.Employing Random Forest, the new classification algorithm in this study has shown high classification accuracy (e.g., overall accuracies of 85-97%) with the ability to discriminate minute details of land cover classification in the complex and heterogeneous landscape represented by such complex, urban land cover.Among the six land cover types, the agricultural class was the dominant land cover type throughout the study period.Its total volume, however, has declined significantly from 921 km 2 in 1972 to 414 km 2 in 2015, with a net decline of 55% over the 43-year study period.This rate was higher in the DCC with 142 km 2 , which accounted for 47% of the DCC area in 1972, and dropped to 33 km 2 , covering only 10% of DCC areas in 2015-a net decline of 77%.Compared to Dhaka City, neighboring cities like Gazipur in the north and Savar in the west experienced higher urban development during that period.The total urban development between 1972 and 1980 in Gazipur was 145% and in Savar was 211%, with annual growth rates of about 12% and 15%, respectively.This rapid urban growth with indiscriminate land reclamation processes has resulted in the degradation of agricultural land, forests, vegetated areas, wetlands, and water bodies across the whole region.The built-up area and fallow land in the same period in DCC increased by 459% and 3561%, respectively, while in the study region the net growth rates of built-up land and fallow land were substantially higher at 1019% and 11,270%, respectively.Major expansion occurred to the north, northwest, and southwest, consuming agricultural land, vegetation cover, and bodies of water.This increase in urban development resulted in a loss of agricultural land, vegetation cover, and bodies of water.In addition, the findings highlight the increasing vulnerability of ecologically sensitive wetlands in urban and peri-urban areas, which have largely been degraded and encroached on by the expansion of brickfield and landfilling.This study has highlighted that time series analysis of remotely sensed data and non-parametric, regression-based machine learning algorithms can be utilized to accurately classify complex, heterogeneous urban land cover with quite high classification accuracy.The methods and techniques employed in this study may serve as guidelines for future applications across varied landscapes.Furthermore, the resulting land cover maps and spatial information generated from this study may assist planners, stakeholders, and academia with gaining a detailed understanding of land cover change, urban growth patterns, and associated ecological risk in Dhaka's mega-urban region of Bangladesh, thus helping to make better planning decisions for a more sustainable urban future.

Sustainability 2018, 10 ,
0010 10.3390/su10010010 4 of 25 the most important means of travel and communication for a long period.Dhaka and its neighboring areas are composed of alluvial terraces of the southern part of the Pleistocene uplands known as the Madhupur tract and low-lying areas around the confluence of the Meghna and Lakhya rivers.

Figure 1 .
Figure 1.Location of the study area and elevation.The boundary highlights the physical delimitation of each growth center.

Figure 1 .
Figure 1.Location of the study area and elevation.The boundary highlights the physical delimitation of each growth center.

Figure 4 .
Figure 4.The proportional distribution (%) of each land cover class in eight periods.(A) Core city area/Dhaka City Corporation area and (B) whole study area.Figure 4. The proportional distribution (%) of each land cover class in eight periods.(A) Core city area/Dhaka City Corporation area and (B) whole study area.

Figure 4 .
Figure 4.The proportional distribution (%) of each land cover class in eight periods.(A) Core city area/Dhaka City Corporation area and (B) whole study area.Figure 4. The proportional distribution (%) of each land cover class in eight periods.(A) Core city area/Dhaka City Corporation area and (B) whole study area.

Figure 4 .
Figure 4.The proportional distribution (%) of each land cover class in eight periods.(A) Core city area/Dhaka City Corporation area and (B) whole study area.

Sustainability 2018, 10 ,
0010 10.3390/su10010010 17 of 25 infrastructure development.As a result, between 1980 and 1990 the built-up area increased to 83%.This time the city is expanding towards the northwest, following the old axis of the Dhaka-Naryanganj highway, but, due to physical constraints such as the river Buriganga to the south and the flood-prone swampy land to the north, urban development away from the road was slow.

25 Figure 10 .
Figure 10.The bar chart on the left (A) compares total built-up area between the core city (DCC) and the peripheral area and the line graph on the right (B) compares annual urban expansion rate (%) with annual urban population growth rate (%) from 1972 to 2015.

Figure 10 .
Figure 10.The bar chart on the left (A) compares total built-up area between the core city (DCC) and the peripheral area and the line graph on the right (B) compares annual urban expansion rate (%) with annual urban population growth rate (%) from 1972 to 2015.

Table 1 .
Descriptions of the major land cover categories.

Table 2 .
List of predictor variables and their descriptions.Topographic variables were derived from the SRTM 30 DEM.Spectral variables and spectral indices were computed using spectral bands of Landsat imagery of TM, ETM+, and OLI-TIRS.

Table 3 .
Accuracy assessment using independent reference points.UA: User Accuracy, PA: Producer Accuracy, OA: Overall Accuracy.

Table 4 .
Area statistics of each land cover category, total and annual increase/decrease rate in core city area (DCC) (Unit: Hectares and %).

Table 5 .
Area statistics of each land cover class, total and annual increase/decrease rate in Savar (Unit: hectares and %).

Table 5 .
Area statistics of each land cover class, total and annual increase/decrease rate in Savar (Unit: hectares and %).

Table 6 .
Area statistics of each land cover category, total and annual increase/decrease rate in Narayanganj (Unit: hectares and %).

Table 6 .
Area statistics of each land cover category, total and annual increase/decrease rate in Narayanganj (Unit: hectares and %).

Table 7 .
Total built-up area (in hectares) in nine connecting cities from 1972 to 2015.The highest urban development in 43 years was documented in the core city area (DCC) and the lowest was in Kaliganj.