Next Article in Journal
The Spatial Dimension of COVID-19: The Potential of Earth Observation Data in Support of Slum Communities with Evidence from Brazil
Previous Article in Journal
Visit Probability in Space–Time Prisms Based on Binomial Random Walk
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Interaction Effect of Population Density Patterns in Sub-Districts of Northeastern Thailand

1
Department of Geography, Faculty of Humanities and Social Sciences, Mahasarakham University, Maha Sarakham Province 44150, Thailand
2
Geography, Geoinformatics and Resources Management Research Unit, Faculty of Humanities and Social Sciences, Mahasarakham University, Maha Sarakham Province 44150, Thailand
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2020, 9(9), 556; https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi9090556
Submission received: 25 August 2020 / Revised: 14 September 2020 / Accepted: 17 September 2020 / Published: 19 September 2020

Abstract

:
The north-eastern region in Thailand is the largest in area and population. Its average income per capita is, however, the lowest in Thailand. This phenomenon leads to migration to big cities, which are considered economic centres. We investigated the effect of spatial interaction on the population density pattern in 20 provinces in north-eastern Thailand. Data was obtained from the compilation and preparation of the demographic data of 2676 sub-districts for 2002–2017. A field survey was conducted through GPS at educational institutions, hospitals, airports, government offices, and shopping malls. The data was analysed using spatial autocorrelation analysis by a global indicator (global Moran’s I) and a local indicator (local Moran’s I and Getis–Ord Gi*). Eight Mueang districts exhibited the high-high (H-H) cluster pattern or hot spot at an increasing yearly rate. In addition, the area with the highest gravity was located near service sources and was found to have the largest population. Moreover, gravity interaction with service sources had a strong positive correlation with migration patterns. Thus, the cluster of areas with the greatest population density is located within the Mueang district in one of the provinces with most service sources, as these places attract people and consequently industrial factories and service trades.

1. Introduction

The north-eastern region of Isan is the largest in terms of area in Thailand and also hosts the largest population. According to the 12th National Economic and Social Development Plan (2017–2021), the north-eastern region is the economic centre of the Greater Mekong sub-region (GMS) with emphasis on the development of the economic area, tourism, and utilisation of inholdings for the accommodation of the new production bases in the industrial sector and the development of a link from Thailand to Lao and Vietnam. A high speed train project was initiated to connect to the People’s Republic of China through Vietnam [1,2]. However, the north-eastern region is still poverty-stricken and features a higher income inequality than other region. This leads to the migration of populations from rural areas to cities, which feature relatively higher economic and industrial growth compared to other regions [3,4], thus, increasing the population density in big cities. The management of urban area expansion and the development of neighbouring towns have, thus, been established to accommodate population growth. Population distribution information is an important variable whose pattern is used at the lowest level of urban development planning of large areas [5,6].
Cluster analysis is one of the methods used to identify patterns in population density. It is also used to drive the creation of economic prosperity, which corresponds to the labour production sector, land, utilities, infrastructure, transportation, and academic institutions [7,8], as support for areas that could potentially be developed with the aim of improving the level of urbanisation, with major cities as bases for economy, education prosperity, trade, services, and industrial factories [9,10,11,12]. This results in an effective management, an equal income distribution, and the reduction of environmental problems. The cluster concept is dependent on data analysis using geographic information systems, spatial automatic correlation statistics, and gravity analysis using the gravity centre model as a means of investigating density and geographical phenomena [13,14].
Data analysis includes spatial automatic correlation statistics based on global Moran’s I, a global spatial-correlation indicator, and local Moran’s I and Getis–Ord Gi* for the local spatial-correlation indicator (LISA) of the population density of a particular region and nearby areas. These techniques were used to identify population density patterns, which were classified as high-high (H-H), high-low (H-L), low-low (L-L), or low-high (L-H) [15,16,17,18]. Moreover, local G-statistics or Getis–Ord Gi* was also used as an indicator of automatic spatial correlation to identify cluster patterns through hot spot analysis [19]. Gravity and network analyses were conducted to indicate the spatial pattern of locations with a gravity of population density, trade factors, services, and access to the city. The spatial cluster technique was mostly used to analyse groups and identify the distribution of variables, such as epidemics, tourism, poor quality housing, residential density, ecosystems, and the economy [20,21,22,23,24,25]. The gravity model was used to indicate the gravity caused by moving into the centre [13,26,27,28]. Many studies have analysed spatial patterns in sub-districts, in other words, at the lowest level of urban development planning, in the long term using spatial autocorrelation analysis (SAA). Moreover, the gravity model has not been applied to SAA for analysis of spatial interaction effects. Herein, the analysis of SAA demonstrates the patterns of population density at the sub-district levels and nearby areas. The spatial interaction based on the gravity index was investigated in the sub-districts of the entire north-eastern region of Thailand, and this index was correlated various factors. Therefore, this study aimed at analysing spatial interaction effects of population density patterns in sub-districts in 20 provinces in the north-eastern part of Thailand. We present spatial population distribution pattern maps and gravity index data that demonstrates the interaction effect, which indicates access of population from sub-districts to the city centre. In addition, we display the correlation between the gravity index and factors related to population. This study provides strong suggestions to evaluate the interactive relationship between population density patterns and the economic sector. This represents a very important basis for government authorities and sub-district administrative organizations in decision making and planning, such as related to the rational distribution of local people and urbanization policies.

2. Materials and Methods

In this study, the method was included into two parts (Figure 1). First, spatial population distribution patterns were analyzed using SAA. Second, gravity centre model was applied with SAA results to measure the spatial interaction effect of density patterns.

2.1. Study Area

The study area was the north-eastern part of Thailand (hereafter referred to as the Northeast), a region located on the Korat plateau that looks like a basin, which drains westward from a high altitude into an eastward lowland. The Northeast is bordered by the Mekong River (along the border with Laos) to the north and east, by Cambodia to the southeast, and Sankamphaeng Range south of Nakhon Ratchasima to the south. It has the Phu Phan range as the internal range. As the largest region in Thailand, the Northeast covers 20 provinces, namely Amnat Charoen, Bueng Kan, Buriram, Chaiyaphum, Kalasin, Khon Kaen, Loei, Maha Sarakham, Mukdahan, Nakhon Phanom, Nakhon Ratchasima, Nong Bua Lam Phu, Nong Khai, Roi Et, Sakhon Nakhon, Si Sa Ket, Surin, Ubon Ratchathani, Udon Thani, and Yasothon (Figure 2). This region covers an area of 168,854 km2. According to national statistics, the Northeast consists of 2676 sub-districts. The Northeast has a population of 21,551,434 people as of 2002, while the population was 21,381,249 in 2007, 21,697,488 in 2012, and 21,989,477 people in 2017. The central and lower parts of this region consist of the vast Thung Kula Ronghai area, where the highest amount of jasmine rice is cultivated in the country. Sugarcane, cassava, and rubber are also grown here as agricultural products for the industrial sector.

2.2. Data Acquisition and Preparation

Data from 2676 sub-districts in the Northeast of Thailand for the period ranging from 2002 to 2017 were obtained from the National Statistical Office and Ministry of Interior for use as the database for the analysis of the relationships among population size, public utilities, and public assistance and excluded the illegal population in the analysis. The National Statistic Office obtained the data online through their network, collected by relevant sub-district administrative organisations. For accuracy, we compared the collected data to the census registration from the sub-district administrative organization. Therefore, the population data were similar when cross-checking the data of the National Statistic Office and the sub-district administrative organisations (SAO) through field surveys and interviews with the authorities. The next step involved the surveyance and data collection of geographic coordinates using global positioning system (GPS) and high-resolution satellite images. The data included five types of public amenities, including schools, hospitals, government offices, airports, and shopping malls. Schools were classified into three categories: large-sized schools with 1500 students or more, medium-sized schools with 500–1499 students, and small-sized schools with less than 500 students. These were used for the analyses of an accessibility and gravity model, which reveals population interactions [29,30,31]. Moreover, data on public utilities, such as main and secondary roads, were also collected to determine factors related to speed, distance, and time according to the standard speed of the Department of Rural Roads. The data were then input into the geographic information system software ArcGIS (ESRI, Canada) for the analysis of road networks and the shortest travel time from sub-districts to the centre of population gravity.

2.3. Spatial Autocorrelation Analysis

This study used SAA at two levels: global and local indicators. Global Moran’s I was used for global indicator spatial autocorrelation analysis, whereas local Moran’s I and Getis–Ord Gi* were used for local indicator spatial autocorrelation. The relationship between population density in each of the sub-districts and in 8–16 nearby sub-districts was examined to analyse the continuity weight using Queen contiguity. Using Euclidean distance, the next step involved a selection of a distance metric via the GeoDa software for calculation. The confidence level for the result was 95% [32,33,34].
To analyse the population density using one variable, global Moran’s I was used to examine the global indicator spatial autocorrelation. For the years 2002, 2007, 2012, and 2017, the data concerning population density were used to identify the overall patterns of the Northeast in one unit at the regional level by comparing the resulting point value in the position under investigation to other positions using a spatial weight matrix [18,35], as shown in Equation (1).
I = n w i j ( x i x ¯ ) ( x j x ¯ ) w ( x j x ¯ ) 2
where I is the spatial correlation based on global Moran’s I, x i is the population per area i , x j is population per area j which is adjacent to area i , x ¯ is the average population density, n is the total number of sub-districts in the Northeast, and W i j is the spatial weight matrix of distance from the sub-districts’ centre point.
Local Moran’s I was used to analyse the local indicator of spatial autocorrelation. Such statistics were used to measure the population density correlation of areas at the sub-district level and nearby sub-districts in the pattern of agglomeration from a total of 2676 sub-districts, as shown in Equation (2).
I i = Z i i W i j Z j
where I i is the population density correlation based on local Moran’s I, Z i is the deviation from the population per area i which is close to area j population per area the deviation from the is Z j , i , and W i j is the spatial weight matrix of distance from the sub-districts’ centre point.
The coefficient of correlation based on Moran’s I is between −1 and +1. A positive spatial autocorrelation indicates a significant spatially aggregated distribution of a high or low population density (cluster). In contrast, a negative value indicates a significant spatial difference between the area of each sub-district in terms of population density (disperse), whereas zero indicates a non-significant spatial autocorrelation (random) [36,37]. Local Moran’s I illustrates the patterns of population density clusters in four classes: high-high (H-H), low-low (L-L), high-low (H-L), and low-high (L-H) [15,17,18].
Getis–Ord Gi* was used to analyse local indicators of spatial autocorrelation for the identification of clusters of data, called hot spot analysis. The areas with the highest and lowest population densities yielded Z-score values that were higher than +1.96, called hot spots, and lower than −1.96, called cold spots, respectively. A positive Z-score indicates a high-value cluster, while a negative one indicates a low-value cluster, determined at a confidence of 95%, Getis–Ord Gi* was obtained from Equation (3).
G i * ( d ) = j = 1 n w i , j ( d ) x j x ¯ j = 1 n w i , j ( d ) s n j = 1 n w i , j 2 ( d ) ( j = 1 n w i , j 2 ( d ) ) 2 n 1
where
x ¯ = j = 1 n x j n ,
S = j = 1 n x j 2 n x ¯ 2
and G i * ( d ) is the distribution of the population density (Z-scores), d is the distance from the central point of one sub-district to that of all neighbouring sub-districts, x i is the population per area i , x j is the population per area j which is close to area i is the spatial weight matrix of the W i j distance from the central point of sub-districts, and n is the total number of sub-districts in the Northeast.

2.4. Gravity Centre Model

The gravity model was used to analyse the interaction between the areas that affect the activities for area development. This model is a well-known method of analysing economic phenomena that affect the population. In this work, the clustered districts from the LISA analysis were specifically selected and created mean centre points of sub-districts for applying the model. This study analysed access to amenities (e.g., schools, hospitals, government offices, airports, and shopping malls), which affects population clusters, and the effectiveness of access to the road network. Applying the gravity centre model, the equation of Huff [38,39] and Lakshmanan and Hansen was set and applied [40]. Two factors affecting the analysis included the population in the area and the potential accessibility to the centre sources [14,31,41].
To analyse the gravity index of living in the area surrounding the centre of services with the spatial interaction, the calculation was performed based on the human resources index and the accessibility index using the gravity model, as expressed below.
G i j = P H i j T H i j
G i j is the index of access to the service area. P H i is the human resources index and T H i j is the accessibility index. The next step involved the classification of the gravity of access to the service area into different levels using natural breaks (Jenks). The groups with G i j at the same level were grouped together; the variance within the same group was lowest but the variance between groups was highest. This study classified G i j into four groups, namely, higher, high, medium, and low.
The human resources index ( P H i ) is an attraction factor and used to analyse the population per area, as expressed in Equation (5) [8].
P H i = ( P i P min ) ( P max P min ) + 1
where P H i is the human resources index, P i is the population in the sub-district i , P max and P min are the maximum and minimum of 2676 sub-districts in the north-eastern part of Thailand.
The accessibility index ( T H i j ) is a factor of travelling to the central sources via the shortest path obtained from the analysis of the spatial network structure. A network analysis tool (closed facility function) was used with time and distance to process the shortest path route for providing the travel time. It was used to obtain data concerning the built road line network to determine limited speed by road type, the centre of the sub-district, and the centre of business quarters, with which the shortest path based on road type for travelling to the centre was determined, as stated in Equation (6).
T H i j F j T i j λ j = 1 n F j T i j λ
where T H i j is the accessibility index (probability of area j being visited by the inhabitant of area i ). T i j is the drive time measured between the gravity point and sub-district central point of the clustering pattern. λ is the exponent of the drive time (equals 1), and F j is the attractiveness of area j (equals 1).
After the survey and the collection of data on the geographic coordinates of public utilities, the centre of attraction for services and the destination point of travel to reach the central source was determined using mean centre analysis to determine the weight value of all positions of public utilities, which equal 1, to indicate the central source of the centre of city clusters.

2.5. Interaction between Gravity Index and Related Factors

The spatial correlation analysis between the gravity index, population density, migration, local Moran’s I, and Getis–Ord Gi* was conducted as expressed in Equation (7).
ρ = i = 1 N ( x i μ x ) ( y i μ y ) i = 1 N ( x i μ x ) 2 i = 1 N ( y i μ y ) 2
where ρ is Pearson’s correlation coefficient. x , y are variables, and μ x , μ y are the mean values of variables x and y .
The correlation value is between −1 and +1: a correlation of 0 is an indication of the lack of correlation between variables x and y . A correlation value that is non-zero indicates a correlation between variables x and y . A positive correlation indicates that the two variables change in the same direction, while a negative one means that the variables change in opposite directions. The approach presented here was examined for appropriateness and approved by the Mahasarakham University Ethics Committee in Human Research No. 131/2019.

3. Results

3.1. Spatial Pattern Analysis Using Spatial Autocorrelation

An analysis of the distribution of population using spatial correlation based on global Moran’s I for the period ranging from 2002 to 2017 yielded a positive global Moran’s I value (p < 0.05) every year. This indicated that the area with a similar density (high or low) was clustered and the value of global Moran’ s I increased every year (Table 1).
The spatial autocorrelation analysis of the population density based on local Moran’s I showed that in the year 2002, most areas had no correlation (2053 sub-districts from a total of 2676 sub-districts). Among the areas that were correlated (p < 0.05), most featured the L-L pattern (485 districts) followed by H-H (97 sub-districts), L-H (17 sub-districts) and H-L (24 sub-districts) with a similar amount, which yielded similar characteristics in the results obtained for years 2007, 2012, and 2017 (Table 2).
From 2002 to 2017, the area with the H-H pattern increased every year, whereas the area with the L-L pattern decreased. This indicated an increase in clustered sub-districts with a high population density, while there was a decrease in clustered sub-districts with a low population density. The area with the H-H pattern was mostly observed in Mueng districts and an area roughly resembling Khorat, Mun River, and Chee River basins, whereas the L-L pattern was mostly observed in the mountain range and forests in the west of the Northeast. This area covers the Phetchabun, Dong Phaya Yen, and Dong Sankampang mountain ranges, which drain into the east of the Phu Phan mountain range. The south was the Khao Phanom Dong Rak mountain range, which runs along the entire extent from east to west, and this area is adjacent to the border with Cambodia (Figure 3).
The data of each province for the period ranging from 2002 to 2007 show that only eight provinces were clustered (from a total of 20 provinces), with the correlation being unequal (p < 0.05) as illustrated in Figure 2. Five provinces were initially found to have the H-H pattern (Buriram, Roi Et, Ubon Ratchathani, Nong Khai, and Surin). Three more provinces were found to have the H-H pattern (Nakhon Ratchasima, Khon Kaen, and Udon Thani), as illustrated in Figure 3 and Figure 4. For the last 15 years, many areas have expanded from the centre of this region, such as Nakhon Ratchasima, being the capital of the lower north-eastern region, the gateway linking the central region and the eastern region. This city is the hub of economic, trade, industrial and agricultural advancements, Khmer cultural tourism, sports, medicine and education, geographically linking the economic advancement with other provinces in the region, including Buriram, Surin, and Ubon Ratchathani. Khonkaen is the capital of the central north-eastern region, which clustered and expanded to nearby districts and created new cities in the area. Udon Thani is the capital of the higher north-eastern region, being the hub of services, agricultural industries, and sugarcane mills and links to Nong Khai, which is the economic hub around the Mekong River, allowing trade across the border. Moreover, the population settled along the main transportation network connecting Nakhon Ratchasima, Khon Kaen, Udon Thani, and Nong Khai.
The analysis through Getis–Ord Gi* indicated that in 2002, clusters of high population density areas or hot spots were in 122 sub-districts. Contrarily, clusters of low population density areas or cold spots were observed in 506 sub-districts. There were 2048 sub-districts that were not clustered (p < 0.05) remaining. These results from 2002 to 2017 also suggested an increase in the number of sub-districts with a hot spot every year, whereas those with a cold spot decreased as illustrated in Table 2. This table shows the result of the analysis based on Getis–Ord Gi* (Figure 4), which corresponded with the result of the analysis based on local Moran’s I.

3.2. Gravity Centre Index

Two factors were used to analyse the gravity-centre index: gravity index and the travel time it took to reach the central source. An analysis of eight cities with cluster patterns indicates that the sub-districts were classified at the highest level (dark red). They were located near the service centre. However, some sub-districts with a gravity index at the highest level were found to be located far from the central source of the city. This yielded different results when compared to those sub-districts with a gravity index at a high level (light red), as these areas had a large population size. Ubon Ratchathani in Rai Noi sub-district, Meung District, is such an example; although it had a population size of 25,190 people and a distance of 11,957 m from the centre, which is quite far, it was still classified at the highest level (dark red). In contrast, while the sub-districts in the cities had a population size of 11,291, and the distance from the centre was around 1649 m, they were classified at a high level (light red). The area within a 30 km radius around the centre of the city was designated as the point of maximum accessibility, where business quarters and services are found, eventually leading to a high population density. The population density, however, decreases when such areas are far from the centre. The centres of regional capitals are service hubs, making them high gravity index centres, including Nakhon Ratchasima, Udon Thani, and Khon Kaen, while Buriram, Roi Et, Surin, and Nong Khai are in the centre of the city, and Ubon Ratchathani is in the educational centre area as illustrated in Figure 5. The gravity model can be used to plan economics and locations, such as regional trade agreements near provinces of geographic indicators. Meanwhile, the effects to economic growth are internal, and it is external economies that create jobs, thus labour movement grows.

3.3. Interaction of Gravity Index and Factors Related to Population

To examine the interaction gravity index and related factors (population density, migration based on local Moran’s I, and Getis–Ord Gi*), this study analysed the area data of the sub-districts with a statistically significant gravity index only. These sub-districts were located in eight provinces, namely Buriram, Khon Kaen, Nakhon Ratchasima, Nong Khai, Roi Et, Surin, Ubon Ratchithani, and Udon Thani.
The results of the analysis of the gravity index and population density interaction suggested that the correlation (r) was positive, between 0.23 and 0.99 (Figure 6), which indicated a gravity index and population density at different levels, except for provinces such as Ubon Ratchathani. An r = 0.23 indicates that the area with the highest gravity was slightly associated with population density.
The gravity index was positively associated with the migration of the population with a value between 0.61 and 0.98. Upon considering all provinces, the r value was higher between the gravity index and other factors (population density, local Moran’s I and Getis–Ord Gi*).
The gravity index was positively associated with local Moran’s I and Getis–Ord Gi*, with a value between 0.54 and 0.98 (local Moran’s I) and 0.46 and 0.93 (Getis–Ord Gi*). Local Moran’s I had a higher value due to the analysis of the standard division value of all the sub-districts in this region, combined with the nearby areas, while Getis–Ord Gi* measured the nearby areas only. Therefore, the gravity index was highly correlated with migration, population density, local Moran’s I, and Getis–Ord Gi*.

4. Discussion

An analysis of spatial data distribution using the spatial autocorrelation method was practically applied to study medium and small-scale areas [13,23,42,43]. Meanwhile, an analysis of large scale areas in different time periods was infrequently performed [44,45]. This study analysed regional distributions at the sub-district level in the years 2002, 2007, 2012, and 2017 using spatial autocorrelation analysis. There are two levels of the spatial autocorrelation analysis: (1) global indicator and (2) local indicator. The global indicator tool is Moran’s I, while the local indicator tools are local Moran’s I and Getis–Ord Gi*. In addition, the analysis of the gravity model of the service centre and the interaction between the gravity index and related factors were performed.
The analysis of the spatial population distribution in the Northeast using global Moran’s I yielded a positive index in every year of study. This indicated the presence of a clustering of high population density sub-districts. Furthermore, the index of analysis increased throughout the 15 years, indicating positive clustering and the presence of relatively stronger positive spatial autocorrelations among provinces over time. The results of development policies in the Northeast over many years, along with the status of the region as an important connecting area for neighbouring countries, has caused the invigoration of the economy for some provinces along with labour related migrations for education also contributed to this trend.
The results of the local Moran’s I and Getis–Ord Gi* analysis were consistent. A dense population was identified in clusters in urban areas near the central areas of eight provinces. This indicated that the neighbouring sub-district population density changed in the same direction. For example, if one sub-district has a high population density, the neighbouring sub-district will also have the same densely populated level. In addition, if one district has a low population density, the neighbouring sub-districts will also have a low population density. Throughout the period ranging from 2002 to 2017, hot spots increased, while cold spots decreased.
One of the reasons for the observed patters is that the population in the Northeast increased rapidly from 21,551,434 people in 2002 to 21,989,477 people in 2017. Public utilities and, agricultural processing plants increased according to the Department of Industrial Works (2014) [46]. There were 41,220 registered factories in 2005; meanwhile, 44,143 factories existed in 2014. The transportation or movement between provinces has become more convenient due to the extension of the main roads, including the fact that in 2015 the airport in Roi Et province was upgraded to accommodate long haul flights. In addition, regional airports were renovated and upgraded to international levels in Khon Kaen, Ubon Ratchathani, and Udon Thani provinces [47]. The economic growth resulted in an average per capita income increase for the population during 2005–2014 from 32,308 to 71,286 Baht per capita [48]. Similarly, Luo [27] found that economic fragmentation was an important factor affecting the patterns of the population distribution in the area of Sichuan-Yunnan-Guizhou province, China.
Most of the areas with high-low and low-high clustering were located around the city areas due to the gap between economic development cities and that in urban areas. Therefore, the rates of economic development and urban population density have also increased. Some reasons for this are migration, whereby people move from the suburbs into the city areas. Similarly, in Nanijung, Jiangsu province, China, a high-low clustering has been reported due to the large population in Nanijung downtown. However, the neighbouring cities has small populations [28].
Eight provinces with H-H clustering were different from the 12 provinces in the Northeast without clustering, especially in the border cities, such as Surin, Nong Khai, Ubon Ratchathani. Meanwhile, Buriram province is considered to be the distinctive point of tourism, sporting destinations for football, and world-class racing tracks. Nakhon Ratchasima, Khon Kaen, Udon Thani provinces are defined as the centres of trade and education among neighbouring countries. A cluster has also been found in the Mueang district area, Roi Et province. Although the province is not considered to be a distinct city like the others, the clustering results revealed a dense population in the centre of the city and areas around the city. The reason for this pattern is likely due to the province being the centre of the Northeast, where it represents a hub for the connection of services and product distribution to nearby provinces.
The results of the gravity analysis in the service centres of cities with high indexes such as Nakhon Ratchasima, Udon Thani, and Khon Kaen, were consistent with the reality of the city areas. As these are a large cities in the region and economic centres of trade, health services, and education, which are utilities for supporting population growth, the cities have expanded rapidly from the past to the present [49]. Similar results from Bunea [50] indicated that migration was related to the attraction of the area, especially the population size and distance of migration. If the population was large, there was also a chance of migration. On the other hand, migration opportunities were reduced in instances of far proximity. Meanwhile, Mafi [51] suggested an analysis of accesses to service areas within the shortest time based on speed and types of roads. The results of the analysis were consistent with the shortest time used to access the service centre, speed limited by road types, and the quality of the road affecting the gravity of the local population travels and transportations.
An analysis of spatial data distribution using spatial autocorrelation and the comparison of gravity indexes based on population-related factors showed that the coherent provinces of the sub-district had a high population density and a large population. The main reasons for the increase in clustering were economic and service sources, for instance, hospitals, multifarious, famous educational institutions, being promoted as a world-class sports city, a border town of international trades, et cetera, which can attract populations of migrants in search of job opportunities [50]. Therefore, the development of public utilities and transportations should play a key role in attracting people, which will also affect industrial properties, trade, and services.
For further suggestions in the study of the phenomenon of populations incoherent with accurate patterns to identify economic conditions, in areas with high growth and urbanisation, the sizes and concentrations will be higher than the urbanised communities. Although this study is limited to the collection of the gross domestic product (GDP) at the district levels of Thailand, population density data per area can show the patterns of population distribution, and correlated to gravity and population-related factors, which can effectively solve the problems of local data storages in large areas, regions, and countries. In addition, the approach can further be developed to quickly identify the concentration and distribution patterns of different economic zones. The results of the study can be utilised to help develop an area in line with the population density of each province in the Northeast with efficiency and accuracy in a wide area to a local level.

5. Conclusions

From the study of spatial interactions based on a population density model of the 20 provinces of the Northeast region of Thailand from 2002 to 2017, the following conclusions are drawn:
1. Eight provinces had a positive clustering, which were all sub-districts in the Muang districts (Buriram, Khon Kaen, Nakhon Ratchasima, Nong Khai, Roi Et, Surin, Ubon Ratchathani and Udon Thani). Clustering of the population also increased in neighbouring sub-districts from 2002 until 2017.
2. The number of adjacent sub-districts with a similarly high population density (H-H or hot spots), which were all urban areas, increased over the period ranging from 2002 to 2017. However, the agglomeration of sub-districts with a low population density (L-L or cold spots), most of which were mountainous and forested areas, decreased. However, most of the areas (approximately 77% of the total number of sub-districts in the Northeast) were not related.
3. The districts located near trade and service centres (educational institutions, hospitals, government offices, and department stores) exhibited more interactions than those in remote areas of the sub-districts. These interactions were also dependent on the population size: they were high among large populations and low among small populations.
4. The gravity index of access to service areas was positively related to migration. This was more intense than the relationships between the gravity index of access to service centres and population density, local Moran’s I, and Getis–Ord Gi*.
5. We used population data from the National Statistic Office in Thailand and excluded data on the illegal population. Future studies should collect and combine the number of households via field surveys or sub-district administrative organizations (SAO) for processing, which will provide more accurate results of spatial population distribution patterns.

Author Contributions

Conceptualization: Savittri Ratanopad Suwanlee and Jaturong Som-ard; data curation: Savittri Ratanopad Suwanlee and Jaturong Som-ard; formal analysis: Jaturong Som-ard and Savittri Ratanopad Suwanlee; investigation Jaturong Som-ard and Savittri Ratanopad Suwanlee; methodology: Jaturong Som-ard and Savittri Ratanopad Suwanlee; writing—original draft preparation: Savittri Ratanopad Suwanlee and Jaturong Som-ard; writing—review and editing: Savittri Ratanopad Suwanlee and Jaturong Som-ard; funding acquisition: Savittri Ratanopad Suwanlee and Jaturong Som-ard. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Faculty of Humanities and Social Sciences, Mahasarakham University, Fiscal Year 2019 under unit number 033/2019. The funding agency made no contribution to the study design, data collection, analysis and interpretation of data, writing of the report, and the decision to submit the article for publication.

Acknowledgments

We are thankful to the National Statistical Office and Ministry of Interior for supply population and migration data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Foran, T.; Ward, J.; Leitch, A.; Smajgl, A. Exploring Mekong Futures: Compilation of Scenarios; Technical Report, June 2012; CSIRO: Canberra, Australia, 2012. [Google Scholar]
  2. National Statistical Office. Number of Population from Registration by Sex, Area, Density, House, Region and Province: 2008–2017. Available online: http://statbbi.nso.go.th/staticreport/page/sector/th/01.aspx (accessed on 26 August 2019).
  3. Chattranond, O. Isan in Asean: Inequality and regional development. J. Public Adm. Pol. 2020, 9, 81–104. [Google Scholar]
  4. Pinyochatchinda, S.; Walsh, J. The role of SEZs in Thailand’s regional economic development. Bus. Manag. Rev. 2015, 6, 199. [Google Scholar]
  5. Turner, B.; Mauer, B.A. Geographical population analysis: Tools for the analysis of biodiversity. J. Ecol. 1994, 82, 981. [Google Scholar] [CrossRef]
  6. Páez, A.; Scott, D.M. Spatial statistics for urban analysis: A review of techniques with examples. GeoJournal 2004, 61, 53–67. [Google Scholar] [CrossRef]
  7. Hu, Z.; Lo, C. Modeling urban growth in Atlanta using logistic regression. Comput. Environ. Urban Syst. 2007, 31, 667–688. [Google Scholar] [CrossRef]
  8. Ramos, R. Gravity models: A tool for migration analysis. IZA World Labor 2016, 239, 1–10. [Google Scholar] [CrossRef] [Green Version]
  9. Carrothers, G.A.P. An historical bedew of the gravity and potential concepts of human interaction. J. Am. Inst. Plan. 1956, 22, 94–102. [Google Scholar] [CrossRef]
  10. Gordon, I. Densities, urban form and travel behaviour. Town Ctry. Plan. 1997, 66, 239–241. [Google Scholar]
  11. Jones, C.; Pimdee, P. Innovative ideas: Thailand 4.0 and the fourth industrial revolution. Asian Int. J. Soc. Sci. 2017, 17, 4–35. [Google Scholar] [CrossRef]
  12. Shao, Y.; Walsh, S.J.; Entwisle, B.; Rindfuss, R.R. Spatial clustering and urban settings of rural migrants in Bangkok, Thailand. Geocarto Int. 2008, 23, 35–52. [Google Scholar] [CrossRef]
  13. Liu, C.; Wu, X.; Wang, L. Analysis on land ecological security change and affect factors using RS and GWR in the Danjiangkou Reservoir area, China. Appl. Geogr. 2019, 105, 1–14. [Google Scholar] [CrossRef]
  14. Orsi, F.; Geneletti, D. Using geotagged photographs and GIS analysis to estimate visitor flows in natural areas. J. Nat. Conserv. 2013, 21, 359–368. [Google Scholar] [CrossRef]
  15. Anselin, L. Local Indicators of spatial association-LISA. Geogr. Anal. 2010, 27, 93–115. [Google Scholar] [CrossRef]
  16. Lei, C.; Zhang, A.; Qi, Q.; Su, H.; Wang, J. Spatial-temporal analysis of human dynamics on urban land use patterns using social media data by gender. ISPRS Int. J. Geo-Inf. 2018, 7, 358. [Google Scholar] [CrossRef] [Green Version]
  17. Mathur, M. Spatial autocorrelation analysis in plant population: An overview. J. Appl. Nat. Sci. 2015, 7, 501–513. [Google Scholar] [CrossRef] [Green Version]
  18. Oden, N. Adjusting Moran’s I for population density. Stat. Med. 1995, 14, 17–26. [Google Scholar] [CrossRef]
  19. Wang, W.C.; Chang, Y.-J.; Wang, H.-C. An application of the spatial autocorrelation method on the change of real estate prices in Taitung city. ISPRS Int. J. Geo-Inf. 2019, 8, 249. [Google Scholar] [CrossRef] [Green Version]
  20. Assuncao, R.M.; Reis, E.A. A new proposal to adjust Moran’s I for population density. Stat. Med. 1999, 18, 2147–2162. [Google Scholar] [CrossRef]
  21. Gutiérrez, A.; Arauzo-Carod, J.-M. Spatial analysis of clustering of foreclosures in the poorest-quality housing urban areas: Evidence from Catalan cities. ISPRS Int. J. Geo-Inf. 2018, 7, 23. [Google Scholar] [CrossRef] [Green Version]
  22. Millward, H.; Bunting, T. Patterning in urban population densities: A spatiotemporal model compared with Toronto 1971–2001. Environ. Plan. A Econ. Space 2008, 40, 283–302. [Google Scholar] [CrossRef]
  23. Rodríguez-Rangel, M.C.; Sánchez-Rivero, M. Spatial imbalance between tourist supply and demand: The identification of spatial clusters in Extremadura, Spain. Sustainability 2020, 12, 1651. [Google Scholar] [CrossRef] [Green Version]
  24. Sarrión-Gavilán, M.D.; Márquez, M.D.B.; Mora-Rangel, E.O. Spatial distribution of tourism supply in Andalusia. Tour. Manag. Perspect. 2015, 15, 29–45. [Google Scholar] [CrossRef]
  25. Schuurman, N.; Peters, P.A.; Oliver, L.N. Are obesity and physical activity clustered? A spatial analysis linked to residential density. Obesity 2009, 17, 2202–2209. [Google Scholar] [CrossRef] [PubMed]
  26. Jeefoo, P. Analyzing spatial clustering and hotspots detection of HIV/AIDS prevalence using GIS technology. Int. J. Geoinf. 2016, 12, 65–73. [Google Scholar]
  27. Luo, Y.; Yu, H.; Liu, S.; Liang, Y.; Liu, S. Spatial heterogeneity and coupling of economy and population gravity centres in the Hengduan Mountains. Sustainability 2019, 11, 1508. [Google Scholar] [CrossRef] [Green Version]
  28. Yang, M.; Ma, J.; Jia, P.; Pu, Y.; Chen, G. The use of spatial autocorrelation to analyze changes in spatial distribution patterns of population density in Jiangsu province, China. In Proceedings of the 2011 19th International Conference on Geoinformatics, Shanghai, China, 24–26 June 2011. [Google Scholar]
  29. Haggett, P. Geography: A Global Synthesis, 4th ed.; Prentice Hall: Harlow, UK, 2001. [Google Scholar]
  30. Haynes, K.E.; Fotheringham, A.S. Gravity and Spatial Interaction Models; Sage: Beverly Hills, CA, USA, 1984; Volume 2. [Google Scholar]
  31. Kurowska, K.; Kryszk, H.; Kietlińska, E. The Use of Gravity Model in Spatial Planning. In Proceedings of the International Conference on Environmental Engineering, ICEE, Vilnius, Lithuania, 27–28 April 2017. [Google Scholar]
  32. Bivand, R.S.; Wong, D.W.S. Comparing implementations of global and local indicators of spatial association. TEST 2018, 27, 716–748. [Google Scholar] [CrossRef]
  33. Mitchell, A. The ESRI guide to GIS analysis. In Spatial Measurements and Statistics; ESRI Guide to GIS Analysis; ESRI Press: Redlands, CA, USA, 2005. [Google Scholar]
  34. Ord, J.K.; Getis, A. Local spatial autocorrelation statistics: Distributional issues and an application. Geogr. Anal. 1995, 27, 286–306. [Google Scholar] [CrossRef]
  35. Junior, G.B.; de Paiva, A.C.; Silva, A.C.; de Oliveira, A.C.M. Classification of breast tissues using Moran’s index and Geary’s coefficient as texture signatures and SVM. Comput. Biol. Med. 2009, 39, 1063–1072. [Google Scholar] [CrossRef]
  36. Rutter, C.M.; Isaaks, E.H.; Srivastava, R.M. An introduction to applied geostatistics. J. Am. Stat. Assoc. 1991, 86, 548. [Google Scholar] [CrossRef]
  37. Fortin, M.-J.; Dale, M.R.T.; Ver Hoef, J. Spatial analysis in ecology. In Encyclopedia of Environmetrics; El-Shaarawi, A.H., Piegorsch, W.W., Eds.; John Wiley & Sons Ltd.: Chichester, England, 2002; Volume 4, pp. 2051–2058. [Google Scholar]
  38. Huff, D.L. A Probabilistic analysis of shopping center trade areas. Land Econ. 1963, 39, 81. [Google Scholar] [CrossRef]
  39. Huff, D.L. Parameter estimation in the Huff model. ArcUser. Ann. GIS 2003, 19, 34–36. [Google Scholar]
  40. Lakshmanan, J.R.; Hansen, W.G. A retail market potential model. J. Am. Inst. Plan. 1965, 31, 134–143. [Google Scholar] [CrossRef]
  41. Lee, S.-H.; Choi, J.-Y.; Yoo, S.-H.; Oh, Y.-G. Evaluating spatial centrality for integrated tourism management in rural areas using GIS and network analysis. Tour. Manag. 2013, 34, 14–24. [Google Scholar] [CrossRef]
  42. Bai, J.; Li, S.; Wang, N.; Shi, J.; Li, X. Spatial spillover effect of new energy development on economic growth in developing areas of China—An empirical test based on the spatial Dubin model. Sustainability 2020, 12, 3249. [Google Scholar] [CrossRef] [Green Version]
  43. Sánchez-Martín, J.-M.; Gallego, J.I.R.; Blas-Morato, R. Hot Spot analysis versus cluster and outlier analysis: An Enquiry into the grouping of rural accommodation in Extremadura (Spain). ISPRS Int. J. Geo-Inf. 2019, 8, 176. [Google Scholar] [CrossRef] [Green Version]
  44. Eichhorn, S. Disaggregating population data and evaluating the accuracy of modeled high-resolution population distribution—The case study of Germany. Sustainability 2020, 12, 3976. [Google Scholar] [CrossRef]
  45. Wei, C.; Wang, Z.; Lan, X.; Zhang, H.; Fan, M. The spatial-temporal characteristics and dilemmas of sustainable urbanization in China: A new perspective based on the Concept of Five-in-One. Sustainability 2018, 10, 4733. [Google Scholar] [CrossRef] [Green Version]
  46. Department of Industrial Works. Statistics on Number of Factories and Permitted to Operate, Northeastern Region; Ministry of Industry: Bangkok, Thailand, 2014. Available online: http://service.nso.go.th/nso/web/statseries/statseries16.html (accessed on 30 August 2019).
  47. Ronghanam, P.; Bejranade, M. Emerging of low-cost carrier and changing of aviation network, airport hierarchy and relationship between the economics of airport located province. Built Environ. 2019, 18, 51–68. [Google Scholar] [CrossRef]
  48. Office of the National Economic and Social Development Board. Per Capita Income of Population by Region and Province; Office of the Prime Minister: Bangkok, Thailand, 2014. Available online: http://service.nso.go.th/nso/web/service/service03.html (accessed on 30 August 2019).
  49. Chotchaiwong, P.; Wijitkosum, S. Predicting urban expansion and urban land use changes in Nakhon Ratchasima city using A CA-Markov model under two different scenarios. Land 2019, 8, 140. [Google Scholar] [CrossRef] [Green Version]
  50. Bunea, D. Modern gravity models of internal migration. The case of Romania. Theor. Appl. Econ. 2012, 4, 127. [Google Scholar]
  51. Mafi, S.; Abdelrazig, Y.; Amirinia, G.; Kocatepe, A.; Ulak, M.B.; Ozguven, E.E. Investigating exposure of the population to crash injury using a spatiotemporal analysis: A case study in Florida. Appl. Geogr. 2019, 104, 42–55. [Google Scholar] [CrossRef]
Figure 1. Framework used to detect spatial interactions effect of population density patterns.
Figure 1. Framework used to detect spatial interactions effect of population density patterns.
Ijgi 09 00556 g001
Figure 2. Study area.
Figure 2. Study area.
Ijgi 09 00556 g002
Figure 3. Local spatial autocorrelation (local Moran’s I) of the Northeast during 2002–2017.
Figure 3. Local spatial autocorrelation (local Moran’s I) of the Northeast during 2002–2017.
Ijgi 09 00556 g003
Figure 4. Local spatial autocorrelation (Getis–Ord Gi*) of the Northeast during 2002–2017.
Figure 4. Local spatial autocorrelation (Getis–Ord Gi*) of the Northeast during 2002–2017.
Ijgi 09 00556 g004
Figure 5. Sub-district location around the central point of service based on gravity index spatial interaction.
Figure 5. Sub-district location around the central point of service based on gravity index spatial interaction.
Ijgi 09 00556 g005
Figure 6. Correlation coefficient of gravity index and related factors.
Figure 6. Correlation coefficient of gravity index and related factors.
Ijgi 09 00556 g006
Table 1. Global Moran’s I and population density pattern for the period ranging from 2002 to 2017.
Table 1. Global Moran’s I and population density pattern for the period ranging from 2002 to 2017.
YearsMoran’s Ip-ValueZ-ScorePattern
20020.0640.0046.8221Clustered
20070.0770.0028.0045Clustered
20120.0750.0037.8022Clustered
20170.0960.0019.6438Clustered
Table 2. Clustering statistics of local Moran’s I and Getis–Ord Gi* statistics for population density of sub-districts (number of sub-district).
Table 2. Clustering statistics of local Moran’s I and Getis–Ord Gi* statistics for population density of sub-districts (number of sub-district).
YearLocal Moran’s IGetis–Ord Gi*
H-HL-LL-HH-LHotCold
2002974851724122506
20071034821517126500
20121064781617129498
20171134721716135490

Share and Cite

MDPI and ACS Style

Suwanlee, S.R.; Som-ard, J. Spatial Interaction Effect of Population Density Patterns in Sub-Districts of Northeastern Thailand. ISPRS Int. J. Geo-Inf. 2020, 9, 556. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi9090556

AMA Style

Suwanlee SR, Som-ard J. Spatial Interaction Effect of Population Density Patterns in Sub-Districts of Northeastern Thailand. ISPRS International Journal of Geo-Information. 2020; 9(9):556. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi9090556

Chicago/Turabian Style

Suwanlee, Savittri Ratanopad, and Jaturong Som-ard. 2020. "Spatial Interaction Effect of Population Density Patterns in Sub-Districts of Northeastern Thailand" ISPRS International Journal of Geo-Information 9, no. 9: 556. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi9090556

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