Next Article in Journal
Design and Implementation of a Robust Decision Support System for Marine Space Resource Utilization
Previous Article in Journal
Occlusion-Free Visualization of Important Geographic Features in 3D Urban Environments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatiotemporal Modeling of Urban Growth Predictions Based on Driving Force Factors in Five Saudi Arabian Cities

1
School of Environmental and Rural Science, University of New England, Armidale, NSW 2351, Australia
2
Department of Geography, Umm Al-Qura University, Makkah 21955, Saudi Arabia
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2016, 5(8), 139; https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi5080139
Submission received: 6 July 2016 / Revised: 1 August 2016 / Accepted: 2 August 2016 / Published: 8 August 2016

Abstract

:
This paper investigates the effect of four driving forces, including elevation, slope, distance to drainage and distance to major roads, on urban expansion in five Saudi Arabian cities: Riyadh, Jeddah, Makkah, Al-Taif and Eastern Area. The prediction of urban probabilities in the selected cities based on the four driving forces is generated using a logistic regression model for two time periods of urban change in 1985 and 2014. The validation of the model was tested using two approaches. The first approach was a quantitative analysis by using the Relative Operating Characteristic (ROC) method. The second approach was a qualitative analysis in which the probable urban growth maps based on urban changes in 1985 is used to test the performance of the model to predict the probable urban growth after 2014 by comparing the probable maps of 1985 and the actual urban growth of 2014. The results indicate that the prediction model of 2014 provides a reliable and consistent prediction based on the performance of 1985. The analysis of driving forces shows variable effects over time. Variables such as elevation, slope and road distance had significant effects on the selected cities. However, distance to major roads was the factor with the most impact to determine the urban form in all five cites in both 1985 and 2014.

1. Introduction

The urban system is widely recognized as a complex process and is very dynamic due to its rapid change over time. Rapid urbanization is linked to population growth and human activities, which have become more frequent in recent years compared to the past [1,2]. Since most of the global population currently lives in cities, the nature of the cities contains a complex human settlement structure compared to rural portions [3]. The population size is expected to increase in the future (i.e., greater than 5 billion by 2030) [4,5]. Most of this growth will occur in developing countries [6], transforming natural ecosystems into urban areas. The increase in urban areas has massive environmental and ecological consequences, including landscape change, hydrological process altering and regional and global climate change [7,8].
Recently, Saudi Arabia, like other developing countries, has faced massive urban growth due to an increase in national economics [9,10], which have been used to support the public and private sectors to increase development across the country. This factor, along with natural demographic growth, has led to massive changes and increased complexity in the structure of cities as well as unsystematic development in the last three decades. According to the Central Department of Statistics and Information 2015, the population in large cities, such as Riyadh and Jeddah, has increased by 4,682,200 and 3,024,368 in Riyadh and Jeddah respectively between 1985 and 2014. With such rapid development, local authorities and planners are challenged to understand the urban growth process and associated consequences [11]. The ability to manage urban growth lies in determining the past and present spatial dimensions and the driving forces behind the distribution of urbanization. Understanding the urban growth process and its driving forces is essential to secure sustainable development for future urban growth and to avoid environmental and ecological consequences.
Urban expansion is controlled by factors or variables that strongly affect urban distribution and forms [12]. These variables include biophysical variables, such as climate and topography; proximity variables, such as distance to water, distance to the central business district and distance to roads or highways; socioeconomic variables, such as population and gross domestic product; and neighborhood variables, such as the density of agricultural land, water bodies and forestland. Several studies [13,14,15,16,17,18,19] have emphasized the effects of these variables in describing urban patterns and to determine the probable urban growth in the future. There has been a considerable effort to identify the effect of driving variables of urban distribution. However, it is important to consider that driving factors may change over time due to the complexity of urban dynamics. Understanding the effect of these variables on urban growth is critical because their effect varies with time.
Geographic Information Systems (GIS) and Remote Sensing provide the required data and tools to spatially detect and model the past, current and future urban growth at the local, regional and global scales. Modeling urban expansion is crucial to understanding the current and future processes based on the driving factors. A wide range of statistical models have been developed and examined to identify the direction of urban distribution and to determine the expansion process [20,21,22]. The statistical methods of urban expansion driving forces mainly include multiple linear regression [23], geographically weighted regression [21,24], bivariate regression [25], an autocorrelation model [21], and a logistic regression model [14,15,16,17,26,27]. The logistic regression model is able to determine the trend of urban expansion in the future based on driving variables that affect the distribution of urbanization. However, most previous studies have used driving force factors to predict urban growth probabilities based on a single date of urban change [16].
This paper aims to investigate the effect of four driving force factors, including elevation, slope, distance to drainages, and distance to major roads, to predict urban growth probabilities in five Saudi Arabian cities (Riyadh, Jeddah, Makkah, Al-Taif and Eastern Area) based on urban changes at two time periods: 1985 and 2014. This paper also examines the relative importance of the selected variables in determination of the shape and pattern of urban expansion in the selected cities. Finally, the paper discusses the environmental and ecological implications of urban growth in the selected cities and the expected impacts.

2. Materials and Methods

2.1. Study Area

The selected study areas include five cities in Saudi Arabia: Riyadh, Jeddah, Makkah, Al-Taif and Eastern Area (Figure 1). These cities are considered the most urbanized and populated in the country, with Riyadh being the capital and largest city. This city is situated in the central part of Saudi Arabia on the large Najd plateau and, according to 2013 census statistics, has a population of 6,079,295. This area includes the Al-Kharj Governorate. Jeddah is the largest sea port on the Red Sea coast, a major urban center of western Saudi Arabia and a nationally important commercial hub. With a population of 3,865,873, according to 2013 census, Jeddah is the largest city in the Makkah Region and the second largest in the country. Makkah is the holy place of the Muslim community and is located in the central part of the region, approximately 70 km inland from Jeddah. In 2013, Makkah had a population of 1,867,886. Al-Taif, which is considered the most important tourist city of Saudi Arabia, is the fourth study area used in our analysis and is located in the south-eastern part of the Makkah region. Al-Taif has a population of 1,083,693. Eastern Area is located in eastern Saudi Arabia on the Arabian Gulf and is home to most of Saudi Arabia’s oil production. Eastern Area includes 11 suburbs and governorates that are spread across the region. Only the most populated six suburbs, which are Dammam, Dhahran, Al-Khobar, Al-Qatif, RasTanura and Al-Jubail, are considered in this study. The population of these towns totaled 2,713,583 in 2013. Dammam city is in the central part and is the administrative city of the Eastern Area.
All five cities lie within a very dry climate with average temperatures of 25 °C, 28 °C, 30 °C, 20 °C and 26 °C for Riyadh, Jeddah, Makkah, Al-Taif and Eastern Area, respectively. All five cities experience very little rainfall. The annual precipitation varies from city to city, but not by much, with 92 mm, 61 mm, 111 mm, 120 mm and 90 mm in Riyadh, Jeddah, Makkah, Al-Taif and Eastern Area, respectively. The topographic structure is also varied among cities. While Jeddah is located in the coastal plain of the Red Sea (<5 m), the eastern part of the city contains an elevated terrain (Hijaz Mountains). Eastern Area is relatively flat, which is located in the coastal plain of the Arabian Gulf. Makkah lies in a corridor between the Hijaz Mountains, and the east portion is the most elevated land in the city. Al-Taif is the highest city among the five. It is located at an elevation of 2000 m on the slopes of the Hijaz Mountains. Riyadh has varied elevations, in which the south-west contains elevated terrains; the eastern part has a lower elevation.
During the last 40 years, the five cities have experienced a high rate of urbanization. The socioeconomic benefits provided in these cities were attractive to people who are seeking a better job or life. The different characteristics of each city have led to increases in both internal and external migration. Riyadh is the capital city and is attractive from internal migration due to job demand in the government sector. The commercial advantages in Jeddah have contributed to increasing urban growth in this city. Oil industry developments in Eastern Area, the growth of tourism sector in Al-Taif and the religious circumstances in Makkah have rapidly urbanized these cities.

2.2. Quantifying the Spatial Distribution of Urban Expansion from 1985 to 2014

2.2.1. Data and Pre-Processing

The data that were used to quantify the rate of urban expansion in the five cities were Landsat images from 1985 and 2014 (Table 1). All of the images were obtained from USGS Global Visualization (GloVis) site. All of the images were georeferenced as level 1 products; however, the TM images path 164 and row 42 and path 165 and row 43 of Eastern Area and Riyadh, respectively, were not correctly georeferenced. An automatic image-to-image registration was applied using Landsat 8 images of the similar path and row as base images and Landsat TM images as wrap images. Mosaic processing was applied to merge two path and row images for Riyadh and Eastern Area. Finally, image subsets were extracted for the five cities, including all of the urban boundaries across the five cities. Subsets of the selected areas of the Makkah and Al-Taif images did not include the cloud cover in the scenes.

2.2.2. Image Classification and Accuracy Assessment

An automated hierarchical classification was applied separately for the ten images using a decision tree classifier. The design of the decision tree classifier was based on classes (land cover classes), feature sets used to derive these classes, and threshold values applied to distinguish certain categories in the Landsat imagery. An object-based image analysis was used in eCognition Developer 8.9 to produce the land cover maps of 1985 and 2014 for the selected five cities. Following the method proposed by Taubenböck, et al. [28], the image classification was divided into four categories, including water, vegetation, bare soil and urban land cover using decision tree method. The feature sets selected in this study to derive land cover classes were the spectral patterns of nine different features. Six features stemmed from the Landsat TM and OLI spectral bands. The remaining three features were Normalized Difference Vegetation Index (NDVI), Soil Adjusted Vegetation Index (SAVI), and Normalized Difference Water Index (NDWI). Since the aim of this paper was to analyze the urban process, the other three categories of water, vegetation and bare soil were reclassified as non-urban areas. Thus, the final results presented two categories, including urban and non-urban areas, which were then used in the statistical analysis (Figure 2). The accuracy assessment of image classification was tested using stratified random sampling method. The results indicated a high overall accuracy (average overall accuracies were 89.15% and 90.84% in 1985 and 2014 respectively) and kappa coefficient (average kappa coefficients were 0.80 and 0.84 in 1985 and 2014 respectively) among the 10 images.

2.3. Potential Driving Factors of Urban Expansion

Typically, four types of driving variables have been considered in the previous studies of urban expansion [17]. These types included the biophysical variables, socioeconomic variables, neighborhood variables, and proximity variables and showed considerable effects on urban expansion in other regions. As no data was available for other socioeconomic and neighborhood variables, four different independent variables were selected to examine the most effective variable in the determination of urban distribution in the selected five cities (Figure 3). Two of these driving factors are biophysical variables and two are proximity variables. The two biophysical variables include the topographic factors of elevation and slope and the biophysical factors are the fundamental and natural drivers of the spatial expansion of urban land [17]. Most previous studies (e.g., [15,16,17]) have emphasized that biophysical factors, including topographical elements, strongly influence the distribution of urban form and control the shapes of cities. In this research, digital elevation model (DEM) layers were obtained from the USGS site for the five cities with a spatial resolution of 30 meters (SRTM data: http://earthexplorer.usgs.gov/). Elevation and slope were evaluated using DEM layers. The slope was calculated as a percentage, dividing the rise by the run and multiplying by 100.
Another two driving factors were proximity variables, which included the distance from drainage network and distance to the major roads. Drainages were delineated from the DEM layers and compared to the topographic map of Saudi Arabia in order to determine the major drainage order. The Euclidean distance algorithm was used to calculate the drainage distance. The second factor is the distance to the major roads. Road networks strongly shape the urban form [16] and play an important role in urban development [26,29]. The proximity of the major roads was measured as the straight-line distance to the major roads in the five cities. Roads were evaluated for the five cities in both 1985 and 2014 using heads-up digitizing in ArcGIS 10.2 for the Landsat TM of 1985 and Landsat OLI panchromatic bands for 2014. A visual comparison to a topographic map of Saudi Arabia was applied for both time periods in order to identify the major roads in the five cities. The Euclidean distance algorithm was used to calculate the street distance.

2.4. Statistical Analysis

2.4.1. Urban and Non-Urban Data Sampling

The data sampling was performed for urban and non-urban features for two reasons. First, the derived layers of urban expansion were very large, difficult to handle and time consuming. For example, it is difficult for computer systems to deal with a very large number of pixels for Riyadh (5290 × 5325) with a 30-meter resolution for dependent and independent variables [14,17,30,31]. Thus, the best way to handle such data in the later statistical analysis is to create a reliable feature set and systematic sampling. The second issue is to avoid spatial autocorrelation between both dependent and independent variables. Thus, we created random points for both urban and non-urban features, which coded 1 (presence) for urban and 0 (absence) for non-urban features. Since the extent of each city is different from each other, the total random points were different from city to city. For a large extent, such as Riyadh, Jeddah and Eastern Area, we created 30,000, 20,000 and 25,000 random samples, respectively. For a relatively small extent, such as Makkah and Al-Taif, 10,000 random sample points were created for both urban and non-urban features. However, samples of non-urban features were much larger than those of urban features. Therefore, we recreated sample points for non-urban features to be equal to the urban sample points.

2.4.2. Logistic Regression

Logistic regression is a mathematical model that can be used to model and describe the relationship between a binary dependent variable (Y) and several independent categorical and continuous variables (Xs) [32,33]. The logistic regression model is the most frequently used model in urban growth studies [34]. The advantage of the logistic regression model over other linear regression models in urban studies is its ability to estimate the urban probability of the occurrence of a discrete response factor [35]. As the computation process of logistic regression is less intensive than other urban prediction methods, it can empirically estimate the influence of independent variables on the dependent variable, and can provide a degree of confidence regarding their contribution [16]. In this research, Equation (1) was used to calculate the binary logistic regression model in the urban growth of the five cities.
log i t ( y ) log ( y 1 y ) = β 0 + β 1 χ 1 + ...... β η χ η + ε
where y is the probability of a cell being converted to urban area 1, x𝜂 is an independent variable or driving force, 𝛽0 is the intercept coefficient, 𝛽𝜂 is the coefficient of variable x𝜂, and ε is the random error term (residuals).
We used a binary logistic regression to associate urban growth with two biophysical factors (elevation and slope) and two proximity driving factors (distance to roads and distance to drainages) and to produce urban growth probability maps for the five cities. The binary logistic regression was used first in 1985 to estimate and compare the urban growth probability with urban maps of 2014 that were generated from satellite data. This regression was then applied for the urban growth of 2014 using Equation (1).

2.4.3. Prediction of Spatial Patterns of Urban Distribution

The probability maps for the five cities were derived using Equation (1). The coefficient values of the driving factors (elevation, slope, drainage distance, and road distance) and the intercept values were calculated in ArcGIS 10.2. The odds ratio of value 1 was used individually for each map of 1985 and 2014 across the five cities. The calculation was first applied to the 1985 data in order to validate the probability maps of 2014. A visual comparison was performed by overlaying the probability map of 1985 and the urban area maps of 2014 that were produced from the image classification of 2014 images. The logistic regression model was then applied to the 2014 data using the new road layer of 2014 for the street distance value.

2.4.4. Model Validation

The logistic regression model performance was validated using the Relative Operating Characteristic (ROC) method. ROC is a widely used method to measure the relationship between simulated change and real change in the field of land use and land cover changes [16,36,37,38]. The statistical approach of the ROC method is the area under the curve that connects the plotted points. The equation that was used to calculate the area under the curve (AUC) was proposed by Pontius Jr and Schneider [39] as:
A U C = i = 1 n [ x i + 1 x i ] [ y i + 1 y i / 2 ]
where xi is the rate of false positive for i, yi is the rate of true positive for i, and n is the number of suitability groups. The rate of false positive xi and the rate of true positive yi were obtained from the logistic regression analysis (Equation (1)). Figure 4 illustrates in brief the methods used in this research.

3. Results

3.1. Spatial Statistical Modeling of Urban Expansion

The analysis of the logistic regression model, including the four independent variables (elevation, slope, drainage distance and street distance), successfully explained the variance of urban distribution in the five cities. Table 2, Table 3, Table 4, Table 5 and Table 6 list the statistical analysis results of the logistic regression model of four independent variables in Riyadh, Jeddah, Makkah, Al-Taif, and Eastern Area, respectively. The percent correct prediction (PCP) values with the cut-off point of 0.5 for Riyadh, Jeddah, Makkah, Al-Taif and Eastern Area were 85.2, 85.2, 71.6, 70.1, and 71.1, respectively, in 1985 and 83.8, 79.7, 73.6, 70.8, and 70.1, respectively, in 2014. The PCP values indicated a high percentage in Riyadh and Jeddah for both 1985 and 2014 (Table 2 and Table 3, respectively), while the PCP values of Makkah, Al-Taif, and Eastern Area were slightly lower. However, PCP values of 70 percent indicated a good prediction of the logistic regression model.
The significant effect of the selected variables showed differences among the five cities. While elevation positively affected the urban distribution in Riyadh in 2014, it had less effect in 1985. Slope also negatively affected the urban distribution in 2014; however, it had no effect in 1985. The distance to a road consistently influenced urban expansion in both 1985 and 2014 in Riyadh (Table 2). Similarity, Jeddah had a variety of influences in the selected variables (Table 3). Elevation positively affected urban distribution in 1985 and 2014. In addition, the distance to a road showed a negative influence for both 1985 and 2014. Slope had an effect in 2014 but showed no influence in 1985. The distance to drainage had an effect in 1985 in Makkah but showed no effect in 2014. The distance to a road also negatively affected the urban distribution in both 1985 and 2014 (Table 4). While all of the variables significantly affected the urban distribution in 1985 in Al-Taif, elevation and drainage distance had less effect in 2014 (Table 5). Similar to other variables, the road distance had a significant effect in Eastern Area for both 1985 and 2014, and slope had an effect in 2014 (Table 6). In summary, the significant effect of road distance on urban distribution was common for all five cities in 1985 and 2014, while other variables showed a variety of influences among the five cities in the selected time period.

3.2. Urban Probability among the Five Cities (1985)

To ideally test the effectiveness of the logistic regression model for predicting urban probability in 2014, we first ran the model using 1985 data and compared the obtained results with the urban areas of 2014 that were produced from satellite images. A visual comparison of the results of urban probability in 1985 and urban areas in 2014 is shown in Figure 5, Figure 6, Figure 7 and Figure 8 and Figure 9A for Riyadh, Jeddah, Makkah, Al-Taif and Eastern Area, respectively. The lighter tones indicate higher probabilities of urban expansion, and the darker tones indicate lower probabilities of urban growth. As presented in the figures, the road networks played a significant role in the urban probability in all five cities in 1985. Most of the urban development was close to major roads. However, the urban probability in Makkah and Al-Taif (Figure 7A and Figure 8A, respectively) was affected by elevation and road distance. The comparison results between 1985 and 2014 show an effective model of urban probabilities in the selected cities.

3.3. Urban Probability among the Five Cities (2014)

Figure 5, Figure 6, Figure 7 and Figure 8 and Figure 9B show the probabilities of urban growth in Riyadh, Jeddah, Makkah, Al-Taif and Eastern Area, respectively in 2014. The lighter tones indicate the higher probabilities of urban expansion, and the darker tones indicate lower probabilities of growth. The effect of biophysical factors in Riyadh and Eastern Area is relatively low, except for a small effect of slope in the southern portion of Riyadh (Figure 5B and Figure 9B, respectively). Thus, the expansion of urban growth in Riyadh is likely to expand to the north and north-west more than to the south. The distribution of urbanization in the relatively flat areas in the northern, eastern and western portions of Riyadh and in Eastern Area coincide with the growth of major roads, which play a more important role in the expansion of growth compared to biophysical factors. The increase in road networks in Eastern Area between 1985 and 2014 permitted growth expansion in the north and south (Figure 9B). The probabilities of urban growth in Jeddah are likely to continue north and south across the costal line. In addition, there is a probable development in the east portions of the city, but growth is likely controlled by elevation and slope (Figure 6B). The probabilities of urban development in Makkah are mostly in all directions. The slope factor can control the growth in the city (Figure 7B). Finally, most probable areas of urban growth in Al-Taif are located in the northern, eastern and north-eastern portions of the city, which have a relatively low topographic structure. In contrast, the probabilities of urban growth in the western and southern portions are controlled by slope and elevation variables (Figure 8B). The future distributions of urban areas in the five cities are easily disenable from these maps.

3.4. Model Validation Using the Relative Operating Characteristics (ROC) Method

Figure 10 shows the ROC curves in both 1985 (red) and 2014 (blue) in (A) Riyadh; (B) Jeddah; (C) Makkah; (D) Al-Taif; and (E) Eastern Area. The curves of the logistic regression model vary from city to city and from one date to another. While this model clearly reaches a stable level in cities such as Riyadh and Jeddah in both 1985 and 2014 (Figure 10A, B, respectively), with AUC values of 0.90 and 0.88, respectively, in 1985 and 0.88 and 0.87, respectively in 2014, it shows unstable performance for the other three cities of Makkah, Al-Taif, and Eastern Area (Figure 10C–E, respectively). However, the AUC values are still fair in these cities, at 0.79, 0.77, and 0.78, respectively, in 1985 and 0.78, 0.77, and 0.76, respectively, in 2014 for Makkah, Al-Taif, and Eastern Area. The performance of the model is reasonable when AUC values are greater than random (0.50) [39].

3.5. The Importance of Predictor Variables

Figure 11 shows the impact of the four variables in the selected five cities. The impact of the selected variables on the urban distribution varied over time. However, it is obvious that the distance to a road was the most important predictor in both 1985 and 2014 in the five cities. This value is close to 1 in Eastern Area in 1985 and 2014, in Riyadh in 1985 and 2014 and in Makkah in 2014. Elevation is the second most important predictor in Jeddah in 1985 and 2014. Additionally, slope is the second most important predictor in Al-Taif in 1985 and 2014 and Makkah in 1985. The least important factor is drainage distance, which shows the least effect among the four variables.

4. Discussion

The results that are presented in this research indicate that road distance is the most important variable among the four selected factors in the five cities for the projection of urban growth. The analysis of data from both 1985 and 2014 showed a domination of the road distance factor because of the extensive transportation systems in Saudi Arabia that have been built in the five decades following 1950 [40]. The transportation systems in the five cities significantly increased the probability of urban expansion and will likely contribute to the expansion of urbanization in the future. The significant effects of road distance indicate the probability of urban areas being close to major roads. Thus, locations close to roads are likely to be developed more than are locations that are far from the roads. The results of this study are consistent with the results presented by Müller, et al. [41] and Li, Zhou and Ouyang [17], who reported the relationship between the distance to roads and the probability of urban expansion.
Biophysical variables, such as elevation and slope, showed variable effects among the five cities. While the elevation effects were lower in Riyadh in 1985, they were significantly positive in 2014. This result may be due to the increase in housing demands in recent years [42], which prompted the government and housing agencies to build new construction without considering practical elevation of the land form. In addition, staying within the urban boundary close to the city center rather than in the margins is on demand, even though the construction cost is high. Jeddah was also affected negatively by elevation in both 1985 and 2014, but the effect in 2014 was greater than that in 1985 because urban growth was distributed in the east where the elevation is high. Thus, the elevation factors control the urban growth in the eastern portion of the city more than in the other parts. While Al-Taif was positively affected by elevation at a level of 0.05 in 1985, elevation had less effect in 2014. Other cities, such as Makkah and Eastern Area, were less affected by elevation in both 1985 and 2014.
Slope also had a significant effect in the five cities. All five cities showed a negative effect of slope, indicating that slope limited the distribution of urban expansion in these cities. However, there were variable effects between 1985 and 2014. While the slope negatively influenced urban expansion in Al-Taif in both 1985 and 2014 due to the complex topographic structure in Al-Taif, the effect of slope in cities such as Riyadh, Jeddah and Makkah was significant in 2014. The urban areas in Riyadh, Jeddah and Makkah expanded after 1985 from the center towards the outer boundary, which potentially contain steep slopes. For example, the analysis of urbanization in 1985 in Riyadh showed that the city expansion was limited and centered near the city center, while in 2014, the urban area expanded to reach the slopes, as is the case in the southern part of the city. Similarly, urban expansion in Jeddah in 2014 was distributed towards the eastern portion, which has relatively steep slopes. The influence of slope will continue to affect these cities in the future because urban expansion is likely to continue in the same direction. Thus, slope is an important factor in determining the distribution of urban areas. Other studies in different places have observed the effect of slope on urban expansion [15,17,43].
While other studies indicated that distance to water has an effect on the urban growth, our results showed the opposite. Distance to drainage showed least effects on urban distribution among the selected variables. This is because of the climate condition in the desert areas. In such areas, water surfaces, except seas, are fed by rain, which is very little in these desert cities. Most of the water surfaces are dried out throughout the year due to rising temperature and evaporation. This condition also applies to the drainage lines. In Saudi Arabia, drainage lines are not permanent waterways. They are normally dry most of the year and fill after a heavy rain.
The analysis also indicates that the selected variables show variable effects over time. However, there is an exception in the distance to roads, which consistently showed an effect in both 1985 and 2014 in all five cities. The accessibility to urban facilities and services and security may be attractive for construction and expansion close to roads. The road distance factor may continue to influence urban growth in the selected cities and to determine the shape and spatial dimension of urban distribution in the future. The influence of biophysical factors, such as elevation and slope, may decrease in the future due to the advancement of technology that can reduce the cost of construction in these locations [17,43] as well as limited low-slope land being available. Thus, driving forces are variable over the time due to the increased urbanization rate.
The logistic regression model has the ability to provide an easily discernible prediction of urban growth [16]. However, the ideal approach to test the reliability of the prediction model is to develop a model based on urban change over one-time period and to apply the model to predict the probabilities of growth after the second time period [44]. Applying the model to predict urban growth in 1985 offered us an opportunity to test the ability of the model to predict urbanization probabilities based on 2014 urban growth. The results of urbanization probability maps based on 1985 urban change and the four driving forces showed higher predictors compared to the actual urban growth in 2014. This result means that the model that was applied to predict the probable areas of urbanization after 2014 provides consistent and accurate probabilities of urban growth. Therefore, multi-temporal datasets offer a technique to compare actual growth with a simulation model of the potential development of urbanization.
While logistic regression model is powerful for analyzing the future urban development based on the trends that were observed in the past [15,45], the limitations of the logistic regression model in the urban change prediction can be summarized in three points. First, the logistic regression method is not able to provide a transition matrix from different land cover areas to urban areas. Moreover, it may not be suitable to provide a detailed description of the temporal dynamics of the future urban growth and to address the correlation between independent variables [46]. Finally, the prediction model of logistic regression ignores the effect of other land cover areas on the expected urban growth which may produce an error of the estimation of the future development. However, logistic regression model is useful to account for the dynamics of the past and future urban growth based on different driving factors.
Modification of the landscape in the five cities to expand construction has led to the loss of several environmental and ecological components, including vegetation cover, biodiversity and soil. The landscape problem in Saudi Arabia lies in that most urban areas have been built near or on sensitive environments, such as mountains, which increases the developmental pressure in these areas and decreases the available green spaces and associated biodiversity. Another issue is the increased karst collapse hazard. The increased karst collapse hazard, particularly in Eastern Area, is due to rapid urban development, which has produced larger buildings and much stronger foundations [47]. The change in the landscape structure and landform of these cities to expand urban areas suggests an increase in the danger of such hazards. In addition to the karst collapse hazard, the higher level of urbanization and impervious surfaces has a significant impact on the hydrological systems and increases the potential hazard of flood risks [48]. Recently, floods have become an annual disaster with a high impact in many places in Saudi Arabia [49]. Thus, the measurement of urban expansion is essential for managing and assessing environmental and ecological impacts in order to avoid an increase in potential risks and to protect such sensitive environments. The analysis of urban change prediction provides a tool to manage and assess the environmental changes in the five cities and provides the needed information that can help planners and local authorities in possible future development.

5. Conclusions

Information on urbanization processes and their deriving factors should help planners to provide adequate strategies for further developments. Urban planners and environmental managers need an answer to the question as to which factors have an influence on the expansion of urban growth. Failure to understand the urbanization process and its driving factors in the past and present growth will result in unexpected expansion in the future. This research tried to address this question by analyzing urban expansion in five Saudi Arabian cities between 1985 and 2014 using satellite data and a logistic regression model.
This study emphasized that incorporation of satellite data and regression statistical analysis (e.g., logistic regression) is effective in improving our understanding of the trend and distribution of the future urban growth in Saudi Arabian cities. It also provides the essential information of the most affected variables on urban expansion in the past and future changes. The results indicate that the roads variable had the most effect on urban developments among the five cities. Urban cover in Saudi Arabian cities is likely to develop near and on the major roads. Slope and elevation considerably controlled urban growth during the past growth in some cities, but these effects are likely to be reduced in the future, especially in the low-slope lands, due to advancement in technology to build on sloping land and the scarcity of flat land.
In this study, we analyzed only four variables, including elevation, slope, drainage distance and street distance, and the main effects of such variables on urban distribution in five Saudi Arabian cities. Other variables, such as population, gross domestic product, distance to the socioeconomic centers, land use policy and master plan may have an important relationship with urban distribution in the selected cities. However, we did not include such variables due to limitations in the availability of these data. Incorporating such variables with the selected factors can provide more refinement of the urbanization process in Saudi Arabian cities. In addition, integrating logistic regression model and other regression models, such as geographically weighted regression (GWR), will be worthwhile to determine the heterogeneity in the estimated relationships among the independent and dependent variables [50,51]. In addition, incorporating logistic regression and cellular automata model may provide more detail regarding the future transition from other land cover areas to urban land cover.

Author Contributions

Abdullah Alqurashi, Lalit Kumar and Khalid Al-Ghamdi conceived and designed the experiments; Abdullah Alqurashi performed the experiments; Abdullah Alqurashi and Lalit Kumar analyzed the data; Khalid Al-Ghamdi contributed reagents/materials/analysis tools; Abdullah Alqurashi wrote the paper. Lalit Kumar and Khalid Al-Ghamdi revised the manuscript. All authors equally contributed to this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Masek, J.G.; Lindsay, F.E.; Goward, S.N. Dynamics of urban growth in the washington dc metropolitan area, 1973–1996, from landsat observations. Int. J. Remote Sens. 2000, 21, 3473–3486. [Google Scholar]
  2. Xu, H. Extraction of urban build-up land features from landsat imagery using a thematicoriented index combination technique. Photogramm. Eng. Remote. Sens. 2007, 73, 1381–1391. [Google Scholar] [CrossRef]
  3. Weng, Q. What is special about global urban remote sensing? In Global Urban. Monitoring and Assessment through Earth Observation; CRC Press: Boca Raton, FL, USA, 2014; pp. 1–12. [Google Scholar]
  4. Seto, K. Global urban issues. In Global Mapping of Human Settlement; CRC Press: Boca Raton, FL, USA, 2009; pp. 3–9. [Google Scholar]
  5. Rahman, A.; Aggarwal, S.P.; Netzband, M.; Fazal, S. Monitoring urban sprawl using remote sensing and gis techniques of a fast growing urban centre, india. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2011, 4, 56–64. [Google Scholar] [CrossRef]
  6. Cohen, B. Urban growth in developing countries: A review of current trends and a caution regarding existing forecasts. World Dev. 2004, 32, 23–51. [Google Scholar] [CrossRef]
  7. Grimmond, S.U.E. Urbanization and global environmental change: Local effects of urban warming. Geogr. J. 2007, 173, 83–88. [Google Scholar] [CrossRef]
  8. Hoymann, J.; Goetzke, R. Simulation and evaluation of urban growth for germany including climate change mitigation and adaptation measures. ISPRS Int. J. Geo-Inf. 2016, 5. [Google Scholar] [CrossRef]
  9. Al-Ahmadi, K.; See, L.; Heppenstall, A.; Hogg, J. Calibration of a fuzzy cellular automata model of urban dynamics in saudi arabia. Ecol. Complex. 2009, 6, 80–101. [Google Scholar] [CrossRef]
  10. Alqurashi, A.; Kumar, L. Land use and land cover change detection in the saudi arabian desert cities of makkah and al-taif using satellite data. Adv. Remote Sens. 2014, 3, 106–119. [Google Scholar] [CrossRef]
  11. Aljoufie, M. Toward integrated land use and transport planning in fast-growing cities: The case of jeddah, saudi arabia. Habitat Int. 2014, 41, 205–215. [Google Scholar] [CrossRef]
  12. Kamusoko, C.; Gamba, J. Simulating urban growth using a random forest-cellular automata (rf-ca) model. ISPRS Int. J. Geo-Inf. 2015, 4. [Google Scholar] [CrossRef]
  13. Batisani, N.; Yarnal, B. Urban expansion in centre county, pennsylvania: Spatial dynamics and landscape transformations. Appl. Geogr. 2009, 29, 235–249. [Google Scholar] [CrossRef]
  14. Cheng, J.; Masser, I. Urban growth pattern modeling: A case study of wuhan city, pr china. Landsc. Urban Plan. 2003, 62, 199–217. [Google Scholar] [CrossRef]
  15. Dubovyk, O.; Sliuzas, R.; Flacke, J. Spatio-temporal modelling of informal settlement development in sancaktepe district, istanbul, turkey. ISPRS J. Photogramm. Remote Sens. 2011, 66, 235–246. [Google Scholar] [CrossRef]
  16. Hu, Z.; Lo, C.P. Modeling urban growth in atlanta using logistic regression. Comput. Environ. Urban Syst. 2007, 31, 667–688. [Google Scholar] [CrossRef]
  17. Li, X.; Zhou, W.; Ouyang, Z. Forty years of urban expansion in beijing: What is the relative importance of physical, socioeconomic, and neighborhood factors? Appl. Geogr. 2013, 38, 1–10. [Google Scholar] [CrossRef]
  18. Liu, H.; Zhou, Q. Developing urban growth predictions from spatial indicators based on multi-temporal images. Comput. Environ. Urban Syst. 2005, 29, 580–594. [Google Scholar] [CrossRef]
  19. Zeng, C.; Zhang, M.; Cui, J.; He, S. Monitoring and modeling urban expansion—A spatially explicit and multi-scale perspective. Cities 2015, 43, 92–103. [Google Scholar] [CrossRef]
  20. Arsanjani, J.J.; Helbich, M.; Mousivand, A.J. A morphological approach to predicting urban expansion. Trans. GIS 2014, 18, 219–233. [Google Scholar] [CrossRef]
  21. Shafizadeh-Moghadam, H.; Helbich, M. Spatiotemporal variability of urban growth factors: A global and local perspective on the megacity of mumbai. Int. J. Appl. Earth Obs. Geoinform. 2015, 35, Part B. 187–198. [Google Scholar] [CrossRef]
  22. Ahmed, B.; Ahmed, R. Modeling urban land cover growth dynamics using multi‑temporal satellite images: A case study of dhaka, bangladesh. ISPRS Int. J. Geo-Inf. 2012, 1. [Google Scholar] [CrossRef]
  23. Maniquiz, M.C.; Lee, S.; Kim, L.-H. Multiple linear regression models of urban runoff pollutant load and event mean concentration considering rainfall variables. J. Environ. Sci. 2010, 22, 946–952. [Google Scholar] [CrossRef]
  24. Crespo, R.; Grêt-Regamey, A. Spatially explicit inverse modeling for urban planning. Appl. Geogr. 2012, 34, 47–56. [Google Scholar] [CrossRef]
  25. Wu, K.-y.; Zhang, H. Land use dynamics, built-up land expansion patterns, and driving forces analysis of the fast-growing hangzhou metropolitan area, eastern china (1978–2008). Appl. Geogr. 2012, 34, 137–145. [Google Scholar] [CrossRef]
  26. Reilly, M.K.; O’Mara, M.P.; Seto, K.C. From bangalore to the bay area: Comparing transportation and activity accessibility as drivers of urban growth. Landsc. Urban Plan. 2009, 92, 24–33. [Google Scholar] [CrossRef]
  27. Tsutsumida, N.; Saizen, I.; Matsuoka, M.; Ishii, R. Addressing urban expansion using feature-oriented spatial data in a peripheral area of ulaanbaatar, mongolia. Habitat Int. 2015, 47, 196–204. [Google Scholar] [CrossRef] [Green Version]
  28. Taubenböck, H.; Esch, T.; Felbier, A.; Wiesner, M.; Roth, A.; Dech, S. Monitoring urbanization in mega cities from space. Remote Sens. Environ. 2012, 117, 162–176. [Google Scholar] [CrossRef]
  29. Meyer, M.D.; Miller, E.J. Urban Transportation Planning: A Decision-Oriented Approach; McGraw-Hill: New York, NY, USA, 2001. [Google Scholar]
  30. Crk, T.; Uriarte, M.; Corsi, F.; Flynn, D. Forest recovery in a tropical landscape: What is the relative importance of biophysical, socioeconomic, and landscape variables? Landsc. Ecol. 2009, 24, 629–642. [Google Scholar] [CrossRef]
  31. Luo, J.; Wei, Y.H.D. Modeling spatial variations of urban growth patterns in Chinese cities: The case of Nanjing. Landsc. Urban Plan. 2009, 91, 51–64. [Google Scholar] [CrossRef]
  32. Kleinbaum, D.G.; Klein, M. Logistic Regression: A Self-Learning Text, 3rd ed.; Springer Science & Business Media: New York, NY, USA, 2010. [Google Scholar]
  33. McCullagh, P.; Nelder, J.A.; McCullagh, P. Generalized Linear Models; Chapman and Hall London: Boca Raton, FL, USA, 1989. [Google Scholar]
  34. Liao, J.; Tang, L.; Shao, G.; Su, X.; Chen, D.; Xu, T. Incorporation of extended neighborhood mechanisms and its impact on urban land-use cellular automata simulations. Environ. Model. Softw. 2016, 75, 163–175. [Google Scholar] [CrossRef]
  35. Rienow, A.; Goetzke, R. Supporting sleuth—Enhancing a cellular automaton with support vector machines for urban growth modeling. Computers, Environ. Urban Syst. 2015, 49, 66–81. [Google Scholar] [CrossRef]
  36. Pontius, R.G. Quantification error versus location error in comparison of categorical maps. Photogramm. Eng. Remote Sens. 2000, 66, 1011–1016. [Google Scholar]
  37. Mas, J.-F.; Soares Filho, B.; Pontius, R.G.; Farfán Gutiérrez, M.; Rodrigues, H. A suite of tools for roc analysis of spatial models. ISPRS Int. J. Geo-Inf. 2013, 2, 869–887. [Google Scholar] [CrossRef]
  38. Pontius, R., Jr.; Parmentier, B. Recommendations for using the relative operating characteristic (roc). Landsc. Ecol. 2014, 29, 367–382. [Google Scholar] [CrossRef]
  39. Pontius, R.G., Jr.; Schneider, L.C. Land-cover change model validation by an roc method for the ipswich watershed, massachusetts, USA. Agric. Ecosyst. Environ. 2001, 85, 239–248. [Google Scholar] [CrossRef]
  40. Dincer, I.; Hussain, M.M.; Al-Zaharnah, I. Energy and exergy utilization in transportation sector of Saudi Arabia. Appl. Therm. Eng. 2004, 24, 525–538. [Google Scholar] [CrossRef]
  41. Müller, K.; Steinmeier, C.; Küchler, M. Urban growth along motorways in switzerland. Landsc. Urban Plan. 2010, 98, 3–12. [Google Scholar] [CrossRef]
  42. Assaf, S.A.; Bubshaitr, A.A.; Al-Muwasheer, F. Factors affecting affordable housing cost in Saudi Arabia. Int. J. Hous. Mark. Anal. 2010, 3, 290–307. [Google Scholar] [CrossRef]
  43. Ye, Y.; Zhang, H.; Liu, K.; Wu, Q. Research on the influence of site factors on the expansion of construction land in the pearl river delta, china: By using gis and remote sensing. Int. J. Appl. Earth Obs. Geoinf. 2013, 21, 366–373. [Google Scholar] [CrossRef]
  44. Aspinall, R. Modelling land use change with generalized linear models—A multi-model analysis of change between 1860 and 2000 in gallatin valley, montana. J. Environ. Manag. 2004, 72, 91–103. [Google Scholar] [CrossRef] [PubMed]
  45. Triantakonstantis, D.; Mountrakis, G. Urban growth prediction: A review of computational models and human perceptions. J. Geogr. Inf. Syst. 2012, 4, 555–587. [Google Scholar] [CrossRef]
  46. Dendoncker, N.; Rounsevell, M.; Bogaert, P. Spatial analysis and modelling of land use distributions in belgium. Comput. Environ. Urban Syst. 2007, 31, 188–205. [Google Scholar] [CrossRef]
  47. Vincent, P. Saudi Arabia: An Environmental Overview; CRC Press: Boca Raton, FL, USA, 2008. [Google Scholar]
  48. Lee, J.; Heaney, J. Estimation of urban imperviousness and its impacts on storm water systems. J. Water Resour. Plan. Manag. 2003, 129, 419–426. [Google Scholar] [CrossRef]
  49. Al Saud, M. Assessment of flood hazard of jeddah area 2009, Saudi Arabia. J. Water Resour. Prot. 2010, 2, 839–847. [Google Scholar] [CrossRef]
  50. Fotheringham, A.S.; Brunsdon, C.; Charlton, M. Geographically Weighted Regression: The Analysis of Spatially Varying Relationships; John Wiley & Sons: Hoboken, NJ, USA, 2003. [Google Scholar]
  51. Long, Y.; Gu, Y.; Han, H. Spatiotemporal heterogeneity of urban planning implementation effectiveness: Evidence from five urban master plans of Beijing. Landsc. Urban Plan. 2012, 108, 103–111. [Google Scholar] [CrossRef]
Figure 1. Geographical locations of the selected cities in Saudi Arabia.
Figure 1. Geographical locations of the selected cities in Saudi Arabia.
Ijgi 05 00139 g001
Figure 2. Urban areas in 1985 in grey and the expansion of urban areas from 1985 to 2014 in black in (A) Riyadh; (B) Jeddah; (C) Makkah; (D) Al-Taif and (E) Eastern Area.
Figure 2. Urban areas in 1985 in grey and the expansion of urban areas from 1985 to 2014 in black in (A) Riyadh; (B) Jeddah; (C) Makkah; (D) Al-Taif and (E) Eastern Area.
Ijgi 05 00139 g002
Figure 3. Driving factors used in this research.
Figure 3. Driving factors used in this research.
Ijgi 05 00139 g003
Figure 4. Flowchart of the research methodology.
Figure 4. Flowchart of the research methodology.
Ijgi 05 00139 g004
Figure 5. Urbanization probability maps in Riyadh (A) a comparison of the urban probability in 1985 and urban areas in 2014; (B) urban probability in 2014.
Figure 5. Urbanization probability maps in Riyadh (A) a comparison of the urban probability in 1985 and urban areas in 2014; (B) urban probability in 2014.
Ijgi 05 00139 g005
Figure 6. Urbanization probability maps in Jeddah (A) a comparison of the urban probability in 1985 and urban areas in 2014; (B) urban probability in 2014.
Figure 6. Urbanization probability maps in Jeddah (A) a comparison of the urban probability in 1985 and urban areas in 2014; (B) urban probability in 2014.
Ijgi 05 00139 g006
Figure 7. Urbanization probability maps in Makkah (A) a comparison of the urban probability in 1985 and urban areas in 2014; (B) urban probability in 2014.
Figure 7. Urbanization probability maps in Makkah (A) a comparison of the urban probability in 1985 and urban areas in 2014; (B) urban probability in 2014.
Ijgi 05 00139 g007
Figure 8. Urbanization probability maps in Al-Taif (A) a comparison of the urban probability in 1985 and urban areas in 2014; (B) urban probability in 2014.
Figure 8. Urbanization probability maps in Al-Taif (A) a comparison of the urban probability in 1985 and urban areas in 2014; (B) urban probability in 2014.
Ijgi 05 00139 g008
Figure 9. Urbanization probability maps in Eastern Area (A) a comparison of the urban probability in 1985 and urban areas in 2014; (B) urban probability in 2014.
Figure 9. Urbanization probability maps in Eastern Area (A) a comparison of the urban probability in 1985 and urban areas in 2014; (B) urban probability in 2014.
Ijgi 05 00139 g009
Figure 10. ROC curve of 1985 (red) and 2014 (blue) in (A) Riyadh; (B) Jeddah; (C) Makkah; (D) Al-Taif and (E) Eastern Area.
Figure 10. ROC curve of 1985 (red) and 2014 (blue) in (A) Riyadh; (B) Jeddah; (C) Makkah; (D) Al-Taif and (E) Eastern Area.
Ijgi 05 00139 g010
Figure 11. The impact of the various predictors of the driving factors across the five cities.
Figure 11. The impact of the various predictors of the driving factors across the five cities.
Ijgi 05 00139 g011
Table 1. Landsat images that were used to produce urban and non-urban areas among the five cities.
Table 1. Landsat images that were used to produce urban and non-urban areas among the five cities.
NoPathRowSensorDateCloud Cover (%)Covered Area
116441TM4-April-19851Eastern Area
2OLI2-August-20140
316442TM4-April-19850
4OLI2-August-140
516543TM13-January-19850Riyadh
6OLI20-January-20140
716643TM14-August-19850
8OLI16-August-20140
916945TM30-March-19858Makah and Al-Taif
10OLI4-May-20143
1117045TM9-June-19850Jeddah
12OLI9-June-20140
Table 2. Estimated results of the logistic regression model in Riyadh of the four independent variables.
Table 2. Estimated results of the logistic regression model in Riyadh of the four independent variables.
Variables1985 2014
βS.E.p (>|z|)Exp (β)βS.E.p (>|z|)Exp (β)
Intercept2.8603 **1.03860.005917.46710.11640.38420.76191.1235
Elevation−0.00030.00180.84900.99970.0032 ***0.00070.00001.0031
Slope−0.00060.02210.97860.9994−0.0332 ***0.00870.00010.9674
Drainages distance−0.0001.0.00000.08840.99990.0000.0.00000.09651.0000
Street distance−0.0014 ***0.00010.00000.9986−0.0012 ***0.00010.00000.9988
PCP85.1911 83.8340
Neg. R20.68 0.61
2 Log Likelihood424.5 1966.0
Significant codes: *** p < 0.001, ** p < 0.01, . p < 0.1.
Table 3. Estimated results of the logistic regression model in Jeddah of the four independent variables.
Table 3. Estimated results of the logistic regression model in Jeddah of the four independent variables.
Variables19852014
βS.E.p (>|z|)Exp (β)βS.E.p (>|z|)Exp (β)
Intercept2.2401 ***0.25350.00009.39462.4729660.1163072.5E−10011.85756
Elevation−0.0212 ***0.00410.00000.9790−0.0153 ***0.0014891.22E−240.984854
Slope0.01490.02260.51111.0150−0.0108 ***0.0091360.2368260.989251
Drainages distance0.00000.00010.86341.0000−0.00013E−05.0.0002070.999889
Street distance−0.0016 ***0.00020.00000.9984−0.0015 ***8.65E−055.01E−710.99846
PCP85.2140 79.69256
Neg. R20.63 0.54
2 Log Likelihood385.5 2150.3
Significant codes: *** p < 0.001, . p < 0.1.
Table 4. Estimated results of the logistic regression model in Makkah of the four independent variables.
Table 4. Estimated results of the logistic regression model in Makkah of the four independent variables.
Variables19852014
βS.E.p (>|z|)Exp (β)βS.E.p (>|z|)Exp (β)
Intercept1.3834 ***0.19390.00003.98851.5837 ***0.27100.00004.8727
Elevation0.00010.00060.86921.0001−0.00020.00090.83660.9998
Slope−0.00440.00430.30070.9956−0.0523 ***0.00740.00000.9491
Drainages distance−0.0002 ***0.00000.00010.99980.00000.00010.69461.0000
Street distance−0.0012 ***0.00010.00000.9988−0.0010 ***0.00010.00000.9990
PCP71.6411 73.6136
Neg. R20.32 0.34
2 Log Likelihood2459.8 1216.6
Significant codes: *** p < 0.001.
Table 5. Estimated results of the logistic regression model in Al-Taif of the four independent variables.
Table 5. Estimated results of the logistic regression model in Al-Taif of the four independent variables.
Variables19852014
βS.E.p (>|z|)Exp (β)βS.E.p (>|z|)Exp (β)
Intercept−1.09691.05570.29880.33392.6620 ***0.73020.000314.3245
Elevation0.0017 *0.00070.01081.0017−0.0008.0.00050.08360.9992
Slope−0.0830 ***0.01350.00000.9204−0.0781 ***0.01010.00000.9249
Drainages distance−0.0004 ***0.00010.00010.99960.00000.00010.99341.0000
Street distance−0.0004 ***0.00000.00000.9996−0.0008 ***0.00010.00000.9992
PCP70.0730 70.8333
Neg. R20.29 0.29
2 Log Likelihood935.4 2024.3
Significant codes: *** p < 0.001, * p < 0.05, . p < 0.1.
Table 6. Estimated results of the logistic regression model in Eastern Area of the four independent variables.
Table 6. Estimated results of the logistic regression model in Eastern Area of the four independent variables.
Variables19852014
βS.E.p (>|z|)Exp (β)βS.E.p (>|z|)Exp (β)
Intercept0.5358 **0.20350.00851.70881.1041 ***0.12700.00003.0166
Elevation0.00330.00900.71471.00330.00080.00520.88031.0008
Slope0.0564.0.03010.06121.0580−0.0400 *0.01740.02160.9608
Drainages distance0.00000.00000.59901.00000.00000.00000.24821.0000
Street distance−0.0007 ***0.00010.00000.9993−0.0011 ***0.00010.00000.9989
PCP71.0983 70.1521
Neg. R20.28 0.31
2 Log Likelihood798.9 2367.2
Significant codes: *** p < 0.001, ** p < 0.01, * p < 0.05, . p < 0.1.

Share and Cite

MDPI and ACS Style

Alqurashi, A.F.; Kumar, L.; Al-Ghamdi, K.A. Spatiotemporal Modeling of Urban Growth Predictions Based on Driving Force Factors in Five Saudi Arabian Cities. ISPRS Int. J. Geo-Inf. 2016, 5, 139. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi5080139

AMA Style

Alqurashi AF, Kumar L, Al-Ghamdi KA. Spatiotemporal Modeling of Urban Growth Predictions Based on Driving Force Factors in Five Saudi Arabian Cities. ISPRS International Journal of Geo-Information. 2016; 5(8):139. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi5080139

Chicago/Turabian Style

Alqurashi, Abdullah F., Lalit Kumar, and Khalid A. Al-Ghamdi. 2016. "Spatiotemporal Modeling of Urban Growth Predictions Based on Driving Force Factors in Five Saudi Arabian Cities" ISPRS International Journal of Geo-Information 5, no. 8: 139. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi5080139

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop