Analysis of Settlement Expansion and Urban Growth Modelling Using Geoinformation for Assessing Potential Impacts of Urbanization on Climate in Abuja City, Nigeria

This study analyzed the spatiotemporal pattern of settlement expansion in Abuja, Nigeria, one of West Africa's fastest developing cities, using geoinformation and ancillary datasets. Three epochs of Land-use Land-cover (LULC) maps for 1986, 2001 and 2014 were derived from Landsat images using support vector machines (SVM). Accuracy assessment (AA) of the LULC maps based on the pixel count resulted in overall accuracy of 82%, 92% and 92%, while the AA derived from the error adjusted area (EAA) method stood at 69%, 91% and 91% for 1986, 2001 and 2014, respectively. Two major techniques for detecting changes in the LULC epochs involved the use of binary maps as well as a post-classification comparison approach. Quantitative spatiotemporal analysis was conducted to detect LULC changes with specific focus on the settlement development pattern of Abuja, the federal capital city (FCC) of Nigeria. Logical transitions to the urban category were modelled for predicting future scenarios for the year 2050 using the embedded land change modeler (LCM) in the IDRISI package. Based on the EAA, the result showed that urban areas increased by more than 11% between 1986 and 2001. In contrast, this value rose to 17% between 2001 and 2014. The LCM model projected LULC changes that showed a growing trend in settlement expansion, which might take over allotted spaces for green areas and agricultural land if stringent development policies and enforcement measures are not implemented. In conclusion, integrating geospatial technologies with ancillary datasets offered improved understanding of how urbanization processes such as increased imperviousness of such a magnitude could influence the urban microclimate through the alteration of natural land surface temperature. Urban expansion could also lead to increased surface runoff as well as changes in drainage geography leading to urban floods.


Introduction
Some of the most dynamic places on planet Earth are urbanized locations developing across multiple dimensions.Globally, the expansion of urban forms is documented to be a principal front in habitat destruction which begins with habitat loss and later results in species extinction [1].Notwithstanding their importance, their growth is associated with huge impacts on contiguous ecosystems [2].For instance, the loss of arable agricultural land to urbanization, especially in developing nations, is flagged to be a result of prevalent anthropogenic activities [3,4].Similarly, the unprecedented transformation of natural landscapes into urban settings significantly affects the natural functioning of ecosystems [5].Hence, urbanization has been the foremost human led land-use anthropogenic activity with huge and irreversible impacts.It is a major force that drives changes such as land-use land-cover change (LULCC), biodiversity loss, the biogeochemical cycle, hydrological systems and climate [6].Another prominent agent that can be linked to the unprecedented growth witnessed in urban expansion is population increase [7].One hundred years ago, of every 10 persons, two resided in urban areas.By 2030, the number of people living in urban areas is likely to hit six, and by 2050 seven out of every 10 [8].Since the 1950s, the number of global urban inhabitants has increased, and by 2050 a two fold increase is anticipated from an approximate value of 3.4 thousands of million as of 2009 to 6.4 thousands of million by 2050 [8].The year 2020 is projected to be when the majority of mega cities in the world will be in developing countries due to differential population growth and anthropogenic activities such as change in LULC [9,10].Therefore, governments in West African countries and Nigeria in particular must act fast to better understand spatial and urban growth patterns for improved municipal planning.
The pathway to understanding the process of urbanization can be established by deliberate monitoring of biophysical and socioeconomic conditions of the existing and transformed urban areas [11].However, reliable information on the biophysical dimensions of urban landscapes, especially the urban LULC of the built environment, can be difficult to obtain.Hence, remotely sensed data and its applications can provide critical information about urbanization in order to advance urban science for improved policy and decision making.However, having access to requisite biophysical, socio-economic characteristics of the urban areas is difficult because gaps exist between these data streams.Moreover, linkages between data collection time lags, administrative and landscape units and spatio-temporal scales pose great challenges.
A paucity of reliable information has coincided with a period when substantial growth in urban areas has been witnessed worldwide, and hopes are high among policy makers and research groups that this information gap can be filled.These groups hope that fine-scale information will be increasingly made available to comprehend the impacts of urbanization on local environments and human security such as temperature variability linked to urban heat and cold islands and local environmental and climate change/variability.This has paved the way lately for more attention being directed towards urbanization science with new directions in data collection and analysis [12], as well as monitoring changes in urban LULC considering its strong influence on ecosystems [13].Remarkably, remote sensing scientists are responding to the call for linked environmental and socio-economic information through the advancement of remotely sensed data application and methods [14][15][16][17].
The support vector machine (SVM) method is one of the latest additions to the existing catalogue of superior and robust image classification techniques for handling multispectral satellite images that support LULC analysis, considering their non-linearity and multidimensionality [18][19][20].The SVM-based approach is a non-parametric machine learning algorithm that uses hyperplanes to separate features of different categories with a maximum distance margin located close to it [21].The best generalization is achieved when the margin distance is farthest from vectors from both classes despite minimal training samples which previously limited numerous remote sensing applications [19].
One of the standard ways of simplifying reality is the use of models and, in spatial analysis, LULCC mapping, monitoring and modelling can be regarded as models used for decision support to assess the root causes and implications.Modelling LULCC leads to improved understanding of human-environmental systems' interaction toward sustainable and spatially framed land-use policy development and planning [22,23].Hence, such spatial models are suitable for assessing arrays of complex biophysical and socioeconomic drivers of spatiotemporal patterns of LULCC and measuring the change consequences [24].Lately, different models such as the cellular automata-based, agent-based, machine learning, and spatially explicit approaches have been applied to study changes in LULC [25][26][27].Similarly, the land change modeler (LCM) has been used to model urban sprawl and growth [28,29].Hua et al [30] ran the SELUTH model which derives its name from input data requirements to run the model slope, land use, exclusion, urban extent, transportation and hill shade over Jimei in Fujian Province, China.Other approaches such as artificial neural networks have also been applied for urban modeling studies [31][32][33].The LCM approach is embedded in the IDRISI package which is an integrated environment suitable for analysis, prediction and validation of LULCC [34].LCM uses categorical gridded images with the same land-cover types sequentially similar in order of arrangement as input for modelling LULCC [35].Land-cover changes are evaluated across multiple time lines, and change results are calculated and presented in the form of graphs and maps.The subsequent step is predicting the future by generating LULC maps based on the transition potential maps [35], trusting the Multi-Layer Perception (MLP) neural networks output [36].For short time frames, the LCM performed better with good prediction accuracies particularly with stable land-covers types against rapid conversion [35].Also, comparing LCM outputs to other LULCC models that predict change based on supervised classifiers such as the weights of evidence (WoE) approach that uses user defined weighting, more accurate change maps are generated.This is because the final change map uses the overall change potential maps which are based on neural network outputs that are capable of expressing changes in various land-cover types much better than single probabilities derived from the WoE approach [36].
The federal capital territory (FCT) of Nigeria was established in 1976 but physical development only began in 1980 [37], and it has been characterized as one of the fastest growing cities in West Africa [38].The territory has experienced rapid LULC changes, urban spatial expansion and transportation infrastructure expansion over the last 30 years and is the major focus of urban spatial analysis in this study.Over time, urban growth significantly changed in Abuja which gave way to complex urban dynamics such as conversion of agricultural land to settlement, road and infrastructure [37], population growth with an estimated annual growth rate of 9.4% [39].Table 1 presents an overview of how significant population growth in the federal capital city (FCC) differs from other parts of the FCT, which indicates that rapid population growth is a major driver to consider in this study.For instance, the decrease in agricultural activities and concomitant loss of cultivated land is likely to contribute to landlessness and food shortages and put the livelihood of inhabitants in jeopardy.To date, obtaining information on the environmental and socio-economic sustainability of Abuja, which is essential for development planning, has received relatively little attention.So far, no available systematic study of the spatiotemporal dynamics of urban growth changes in Abuja in the context of climate impacts has been conducted.In terms of short, medium and long-term development, this current study is appropriate to bridge the knowledge lacuna between these urban events.One of the few studies that has previously been conducted was detecting Land-use Land-cover change (LULCC) in Abuja and that was based on the maximum likelihood classifier [40].This study aims to integrate remotely sensed and ancillary data such as a master plan, digital elevation model and the population to detect LULCC from 1984 to 2014, and spatiotemporally analyze the settlement expansion pattern and model changes to project the LULC in 2050 using LCM in the context of climate impact.Major themes considered to achieve the goal of this study included biophysical information extracted from remotely sensed data as a baseline for subsequent applications such as analyzing settlement expansion which is the focus of this work.Other aspects investigated are the development of a land use change index and modelling urban growth (from the past into the future using LCM which can be contextualized for potential climate change impacts).The subsequent subsection presents and describes the study area considering the anthropogenic land-use situation and potential climate change impacts.

Study Area
The study area (Figure 1), also known as the FCC, falls within the FCT occupying approximately 400 km 2 .It is located in the center of Nigeria between 7 ˝20 1 and 9 ˝15 1 North of the Equator and longitudes 6 ˝45 1 and 7 ˝39 1 East.In the tropics under the Köppen climate classification, Abuja features a tropical wet and dry climate thereby experiencing three weather conditions annually.This includes a warm, humid rainy season and an intense dry season.The rainy season begins in April and ends in October.There is a brief interlude of harmattan occasioned by the northeast trade wind, with the main features being dust haze and dryness.This begins in November and lasts until January/February.The topography, high altitudes and undulating terrain of the FCT act as a moderating influence on the weather of the territory.Rainfall in the FCT reflects the territory's location on the windward side of the Jos Plateau and as a zone of rising air masses.The annual total rainfall is in the range of 1100-1600 mm.This study aims to integrate remotely sensed and ancillary data such as a master plan, digital elevation model and the population to detect LULCC from 1984 to 2014, and spatiotemporally analyze the settlement expansion pattern and model changes to project the LULC in 2050 using LCM in the context of climate impact.Major themes considered to achieve the goal of this study included biophysical information extracted from remotely sensed data as a baseline for subsequent applications such as analyzing settlement expansion which is the focus of this work.Other aspects investigated are the development of a land use change index and modelling urban growth (from the past into the future using LCM which can be contextualized for potential climate change impacts).The subsequent subsection presents and describes the study area considering the anthropogenic land-use situation and potential climate change impacts.

Study Area
The study area (Figure 1), also known as the FCC, falls within the FCT occupying approximately 400 km 2 .It is located in the center of Nigeria between 7°20′ and 9°15′ North of the Equator and longitudes 6°45′ and 7°39′ East.In the tropics under the Köppen climate classification, Abuja features a tropical wet and dry climate thereby experiencing three weather conditions annually.This includes a warm, humid rainy season and an intense dry season.The rainy season begins in April and ends in October.There is a brief interlude of harmattan occasioned by the northeast trade wind, with the main features being dust haze and dryness.This begins in November and lasts until January/February.The topography, high altitudes and undulating terrain of the FCT act as a moderating influence on the weather of the territory.Rainfall in the FCT reflects the territory's location on the windward side of the Jos Plateau and as a zone of rising air masses.The annual total rainfall is in the range of 1100-1600 mm.

Data
In this current study, three epochs of cloud free Landsat series images were downloaded through the Earth Explorer web portal [41] for 1986, 2001 and 2014 at no cost (Table 2).Other ancillary

Data
In this current study, three epochs of cloud free Landsat series images were downloaded through the Earth Explorer web portal [41] for 1986, 2001 and 2014 at no cost (Table 2).Other ancillary datasets Remote Sens. 2016, 8, 220 5 of 24 used comprised of a hard copy of the Abuja Master Plan (AMP) obtained from the survey and mapping department of the federal capital development authority (FCDA) [42], a 20 m digital elevation model (DEM) purchased from the office of the surveyor general of the federation (OSGOF) [43], used for slope generation and GoogleEarth Maps [44].The municipality administrative vector data was acquired from the National Space Research and Development Agency (NASRDA).Socioeconomic data such as census population was collected from the National Population Commission.Major road networks derived from a global position system field campaign in 2008 and a 2010 WorldView-2 very high resolution (VHR) satellite image were obtained from OSGOF.

Overview of Methodology
The downloaded satellite images for this study were ortho-rectified/georeferenced L1T (terrain corrected) product from source.However, the geometric accuracy was verified by overlaying and comparing with existing maps.Here, coordinate system verification and projection to UTM zone 32, WGS1984, Minna Datum was ascertained.Radiometric correction was performed by converting the digital numbers (DN) to at-sensor radiance using the radiometric calibration module in ENVI-5 in conjunction with the provided metadata in the header file.The Flash Line-of-sight Atmospheric Analysis of Spectral Hypercubes (FLAASH) atmospheric correction in ENVI-5 was applied to the at-sensor radiance corrected images of 1986, 2001 and 2014 with the appropriate atmospheric and aerosol models (e.g., Kaufman Tenre aerosol retrieval, tropical setting and urban aerosol model) to produce atmospherically corrected at sensor reflectance images.
Basically, the SVM classifier was applied to the preprocessed multi-temporal epochs of Abuja for 1986, 2001 and 2014 to produce retrospective and recent LULC maps of the study area.The classified images were validated using the accuracy assessment measure based on reference points from multiple sources (e.g., pre-existing points, historical GoogleEarth maps and field campaigns.Subsequently, to determine and quantify LULCC amounts for the three time stamps considered in this study, a robust change detection approach for time t1 to detect what has been converted to another category in time t2 was explored.The transitions between t1 and t2 are presented as a change matrix which serves as the input to the proceeding step; mainly, the calibration and transition modelling of target LULC types.To derive relevant LULC information for city or municipality managers that could be used for settlement expansion analysis, qualitative and quantitative settlement expansion footprints and annual land-use change rate (ALUCR) were computed.In addition, driving forces analysis was conducted for modelling the transition between 1986 and 2001 and was subsequently used to test the predictive power of the LCM by producing a predicted LULC map for 2014.The real LULC map of 2014 was compared to the modelled LULC map predicted from the t1 and t2 categorical maps as a way of validating the LCM LULCC model.After validating the modelling step between 2001 and 2014, the model was considered fit for using a similar prediction timeframe (approximately 15 years) to predict the LULC scenario for 2050.

Production of LULC Maps and LULCC Detection
The LULC maps of 1986, 2001 and 2014 produced in this study were generated using the SVM-based classifier, and five major LULC classes were extracted (Table 3).To avoid the "salt and pepper or speckle effect" often associated with pixel based image classification, a post-classification step based on the majority/minority analysis approach was applied.The option of the 3 by 3 filtering kernel tool was chosen which is suitable for enhancing the quality of the produced maps as well as retaining the Minimum Map Unit (MMU) in urban studies such as the smaller features representative of the urban landscapes.
The LULCC detection was based on two main approaches: the first is the binary method that relays areas/information of change and no change.The output of this technique has only two options, it shows whether changes in pixels have occurred or not between two epochs [45].A more detailed change detection technique is the post-classification comparison approach.This second method yields a "from-to" trajectory with an accompanying complete matrix containing records of exact count of transformed pixels from an initial class to another category [45].However, the change or no-change detection method is often associated with the problem of precise threshold identification.The thresholding methods usually have external influences on the resulting differences and can be linked to atmospheric conditions (e.g., sun angle, variable soil moisture, phonological differences) coupled with the threshold which is subjective to the image scene and is somewhat dependent on analysis of the study area [45].In the case of the post-classification method, autonomous image classification output is proceeded by applying a thematic overlay concept.This results in a complete "from-to" change matrix informative of transformations between categories of the multi-date maps [46].This technique can be associated with high error propagation effects emanating from independent LULC maps.For instance, images of two dates with a classification accuracy of 80% is likely to have 0.80 ˆ0.80 ˆ100%= 64% combined classification accuracy [46].The combination of the binary and post-classification change detection techniques yields a binary mask containing the extent of changed and unchanged areas at two time stamps [47].The resultant change mask was subsequently overlaid on the second classified map which facilitates the identification of changed pixels at the t2 classified imagery [47].In this current study, a hybrid approach which combines the change and no-change area binary maps of 1986-2001 and 2001-2014 derived from the Normalized Difference Vegetation Index (NDVI) was considered.Applying a NDVI thresholding approach directly for change detection underperformed in detecting many settlement and urban features, however, a direct image differencing approach by combining the near infrared difference with red difference performed better for urban change detection studies [48].Consequently, an overlay of the areas of change was implemented to construct a change matrix of the study periods 1986-2001 and 2001-2014.

Accuracy Assessment
Accuracy assessment was based on field reference data and an Orthophoto image from the survey and mapping department of FCDA for the 1986 image.The "Historical View" tool in Google Earth engine was used to validate the image from 2001.Field data, 2013 VHR and Google Earth images were used for the 2014 classification as reference.Given that the stratified random sampling approach employs the probability sampling design techniques which is key in demonstrating statistically rigorous assessment, for each of the produced LULC maps, a set stratified random points (200, 250 and 250 for 1986, 2001 and 2014, respectively) were generated using the random point tool in ArcGIS 10.1 toolbox (Table 4).The smoothed LULC maps were converted to polygons using the data conversion tool in ArcGIS 10.1.Subsequently, a column/field titled "Actual" was created for validating the classified map against the actual LULC.The randomly created points were used to verify the smoothed classified maps based on the reference data and actual LULC visualized from the reference data/VHR images (see italics row and columns in Table 4).The output matrices were further used to calculate overall, user's and producer's accuracies.The Kappa statistic, one of the commonly applied measures of difference between actual and change agreement, was computed [49].

Error Adjusted Area Assessment
In order to derive an error free area of LULCC information, the adjusted error matrix was computed.This involves three basic steps namely: (i) sampling design; (ii) response design and error adjusted area assessment (EAAA) analysis.The first step in EAAA is the sampling design stage.This is a practice used for determining the spatial unit subset that is subsequently used at the final EAAA step [50,51].For these EAAA approaches, a collection of reference points were also based on the stratified random points to facilitate the computation of error adjusted change area detail as well as validate the SVM-based LULCC maps which requires a probability sampling approach [52].This procedure allows the scientist to compute confidence interval estimates with statistical inferences and uncertainty measures.
Methodically, the five major LULC categories were used as sample strata.The LULC categories were suitable for stratification since they are the region of interest for communicating accuracy assessment results as well as the error adjusted class area estimates [50].The stratified sampling is achievable by either equalizing the sample size of each class or making the sample size proportional to the target category spatial extent on the map.Equalizing the samples is more advantageous to the user's accuracy than the producer's and overall accuracies [53], whereas using the proportional sample size approach makes the standard errors of the estimated overall accuracies and producer's appear smaller in comparison with the equal sampling [50].To robustly compute the EAAA, both sampling approaches were adopted, such as increasing the sample size of the target and s rear category such as built-up area in 1986 which conforms to Olofsson et al.'s [50] recommendation for obtaining good EAAA results.
The response design stage is the proceeding step used to compute the EAAA.It is a technique suitable for determining the sampling unit of reference data for LULC classification [53].Accordingly, Table 4 was used for collecting the reference data showing the sample sizes and strata of LULC maps of 1986, 2001 and 2014.The final accuracy assessment analysis involves using the classified map containing m classes for constructing an error matrix [54].Adopting the recommendation proposed by Olofsson et al. [50] for computing LULCC maps accuracies using stratified random sampling techniques in pixel-based classification maps, the area error inherent in the confusion matrices of the produced LULC classification and LULCC maps was adjusted for each category on the maps, while their area of proportions p Pij q was computed: where the area of proportion belonging to a category i in the map is W i , while n ij refers to the sample count labelled as i which belongs to class j contained in the reference dataset, the sample count mapped in category i is known as n i`o n the map.The adjusted error matrix of each cell element Pij indicating the probability of a randomly determined area is classified to category i on the image data to class j in the reference data [55].
Using the newly computed error matrix, overall, user's and producer's accuracies were computed.The overall accuracy p Ôq relays information on the overall proportion of correctly classified areas based on the sum of the diagonal adjusted error matrix represented as Pii .
Ûi " Pii Pi`( 3) The user's accuracy for class i expressedãs p Ûi q which is the proportion of area labelled as i having a class reference i, while producer's accuracy denoted as p Pj q of category j is proportional to the area of reference labelled j mapped as category j.
The area estimate and uncertainty were systematically computed using the approach proposed by Olofsson et al. [50], the constructed adjusted error matrix facilitated the computation of the area estimator using the area of proportion for category j.Here, the area of category j denoted as p Âj q was derived from Equation (5): A tot refers to the total area, while Equation ( 6) serves as the area estimator embedded in an error adjusted estimator.P`j " The area error adjusted estimator takes into account the area containing the error of omission present in the map category j and it eliminates the area of the map containing the error of commission [50].The computed standard error Sp Âj q was achieved based on Equation (7).
Sp Âj q " A tot ˆS `P `j˘( 7) Accordingly, the standard error computed for the stratified estimator based on the proportional area Sp P`j q was derived from Equation (8).
In order to assess uncertainty of the area estimates with sampling variability for computing confidence interval, the use of P`j derived from the reference samples was considered rather than the Pi`f rom the map areas.Hence, Âj which is the approximate of 95% confidence interval (CI) was computed based on Equation (9).
CI " Âj ˘z ˆSp Âj q (9) Remote Sens. 2016, 8, 220 9 of 24 where z denotes the percentile of the standard normal distribution curve and 95% confidence correspond to z score of 1.96.

Land Use Change Analysis and Urban Growth Indicators
Information on land use change is vital, not only in spatio-temporal urban growth and transport analysis, but also in different global, regional, local and urban context analyses.The dynamics of urban areas such as the driving forces of urban development is reflected by land use change [56][57][58].Hence, relevant indicators can be used as an effective means of quantifying and analyzing the spatio-temporal relationships between land use change and urban growth.In this study, one indicator was implemented to annualize the quantification of the spatio-temporal LULC change and urban expansion in Abuja.First, the so-called Annual Land Use Change Rate (ALUCR) was calculated, which was previously defined by Tian et al. [59] as follows: ALUCR a,t " pLU a,t ´LU a,t´1 q {LU a,t´1 pN t ´Nt´1 q ˆ100 (10) where ALUCR t (%) is the land use change rate; LU a,t and LU a,t´1 are the total land area of land use class a in hectares at the time t (current year) and time t ´1 (former year); N is the total number of years from time t (current year to time t ´1 (former year)).In this way, ALUCR enables class-wise change analysis.

LULCC Model Implementation and Validation
The LCM analysis was implemented using the IDRISI Selva software.The modelling procedure involves change analysis, modelling transition potential and determination of driving forces, predicting change and validation of the model.The generated LULC maps of Abuja city for 1986, 2001 and 2014 met the minimum requirement for implementing the LCM change analysis and prediction in 2050.The LCM is a rigid model structure which is based on a fixed flow of procedure with a minimum requirement of two LULC maps and the initial map serves as t1 and later serves as t2.The 1986 and 2001 LULC maps were used for the change analysis step to assess changes between the two time stamps.The two dates had the same specification (matching backgrounds, legends and spatial characteristics) and were used as the basis for understanding the nature of change in Abuja.These maps were used to generate transition potential maps and the probability matrix which were suited for identifying prevalent transitions to the urban category and this, in turn, was used as an input for the modelling of land cover changes.In LCM analysis, the modelling transition potential is vital for change location determination [34].The output of this step generates a series of transition potential maps that corresponds to significant land-cover transition into urban based on the change analysis phase implemented [34].The transition maps must consider the suitability of image pixels that have transformed into urban pixels considering a number of driving forces or factors useful for modeling processes of historical change.In this study, elevation, slope and distance to existing urban forms in 1986 were used as predictive variables.The distance to urban areas was used in the modeling phase, while elevation and slope were used independently as have been considered in similar recent studies [29,60] which were documented to have influence on settlement expansion, urban growth or sprawl and were evident in the LULCC between 1986 and 2001.In 1986, the distance to urban areas was set as a dynamic factor for recalculating the prediction period of (14 years), from 2001 to 2014.As settlements expands, so does the distance to urban forms and this changes over time from 1986.To assess the influence of topography as a driving force, a spatial overlay of the urban footprint was performed on the slope layer.Likewise, randomly generated elevation points corresponding to settlement were extracted from the digital elevation model (DEM) and superimposed on the computed hill-shade as a backdrop layer.
The LCM also allows users to conduct a quick optional test of the potential exploratory model power for driving force which is presented in the Cramer's V.This depicts the correlation coefficient indicating either no correlation (meaning a redundant variable) or perfect correlation referring to a significant potential variable ranging from 0.0 to 1.0, respectively [34].Notably, the Cramer's V does not assure good model performance due to its lack of consideration of mathematical variables and their complex relationship [28], but it helps to determine the usefulness of a driving force [34].In LCM, the MLP neural network is a feedforward neural one direction network that flows from an input to output with hidden layer(s) in between [61].The computation of nodes can be grouped as layers where each node receives an input signal to the destination node [33].For every input signal to another node, the subsequent layer contains the original input multiplied by a weighting and combined with a threshold which is channeled through the activation function in hidden linear or non-linear units [33].The weights are defined in the training stage preceding the network prediction stage which uses a proportion of the data to change the weight in order to minimize errors in the observed and predicted outputs [33].Generally, the MLP neural network is suitable for modelling multiple transitions at once [34], and was applied in this study.
The LCM prediction steps calculates the change rate from both the initial step and the transition potential maps generated from the second step in order to predict the future scenario of 2014.Based on the Markov-Chain analysis, implementing this step helps to ascertain the extent of transformation of other land-cover to urbanized areas in every transition that occurred in 2014 [34].From this step, the hard and soft predictions are the two basic types of predictions derived as outputs [34].The hard prediction result is the projected 2014 map showing pixels of the specific land-cover category which is the most likely class to be converted into.The soft prediction is to some extent a vulnerability map depicting pixels assigned to values ranging from 0.0 to 1.0, which is indicative of a pixel's probability to transform into an urban pixel in 2014 [34].
Model validation is an important step.In this study, validation was used to ascertain the quality of the predicted map of 2014 in comparison to actual 2014 LULC maps which is reality.Two notable approaches used in validating such models are the visual and statistical procedure [62].Using the visual approach yields a crisscrossed tabulation between the 2001 LULC map, the predicted 2014 map and the 2014 actual map used for assessing the accuracy of the map model.The result of this process is a map containing four main categories [62], namely: (i) Hits which signify model prediction of change and that the situation actually occurred; (ii) False alarms that refer to the predicted change from one category to urban area with persistence in urban class; (iii) Misses suggest persistent model prediction that eventually changed to urban class in actual sense and (iv) Null success depicts areas where the model predicted no change and, indeed, no change occurred.The hits and null success are indicative of model correctness, while the false alarms and misses signify errors due to disagreement between the model simulation map and the reference map [62].The figure of merit (FOM) is a ratio between hits, simulation hits, misses and false alarms and was computed to ascertain the overall agreement between the observed and predicted maps.FOM ranges from 0% signifying absence of overlap between the predicted and observed change, to 100% suggesting a perfect agreement between the predicted and observed change [63].The statistical procedure measures the agreement between paired maps showing variable numbers of LULC categories [34].The actual LULC map of 2014 is used as the reference map while the model simulation map was utilized for comparison.The kappa variation techniques were applied as the statistical validation procedure in this study.Here, the Kappa for no information (K no ) signifies overall accuracy obtained in the simulation run, while Kappa for location (K location ) means agreement level in a location [64].Considering the difficulty in interpretation encountered with the Kappa for correctly assigned proportion against the proportion of incorrectly assigned by change (K standard ) [65], the K standard was not used in this study.However, K location was somewhat useful for the validation in the absence of K standard [66].After assessing the model's predictive power, the model was applied to predict the LULC map of the 2050 scenario with the assumption of there being similar driving forces based on the changes between 1986 and 2014 LULC maps.Driving forces are not static, especially since the world is a dynamic system, so different triggers in variable rate and velocity are likely to evolve and may influence the urban evolution process as well as model results.

Production of LULC Maps and Accuracy Assessment
Figure 2 presents the SVM-based LULC maps of Abuja city for 1986, 2001 and 2014.Five major classes, namely water, built-up, vegetation, bare/arable land and complex-landscapes were derived from the multi-date satellite images.It is apparent from Figure 2 that changes occurred in proportions of area covered by the retrospective LULC map categories, especially for built-up land which is the major focus of this study, and analyzing these patterns is useful for gaining insights about the composition of the total map area and changes in patterns during the considered years of analysis.
Table 5 presents the accuracy assessment obtained based on the pixel count error matrices from the LULC maps.The overall accuracies derived from the pixel count matrices of the generated retrospective LULC maps for 1986, 2001 and 2014 were 83%, 92% and 94%, respectively, while kappa accuracies were 81%, 93% and 94% in 1986, 2001 and 2014, respectively (Table 5).
To derive accurate proportions of the targeted LULC categories, the error adjusted area matrix (EAAM) computed based on the pixel count matrix was produced (Table 6).The EAAM gives more informative and accurate detail in terms of omission and commission errors leading to confusion and misclassification in the three LULC maps (Table 6).The robustness of the EAAM can be confirmed by comparing the overall accuracies derived from the pixel count with the one recorded from the EAAM.
triggers in variable rate and velocity are likely to evolve and may influence the urban evolution process as well as model results.

Production of LULC Maps and Accuracy Assessment
Figure 2 presents the SVM-based LULC maps of Abuja city for 1986, 2001 and 2014.Five major classes, namely water, built-up, vegetation, bare/arable land and complex-landscapes were derived from the multi-date satellite images.It is apparent from Figure 2 that changes occurred in proportions of area covered by the retrospective LULC map categories, especially for built-up land which is the major focus of this study, and analyzing these patterns is useful for gaining insights about the composition of the total map area and changes in patterns during the considered years of analysis.
Table 5 presents the accuracy assessment obtained based on the pixel count error matrices from the LULC maps.The overall accuracies derived from the pixel count matrices of the generated retrospective LULC maps for 1986, 2001 and 2014 were 83%, 92% and 94%, respectively, while kappa accuracies were 81%, 93% and 94% in 1986, 2001 and 2014, respectively (Table 5).
To derive accurate proportions of the targeted LULC categories, the error adjusted area matrix (EAAM) computed based on the pixel count matrix was produced (Table 6).The EAAM gives more informative and accurate detail in terms of omission and commission errors leading to confusion and misclassification in the three LULC maps (Table 6).The robustness of the EAAM can be confirmed by comparing the overall accuracies derived from the pixel count with the one recorded from the EAAM.The drop in the 1986 LULC map overall accuracy threshold of 83% based on pixel count to 69% using the EAAM is due to the inclusion of error of commission in estimating the proportion of area.
Similarly, a decline in the resultant pixel count overall accuracies was observed in the 2001 and 2014 LULC maps from 92% and 94% to 91% for both based on the EAAM (Table 6).Some confusion was observed between categories (e.g., BUP, VEG, BAL and CL).Major confusion occurred in the form of errors of commission given that incorrectly classified pixels are prominent in the rows, especially between BAL and CL in 1986 and 2001 (Table 6).Similarly, errors of omission are evident between VEG, BAL and CL in 1986 and 2001 (Table 6).

LULCC Detection and Spatio-Temporal Analysis of LULC Distribution
The area distribution of the various LULC is reported in Tables 7-9 respectively.Comparing the pixel count matrices and the error adjusted matrix revealed a slight difference between the mapped and estimated area for the LULC categories.For all the categories, the produced maps fell within a 95% confidence interval of the estimated area which is indicative of the robustness and reliability of the produced maps [50].
The spatial arrangement of physiographic features was obtained by subjecting the classified images to spatial analysis, the information about changes in land cover proportion was derived.Subsequently, changes in the spatial composition of geographic features was obtained by comparing the classified images of 1986, 2001 and 2014 against each other (Table 10).
As displayed in Table 11, noticeable LULC changes appeared during the 28-year study timesteps.In this study, significant transitions from other classes into the water category were from vegetation and bare/arable land.During the 28 year time period, 0.03% of the vegetation transformed into water, as did 0.2% of bare/arable land.Approximately, 24.9% of bare/arable-land had been converted to built-up while 4.8% of area mapped as bare/arable-land in 1986 is now vegetation.Generally, persistence in all the classes stood at 4.3%.

Analysis of Settlement Expansion, Floodplain Encroachment and Driving Forces
The evolution of the built-up land category in Abuja city can also be seen in the LULC maps presented in Figure 2. From the time stamps, it was apparent built-up land increased by more than 50% from 1986 to 2001, and between 2001 and 2014, the expansion was approximately 48%, which is indicative of unprecedented urban growth and which conforms to the quantitative measures in Table 10.
For further spatio-temporal analysis, the overlay analysis of the extracted 2014 built-up time stamp onto the AMP map facilitated the evaluation of how settlement development did not conform to the plan.The AMP is a comprehensive document prepared to guide development in Abuja city to avoid chaotic urban development.From the overlay, a lateral settlement development situation was observed in Figure 3 and the uncontrolled lateral settlement expansion in Abuja suggests possible environmental challenges, such as floods due to floodplain encroachment, as highlighted with black boxes (Figure 3).
The computed annual land use change rate (ALUCR) analysis further established the LULC transformation rate in Abuja between 1986 and 2014 (Table 12).Built-up land being the core interest for Abuja city dramatically increased by 8.51% between 1986 and 2001 and by 6.31% from 2001 to 2014 (Figure 2 and Table 12).

Analysis of Settlement Expansion, Floodplain Encroachment and Driving Forces
The evolution of the built-up land category in Abuja city can also be seen in the LULC maps presented in Figure 2. From the time stamps, it was apparent built-up land increased by more than 50% from 1986 to 2001, and between 2001 and 2014, the expansion was approximately 48%, which is indicative of unprecedented urban growth and which conforms to the quantitative measures in Table 10.
For further spatio-temporal analysis, the overlay analysis of the extracted 2014 built-up time stamp onto the AMP map facilitated the evaluation of how settlement development did not conform to the plan.The AMP is a comprehensive document prepared to guide development in Abuja city to avoid chaotic urban development.From the overlay, a lateral settlement development situation was observed in Figure 3 and the uncontrolled lateral settlement expansion in Abuja suggests possible environmental challenges, such as floods due to floodplain encroachment, as highlighted with black boxes (Figure 3).
The computed annual land use change rate (ALUCR) analysis further established the LULC transformation rate in Abuja between 1986 and 2014 (Table 12).Built-up land being the core interest for Abuja city dramatically increased by 8.51% between 1986 and 2001 and by 6.31% from 2001 to 2014 (Figure 2 and Table 12).The assessment of driving forces such as transportation infrastructure and suitable topography is presented in Figures 4 and 5  The assessment of driving forces such as transportation infrastructure and suitable topography is presented in Figures 4 and 5, but is later discussed in Section 4.2.

Land Use Change Model Implementation and Validation
Based on the LULCC maps, significant transitions occurred for all four categories: from water, vegetation, bare/arable-land and complex-landscape to built-up, hence, all were used as the major transitions in the model.
The transitions displayed varying urban spatial trends indicating that a divergent predictor variable will affect them.Table 13 presents the Cramer's V of the relevant and considered factors.It is a representation of the potential explanatory power of the major LULC category as a driving force in the settlement expansion process.The Cramer's V of 0.15 and above for any variable is useful, while values ≥ 0.4 are good [34].The proceeding step was to run the sub-model of relevant transitions and produce transition maps.

Driving Force
Cramer's V Distance to Built-up in 1986 0.51  The assessment of driving forces such as transportation infrastructure and suitable topography is presented in Figures 4 and 5, but is later discussed in Section 4.2.

Land Use Change Model Implementation and Validation
Based on the LULCC maps, significant transitions occurred for all four categories: from water, vegetation, bare/arable-land and complex-landscape to built-up, hence, all were used as the major transitions in the model.
The transitions displayed varying urban spatial trends indicating that a divergent predictor variable will affect them.Table 13 presents the Cramer's V of the relevant and considered factors.It is a representation of the potential explanatory power of the major LULC category as a driving force in the settlement expansion process.The Cramer's V of 0.15 and above for any variable is useful, while values ≥ 0.4 are good [34].The proceeding step was to run the sub-model of relevant transitions and produce transition maps.

Driving Force
Cramer's V Distance to Built-up in 1986 0.51

Land Use Change Model Implementation and Validation
Based on the LULCC maps, significant transitions occurred for all four categories: from water, vegetation, bare/arable-land and complex-landscape to built-up, hence, all were used as the major transitions in the model.
The transitions displayed varying urban spatial trends indicating that a divergent predictor variable will affect them.Table 13 presents the Cramer's V of the relevant and considered factors.It is a representation of the potential explanatory power of the major LULC category as a driving force in the settlement expansion process.The Cramer's V of 0.15 and above for any variable is useful, while values ě 0.4 are good [34].The proceeding step was to run the sub-model of relevant transitions and produce transition maps.Figure 6 is a map illustration of the ocular validation of the projected urban LULCC of 2014 showing the precision achieved by the LCM model.Figure 7 depicts the validation map based on the crisscross which consists of hits, misses, false alarms and null success, with 1.08% hits achieved, 6.04% misses recorded, false alarms being 5.89% and null successes 86.99%.6 is a map illustration of the ocular validation of the projected urban LULCC of 2014 showing the precision achieved by the LCM model.Figure 7 depicts the validation map based on the crisscross which consists of hits, misses, false alarms and null success, with 1.08% hits achieved, 6.04% misses recorded, false alarms being 5.89% and null successes 86.99%.Bare/Arable Land 0.73 Vegetation 0.21 Figure 6 is a map illustration of the ocular validation of the projected urban LULCC of 2014 showing the precision achieved by the LCM model.Figure 7 depicts the validation map based on the crisscross which consists of hits, misses, false alarms and null success, with 1.08% hits achieved, 6.04% misses recorded, false alarms being 5.89% and null successes 86.99%.

Relevance of Data Adequacy for LULCC Mapping and Spatio-Temporal Analysis
Depending on the overall goal and context, target parameters, such as LULCC, could either be mapped or monitored, which is mostly dynamic and is best studied within a continuity framework.Given that LULCC can be linked to urbanization and global warming trends, adequate image time series are needed to identify changes and track trends from a historical perspective which is the basis for modelling the future.However, in tropical regions such as West Africa, especially Nigeria, cloud cover hampers the opportunity presented by the Landsat image archives to select the desired image time series from a particular Landsat sensor and, thus, leads to image scarcity.This challenge hindered LULCC analysis, such as of the settlement expansion and urbanization phenomena this study investigated.However, mapping and monitoring urbanization for analyzing imminent problems with limited data is still valuable for apt policy and decision-making, without which city managers would not be able to compare such empirical findings of urban development processes with other socio-economic indicators.Hence, due to the scarce availability of cloud free Landsat images for the same season over the study area, the authors adopted the multi-sensor and multitemporal (MSMT) concept demonstrated in this study.However, this approach is associated with challenges that require thorough image pre-processing steps that must be fulfilled (see Section 2.3) since the direct application of Digital Numbers (DNs) will not suffice anymore in a MSMT context [67,68].Therefore, conversion of DN to at-surface reflectance was performed to produce high quality scientific data for reliable downstream products such as the LULC maps.This procedure yields data that is illumination and atmospheric artefact independent and, therefore, plausible for making comparisons, such as in the change analysis performed in this study.
From the LULC change analysis of Abuja, it was found that water increased by 0.1% between 1986 and 2001 and, by 2014, another 0.1% increase was recorded.Generally, it is apparent that water surfaces expanded by 2% between 1986 and 2014 (Table 10).Substantial increase in the built-up class occurred between 1986 and 2001, by 403,598.3ha(20.8%), which is an average of more than 26,906 ha/year for the timeframe.Built-up area significantly expanded, by 590,274.1 ha (37.8%) from 2001 to 2014, which is more than 42,162.4ha/year, with a net growth of approximately 993,872.4 ha in urban area in the study period (Table 10).A vegetation increase of 51,703.8ha could be measured in 1986-2001 which is also apparent in Figure 2.However, between 2001 and 2014, vegetation declined drastically by 510,329.3 ha which can be linked to rapid deforestation, urban growth and a lack of parks and

Relevance of Data Adequacy for LULCC Mapping and Spatio-Temporal Analysis
Depending on the overall goal and context, target parameters, such as LULCC, could either be mapped or monitored, which is mostly dynamic and is best studied within a continuity framework.Given that LULCC can be linked to urbanization and global warming trends, adequate image time series are needed to identify changes and track trends from a historical perspective which is the basis for modelling the future.However, in tropical regions such as West Africa, especially Nigeria, cloud cover hampers the opportunity presented by the Landsat image archives to select the desired image time series from a particular Landsat sensor and, thus, leads to image scarcity.This challenge hindered LULCC analysis, such as of the settlement expansion and urbanization phenomena this study investigated.However, mapping and monitoring urbanization for analyzing imminent problems with limited data is still valuable for apt policy and decision-making, without which city managers would not be able to compare such empirical findings of urban development processes with other socio-economic indicators.Hence, due to the scarce availability of cloud free Landsat images for the same season over the study area, the authors adopted the multi-sensor and multi-temporal (MSMT) concept demonstrated in this study.However, this approach is associated with challenges that require thorough image pre-processing steps that must be fulfilled (see Section 2.3) since the direct application of Digital Numbers (DNs) will not suffice anymore in a MSMT context [67,68].Therefore, conversion of DN to at-surface reflectance was performed to produce high quality scientific data for reliable downstream products such as the LULC maps.This procedure yields data that is illumination and atmospheric artefact independent and, therefore, plausible for making comparisons, such as in the change analysis performed in this study.
From the LULC change analysis of Abuja, it was found that water increased by 0.1% between 1986 and 2001 and, by 2014, another 0.1% increase was recorded.Generally, it is apparent that water surfaces expanded by 2% between 1986 and 2014 (Table 10).Substantial increase in the built-up class occurred between 1986 and 2001, by 403,598.3ha(20.8%), which is an average of more than 26,906 ha/year for the timeframe.Built-up area significantly expanded, by 590,274.1 ha (37.8%) from 2001 to 2014, which is more than 42,162.4ha/year, with a net growth of approximately 993,872.4 ha in urban area in the study period (Table 10).A vegetation increase of 51,703.8ha could be measured in 1986-2001 which is also apparent in Figure 2.However, between 2001 and 2014, vegetation declined drastically by 510,329.3 ha which can be linked to rapid deforestation, urban growth and a lack of parks and gardens.
The increase in vegetation in 1986-2001 can be attributed to annual variation (in 2001 comparatively high precipitation values were recorded) and increasing precipitation trends after the 1980 drought period.In addition to variability of natural conditions, the protection of green areas received more focus, which can be seen in the establishment of the millennium park and the developments of parks and gardens in the city center of Abuja.
Of the total change in built-up area, most of the areas were previously bare/arable land which formed the newly urbanized areas in the study period.Apparently, the dominant land use types in 1986 were natural vegetation cover and bare/arable land for cultivation, and the urban area (referred to as the built-up category herein) was found in the city center (Figure 2).The increments seen in the built-up category reveal rapid urban development was a result of the demand for residential, public offices, academic and business/commercial facilities.This is in line with the substantial growth witnessed in Abuja between 1986 and 2001, as highlighted in Table 10.From the classification results, it can be seen that urban expansion started to evolve in all directions (Figure 2).Noticeably, between 1986 and 2001 the observed decline in bare/arable land category by 17.8% is a significant change which can be attributed to the drought in the 1980s.Vegetation and complex landscape occupied 23.9% and 14.1%, respectively, of the study area in 2001, but this declined to 9.2% and 6.3% by 2014, respectively, due to both human and natural processes as well as misclassification errors.For instance, the complex landscape is a complicated category because, in this class, several activities are taking place, such as farming as well as quarrying in some locations.Similarly, landscapes such as murky water, dark soil and wetlands have identical spectral signatures and remain a source of misclassification in this study.Remarkably, the built-up category grew to 17% in the 2001-2014 study period.An explanation for the slight increase of (5.4%) bare/arable land can be linked to farmland expansion while the decline in vegetation may be attributed to intensified urban expansion in Abuja city.

Settlement Expansion, Floodplain Encroachment and Driving Force Analysis
With climate change impacts such as increased rainfall and decreased flood return period, encroached areas make buildings and people in such locations susceptible to flood impacts.Uncontrolled lateral expansion in Abuja also poses direct and indirect threats to essential ecosystems from a future settlement development perspective.According to Seto et al [69], in the coming decade, rapid land conversion is likely to occur in biodiversity hotspots, which in year 2000 were undisturbed.
In the case of Abuja, urban growth was expected as there was provision for expansion in the master plan.However, according to FCDA [37], the actual land budgeted for specific purposes had been exceeded, and proximity to existing buildings or infrastructure is a factor that puts pressure on the available land irrespective of the land use specified for such parcels.Hence, for integrated spatial and socioeconomic analysis, the ALUCR can be compared with the population growth rate and other relevant indices for data driven decision making.The rapid growth in Abuja can be tied to population growth, construction of transportation infrastructure, business development and the government's development policies, such as massive housing estate development projects.Similarly, other potential social problems in Abuja due to rapid urban growth include the accelerated loss of highly productive farmlands, alteration of energy budgets with modifications to climatic, hydrologic, and biogeochemical cycles, habitat fragmentation or reduction in biodiversity [70].Looking at the negative ALUCR of ´4.73 between 2001 and 2014 for vegetation and the positive ALUCR of 1.02 between 2001 and 2014, this situation can be linked to micro-climate modifications such as urban heat island development in Abuja.
The observed LULC changes and settlement expansion in Abuja are driven by intertwined socioeconomic and environmental factors.Socioeconomically, some of the growth in Abuja can be associated with the oil boom and the relocation of government department headquarters to the center.During this period, a dramatic settlement expansion pattern was noticed which can be linked to the upward changing pattern of both the economy and the population.This was a period of rapid transportation infrastructure development in both the city and satellite towns (Figure 4).Two notable patterns of expansion are outward growth of the city based on the master plan and sprawl development towards the satellite towns with a pattern discerned along the outer northern and southern expressways of the city.
Geographically, the spatial overlay of extracted built-up land for 2014 onto the slope map generated from a 20 m DEM proved that the built-up development pattern of Abuja city is supported by the suitable topography and slope of the converted landscape (Figure 5a).Results obtained from the overlaid 2014 LULC built-up footprint on DEM-derived hill-shade as the backdrop to Abuja city are presented in Figure 5b.The constrained analysis performed on built-up samples shows that significant urban expansion occurred in topographically flat and suitable landscapes of ď539 meters (Figure 5b).

Land Use Change Model Implementation and Validation
The visual comparison presented in Figure 6 showed that the simulated 2014 built-up footprint performed satisfactorily.However, over prediction is apparent in the upper right corner of Figure 6a when compared to the actual data in Figure 6b which could be due to inadequacy of driving force restriction and model uncertainty.A 37.1% built-up extent in 2014 was achieved by the model prediction for the 2014 urban landscape, while the actual built-up extent in 2014 was 40.4%.A difference of 3.3% was observed, and this can be a source of the error found in the model.In this study, a FOM of 12.6% was achieved which can be regarded as a poor performance, however, it is still comparable with some recent studies such as by Olmedo [71] and Megahed [29] where FOM was 10.4% and 12.76%, respectively.This also suggests that, in modelling large areas such as Abuja, consideration of more complex predictive variables as driving forces will be needed to achieve better results.While visual assessment remains the swiftest approach to analyze spatial patterns, which is usually a weakness associated with the current statistical technique, ocular examination remains subjective and somewhat misleading.Hence, the current statistical technique is still vital in model validation [62].The Kappa variation was used to compare the projected 2014 LULC map to the real LULC map of 2014 and realized a K no of 85% and a K location of 95%.These high accuracy measures suggest that changes in the categorical maps are systematic for the 28 year (1986-2014) epoch assessed.It is also critical to carry out visual assessment throughout the FOM analysis because it helps the analyst identify localized areas of change rather than regions of no change [29].
The LCM prediction Figure 8a shows how built-up area will significantly grow and take over adjacent land categories in the absence of policies that ensure strict adherence to the master plan or sustainable spatial planning, in comparison to Figure 8b in the middle.This is indicative of continuous removal of natural cover as well as ecosystem destruction.Modification of natural cover such as vegetation will certainly reconfigure the energy balance (temperature development) such as the formation of an urban heat island (UHI) due to impervious surfaces, which has been established and well documented based on building configurations including shape, size, and geometry [11].Therefore, the model result provides an insight into the need for a carbon mitigation strategy and green urban design in Abuja due to potential disappearance of vegetation cover.Also, Figure 8c illustrates the spatial pattern of built-up land in Abuja suggesting numerous environmental challenges (e.g., UHI potential, distortion of urban hydrology).Complementary to the UHI, urbanization can also influence precipitation patterns.In literature, two factors noted are: (i) the influence of urbanization in altering convective forces driving precipitation and; (ii) the variability of aerosol abundance in urbanized spaces [72].Also, an observation of dedicated urban climate studies has been the occurrence of locational shifts and regional increases in convective precipitation linked to urbanization which thus results in an increase in surface temperature [73].Thus, with the dozens of studies pointing to the impact of urbanization on precipitation, Abuja is unlikely to be an exception even though research outputs might differ if it is further explored.

Framework for Future Research
The absence of previous research on a city such as Abuja hampers opportunities for a comparison of the findings of this study regarding data selection, methodology used and results obtained.This study is just the tip of the iceberg.However, the research provides a sufficient basis for further studies, and, therefore, can be classified as foundation research that urgently needs to be followed-up with more integrated approaches that will address aspects this study did not consider.For instance, in accordance with numerous projections, rapid settlement expansion can be seen as unavoidable, especially considering the total world urban population is currently at approximately 7.4 billion, and by 2050, it is anticipated to range from 8.3 to 10.9 billion.In developing countries, 70%-80% of the population will be urban dwellers.As a result of this, urban planners proposed two coping strategies to address the ever-increasing population pressure exerted on urban areas.The first is through population control while the second put forwards the need for an efficient urban management agenda.The initial proposal has proven to be very challenging, especially in developing regions like West Africa, including Nigeria.Sustainable urban management is founded on two major planning processes, namely strategic and operational planning.In the case of urban planning, one of the principal datasets needed is land use.Fundamentally, land use data is useful for strategic and operational planning, particularly for detecting and monitoring changes, trends analysis and for future prediction based on scenarios/consequences and review of master plans.Therefore, a relevant aspect to be further investigated that is related to the generated dataset includes assessing the influence of settlement expansion on land surface temperature and UHI formation, flood risk zonation of the city and its integration with urban density information/urban structural types.This is integral to assessing risk and vulnerability to environmental phenomena such as climate change.It will also be worthwhile to compare the SVM based approached used in this study with other similar methods that have an integrated GIS approach with artificial neural network and fuzzy logic by incorporating additional driving force constraints based on additional data [32,74,75], and other tested integrated urban remote sensing and GIS approaches [76].

Conclusions
The FCC of Nigeria, Abuja, was used as the test site for this study using integrated remote sensing datasets and GIS modelling approaches.It was established that significant LULC change and settlement expansion has occurred.Due to the study timeframe and there being limited cloud free datasets available from a single Landsat sensor, the multi-sensor and multi-temporal images of Landsat (TM, ETM+ and OLI for 1986, 2001 and 2014) were used for this study.The utility of these datasets was realized through rigorous image pre-processing to ensure that real change over time was measured and the results obtained were reliable and fit for further use (e.g., informed policy and decision-making).Also, the application of a robust information extraction algorithm such as SVM required limited training samples and yielded good image classification.In future, the SVM and other advanced methods such as object oriented image analysis and random forest can either be integrated or compared to assess their performance in problematic or heterogeneous areas.
The core objective was to spatio-temporally analyze LULC for 1986, 2001 and 2014.Based on the pixel count error matrix, overall accuracies of the three LULC maps ranged from 82% to 94%.Contrary to these statistics, the adjusted area error gave a somewhat different accuracy measure that ranged from 69% to 91%.The study found that computing area information directly from pixel count for accuracy assessment can be misleading.Hence, in producing local maps that can be useful for evaluating the accuracy of global maps, climate modelling and other relevant applications, assessing error propagation in regional maps by applying the error adjusted area estimator yields a more reliable and informative result including confidence intervals and uncertainty measures.The different components of LULC change analysis allowed better understanding of the transformation processes, especially for the transition from the bare/arable land and vegetation categories to the built-up area class.The computed indices provided empirical insight into a realistic spatio-temporal situation, and detailed annual land use and urban spatial expansion change rates in the study timeframe with major impacts on the landscape and potential influence on the local climate of Abuja.The forces governing such a momentous expansion might be numerous, however, apparent drivers in the case of Abuja included suitable topography, availability of public infrastructure, such as social amenities, and population growth, particularly in the past two decades.
Two salient aspects ascertained by this research are that massive impervious surface development is occurring due to urbanization, which may lead to elevated urban temperatures known as the UHI phenomenon.This phenomenon, is one of the likely climate control factors at a local scale highlighted by the Fifth Assessment Report of the IPCC [77], and COP 21 [78].Complementary to the UHI is the change in drainage geography that may also increase surface runoff which translates to flash flood events in cities.With rapid urban expansion in Abuja, when rainfall increases, surface runoff is expected to increase and can trigger flash flood events in areas deficient of adequate drainage, especially looking at the 2050 built-up area projection in the context of climate change.From this study, further research is needed to empirically verify present and future impacts of urbanization on Abuja city.

Figure 1 .
Figure 1.Location of study area.(a) West Africa inset Nigeria; (b) Nigeria inset Federal Capital City (FCT) and Abuja city in red rectangle and (c) Zoom of Abuja City showing major districts.

Figure 1 .
Figure 1.Location of study area.(a) West Africa inset Nigeria; (b) Nigeria inset Federal Capital City (FCT) and Abuja city in red rectangle and (c) Zoom of Abuja City showing major districts.

Figure 3 .
Figure 3. (a) Overlay of digitized open space and undevelopable area from the Abuja Master Plan (AMP) onto extracted 2014 built-up layer and AMP; (b) zoom-out of the overlay of digitized open space and undevelopable area from the AMP of Abuja for better visualization and (c) zoom of the AMP legend.

Figure 3 .
Figure 3. (a) Overlay of digitized open space and undevelopable area from the Abuja Master Plan (AMP) onto extracted 2014 built-up layer and AMP; (b) zoom-out of the overlay of digitized open space and undevelopable area from the AMP of Abuja for better visualization and (c) zoom of the AMP legend.

Figure 4 .
Figure 4. Example of the road expansion between 2003 (a) and 2015 (b) in the outer Northern Express Way of Abuja (image source Google Earth).

Figure 5 .
Figure 5. Assessment of the influence of topography on urban expansion.(a) Shows the GIS overlay of built-up on DEM derived slope.While (b) presents a constraint analysis showing built-up dependency on suitable heights as overlaid onto hillshade.

Figure 4 .
Figure 4. Example of the road expansion between 2003 (a) and 2015 (b) in the outer Northern Express Way of Abuja (image source Google Earth).

Figure 4 .
Figure 4. Example of the road expansion between 2003 (a) and 2015 (b) in the outer Northern Express Way of Abuja (image source Google Earth).

Figure 5 .
Figure 5. Assessment of the influence of topography on urban expansion.(a) Shows the GIS overlay of built-up on DEM derived slope.While (b) presents a constraint analysis showing built-up dependency on suitable heights as overlaid onto hillshade.

Figure 5 .
Figure 5. Assessment of the influence of topography on urban expansion.(a) Shows the GIS overlay of built-up on DEM derived slope.While (b) presents a constraint analysis showing built-up dependency on suitable heights as overlaid onto hillshade.

Figure 7 .
Figure 7. Ocular map validation of correctness using 2001 and 2014 as references and the projected LULC map for 2014.

Figure 8
Figure 8 presents the LULC projection into the year 2050 that was performed, similar to the model testing implemented between the 2001 and 2014 time stamps, the actual 2014 LULC map and built-up area expansion from 1986 to 2050.

Figure 7 .
Figure 7. Ocular map validation of correctness using 2001 and 2014 as references and the projected LULC map for 2014.

Figure 8
Figure 8 presents the LULC projection into the year 2050 that was performed, similar to the model testing implemented between the 2001 and 2014 time stamps, the actual 2014 LULC map and built-up area expansion from 1986 to 2050.

Figure 7 .
Figure 7. Ocular map validation of correctness using 2001 and 2014 as references and the projected LULC map for 2014.

Figure 8
Figure 8 presents the LULC projection into the year 2050 that was performed, similar to the model testing implemented between the 2001 and 2014 time stamps, the actual 2014 LULC map and built-up area expansion from 1986 to 2050.

Table 1 .
Population statistics of the municipalities in FCT for the year 2006.

Table 2 .
Satellite Image description.

Table 3 .
Land cover classes used in the study.

Land Use/Cover Types Description Water
River, permanent open water, lakes, ponds and reservoirs Built-Up Residential, commercial and services, industrial, transportation, roads, mixed and other urban features Vegetation Coniferous and deciduous forest, mixed forest lands, shrubs and grassland (green areas) Bare/Arable-land Exposed soils, area of active excavation, Agricultural areas crop fields and fallow land Complex Landscapes Rock-outcrops metamorphosed rock dark soils, wetland and mucky water

Table 4 .
Sample sizes allotted to the targeted LULC classes in 1986, 2001 and 2014.

Table 5 .
Summary of Land Use Land Cover (LULC) map accuracies (%) for 1986, 2001, and 2014 based on the pixel count matrices.

Table 5 .
Summary of Land Use Land Cover (LULC) map accuracies (%) for 1986, 2001, and 2014 based on the pixel count matrices.

Table 7 .
Proportion of LULC category for 1986.

Table 8 .
Proportion of LULC category for 2001.

Table 9 .
Proportion of LULC category for 2014.

Table 10 .
Spatial analysis result of error adjusted matrices of the 1986, 2001 and 2014 LULC maps from showing map category, class area in hectares and percentage and area changed in hectares.

Table 12 .
Annual Land Use Change Rate (ALUCR) and use change in Abuja from 1986 to 2014.

Table 12 .
Annual Land Use Change Rate (ALUCR) and use change in Abuja from 1986 to 2014.
but is later discussed in Section 4.2.

Table 13 .
Cramer's V driving force threshold for potential LULCC to built-up.

Table 13 .
Cramer's V driving force threshold for potential LULCC to built-up.

Table 13 .
Cramer's V driving force threshold for potential LULCC to built-up.