Characteristics and Influencing Factors of Spatial Differentiation of Market Service Industries in Rural Areas around Metropolises—A Case Study of Wuhan City’s New Urban Districts

Abstract

Market services industries are closely related to residents’ lives, and its spatial distribution has an important impact on satisfying residents’ consumption needs and promoting economic development. In recent years, with the rapid development of urban–rural integration and the gradual implementation of a rural revitalization strategy, rural areas around metropolises have become the “frontier” of urban geographic expansion, causing the rural market services industries to specialize, commercialize, and modernize. Taking Wuhan as a case study, the spatial distribution characteristics and influencing factors of the market services industries were measured using the average nearest neighbor index, the kernel density estimate, the spatial correlation analysis, and the Geodetector method. The results are as follows. (1) The market services industries in Wuhan’s new urban districts, as a whole and individually, showed characteristics of agglomeration. The market services industries overall formed two high-density areas and multiple agglomeration areas, and the high-density areas of different types of market services industries showed characteristics of partial spatial overlap. (2) There was a significant positive spatial correlation in market services industries, as a whole and individually. Among them, the life services industry had the strongest spatial correlation, while the accommodation services industry had the weakest. (3) Market demand and traffic conditions were the core factors influencing the distribution of the market services industries in Wuhan’s new urban districts, followed by urbanization and economic levels; and tourism conditions had a lower impact. Each detector had a different impact on the spatial distribution of different market services industries, and the interaction research showed that the spatial distribution of the market services industries was the result of a combination of multiple factors. This research provides a future development direction for market service industries in rural areas.

1. Introduction

With the rapid development of globalization and technology information, the industrial structure of China has gradually transformed from “secondary-industry-oriented” to “service-industry-oriented”, and the service industry has increasingly become the engine of regional economic growth. As a branch of the service industry, the distribution and location selection of the market services industries, which are closely related to the daily life of residents, have always been important research fields in geography. ‘Market services industries’ refers to the industries that directly provide residents with material and spiritual services, and they have the characteristics of wide distribution, fierce competition, and complex formats. The spatial distribution not only reflects the law of industry location selection, but also affects the living standards of residents. In recent years, as the relationship between urban and rural areas in China has shifted from rural areas dominated by urban areas to urban–rural integration, the rural areas around metropolises have become node areas with specific complex values that extend from the urban center. Affected by this, the rural economy has begun to transform to specialization, commercialization, and modernization, and various types of market services industries have gradually emerged in rural areas. Studying the spatial distribution of rural market services industries is of vital for optimizing commercial land use, rationally allocating resources, stimulating economic growth, and increasing employment opportunities.

At present, there is a substantial amount of research concerning various types of market services, mainly focusing on the following aspects: location selection and optimization [1,2,3], the recognition and evolution of spatial characteristics [4,5,6,7,8,9,10,11,12], accessibility and spatial equality [7,13,14,15,16,17], influencing factors of spatial distribution [2,18,19,20], policy supply and present allocation, etc. [21,22,23,24,25]. In terms of research scale, the current research covers macro-, meso-, and microgeographic scales, e.g., country, city cluster, province, and city scales, and is carried out based on different geographic statistical units, but most studies are concentrated in economically developed urban areas, and the research on underdeveloped areas and rural areas is relatively lacking. In terms of research methods, spatial analysis methods, neural network analysis models, weighted evaluation methods, and 2SFCA have mainly been used to explore the spatial structure, optimization, and location selection of the market services industries [15,26,27,28,29]; geographic weighted regression, the Geodetector method, the Pearson correlation coefficient matrix, ordinal logistic regression, etc., are used to explore the influencing factors [19,30,31,32,33,34,35,36]. Based on the above studies, we found the following: (1) Research on the market services industries is mostly concentrated on a single type of market services industry, such as retail, catering, finance, leisure and entertainment, and there is relatively little research on the spatial distribution of all market service industries. (2) Most existing research is concentrated in urban areas; research on towns and rural areas is scarce. (3) The data used in the analysis of influencing factors are relatively simple, and most of them are economic data published by statistical departments. It is difficult to accurately characterize the distribution characteristics and laws of the market services industries.

In recent years, big data and geographical conditions census data have provided new opportunities for refined geography research. POI data (point of interest), as a type of big data, represent geographical entities in daily life (shopping malls, convenience stores, hospitals, schools, etc.) and have the characteristics of diversification, high accuracy, fast update speed, and large quantity. Moreover, they can accurately reflect the social and economic activities and their relationship with geographic locations, which is of great significance to the study of human–land relationships [37]. Nowadays, POI data are applied to the division of urban functional areas [38,39,40], city center and boundary identification [41,42], industry cluster analysis [43,44], points of interest recommendations [45,46], etc., but the refined research applied to the rural market services industries is rare. Geographical condition census data include not only detailed spatial data, but also a large amount of nonspatial data. The data are authoritative and accurate, providing a reliable basis for the research on the spatial distribution and influencing factors of the rural market services industries in terms of microlocation and microscale.

Regarding rural areas surrounding metropolises, the development level of the market services industries in Wuhan’s new urban districts is relatively lagging. A scientific understanding of the spatial distribution of the market services industries and their influencing factors plays an important role in optimizing spatial distribution and improving the comprehensive level of the market services industries. Based on POI data and geographical condition census data, this article uses nearest neighbor index, kernel density, and spatial autocorrelation methods to analyze the characteristics of the market services industries in Wuhan’s new urban districts. Based on this, the Geodetector method is used to explore the factors affecting the spatial distribution of the market services industries in Wuhan’s new urban districts. This research provides a future development direction for the market services industries in rural areas. Moreover, it has important theoretical and practical significance for understanding the spatial distribution rules and influence mechanism of the market services industries in rural areas in China.

This article is organized as follows (Figure 1). The first section explains the significance and goal of studying the rural market services industries. The second section introduces the research area, the data source, and the methodology. The third section analyzes the characteristics and influencing factors of the spatial distribution of the rural market services industries in Wuhan’s new urban districts. Based on the previous research, the fourth section discusses the practical meaning of these findings and states the limitations of our research. The last section contains the conclusions and discussion.

2. Materials and Methods

2.1. Research Area

Wuhan is located in center of the hinterland of China, the eastern part of Hubei Province, and at the intersection of the Yangtze River and Hanjiang River. It is the central city of central China, the capital of Hubei province, a national historic and cultural city, and one of the 14 megacities in China. Wuhan has 13 municipal districts, including seven central urban districts and six new urban districts. This article studies the current situation of the market services industries in rural area, and the central urban districts are mainly economically developed urban built-up areas; therefore, the research area comprises new urban districts dominated by rural areas, including Huangpi District, Xinzhou District, and Jiangxia District, Dongxihu District, Hannan District, and Caidian District (Figure 2). The new urban districts consist of 2206 administrative villages with a total area of 7614 km². As of 2019, the permanent population of the new urban districts was 4.177 million, of which 2.166 million was the agriculture population, which accounts for 55.95% of the total population. The GDP of Wuhan city was 636.72 billion yuan, of which the primary industry was 36.718 billion yuan, the secondary industry was 342.543 billion yuan, and the tertiary industry was 257.459 billion yuan. In the whole year, the total retail sales of consumer goods in Wuhan’s new urban districts amounted to 239.048 billion yuan, and the disposable income was 21,462 yuan. Wuhan city has a solid foundation in the market services industries, especially after entering the 21st century. Wuhan’s market services industries have made great progress and have become an important economic pillar. The rural area of new urban districts is vast, and it is obviously driven by the radiation of the central urban area. Therefore, selecting the new urban districts as a typical case study has representative significance for studying the development of the market services industries in the rural areas surrounding the metropolis.

2.2. Data Source

The data used in this study mainly included POI data on the market services industries in Wuhan, data on basic geographic information in Wuhan, data on mobile phone signaling and data on socioeconomic development. Data details are as follows.

2.2.1. POI Data on Market Services Industries in Wuhan

This dataset was derived from the Baidu Maps website (http://lbsyun.baidu.com/, accessed on 31 March 2019). The data of the market services industries in Wuhan city’s new urban districts mainly included the ID, name, address, type, administrative division, and coordinates. Referring to the “Industrial Classification for National Economic Activities” established by the People’s Republic of China, after a series of preprocessing, e.g., cleaning, double checking, and reclassification of the original data, as of 31 March 2019, we obtained a total of 40,635 market services industries’ POI data and divided them into six types of market services industries: shopping, life, catering, leisure, financial, and accommodation services (Figure 3,

Table 1. The classification and statistics of market services industries’ POI data in Wuhan city’s new urban districts.

Type	POI Data Content	Number	Percentage
Shopping services	Convenience store, supermarket, food market, specialty store, department store, building materials market	7922	19.49%
Life services	Beauty salon, photography shop, repair shop, laundry shop, mobile business hall, post office, express station	11,644	28.65%
Catering services	Chinese restaurant, breakfast restaurant, casual restaurant, western restaurant, fast food restaurant, dessert shop	14,067	34.62%
Leisure services	KTV, sanatorium, gymnasium, chess room, billiards room, Internet cafe, theater, game center, scenic spot	2813	6.93%
Financial services	Insurance company, bank, securities company	1756	4.33%
Accommodation services	Apartment hotel, economy hotel, hostel, guest house, starred hotel	2433	5.98%

).

2.2.2. Data on Basic Geographic Conditions in Wuhan

The data on basic geographic conditions came from the Wuhan Geomatics Institute. According to the basic geographic conditions data, the total number of the administrative villages identified in Wuhan city’s new urban districts was 2206 in 2019, and details are stored in text and table formats. In addition, the data also included land, population, traffic conditions, etc. There are technical and procedural regulations concerning every step of data collection, including filed surveys, data checking, data entry into the database, and statistical analysis. Data of the basic geographic information of Wuhan are kept in the format of ArcGIS vector data. The obtained data of basic geographic conditions were geospatially matched with the data of market services industries in Wuhan.

Table 2. List of Wuhan geographical conditions census data used in this study.

Category	Feature	Attributes Used
Land coverage	Built-up area data	Location, area
Road networks	Centerline of city road Centerline of country road and rural road Bus station	Type, length, width, road grade, location
Geographic units	District (county) boundaries Village boundaries	Location, area
Thematic data	Survey of population	population

presents the basic geographical condition data used in this study.

2.2.3. Mobile Phone Signaling Data

The mobile phone signaling data were provided by China Unicom Company in June 2019, including the MSIID code, date, location, mobile phone number attribution, and stay time. The data regarding the employed population, local tourists, and nonlocal tourists were obtained through mobile phone signaling data. If mobile phone users were connected to the same base station for more than 4 h between 08:00 and 17:00 for 3 consecutive working days, they were identified as an employed population [47]; if mobile users received signals of more than 1 h from different base stations within the same scenic spot three times in a day, they were identified as tourists [48,49], and local and nonlocal tourists were distinguished by mobile phone number attribution.

2.2.4. Data on Social and Economic Development

This dataset was derived from the Wuhan Statistical Yearbook 2019, Wuhan Census of the Geographical Conditions Report 2019. Permanent resident population data are from the survey of “Actual Housing and Actual Population” in Wuhan in 2019.

2.3. Research Methods

2.3.1. Average Nearest Neighbor Index

The average nearest neighbor index is an index that indicates the mutual proximity of point elements in geographic space, and can well-reflect the spatial distribution characteristics [50]. In this study, the average nearest neighbor index was used to determine the distribution of types of market services industry in Wuhan’s new urban districts. The index value was calculated as follows:

$$R=2\left[{\displaystyle \sum _{i=1}^N\frac{MIN\left(d_{ij}\right)}{N}}\right]/\sqrt{\frac{A}{N}}$$

(1)

where $R$ is the nearest neighbor index of market services industries, $MIN\left(d_{ij}\right)$ is the theoretical nearest distance of market services industry points, $A$ is the area of administrative region in Wuhan’s new urban districts, and $N$ is the number of market services industry points. $R=1$ indicates that market services industry points generally show a random distribution, whereas $R<1$ indicates an agglomerated distribution, and $R>1$ indicates a uniform distribution.

2.3.2. Kernel Density Estimation

Although the nearest neighbor index can reflect the spatial distribution type of market services industries, it fails to observe their spatial dispersion or agglomeration characteristics. Kernel density estimation is a method to study the probability of occurrence of points at different locations in space [51]. In this study, kernel density estimation was used to identify the spatial distribution density of the market services industry points as follows:

$$f_n(x)=\frac{1}{n{r}^d}{\displaystyle \sum _{i=1}^{n}K\left(\frac{x-x_i}{r}\right)}$$

(2)

where $K\left(\frac{x-x_i}{r}\right)$ is the kernel density function of market services industries, $x$ is the location of market services industry points, $x_i$ is the market services industry points set to be estimated, $r$ is the bandwidth for the search radius distance ($r>0$), $n$ is the number of market services industry points in the range of bandwidth, and $d$ is the dimensionality of market services industries.

2.3.3. Spatial Autocorrelation Analysis

Spatial autocorrelation analysis is used to explore the spatial agglomeration and correlation of the number of market services industries in different administrative villages [52]. In this study, Moran’s I index is used to analyze the global spatial correlation characteristics, and local Moran’s I index is used to analyze the spatial correlation degree between area $i$ and surrounding areas.

Moran’s I index is used to measure the spatial correlation model, which is obtained by the following formula [53]:

$$I={\displaystyle {\sum }_{i=1}^n}{\displaystyle {\sum }_{j=1}^n{w}_{ij}\left(x_i-\overline{x}\right)/S^2\cdot {\displaystyle {\sum }_{i=1}^n}}{\displaystyle {\sum }_{j=1}^n{w}_{ij}}$$

(3)

$$S^2=\frac{1}{n}\cdot {\displaystyle {\sum }_{i=1}^n\left(x_i-\overline{x}\right)}$$

(4)

where $x_i$ and $x_j$ are the number of rural market services industries in the geographic units of area $i$ and $j$, and $w_{ij}$ is the spatial weight matrix. Moran’s I > 0 indicates a positive spatial correlation, which means that market services industries in Wuhan’s new urban districts are geospatially agglomerated, and the larger the value, the more obvious agglomeration of market services industries; Moran’s I < 0 indicates a negative spatial correlation, which means that rural market services industries in Wuhan’s new urban districts exhibit significant spatial differences in geospace; and Moran’s I = 0 indicates that the distribution of market services industries in Wuhan’s new urban districts is random.

The local Moran’s I index for an area $i$ measures the association between a value at $i$ and values of its nearby areas, defined as [54]:

$$I_i=\frac{\left(x_i-\overline{x}\right)}{s_x^2}{\displaystyle {\sum }_j\left[w_{ij}\left(x_j-\overline{x}\right)\right]}$$

(5)

where $s_x^2={\displaystyle {\sum }_j{\left(x_j-\overline{x}\right)}^2/n}$ is the variance, and other notations are the same as in Formula (3). Note that the summation over $j$ does not include the area $i$ itself, that is, $j\ne i$. A positive $I_i$ means either a high value surrounded by high values (high-high) or a low value surrounded by low values (low-low). A negative $I_i$ means either a low value surrounded by high values (low–high) or a high value surrounded by low values (high–low).

2.3.4. The Geodetector Method

Geodetector is a new statistical method for detecting spatial differentiation and revealing the driving factors [55].

• Factor Detection

A factor detector detects the spatial heterogeneity of the dependent variable and the power of determinant of the independent variables on the dependent variable. The value is measured by the $q$ value.

$$q=1-\frac{{\displaystyle {\sum }_{h=1}^L{N}_h{\sigma }_h^2}}{N{\sigma }^2}=1-\frac{SSW}{SST}$$

(6)

where $h=1,\dots ,L$ is the layer of independent variable $X$, $N_h$ and $N$ are the number of sample units in layer $h$ and the total region $Y$, respectively. Moreover, ${\sigma }_h^2$ and ${\sigma }_{}^2$ are the variance in the $h$ layer and the variance in the region. The $q$ value lies in [0, 1]. The larger the value of $q$, the stronger the explanatory power of this factor on the differentiation of market services industries, and vice versa. If factor $X$ completely controls the distribution of market services industries, the $q$ value equals 1; if the factor $X$ is completely unrelated to the distribution of market services industries, the $q$ value equals 0.

2.

Interaction Detection

An interaction detector is mainly used to identify the interaction between different factors, that is, to evaluate whether the two factors will increase or decrease the explanatory power of market services industries or whether the influences of such factors on market services industries are independent of each other. The classification of related interaction types is shown in

Table 3. Types of interaction between two covariates.

Description	Interaction
$q\left(X_1\cap X_2\right)<Min\left(q\left(X_1\right),q\left(X_2\right)\right)$	Weaken, nonlinear
$Min\left(q\left(X_1\right),q\left(X_2\right)\right)<q\left(X_1\cap X_2\right)<Max\left(q\left(X_1\right),q\left(X_2\right)\right)$	Weaken, single factor nonlinear
$q\left(X_1\cap X_2\right)>Max\left(q\left(X_1\right),q\left(X_2\right)\right)$	Enhanced, double factors
$q\left(X_1\cap X_2\right)=q\left(X_1\right)+q\left(X_2\right)$	Independent
$q\left(X_1\cap X_2\right)>q\left(X_1\right)+q\left(X_2\right)$	Enhanced, nonlinear

.

3.

Construction of the Detection Indicator System

The spatial distribution of the market services industries is affected by many factors such as economy, society, infrastructure, and policy, and the impact power of different factors is also different. Most existing studies consider the built-up area, population size, GDP, urbanization level, etc., all to be factors that affect the spatial distribution of market services industries [33,34,35,36]. Taking into account the actual conditions of Wuhan city’s new urban districts and the data available based on the degree of difficulty, 10 indicators were selected as factors influencing the spatial distribution characteristics of the market services industries based on the principles of scientificity, representativeness, and accessibility (

Table 4. The detection indicator system for influencing factors.

System	Detection Factor	Detection Indicator	Remark
Economic level	X1	GDP per capita (yuan)	GDP/permanent population of administrative villages
	X2	Disposable income (yuan)	/
Urbanization level	X3	Built-up area ratio (%)	Built-up area/administrative village area
Market demand	X4	Permanent population (person)	Permanent population of administrative villages
	X5	Employed population (person)	Mobile phone signaling data in June 2019
	X6	Population density (person/m2)	Permanent population/administrative village area
Transportation conditions	X7	Road density (m/km2)	Total mileage of administrative village Road/area of administrative village
	X8	Coverage percentage of the public transport (%)	Buffer area of 500 m from the bus station/area of administrative village
Tourism conditions	X9	Number of local visitors (person)	Mobile phone signaling data in June 2019
	X10	Number of nonlocal visitors (person)	Mobile phone signaling data in June 2019

).

3. Results

3.1. Spatial Differentiation Characteristics

3.1.1. Spatial Differentiation Characteristics of All Market Services Industries

The nearest neighbor index of market services industries in Wuhan’s new urban districts was 0.1998 and passed the significance test, indicating an agglomerated distribution. The spatial agglomeration characteristics of market services industries are discussed through kernel density mapping. Figure 4 illustrates that the kernel density of market services industries in Wuhan’s new urban districts had two high-density areas and multiple agglomeration areas. The high-density values were mainly distributed in the southeast of Dongxihu District and the north of Jiangxia District, such as Wenjiawan Community and Kaipingli Community. The kernel density of the center of this area was 233.95–481.09/km²; the agglomeration areas were located in central Huangpi District, northern Jiangxia and Caidian District, and southern and central Xinzhou District, such as Yangyuan Community and Panjiatian Community. The kernel density of the center of this area was 90.57–233.94/km². The market services industries in most other areas were relatively sparse. In general, the high-density areas and agglomeration areas were either adjacent to the central urban districts or located in the administrative and economic centers of new districts. The construction of infrastructure in these areas was relatively complete, and these areas gathered a large population and were greatly affected by the radiation of the central urban districts.

Moran’s I index of the market services industries was 0.5646 > 0, indicating that market services industries were spatially positively correlated. Taking the administrative village as a unit, the number of market services industries in each administrative village was counted, and the local Moran’s I was calculated. The results show that the number of H-H areas was 183, accounting for 8.3% of the total; these 183 H-H areas were basically distributed in the east of Dongxihu District, the south of Huangpi District, the north of Jiangxia District, and the northeast of Hannan District. The number of L-L areas was 219, accounting for 9.93% of the total (the highest number); these 219 L-L areas were widely distributed along the edge of the new urban districts. The number of L-H areas was 42, and they were mainly distributed around the H-H areas, accounting for 1.90% of the total. The number of H-L areas was five (the lowest number), and they were scattered across Huangpi District, Xinzhou District and Jiangxia District, only accounting for 0.23% of the total. Furthermore, 79.65% of the areas were not significant, thus showing a trend of random distribution (Figure 5).

3.1.2. Spatial Differentiation Characteristics for Different Types of Market Services Industries

The nearest neighbor index values of different types of market services industries are summarized in

Table 5. The nearest neighbor index of different types of market services industries in Wuhan city’s new urban districts.

Type	Average Observation Distance (m)	Expected Mean Distance (m)	Nearest Neighbor Index	p Values	Type of Spatial Distribution
Shopping service	101.844	595.808	0.171	0.00	agglomeration
Life service	68.832	488.148	0.141	0.00	agglomeration
Catering service	58.171	459.814	0.126	0.00	agglomeration
Leisure service	308.179	1023.137	0.301	0.00	agglomeration
Financial service	194.349	1238.935	0.157	0.00	agglomeration
Accommodation service	212.453	1068.461	0.199	0.00	agglomeration

. The ascending order of the $R$ values of market services industries was catering services (0.126), life services (0.141), financial services (0.157), shopping services (0.171), accommodation services (0.199), and leisure services (0.301). The results were all less than 1, with different levels of agglomeration. Among them, the catering services industry had the smallest $R$ value, with the highest degree of agglomeration, and the leisure services industry had the lowest degree of spatial agglomeration.

The kernel density result showed that different types of market services industries in Wuhan’s new urban districts had formed kernel centers with different densities (Figure 6). Further analysis showed the following: (1) The high-density areas of each type of market services industry presented the characteristics of partial spatial overlap, and they were concentrated in the east of Dongxihu District, the north of Jiangxia District, and the center of Huangpi District. (2) From the perspective of the kernel density levels of different market services industries, the kernel density levels of the life and catering services industries had comparative advantages, while the density levels of the leisure, financial, and accommodation services industries were much lower than those of other market services industries. The kernel density of the shopping market services industry was at an intermediate level. (3) The areas with higher kernel density values for shopping, life, and financial services were more widely distributed than the other three types, which is related to the nature of the services industries and China’s policy of promoting financial service in rural areas. (4) In terms of regional differences, no administrative district formed a complete kernel density grading system in all types of market services industries. There were three types of market services industries that had formed a complete kernel density rating system in Dongxihu District; namely, the catering, leisure, and financial services industries. Four types of market services industries formed a complete kernel density rating system in Jiangxia District. However, the kernel density grading system of the remaining four districts is incomplete.

Moran’s I index of all types of market services industries had positive index values in the range of 0.1861 to 0.6036, and passed the significant test at the 5% level, indicating positive spatial correlation. The largest value of Moran’s I index was 0.6036 for the life services industry, which had the most obvious spatial correlation, while the smallest value of Moran’s I index was 0.1861 for the accommodation services industry, which had the weakest spatial correlation (

Table 6. The Moran’s I index of different types of market services industries in Wuhan city’s new urban districts.

Type	Moran’s I Index	Expectations Index	Variance	Z Score
Shopping services	0.5292	−0.0004	0.0001	62.7768
Life services	0.6306	−0.0004	0.0001	74.5563
Catering services	0.5015	−0.0005	0.0001	59.7953
Leisure services	0.3377	−0.0005	0.0001	40.8470
Financial services	0.5423	−0.0005	0.0001	64.3360
Accommodation services	0.1861	−0.0005	0.0001	23.0164

).

Combining the LISA map (Figure 7), we found the following: (1) The spatial correlation patterns of different types of market services industries were mainly positively correlated. The largest proportion of the H-H areas was the life services industry (10.11%), and the smallest proportion was the accommodation services industry (5.3%). The largest proportion of the L-L areas was the financial services industry (4.4%), and the smallest proportion was the catering services industry (2.63%). (2) The distributions of spatial clustering in the different types were basically the same; i.e., the H-H and L-L areas were agglomerated, while the H-L and L-H areas were scattered. (3) The spatial clustering of shopping, life, and catering services were distributed across all districts, while the spatial clustering of leisure, finance, and accommodation services were mainly concentrated in the central and southern parts of the new districts, and the northern part was not significant.

3.2. Factors Influencing Spatial Characteristics

3.2.1. Factor Detection Results

(1)

Factors Influencing the Distribution of All Market Service Industries

The distribution maps of influencing factors are displayed in Figure 8. According to the calculation using the Geodetector model, each detection factor had a positive impact on the spatial distribution of the market services industries, and a descending order of influence power was employed population (0.492), permanent population (0.467), population density (0.461), coverage percentage of the public transport (0.423), road density (0.338), proportion of built-up area (0.291), GDP per capita (0.183), disposable income (0.181), local visitors (0.164), and nonlocal visitors (0.155). The influence force was ranked by intensity in the order of market demand > transportation conditions > urbanization level > economic level > tourism conditions.

Market demand had the strongest explanatory force for the distribution of the market services industries, indicating that this was the core factor influencing their spatial distribution. The market services industries were generally located around the residential area and business area, forming a 15-min life circle and providing people with products and services to meet their material and spiritual needs. At the same time, the market services industries also provided employment. For example, the core business district of Wenyuan Community in Jiangxia District has concentrated large-scale corporate companies such as COFCO and office buildings, which incubates important market service industries, e.g., Zhongbai Plaza, Wanhui Mall, Marriott Century Day Shopping Mall, and Walmart. In summary, market demand was the most important factor affecting the spatial distribution of the market services industries. Whether it was a residential-oriented area or an employment-oriented area, the market services industries were denser, and the development was more mature.

The explanatory power of transportation conditions was relatively high, indicating that the higher the traffic accessibility, the denser the market services industries. To further explain the spatial relationship between transportation and the distribution of the market services industries, the ArcGIS 10.2 proximity analysis tool was used to draw a 300-m buffer zone of bus stations and main roads (including township roads, county roads, and national roads), and the spatial join tool was used to count the number of market services industries in the buffer zone (Figure 9). The results showed that 38,528 market services industries were located in the 300-m buffer zone of bus stations, accounting for 94.81% of the total, and 39,958 market services industries were located in the 300-m buffer zone of main roads, accounting for 98.33% of the total. Good transportation conditions promote the circulation and connection of various factors, and reduce time costs for consumers. At the same time, market services industries can also expand their service radius with good traffic conditions, increasing the number of consumers and obtaining greater economic benefits. Additionally, the explanatory power of public transportation coverage was slightly higher than that of road density. This was due to the limited number of cars owned by households in Wuhan’s rural areas and because public transportation is the main mode of travel. With the implementation of the rural revitalization strategy and “the 13th Five-Year Plan for Modern Comprehensive Transportation System Development”, Wuhan’s new urban districts had fully completed the project goal of “Public Transport Covers All Villages”, which can meet the basic consumption and travel needs of rural residents.

The economic level and urbanization level affected the spatial distribution of the market services industries to a certain extent. The market services industries, as service-oriented industries that directly face service targets, have a higher need for agglomerated customers. A high level of urbanization indicates a high quality of urban life, which guides the rapid flow of population across districts and attracts more immigrants. The low dependence degree of the market services industries on the economic level showed that it did not overly affect their spatial distribution. The economic level mostly affects the level of the market services industries, but the impact on the quantity is not significant.

Tourism conditions had the lowest explanatory power. This is because the research area mainly comprises the rural population, whose tourism awareness is relatively weak and whose daily life consists of agricultural production activities. In addition, except for Huangpi District, the districts lack resources to attract tourists, so their influence on the spatial distribution of market service industries is limited.

(2)

Factors Influencing the Distribution of Different Types of Market Services Industries

Geodetector was used to measure the explanation power of detection factors on the spatial distribution of different types of market services industries, and the results are shown in

Table 7. Factor detection results of factors influencing the spatial characteristics of different types of market services industries in Wuhan city’s new urban districts.

Type	${\mathit{X}}_{\mathbf{1}}$	${\mathit{X}}_{\mathbf{2}}$	${\mathit{X}}_{\mathbf{3}}$	${\mathit{X}}_{\mathbf{4}}$	${\mathit{X}}_{\mathbf{5}}$	${\mathit{X}}_{\mathbf{6}}$	${\mathit{X}}_{\mathbf{7}}$	${\mathit{X}}_{\mathbf{8}}$	${\mathit{X}}_{\mathbf{9}}$	${\mathit{X}}_{\mathbf{10}}$
Shopping services	0.1623 ***	0.1592 *	0.2263 **	0.4276 ***	0.4375 ***	0.4218 ***	0.3957 **	0.4185 ***	0.1858	0.1575 ***
Life services	0.1641 ***	0.1679 **	0.1551 *	0.5112 ***	0.5341 ***	0.4328 ***	0.3844 *	0.4169 **	0.1732 **	0.1424 **
Catering services	0.1698 ***	0.1780 ***	0.2891 *	0.3935 ***	0.4254 ***	0.4125 ***	0.3168	0.4032	0.1782 **	0.1468 *
Leisure services	0.1615 ***	0.1782 **	0.1639 *	0.3664 ***	0.4167 ***	0.4073 ***	0.3180	0.4012	0.1920 **	0.1632 **
Financial services	0.2023 ***	0.1821 ***	0.1653 *	0.4249 ***	0.4389 ***	0.4204 ***	0.3841 **	0.4079 **	0.1709 ***	0.1217
Accommodation services	0.1525 ***	0.1480 *	0.2920 ***	0.3590 ***	0.4029 ***	0.4050 ***	0.3188	0.4037 *	0.1841 *	0.1600 **

Note: ***, ** and * indicate statistical significance at the 1%, 5%, and 10% level, respectively.

. The influence of each factor on the spatial differentiation of various market services industries is different. The permanent population and population density were the main reasons for the spatial distribution of the life services industry, which mainly provided residents with services closely related to their daily lives. The permanent population and population density reflect the market scale and determines the economic benefits of the life services industry. The employed population had the greatest explanatory power for life and financial services, indicating that employed population had a greater demand for these types of market services. At the same time, the life and financial service industries had also created relatively more jobs that attracted an influx of the employed population, which led to the agglomeration of market services industries.

Per capita GDP and disposable income had the greatest explanatory power for the spatial distribution of the financial services industry. Statistical results showed that 68.28% of the financial services industries were located in the administrative villages with a per capita GDP of over 100,000 yuan. The bank loans provided by the financial services industry can meet the liquidity needs of companies, which, in turn, can be fed back to the financial services industry. In total, 39.97% of the financial services industries were distributed across the administrative villages with a disposable income of more than 15,000 yuan. The income level determines the ability of rural residents to pursue higher-level needs. When people’s income reaches a certain level, they are willing to purchase commercial insurance and manage their finances to meet safety requirements.

Road density and public station coverage had the greatest explanatory power for the spatial distribution of the shopping services industry, followed by the life services industry. Compared with the shopping services industry, the life services industry generally has a smaller service radius. Rural residents can arrive at life services points by walking or by electric bicycle. For higher-level shopping needs, travel distances are generally longer, and shopping activities mainly rely on public transport and roads for their service functions.

The built-up area ratio had the greatest explanatory power for the spatial distribution of the accommodation service industry. This is because, whether it is run in other building’s functional areas or independently operated in detached building, the development of the accommodation service industry depends on regional infrastructure facilities, i.e., the water and power supply, sewage treatment, fire protection, and other facilities. In areas where the built-up area ratio is relatively high, a complete infrastructure can greatly reduce the operating costs of the accommodation services industry and can lead to greater operating benefits.

Both local tourists and nonlocal tourists had the greatest explanatory power for the spatial distribution of the leisure services industry. As tourism has become a part of people’s lifestyles, tourist destinations are not limited to scenic spots; they include the daily leisure and entertainment spaces of residents. The leisure activities of tourists and local residents have been integrated, which has promoted the development of the leisure service industry. The explanatory power of nonlocal tourists was higher than that of local tourists, because most of the local tourists are short-distance, single-day travelers, and very few are rural leisure tourists that require accommodation during the holidays.

3.2.2. Interaction Detection Results

The interaction detection results are shown in

Table 8. Interaction detection results of factors influencing the spatial characteristics of different types of market services industries in Wuhan city’s new urban districts.

Type	Rank 1			Rank 2			Rank 3
	${\mathit{X}}_{\mathit{i}}\mathbf{\cap }{\mathit{X}}_{\mathit{j}}$	$\mathit{q}\mathbf{\left(}{\mathit{X}}_{\mathit{i}}\mathbf{\cap }{\mathit{X}}_{\mathit{j}}\mathbf{\right)}$	Interaction Type	${\mathit{X}}_{\mathit{i}}\mathbf{\cap }{\mathit{X}}_{\mathit{j}}$	$\mathit{q}\mathbf{\left(}{\mathit{X}}_{\mathit{i}}\mathbf{\cap }{\mathit{X}}_{\mathit{j}}\mathbf{\right)}$	Interaction Type	${\mathit{X}}_{\mathit{i}}\mathbf{\cap }{\mathit{X}}_{\mathit{j}}$	$\mathit{q}\mathbf{\left(}{\mathit{X}}_{\mathit{i}}\mathbf{\cap }{\mathit{X}}_{\mathit{j}}\mathbf{\right)}$	Interaction Type
Market services	$X_5\cap X_6$	0.9582	Enhanced, double factors	$X_5\cap X_9$	0.8559	Enhanced, nonlinear	$X_3\cap X_5$	0.8457	Enhanced, nonlinear
Shopping services	$X_5\cap X_6$	0.8592	Enhanced, double factors	$X_5\cap X_7$	0.8579	Enhanced, nonlinear	$X_5\cap X_8$	0.8572	Enhanced, nonlinear
Life services	$X_5\cap X_6$	0.8598	Enhanced, double factors	$X_3\cap X_5$	0.7589	Enhanced, nonlinear	$X_5\cap X_9$	0.7584	Enhanced, nonlinear
Catering services	$X_5\cap X_6$	0.8353	Enhanced, double factors	$X_5\cap X_9$	0.8251	Enhanced, nonlinear	$X_3\cap X_5$	0.8050	Enhanced, nonlinear
Leisure services	$X_4\cap X_6$	0.8346	Enhanced, nonlinear	$X_5\cap X_6$	0.8338	Enhanced, nonlinear	$X_4\cap X_5$	0.7328	Enhanced, double factors
Financial services	$X_3\cap X_5$	0.8510	Enhanced, nonlinear	$X_5\cap X_6$	0.8506	Enhanced, double factors	$X_5\cap X_9$	0.8497	Enhanced, nonlinear
Accommodation services	$X_5\cap X_6$	0.8301	Enhanced, nonlinear	$X_4\cap X_6$	0.8280	Enhanced, nonlinear	$X_3\cap X_5$	0.8269	Enhanced, nonlinear

(due to limited space, only the top three interaction values are listed). The following was found: (1) The interaction detection results of the market services industries, as a whole and individually, were nonlinear-enhanced or double-factor-enhanced, indicating that the influence of interactions on the spatial distribution of market services industries is greater than that of single factors. (2) The interaction factors that ranked first in each market services industry further determined that the market demand and urbanization level were core factors influencing its spatial distribution. (3) The explanatory power of factors with a smaller $q$ value was greatly improved after the interaction, e.g., $X_5\cap X_7$ in the shopping services industry and $X_5\cap X_9$ in the catering services industry.

The interaction detection results of the market services industries, either as a whole or individually, did not show a nonlinear weakening, a single-factor nonlinear weakening, or independence, indicating that the spatial distribution of the market services industries in Wuhan’s new urban districts was not caused by a single factor, and was the result of the combined effect of different influencing factors.

4. Discussion

4.1. The Spatial Distribution Characteristics of Market Services Industries in the Rural Area Surrounding the Metropolis and Nonmetropolis Are Different

Due to the relatively slow development of nonmetropolitan areas, the population and market services industries are still concentrated in central urban districts. However, because of rising rents and land prices, as well as increasing traffic and environmental pressures in metropolises, surrounding rural areas have become important spillover areas for functions, industries, and populations of central urban districts, relying on their superior geographic locations. As one of the metropolises in China, Wuhan city is in an advanced stage of urban development, with a strong urbanization drive and a strong land demand. In the process of spatial structural adjustments of central urban districts, some industries, populations, and functions have been transferred to surrounding rural areas, leading to upgrades to the industrial structure, optimization of the spatial structure, and transformation of land use in the surrounding rural areas, thereby promoting the agglomeration and development of market services industries.

4.2. Spatial Characteristics Promote Top-Level System Design of the Market Services Industries

Kernel density estimation and spatial autocorrelation can discern the development status and agglomeration patterns of different types of market services industries, and analyze the spatial distribution characteristics of market services industries in Wuhan’s new urban districts. At present, the spatial distribution of the market services industries presents a distinct center–periphery structure, but due to the limited diffusion, the spatial distribution of the market service industry is uneven. From the top-level system design, the Wuhan Municipal Government needs to cultivate new market services industry agglomeration centers, especially in the edge areas; be the fulcrum of service-industry clusters to drive the development of surrounding areas; achieve the strategic goal of coordinated and balanced development of the market services industry; and realize the shaping of market services industries on rural spatial structure.

4.3. Influencing Factors Are an Important Basis for Formulating the Development Policy of Market Services Industries

As a surrounding area of a metropolis, Wuhan’s new urban districts are in a period of rapid urbanization and economic development. Providing high-quality and convenient market services for rural residents is an important issue in the process of rural revitalization. This research studies the factors influencing the spatial distribution of market services industries in Wuhan’s new urban districts, which provides an important basis for the government to formulate development policies for the market services industries. According to the results of this research, the Wuhan Municipal Government can optimize the spatial distribution of market services industries in three ways. ① Strengthen the radiation scope of the subcities of Optics Valley, Chegu, Airport Economic Zone, and Yangtze River New Area. Based on the degree of population agglomeration, a “15-min life circle” equipped with daily life needs such as shopping, life, catering, leisure, and other services should be created. At the same time, the convenience of market services can also be used as an attraction, as it will benefit most of the population. ② Construct an efficient and convenient rural road network that is compatible with the land-use layout, accelerate the construction of the subway rail network, and improve the coverage rate of rail transit. ③ Strengthen the construction of rural water conservancy, energy, information, and other municipal infrastructures, improve planning and construction standards, and address the shortcomings of rural facilities.

4.4. The Application of Big Data Removes the Need for Reliance on Traditional Statistical Data

In this research, big data, such as POI data and mobile phone signaling data, were applied to the study of spatial distribution and influencing factors of the rural market services industries, removing the need for reliance on traditional statistical data, and providing new opportunities for the spatial analysis of socioeconomic phenomena in the field of humanities and economic geography. The combination of big data and statistical yearbook data reflects the cross-application of big data and traditional data. It not only provides new solutions for the analysis of socio-economic phenomena and laws in the field of human–economic geography, but also improves the accuracy of research on rural geography.

4.5. Limitations of the Study

Some limitations exist in our study. For example, this research only uses 2019 data to analyze the spatial pattern and influencing factors of the rural market services industries in Wuhan’s new urban districts, and lacks research on its temporal and spatial laws. With the further implementation of the rural revitalization strategy and the reform of the household registration system, the vitality and attractiveness of rural areas have been stimulated, which will inevitably bring new opportunities and challenges to rural development. The existing rural market services industries are likely to conflict with the growing needs of rural residents. Future research needs to obtain dynamic data to predict and simulate the future development of rural market services.

Moreover, due to the lack of rural statistical data, many influencing factors have not been adopted. For example, market demand is the dominant factor influencing the spatial distribution of market services industries. In the future, educational level, gender structure, occupation, etc., can be included in the index system to explore the specific impact mechanism on the market services industries.

5. Conclusions

This study analyzed the characteristics and influencing mechanism of the spatial distribution of rural market services industries in Wuhan’s new urban districts by using different methods such as ArcGIS spatial analysis and Geodetector. The main conclusions are as follows:

(1)

The market services industries in Wuhan’s new urban districts, as a whole and individually, showed characteristics of agglomeration, but the degree of agglomeration was different. In terms of industry type, the catering services industry had the highest degree of spatial agglomeration, and the leisure services industry had the lowest. The market services industries, as a whole and individually, had formed kernel centers with different levels of density. The high-density areas of each type of market services industry presented the characteristics of partial spatial overlap, and the high-density areas of each market services industry were concentrated in the east of Dongxihu District, the north of Jiangxia District, and the center of Huangpi District. No administrative district had formed a complete kernel density grading system in all types of market services industries.

(2)

There was a significant positive spatial correlation in the market services industries, as a whole and individually. Among them, the life services industry had the strongest spatial correlation, while the accommodation services industry had the weakest spatial correlation. The distributions of spatial clustering in different types were basically the same, i.e., the H-H and L-L areas were agglomerated, while the H-L and L-H areas were scattered. The spatial clustering of shopping, life, and catering services was distributed across all districts, while the spatial clustering of leisure, finance, and accommodation services was mainly concentrated in the central and southern parts of the new districts, and the northern part was not significant.

(3)

Market demand and traffic conditions were the core factors influencing the distribution of the market services industries in Wuhan’s new urban districts, followed by urbanization and economic levels, and tourism conditions had a lower impact. Each factor had a different impact on the spatial distribution of different market services industries. The interaction detection results of the market services industries, as a whole and individually, were nonlinear-enhanced or double-factor-enhanced, indicating that the influence of interactions on the spatial distribution of market services industries was greater than that of a single factor.