1. Introduction
Recent developments in remote sensing and GIS make it possible to detect changes and devise strategies based on these changes. GIS-based geostatistical tools for urban planning have been increasingly powerful and versatile. GIS can be used to map and analyze the geographical distribution of population for the purpose of making future decisions and increase planning for agriculture, economic, education and health development (Dionissios et al., 2013). Geostatistics is a methodology for integrating spatial and temporal coordinates of observations into data processing (Simon, 2000). The precision of various triangular irregular approaches to spatial interpolation of spatial distribution patterns has been tested by many scientists, (Goovaerts, 2006) provides a generalization of Poisson kriging, according to which the size and shape of administrative units and the population density (Xiao and Alexis, 2019).shows the effectiveness of parcel and decennial census data in highly developed urban cores for areal geostatistical interpolation, and calls for future research on population characteristics (Ferry et al., 2016). reviewed and contrasted three commonly used spatial interpolation techniques: ordinary kriging, inverse distance weighting and triangular irregular networks, and developed a new distribution-based distance weighting (DDW) spatial interpolation method to be used with temporal dimensions for large-scale spatiotemporal interpolation in environmental modeling. The geographical distribution of the population of any region reflects natural, human and financial circumstances, because these variables have a direct effect on the population and thus on the population concentration in a specific region or on the random allocation in another region. Natural factors may be the direct impact on the spatial distribution of the population in a given region, but economic and social factors may also make the area densely populated despite the lack of relevance in terms of natural human habitat (El Nahry, 2007).
Generally, by knowing certain patterns and understanding which factors have the greatest effect on population density and overall population, it’s possible to forecast future growth or possible decrease in population at any area of the world, as well as spatially locate them all over the world. It will then be possible to create policies that can safeguard the environment, enable the planet’s sustainability and enable the worldwide population to continue the alter (Minmin et al., 2018).
The present study is focus on applying GIS for the purpose of geostatistical analyzing and exploring data variability and spatial relationships of the geographical distribution pattern with the size of population according to the area in the residential neighborhoods of in the both sides of Mosul city. It’s attempted to find the relationship between the statistical variables and the spatially population size which will be provide a new thinking for the future urbanization of Mosul city.
2. Study Area
1) Location
Mosul city is the center of Nineveh governorate which is located in the northern part of Iraq. The geographical location the city is shown in Fig. 1 and can be noted to be found at longitudes 43°02′59.65″E and 43°13′57.89″E. and at latitudes 36°17′23.86″N and 36°25′45.05″N. Mosul is located on both sides of the Tigris river, dividing it into two parts locally known as left and right sides. The morphological characteristics of both sides of the city differ. The right side is characterized by its height as compared to the left side, where the heights exceed (289 m) in some places. Due to the existence of the Tigris flood plain, which represents the largest percentage of the left side and reaches the height of (220 m), the left side, it appears flatter than the right side.
Fig. 1. Map of the study area.
2) Climate
According to the Koppen classification, Mosul city has a hot semi-arid climate (BSh), verging on the Mediterranean climate (Csa), with extremely hot dry summers and moderately wet, relatively cool winters (Deliang and Hans, 2013).
3) Population Size
Population distribution refers to the distribution of a population at a given time in a certain geographical area or spatial form of a population (Chen, 2015), that’s mean, a certain time connected with the total summation must be noted when speaking about the population and their distribution, because there is no a fixed population group but it’s always subject to change by influencing variables (births rate, deaths rate, all kinds of migration and administrative changes).As her process of urbanization accelerates, the population is increasingly concentrated in urban areas, as is the case in Mosul city according to the census carried out. Table 1 lists the total population of Mosul city according to population censuses(Special Development Plan-Nineveh (2010-2020)).
Table 1. The total population versus population censuses of Mosul city
From Table 1, it concludes very important data to enable us to draw the results of the population increases of the Mosul city over time periods, of course these population increase sled to the establishment of modern residential neighborhoods within the city. However, the increase in population expansion on the left side is more than the right side of the city due to the modern urban planning and the increase of commercial and economic centers as well as the presence of most government departments in this side, which led to the movement of many families from the right side and even from the rural areas in the districts of the govern orate to the left side of city.
3. Theoretical Consideration
Semi-variogram is a graphical demonstration for presenting an image of spatial correlation of data.AGIS geostatistical analyst was used to determine Mosul’s spatial distribution of population variable through Kriging approach. Kriging uses the covariance and semi- variogram to make optimal prediction (Omer et al., 2012) Semi-variogram is a graphical illustration of the spatial correlation of data and a high value indicates a high level of co-movement.Thus, it takes two sample location and find the distance between them as it showed in equation (1) (Johnston et al., 2003):
\(\gamma(h)=\frac{1}{2 \pi(h)} \sum_{i=1}^{n(h)}[Z(x i)-Z(x i+h)]^{2}\) (1)
Where n(h) is the number of pairs of data that are approximately separated by lag (h) within a given distance and direction class. If the values at z (xi) and z (xi + h) are automatically correlated, the result will be small relative to an uncorrelated pair. The semi variogram modeling also can determine the best fit for the model that passes through semi-variogram points by calculating the variance between the input pairs variables. The Cross-validation was used to assess performance of predictions. Possible and the root mean standardized square error should be near 1 (Bilgehan and Ali, 2010). The correlation between the variables is an indication of the changes in the variance of the variables and is denoted as γ(h) using the formula in equation (2) (Kazem et al., 2015):
2γ(h) = 1/n∑in = 1[Z(xi + h) – Z(xi)] (2)
The distance is indicated by h, Z(xi + h) and Z(xi) are indicated by point xi + h and xi. Three parameters; nugget, sill and range are important in the interpretation of the semi-variogram as shown in Fig. 2. The sill corresponds to the overall variance in the dataset and the range is the maximum distance of spatial auto correlation.The nugget variance is the positive intercept of the semi-varicogram and can be caused by measurement errors or spatial sources of variation at distances smaller than the sampling interval or both (Mauricio et al.,2014). Finally, the model’s efficiency can be evaluated using the root mean square error (RMSE) that is a Z*(xi) feature and can use the following expression:
\(\mathrm{RMSE}=\sum_{i=1}^{n}\left(\mathrm{Z}^{*}(\mathrm{xi})-\mathrm{Z}(\mathrm{xi})\right)^{2}\) (3)
Where, Z(xi) and Z*(xi) are the observed and predicted values, respectively.
Fig. 2. Typical semi-variogram and its parameters
4. Methodology
1) Data
The data adopted for this research were taken from Nineveh Planning Directorate, its included spatial and descriptive data about population censuses of Mosul city for the year 2013. The dataset also included a Shape files for the Mosul district and its neighborhoods as well as a shape file for the Tigris River.The dataset also included a Panchromatic image satellite (Worldview 2) with a spatial resolution of (46cm) covered the study area (https://www.geoimage.com.au/satellite/worldview-2). The sensor bands and resolution are listed in Table 2.
Table 2. Worldview-2 bands and resolution
2) Method
Geostatistics is a combination of techniques and accidental models that analyze the characteristics of spatial data. This technique has the capacity and power of estimation and can also express an error in estimation. The estimation of spatial structure and distribution variable values and the evaluation of the error associated with this estimate is called kriging geostatistical approach (Andreas, 2013).
In current study, the ArcGIS10.3 / ordinary kriging approach was used to predict the relationship between the area and population attribute data fields provide by the Nineveh Planning Directorate (NPD) for the year (2013) for both sides of the Mosul city as shown in Fig. 3. Mosul city is divided into nine sectors (five in the left side and four in the right side) for the purpose for facilitating the management of municipal services in the city (Mosul Municipality Directorate, 2019), all these sectors are including one hundred forty-one neighborhoods and covering an area of about 240km2. Based on the population size and area attribute field, the total sample points taken by the study are 268 samples (including those who transgressed on agricultural land).
Fig. 3. The two sides of the study area.
The spatial distributions of the samples is shown in the Fig. 4 Below for the years (2013), while Fig. 5 shows the classification map of the neighborhoods with respect to population size of the study area. Fig. 4 and Fig. 5 gives an important indicator of population distribution on both sides of the city of Mosul, as we note from the two figures, that the density of population distribution in areas with proper municipal planning is greater than the areas spread across the borders of the city of Mosul, which have been distributed by cooperative societies and are mainly agricultural land and are not subject to municipal planning. Also, the high-density areas Shown in Fig. 4 are characterized by offering health and municipal services as well as the available of schools, water and electricity.
Fig. 4. The spatial distribution map of the adopted samples.
Fig. 5. The spatial distribution map of the adopted samples
The samples dataset were used for prediction the spatial distribution of the population size with the area of the neighborhoods by using Ordinary Kriging interpolation. Before the actual interpolation was done, a semi-variogram analysis was conducted to determine the spatial similarity of the data samples using the ArcGIS10.3 geostatistical analysis tool. After defining and analyzing the spatial pattern distribution of the samples using the semi-variogram function, cross-validation have been used to decrease the errors of sample points distribution by implementing the regression function.The prediction map was calculated to estimate the unsampled points.
5. Results and Discussion
Ordinary kriging is the most widely used method of kriging interpolation. This helps to calculate a value at a point in a region known to be a semi-variogram by using estimation location data in the neighborhood. Ordinary kriging implicitly test the mean in a moving neighborhood with local second-order stationarity. Many models of semi-variogram are checked to get the one that fits the dataset of samples, Spherical model was chosen as a best one to illustrate the spatial autocorrelation (covariance scales of input data) of the studied samples dataset and gives the best interpolation map.
Fig. 6 shows the semi-variogram model for the data of population samples in both sides of Mosul city. The semi-variogram displays the variances within groups of observations (red dots) plotted as a function of distance between the observations to display the tendency for nearby observations to be more alike than distant observations.
Fig. 6. Semi-variogram model for the samples.
From Fig. 6, it is noted that the semi-variogram model is flattened at a range of (2351.7) in the x-axis which’s corresponding to sill of (7*107) at the y-axis, then all the sample locations separated by distances closer to the range are spatially auto correlated. That is to say, the covariance function decreased between points and obtain similar function, whereas locations farther apart than the range are not auto correlated (increase the distance between points).
Table 3 list the resulted values of the semi-variogram model, the sill was calculated by adding the partial sill to the nugget. To evaluate the accuracy and precision of the resulted spatial interpolation, the cross-validation graph must be determined.
Table 3. The semi-variogram modeling results of the samples
Fig. 7 shows the resulted cross-validation graph which's illustrate the relation between the measured and predicted error values using spherical model.
According to the semi-variogram model and the cross validation approach, the anticipated error values are presented in Table 4.
Fig. 7. The predicted versus measured errors.
Table 4. Performance of the implemented criteria
After analysis of spatial dependence, we evaluated sites that had not been sampled in order to produce surface maps by ordinary kriging. This has been done by defining a circle with four sector to enclose the dataset that are used to predict a value at an unmeasured location as shown in Fig. 8.
Fig. 8. Sampled and unsampled points for the Kriging interpolation
The final spatial prediction map by ordinary kriging method are presented in Fig. 9 for the population size at the year of 2013 in the both sides of the Mosul city.
Fig. 9. Prediction map of the population size.
The prediction map of the population size shown in Fig. 9 Confirms the validity of the results shown in Fig. 6 and Fig. 8 that prediction standard error of that location near sample points generally have a lower error and spatially autocorrelation unlike the far sample points. Fig. 5 and Fig. 9 gives the impression of an irregular distribution of the population size in relation to the area of the neighborhood. There are small areas densely populated, as in many neighborhoods on the right side (spatially in old Mosul city), while there are large areas with little population, and this is evident in the new neighborhoods in the city of Mosul, whether on the right side or the left side of them, noting that there is a large urban expansion on the left side.
6. Conclusion
This study has attempted to predict the spatial distribution of the population size in the Mosul city by using 268 sample dataset using GIS geostatistical analyst / ordinary kriging interpolation method. Generally, results showed that, the semi-variogram/ spherical model was selected as the best to illustrate the spatial autocorrelation of the studied samples dataset with standardized mean square error of (1.014) and addresses spatial autocorrelation and cross autocorrelation related to the distribution of people in urban areas. Rather than ignoring spatial dependency, it models the spatial autocorrelation of population and the impermeable surface fraction by means of variograms and applies them in population interpolation. The resulted prediction map illustrates that the population density is concentrated in the old neighborhoods on the right side and then the left side, while the distribution is random in many neighborhoods that are not subject to city planning. The limitations of the present study were the effort involving of obtaining data from several official authorities in the city, after which it was re-revised and corrected before carrying out the geostatistical analysis and preparing the required maps for the study.
Future work should be done on implementation of GIS classification model to identified the land use of the Mosul city and then find the relation between the population spatial distribution and the land use maps to obtain new map for classifying population distribution according to land use of the both side. Another study can be done to study the distribution of health centers and the services they offer according to the spatial distribution of the population in every neighborhoods of the city.
The last point is that estimates of obtained population spatial distribution are essential for urban planning applications. For example, the population density in residential neighborhoods is the important indicator in sustainability studies of self-dependent regions and in future population census process.
Acknowledgements
The authors thank the Nineveh Planning Directorate for the facilities provided by them in preparing the data adopted in the study. We also thank the University of Mosul for the continuous scientific support in carrying out scientific research that serves the community.
참고문헌
- Andreas, L., 2013. Kriging methods in spatial statistics, Bachelor's Thesis, Technische Universitat Munchen, Munchen, Germany.
- Bilgehan, N. and B. Ali, 2010. Groundwater Quality Mapping in Urban Groundwater Using GIS, Environmental Monitoring and Assessment, 160(1-4): 215-227. https://doi.org/10.1007/s10661-008-0689-4
- Chen, L., 2015. A framework of evolution and optimization of regional population distribution for sustainable development in China, The Hong Kong Polytechnic University, Hung Hom, Hong Kong, pp. 88-96, http://hdl.handle.net/10397/35467.
- Deliang, C. and W. C. Hans, 1913. Using the Koppen classification to quantify climate variation and change: An example for 1901-2010, Environmental Development, 6: 69-79. https://doi.org/10.1016/j.envdev.2013.03.007
- Dionissios, P. K., J. K. Vassiliki, and H. A. Ektoras, 2013. Evaluation of three spatial interpolation methods to estimate forest volume in the municipal forest of the Greek island Skyros, Geo-Spatial Information Science, 16(2): 100-112. https://doi.org/10.1080/10095020.2013.766398
- El Nahry, A.H., 2007. Using geostatistical analysis and spectral signature for identifying shale distribution and type in El-Salhia project, Ismaillia, Egypt, International Journal of Geology, 1: 22-26.
- Geoimage Worldview-2, https://www.geoimage.com.au/satellite/worldview-2, Accessed on Feb. 15, 2020.
- Goovaerts, P., 2006. Geostatistical analysis of disease data: accounting for spatial support and population density in the isopleth mapping of cancer mortality risk using area-to-point Poisson kriging, International Journal of Health Geographic, 5(1): 52. https://doi.org/10.1186/1476-072X-5-52
- Houlding, S., 2000. Practical geostatistics: modeling and spatial analysis manual, Springer Science & Business Media, Berlin, Germany, 159p.
- Iraqi Agricultural Meteorological Network, https://agromet.gov.iq, Accessed on Sep. 26, 2019.
- Johnston, K., J.M. Ver Hoef, K. Krivoruchko, and N. Lukas, 2003. Using ArcGIS Geostatistical Analyst, ESRI (Environmental Systems Research Institute), Redlands, CA, USA.
- Kazem, M.A., J. Reza, and V. Javad, 2015. Analyzing Kriging and CoKriging Methods by Using ArcGIS Software in Preparing SP Map of Farahan Plain Soil, Biological Forum - An International Journal, 7(1): 836-846.
- Mauricio, P. B., D. L. Alessandro, and C. F. Alberto, 2014. Semi-variogram models for estimating fig fly population density throughout the year, Pesquisa Agropecuaria Brasileira, 49(7): 493-505. https://doi.org/10.1590/S0100-204X2014000700001
- Minmin, L., H. Biao, G. Renzhong, L. You, C. Yu, and F. Yong, 2018. Study on Population Distribution Pattern at the County Level of China, Sustainability, 10(10): 3598. https://doi.org/10.3390/su10103598
- Mosul Municipality Directorate, http://www.mosulbaladiya.gov.iq, Accessed on Nov. 4, 2019.
- Omer, K., O. Husaini, S. Mohammed, R.M. Ahmad, and M. Shattri, 2012. Geostatistical Analysis of Karst Landscapes, Electronic Journal of Geotechnical Engineering, 15: 913-934.
- Special Development Plan-Nenavah (2010-2020), Ministry of Planning, Report from Nineveh Planning Directorate, 563p.
- Susanto, F., P. De Souza, and J. He, 2016. Spatiotemporal Interpolation for Environmental Modelling, Sensors, 16(8): 1-20. https://doi.org/10.1109/JSEN.2016.2616227
- Xiao, H. L. and M. Alexis, 2019. Areal Interpolation Using Parcel and Census Data in Highly Developed Urban Environments, International Journal of Geo-Information, 8(302): 1-13.