Spatial Differentiation and Influencing Factors of Tertiary Industry in the Pearl River Delta Urban Agglomeration

: The tertiary industry has become the main driving force for China’s economic development, and the adjustment and optimization of its structure are important prerequisites for achieving high-quality economic development. Existing studies have mostly focused on the spatial layout and influencing factors of the tertiary industry, with insufficient exploration of its internal structure. In this study, the PRD urban agglomeration is selected as the study area. On the basis of classifying the tertiary industry, the Dagum Gini coefficient, kernel density estimation, and local spatial autocorrelation are used to explore the spatial differentiation of various tertiary industries. The influencing factors are analyzed using geographical detectors, and suggestions for future development strategies are proposed. The results show that in terms of regional differentiation, the agglomeration of various tertiary industries in Guangzhou and Shenzhen is the most significant, but there is insufficient spillover to surrounding cities. In terms of development structure, the level of agglomeration of the consumptive tertiary industry is higher, the public tertiary industry tends to be more evenly distributed, and the productive tertiary industry is relatively dispersed. In terms of influencing factors, the interaction between population and employment dominates the spatial differentiation and evolution of the tertiary industry in the PRD urban agglomeration. Therefore, in the future, the tertiary industry in PRD urban agglomeration should promote the optimization of industrial structure and regional coordinated development under the guidance of the government.


Introduction
With the development of economic globalization, the international macro-economy has shifted from an "industrial-oriented economy" to a "service-oriented economy" [1,2], and the agglomeration of the tertiary industry has become an important path to improving the high-quality development of regional economies [3].Driven by factors such as industrialization, urbanization, opening-up, and the continuous growth of rural and urban residents' incomes, China's tertiary industry has continued to increase its share in the national economy and has become an important driving force for national and regional development [4].Under the influence of the 21st-century global economic transformation, the economic development model of China's coastal central cities first began to shift from an industrial-dominant type to a service-dominant type [5].After decades of development, the scale of China's tertiary industry has continued to expand, but there are still deficiencies in its structure, which to some extent affect the quality of economic development and the pace of economic transformation.Entering a new stage of development, China's economic development model starts to shift from a scale expansion-oriented approach to a quality-and efficiency-oriented one.Adjusting and optimizing the structure of the tertiary industry to enhance its level of development is an important strategic choice for China to achieve sustainable economic development [6,7].Consequently, the Chinese government has proposed accelerating the development of modern tertiary industries and upgrading the structure and optimization of the tertiary industry.Guided by national policies, economically developed regions in China should respond first based on their development advantages, with the Pearl River Delta (PRD) urban agglomeration being one of the most vibrant areas of economic development in China and a frontier position for achieving high-quality development.From the current perspective of the development trend of the PRD urban agglomeration's economic structure, its tertiary industry is still in urgent need of accelerated development [8,9].
The study of the spatial layout of the tertiary industry and the spatial differentiation of different types of tertiary industries can explore the development of the tertiary industry from a spatial perspective, which is of great significance for the structural adjustment, optimization, and upgrading of the tertiary industry.As a complex industrial system, the tertiary industry has multiple types, and the spatial layout of different types of tertiary industries varies, but there is currently no unified academic standard for classifying them.This study refers to the classification methods of relevant studies [10][11][12] and divides the tertiary industry into three categories: the productive tertiary industry, the consumptive tertiary industry, and the public tertiary industry.The consumptive tertiary industry refers to services purchased by consumers in the private market, while public tertiary industries are mainly provided by the government for consumer use, and productive tertiary industries refer to services purchased by producers in the private market for further production of goods and services.
Currently, many studies have been conducted on the layout and influencing factors of various types of tertiary industries.The productive tertiary industry is gradually differentiated from the development of the manufacturing industry, and its agglomeration is conducive to forming economies of scale [13], deepening labor division, and thus promoting urbanization [14].Since its development will be significantly affected by the manufacturing industry [15], productive tertiary industry tends to agglomerate around the manufacturing industry in terms of spatial distribution [16].In addition, well-developed transportation infrastructure, a good market competition environment, complete information and communication facilities, high-quality talent resources, and convenient urban facilities also have a significant impact on the agglomeration of the productive tertiary industry [17].On the other hand, the consumptive tertiary industry involves a wide range of fields and a high degree of labor intensity, with diverse forms of specialized development and a huge employment capacity [18,19].As an industry that provides final services to households or individuals, the layout of the consumptive tertiary industry is strongly influenced by population and commercial development levels, and its development can effectively promote the growth of urban centers [20][21][22][23].With the development of social and economic levels, there is a trend for the final demand for services to shift from the private to the public sector [24].Public tertiary industries, as quasi-public products, are greatly influenced by policies in terms of spatial distribution, which is reflected in the fact that public tertiary industries are usually constructed according to population density and have different layout principles during different stages of urbanization [25].The spatial agglomeration of public tertiary industries is not only influenced by population distribution but also by factors such as land prices, transportation accessibility, service targets, and service radius.The spatial distribution of the three types of tertiary industry mentioned above varies due to different service objects and service providers, which are influenced by different factors.The spatial layout and agglomeration of the productive tertiary industry have significant manufacturing orientation and are more sensitive to regional transportation conditions.The consumptive tertiary industry is more affected by market demand and industrial supply; that is, it is more sensitive to population distribution and purchasing power.The public tertiary industry, due to its universal and shared nature, is mostly planned and laid out based on population density and has certain policy orientations.Analyzing the spatial pattern and influencing factors of the tertiary industry is helpful in exploring the impact of different types of tertiary industry on regional economic development, thus providing ideas for the adjustment, optimization, and upgrading of the tertiary industry to a certain extent.
In previous related studies, social and economic data were mostly used for quantitative or qualitative analysis, which has a certain rationality but also leads to problems such as lagging data, poorer consistency, low efficiency, difficulties in the spatialization of these data, and so on.In recent years, geographic information technology has developed rapidly, and more and more types of data can be utilized with the exploration of geospatial and temporal features.Regional industrial space is closely related to population and economic factors, and its spatial characteristics can be explored using a variety of data.However, different data also have their own advantages and disadvantages, such as nighttime lighting data, which can monitor human economic and industrial activities through lighting information, but often has the limitations of the "oversaturation" phenomenon and "spillover effect", which reduces the precision of the research results.Enterprise databases refer to enterprise platform databases that allow third parties to access their data, and can easily access a wide range of enterprise data through enterprise information websites, which is cost-effective and promotes data sharing, but also has limitations in terms of data security, data consistency, and legal and ethical aspects.With the rapid emergence and popularization of web map services, Point of Interest (POI), a data type with geographic coordinate information, has gradually become an important data type for urban research.For the analysis of industrial space, POI data can comprehensively reflect the interaction of population, land, economy, society, and other major elements of the city, integrating geographic location information and functional classification information, and reflecting the spatial distribution characteristics of the object of study.POI data are characterized by easy accessibility, strong current status, rich data volume, high positioning accuracy, better reflection of micro-detailed information, etc. [26].POI data also have the advantages of large scale, wide coverage, multiple categories, easy accessibility, fast updating speed, etc., which have been more and more widely used in the study of the urban spatial field compared with the traditional data [27,28].The spatial distribution of different types of tertiary industry may be driven by different factors such as economic development level and natural environment.Previous studies have mostly analyzed the influencing factors based on socio-economic statistics or questionnaire data [29,30].This study utilizes multi-source big data to conduct a more objective and quantitative analysis of the influencing factors of the spatial differentiation of various types of tertiary industry.For example, this study uses geographic detectors to analyze the main factors affecting the distribution of urban public health risks, breaking the limitations of administrative boundaries in analyzing influencing factors [31,32].
Generally speaking, although studies on the spatial distribution and influencing factors of tertiary industry in urban agglomerations are relatively abundant, most of them focus on the spatial layout and influencing factors of a single type of tertiary industry, ignoring the analysis of the various types included in the tertiary industry, missing a more detailed discussion on the kilometer scale, and failing to include an analysis of the strategic discussion regarding tertiary industry structural adjustment, optimization, and upgrading.Therefore, this study divides the tertiary industry into three types: the productive tertiary industry, the consumptive tertiary industry, and the public tertiary industry.Based on social and economic data, the Dagum Gini coefficient was used to analyze the intra-regional and inter-regional differences in the development of tertiary industry in the Pearl River Delta.Based on POI data, kernel density estimation and local spatial autocorrelation were used to analyze the spatial distribution and agglomeration characteristics of the whole and various types of tertiary industries in the Pearl River Delta.We further explore the influencing factors of their spatial distribution through geographic detectors, using various data sources such as POI, LandScan, the annual China Land Cover Dataset (CLCD), Open Street Map (OSM), the digital elevation model (DEM), and socio-economic data.Based on this, this study proposes strategies for the adjustment and optimization of the Land 2024, 13, 172 4 of 23 tertiary industry structure.This study extends the study content of the tertiary industry in urban agglomerations and enriches the practical significance of spatial study of the tertiary industry.The study results can provide a basis for judging the rationality of the layout of tertiary industry structure, promote the adjustment and optimization of tertiary industry structure in the PRD urban agglomeration, and contribute to the high-quality development of regional economy, which is of great theoretical and practical significance for exploring the formation of spatial characteristics of tertiary industry at the scale of Chinese urban agglomerations and the future development directions of tertiary industry.Moreover, it also has important theoretical and practical significance for regional coordinated development and realizing common prosperity.

Research Area
Located in the south-central part of Guangdong Province and at the estuary of the Pearl River, the PRD urban agglomeration includes nine cities, namely, Guangzhou, Shenzhen, Foshan, Dongguan, Huizhou, Zhongshan, Zhuhai, Jiangmen, and Zhaoqing (Figure 1), covering an area of 55,400 km 2 and a population of 39,249,500 by 2022.The PRD urban agglomeration is one of the most urbanized areas in China and one of the most dynamic and influential urban agglomerations in the country.In 2016, with an urbanization rate of 84.9%, the PRD urban agglomeration surpassed Tokyo, Japan, to become the world's largest urban agglomeration in terms of population and area.By 2022, its urbanization rate had increased to 87.5%.At the same time, the industrial structure of the PRD urban agglomeration is also gradually transforming and upgrading.In 2009, the tertiary industry accounted for 49.9% of GDP, exceeding the secondary industry for the first time.As of 2021, its share has reached 57.57% [33].Although the proportion of the tertiary industry in the PRD urban agglomeration has been rising in recent years, its structure still needs to be further optimized and adjusted [34,35].Therefore, clarifying the spatial distribution characteristics of the tertiary industry in PRD urban agglomeration and the distribution differences of various types of tertiary industry and exploring their influencing factors are of great practical significance for promoting the sustainable development of the urban agglomeration.
Land 2024, 13, x FOR PEER REVIEW 4 of 23 data.Based on this, this study proposes strategies for the adjustment and optimization of the tertiary industry structure.This study extends the study content of the tertiary industry in urban agglomerations and enriches the practical significance of spatial study of the tertiary industry.The study results can provide a basis for judging the rationality of the layout of tertiary industry structure, promote the adjustment and optimization of tertiary industry structure in the PRD urban agglomeration, and contribute to the high-quality development of regional economy, which is of great theoretical and practical significance for exploring the formation of spatial characteristics of tertiary industry at the scale of Chinese urban agglomerations and the future development directions of tertiary industry.Moreover, it also has important theoretical and practical significance for regional coordinated development and realizing common prosperity.

Research Area
Located in the south-central part of Guangdong Province and at the estuary of the Pearl River, the PRD urban agglomeration includes nine cities, namely, Guangzhou, Shenzhen, Foshan, Dongguan, Huizhou, Zhongshan, Zhuhai, Jiangmen, and Zhaoqing (Figure 1), covering an area of 55,400 km 2 and a population of 39,249,500 by 2022.The PRD urban agglomeration is one of the most urbanized areas in China and one of the most dynamic and influential urban agglomerations in the country.In 2016, with an urbanization rate of 84.9%, the PRD urban agglomeration surpassed Tokyo, Japan, to become the world's largest urban agglomeration in terms of population and area.By 2022, its urbanization rate had increased to 87.5%.At the same time, the industrial structure of the PRD urban agglomeration is also gradually transforming and upgrading.In 2009, the tertiary industry accounted for 49.9% of GDP, exceeding the secondary industry for the first time.As of 2021, its share has reached 57.57% [33].Although the proportion of the tertiary industry in the PRD urban agglomeration has been rising in recent years, its structure still needs to be further optimized and adjusted [34,35].Therefore, clarifying the spatial distribution characteristics of the tertiary industry in PRD urban agglomeration and the distribution differences of various types of tertiary industry and exploring their influencing factors are of great practical significance for promoting the sustainable development of the urban agglomeration.

Data Source 2.2.1. POI
This study uses POI data, which is widely used in spatial studies of urban agglomerations.These data are derived from the 2022 national POI dataset and the industry-related POI data of each prefecture-level city within the scope of the PRD urban agglomeration are extracted.Subsequently, the categories related to tertiary industry are selected according to the classification of tertiary industry in existing studies, and some subcategories with ambiguous industry characteristics are examined as well as re-categorized.A total of 2,329,680 POI data points are obtained after spatial matching and de-duplication of the acquired data, and the POI data are then classified into types according to the industry sectors included in each type of tertiary industry [36] (Table 1).This study uses POI data to investigate the spatial differentiation of tertiary industry in the PRD urban agglomeration, and uses the increment of POI data from 2012 to 2022 to explore the influence of policies on the layout of tertiary industry [37].[38][39][40], consisting of nine land cover types, including impervious, cropland, forest, water, shrub, etc., which make it suitable for extracting urban built-up areas.In this study, the CLCD data are used to extract the built-up area of the PRD urban agglomeration and calculate its size.

LandScan
The population data used in this study are sourced from the LandScan dataset from Oak Ridge National Laboratory, U.S. The LandScan population database utilizes the best available census data, and establishes a weight model based on geographic information systems and partition density models, with a spatial resolution of 1 km, making it a highresolution global population distribution dataset that can better reflect population spatial distribution [41].In this study, the LandScan data for PRD urban agglomeration in 2021 are used to illustrate urban population density.

OSM
The road network data used in this study are obtained from the Open Steet Map website (https://www.openstreetmap.orgaccessed on 1 February 2023).Open Street Map (OSM) is an online multi-source map database that provides geospatial data including highways, railroads, water systems, waters, land use, buildings, etc. OSM aims to provide users with free and easily accessible digital map resources, making it the most popular Volunteered Geographic Information (VIG) source currently available [42].OSM roads are of high quality in terms of fixing accuracy and topological relationships, containing basic spatial information such as latitude and longitude as well as attribute information such as road name, road type, maximum travel speed, and one-way streets [43,44].In this study, OSM data of the PRD urban agglomeration in 2021 are used to reflect the traffic level of the urban agglomeration by calculating the road network density.

DEM
The Digital Elevation Model (DEM) is an important foundation for terrain analysis, such as slope, aspect, and hydrological analysis [45], among others.This study uses DEM data from the geospatial Data Cloud website (https://www.gscloud.cn/searchaccessed on 1 February 2023) to explore the terrain factors that contribute to spatial differentiation in the tertiary industry in the PRD urban agglomeration.

Socio-Economic Data (1) Revenue Data of Tertiary Industry Enterprises
To evaluate the current development status of the tertiary industry in PRD urban agglomeration, this study selects data from tertiary industry enterprises above designated size in Guangdong Province from 2016 to 2022, as reported in the "Guangdong Provincial Statistical Yearbook".The proportion of enterprise revenue (the proportion of revenue of tertiary industry enterprises above designated size = the total revenue of tertiary industry enterprises above designated size in different cities/the total revenue of tertiary industry enterprises above designated size in PRD urban agglomeration) is used to measure the level of tertiary industry development.The revenue proportion data ranges from 0 to 1, and the PRD urban agglomeration is further divided into the Shenzhen-Dongguan-Huizhou economic circle, the Guangzhou-Foshan-Zhaoqing economic circle, and the Zhuhai-Zhongshan-Jiangmen economic circle, which are, respectively, denoted by S-D-H Economic Circle, G-F-Z Economic Circle, and Z-Z-J Economic Circle in the calculation.
(2) Employment Data These data are obtained from the seventh national census data released by the National Bureau of Statistics in early 2023, which is connected to the kilometer grid through ArcGIS and combined with LandScan data for optimization and processing, and finally spatialized employment raster data are obtained.In this study, the employment data of PRD urban agglomeration are used to analyze the human resource factors that contribute to the differentiation of the tertiary industry in the urban agglomeration.
(3) GDP spatial distribution grid data These data are the GDP spatial distribution kilometer grid data of 2019 in China, sourced from the Resource and Environment Science and Data Center of the Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences [46].These data are 1 km × 1 km spatial grid data generated by spatial interpolation on the basis of national sub-county GDP statistics, considering the spatial interaction pattern between human activity-related land use types, night-time light brightness, residential density data, and GDP.In this study, the GDP spatial distribution kilometer grid data of PRD urban agglomeration is used to analyze the economic factors contributing to the differentiation of tertiary industry within the urban agglomeration.In summary, all types of study data contain relatively comprehensive information, which can comprehensively discuss the factors affecting the spatial distribution of the tertiary industry in PRD urban agglomeration.All kinds of data will be established through ArcGIS 1 km × 1 km fishing nets to unify their resolution.The basic information of various research data is shown in Figure 2 and Table 2.

REVIEW 7 of 23
In summary, all types of study data contain relatively comprehensive information, which can comprehensively discuss the factors affecting the spatial distribution of the tertiary industry in PRD urban agglomeration.All kinds of data will be established through ArcGIS 1 km × 1 km fishing nets to unify their resolution.The basic information of various research data is shown in Figure 2 and Table 2.

Research Methods
On the basis of the classification of tertiary industry, this study uses the Dagum Gini coefficient, kernel density estimation, and Anselin Local Moran's I to explore spatial differentiation of tertiary industry.Based on this, this paper quantitatively analyzes the factors influencing the spatial differentiation of the tertiary industry in PRD urban agglomeration, so as to better understand the problems existing in the tertiary industry structure and put forward some strategic suggestions.The specific research framework is shown in Figure 3.

Research Methods
On the basis of the classification of tertiary industry, this study uses the Dagum Gini coefficient, kernel density estimation, and Anselin Local Moran's I to explore spatial differentiation of tertiary industry.Based on this, this paper quantitatively analyzes the factors influencing the spatial differentiation of the tertiary industry in PRD urban agglomeration, so as to better understand the problems existing in the tertiary industry structure and put forward some strategic suggestions.The specific research framework is shown in Figure 3.

Dagum Gini Coefficient and Its Decomposition
The classical Gini coefficient is constructed based on the assumptions of normal distribution and homoscedasticity.Although it imposes strict constraints on the absence of overlapping parts between grouped samples, it is difficult to decompose several sub-indices with reasonable economic meanings.To address these limitations, Dagum proposed Dagum Gini coefficient, which can decompose the overall Gini coefficient G of the sample into three parts, namely, intra-group differences  , inter-group net differences  and inter-group hypervariance density  [47,48].The intra-group  reflects the level gap within each region, the inter-group  reflects the level gap between regions, and the hypervariable density  reflects the cross-overlapping phenomenon of each region, reflecting the relative gap.
The Dagum Gini coefficient was originally used to measure regional income gaps.Its decomposition method by subgroups can effectively solve the source of regional differences.Therefore, it is widely used in many fields to describe the problem of uneven regional development.The development level of the tertiary industry in the Pearl River Delta urban agglomeration is uneven in regional space.This paper uses the Dagum Gini coefficient and its decomposition measure to measure the spatial differences in the development level of the tertiary industry in the Pearl River Delta urban agglomeration and explore the sources of the differences.The specific calculation process is as follows: where  is the number of regional divisions,  and ℎ are different ranges within the  regions,  is the sum of the number of cities in each region,  ( ) is the development level of the tertiary industry in the city () within the (ℎ),  ( ) is the number of cities

Dagum Gini Coefficient and Its Decomposition
The classical Gini coefficient is constructed based on the assumptions of normal distribution and homoscedasticity.Although it imposes strict constraints on the absence of overlapping parts between grouped samples, it is difficult to decompose several subindices with reasonable economic meanings.To address these limitations, Dagum proposed Dagum Gini coefficient, which can decompose the overall Gini coefficient G of the sample into three parts, namely, intra-group differences G w , inter-group net differences G b and inter-group hypervariance density G t [47,48].The intra-group G w reflects the level gap within each region, the inter-group G b reflects the level gap between regions, and the hypervariable density G t reflects the cross-overlapping phenomenon of each region, reflecting the relative gap.
The Dagum Gini coefficient was originally used to measure regional income gaps.Its decomposition method by subgroups can effectively solve the source of regional differences.Therefore, it is widely used in many fields to describe the problem of uneven regional development.The development level of the tertiary industry in the Pearl River Delta urban agglomeration is uneven in regional space.This paper uses the Dagum Gini coefficient and its decomposition measure to measure the spatial differences in the development level of the tertiary industry in the Pearl River Delta urban agglomeration and explore the sources of the differences.The specific calculation process is as follows: where k is the number of regional divisions, j and h are different ranges within the k regions, n is the sum of the number of cities in each region, y ij (y hr ) is the development level of the tertiary industry in the city i(r) within the j(h), n j (n h ) is the number of cities in the region Land 2024, 13, 172 9 of 23 j(h), and µ j (µ h ) is the average development level of the tertiary industry in the region j(h).
The Gini coefficient G jj for region j and the intra-group difference G w can be expressed as Equations ( 2) and ( 4), respectively.The Gini coefficient G jh and the inter-group net value difference G nb between regions j and h can be represented by Equations ( 3) and ( 5), and the hypervariance density G m can be expressed as Equation (6).
In Equation ( 5), p j = n j /n, s j = n j µ j /nµ, and the same applies to p h , s h .In Equation ( 6), D jh measures the interactive influence of tertiary industry development levels between regions j and h, and it is calculated specifically in Equation (7).In Equation ( 8), q jh refers to the difference in tertiary industry development levels between the three major regions, representing the mathematical expectation of y ji − y hr > 0 between regions i and h, and E jh is the hypervariable first-order moment, indicating the mathematical expectation of y ji − y hr < 0 between regions j and h.

Kernel Density Estimation
Kernel density estimation is a method of studying the probability of point occurrence in different locations in space, based on point data, by calculating the density of elements around them, which dates back to Rosenblatt [49] and Parzen [50].Different kernel density values are used to represent the spatial distribution characteristics of features, with the maximum density at the center point gradually decreasing with increasing distance.Kernel density estimation is widely used in many fields, such as geographic information and urban research.Its classic applications include various prediction and clustering tasks.This study explores the spatial differentiation characteristics of the tertiary industry in the Pearl River Delta (PRD) based on POI point data.Kernel density estimation is able to spatially cluster the point data better, so this study utilizes the selection of kernel density estimation to identify the spatial structural characteristics of the tertiary industry in the PRD urban agglomerations and conducts a comparative analysis of the spatial differentiation of each type of tertiary industry.The calculation formula is as follows: Land 2024, 13, 172 10 of 23 where f (s) is the kernel density estimate value at feature s, n is the number of features whose distance to s is less than or equal to h, h is the search radius, and k is the spatial distance weight.
In addition, different search radii can lead to different kernel density analysis results.Therefore, the result of the kernel density analysis can be determined by a validated search radius, and the formula for the search radius is as follows: where SD is the standard distance, D m is the median distance, and n is the number of event points.

Anselin Local Moran's I
Tobler's First Law of Geography states that geographic phenomena or attributes are spatially interrelated and exhibit clustering, randomness, or regularity in their distribution [51].In order to explore the spatial patterns or distribution characteristics of natural or social phenomena, spatial autocorrelation indices are often used to measure the degree of correlation among natural or social attributes in space, since the size of the correlation can indicate the spatial pattern and distribution characteristics of the attribute [52].Anselin Local Moran's I can be used to test whether there is a clustering of variables in a local area, as well as the degree and significance of spatial differences between adjacent units [53].Local indicators of spatial association (LISA) are commonly used to analyze spatial correlation and decompose Global Moran's I into each spatial unit [54].In urban research, spatial autocorrelation analysis is used to determine whether a variable is correlated in urban space, and local spatial autocorrelation analysis can further reflect the correlation between each parcel and its neighboring parcels in terms of a certain attribute.Since the spatial layout of the service sector in the PRD city cluster has significant differences between the central city and other regions, this work uses a 1 km 2 grid as a geographic unit, and Anselin Local Moran's I is used to reflect the degree of difference between individual geographic units of the service sector in the PRD city cluster and other geographic units.
The formula is as follows: Anselin Local Moran's I can be used in conjunction with statistical Z-tests to identify spatial correlation patterns with a 95% confidence level of HH (high-high agglomeration), LL (low-low agglomeration), LH (high-value outliers surrounded by low-values), and HL (low-value outliers surrounded by high-values).

Geodetector
Geodetector is a statistical method that investigates the differentiation patterns of various geographic phenomena in space and attempts to reveal their driving factors [55].Compared with traditional statistical analysis methods, geodetectors have obvious advantages in exploring the influence mechanism of spatial heterogeneity of elements without too many assumptions, and therefore are widely used in urban studies [32,56].
Currently, the main research method for studying influencing factors is to use various quantitative models for quantitative research.However, the analysis of influencing factors on the regional layout of tertiary industry is generally based on a certain level of administrative districts.For example, the study of the spatial distribution characteristics and influencing factors of the tertiary industry in the Wuhan market [25] is of great significance for the development of tertiary industry and the economic development of surrounding areas.However, the influencing factors of tertiary industry layout occur based on geospatial proximity rather than being limited to administrative boundaries, and the influence of administrative boundaries is diminishing under the continuous advancement of the policy of regional integration and development of the PRD city cluster.Therefore, this study divides the study area into kilometer grid scales and explores the influencing factors of the spatial distribution characteristics of the tertiary industry in the PRD city cluster and the magnitude of their explanatory power by using factor probes and interaction probes in a geodetector.
(1) Factor Detection The factor detection can detect the spatial heterogeneity of the dependent variable and the explanatory power of the independent variable on the dependent variable.This value is measured by the q-value.
The formula is: where q has a value range from 0 to 1, with a larger value indicating stronger spatial heterogeneity.If the stratification is generated by the independent variable X, a larger q-value indicates that the independent variable has a stronger explanatory power for the dependent variable, whereas a smaller q-value indicates a weaker explanatory power. (

2) Interaction Detection
Interaction detection is mainly used to identify interactions between different factors, i.e., to evaluate whether two factors increase or decrease the explanatory power of tertiary industry, or whether these factors have independent impacts on tertiary industry.
(3) Factor system construction The spatial distribution of the tertiary industry is influenced by various factors such as the economy, society, and infrastructure, and the degree of influence of different factors also varies.Existing studies mostly believe that factors such as built-up area, population density, GDP, urbanization level, and policies all affect the spatial distribution of the tertiary industry.Considering the actual situation in PRD urban agglomeration and the availability of data, nine indicators are selected as the factors influencing the spatial distribution of the tertiary industry according to the principles of scientificity, representativeness, and accessibility.The dependent variable Y is the distribution density of various tertiary enterprises in each city by January 2022, with the independent variable being the 7 indicators in Table 3. Notes: The abbreviations of topography, GDP, built-up area ratio, population density, employment-population ratio, and government are written as X 1 , X 2 , X 3 , X 4 , X 5 , X 6 , and X 7 .
"Topography" has a basic influence on the centralized and continuous layout of industries; "GDP" represents the economic development of the city; "Built-up area ratio" can reflect the development of urbanization; and "Population density" can reflect the tertiary industry in a certain geographic space."Population density" reflects the size of the consumer market in a certain geographic space; "Road density" reflects the city's transportation situation, and to some extent, the distance to points of interest and facilities; "Employment population ratio" reflects the city's employment situation and to some extent the distance to facilities; "Employment population ratio" can reflect the employment situation of the city, and to a certain extent, the ability of the city to retain young people and qualified population; and "Government" can represent a series of top-level designs of the government, such as promoting the high-quality development of the tertiary industry through land-use planning, industrial development planning, and various preferential policies.The abbreviations of topography, GDP, built-up area ratio, population density, employment-population ratio, and government are written as X 1 , X 2 , X 3 , X 4 , X 5 , X 6, and X 7 in Table 3.

Differences and Decomposition of Tertiary Industry Development Levels in the PRD Urban Agglomeration
This study adopts the Dagum Gini coefficient and decomposition method to measure and decompose the relative level of development of the tertiary industry in the three major economic regions of the PRD (the Shenzhen-Dongguan-Huizhou Economic Circle, Guangzhou-Foshan-Zhaoqing Economic Circle, and Zhuhai-Zhongshan-Jiangmen Economic Circle), as well as their evolutionary trends.This part mainly covers two aspects.Firstly, from an overall perspective, the overall relative level and evolving trend of tertiary industry development are measured to demonstrate the development gaps and changes in the tertiary industry among the Shenzhen-Dongguan-Huizhou economic circle, the Guangzhou-Foshan-Zhaoqing economic circle, and the Zhuhai-Zhongshan-Jiangmen economic circle.Secondly, from a regional perspective, the relative levels and evolving trends of industrial development in various cities within each region are measured to showcase the development gaps and changes in the tertiary industry among different cities within each region.
(1) Overall Relative Level of Tertiary Industry Development Since the Statistical Yearbook of Guangdong Province did not collect statistics on the revenue of tertiary industries above designated size before 2015, relevant data from 2015 to 2022 are selected for analysis.The overall relative level differences and evolving trends of tertiary industry development in PRD urban agglomeration from 2015 to 2022 are shown in Table 4.The overall Gini coefficient of the Pearl River Delta urban agglomeration has continued to increase in recent years, exceeding the 2015 level.The intra-group Gini coefficient G w has increased in fluctuation and has remained relatively stable in recent years; that is, the differences in the development levels of the service industry among the three major economic circles in the Pearl River Delta urban agglomeration have not changed much.The Gini coefficient G b between groups is generally high and the fluctuations increase significantly, indicating that the development levels of the service industry among the economic circles of the Pearl River Delta urban agglomeration are quite different.The super-variable density Gini coefficient G t shows a fluctuating downward trend, indicating that the degree of inequality between and within the economic circles of the Pearl River Delta urban agglomeration is decreasing.(2) Intra-regional and Inter-regional Differences in Tertiary Industry Development The regional differences in the development level of the tertiary industry mainly include intra-regional differences in the level of development, inter-regional differences in the level of development, and the contribution rate of the overall differences.The contribution rate of overall differences includes intra-regional differences, inter-regional differences, and hypervariable density.Among them, the hypervariable density reflects the influence of cross-overlapping between regions on the overall difference.Table 5 presents the decomposition results of the Dagum Gini coefficient differences, and Figure 4 depicts the evolving trends of regional differences in tertiary industry development and their decomposition.Figure 4a illustrates the relative differences in industrial development within each region, Figure 4b shows the relative differences in industrial development between the three major regions, and Figure 4c displays the contribution sources to the overall difference in tertiary industry development.

Intra-Group Gini Coefficient
Inter-Group Gini Coefficient (2) Intra-regional and Inter-regional Differences in Tertiary Industry Development

S-D-H G-F-Z Z-Z-J S-D-H & G-F-Z S-D-H & Z-Z-J G-F-Z & Z-Z-
The regional differences in the development level of the tertiary industry mainly clude intra-regional differences in the level of development, inter-regional differences the level of development, and the contribution rate of the overall differences.The con bution rate of overall differences includes intra-regional differences, inter-regional diff ences, and hypervariable density.Among them, the hypervariable density reflects the fluence of cross-overlapping between regions on the overall difference.Table 5 prese the decomposition results of the Dagum Gini coefficient differences, and Figure 4 depi the evolving trends of regional differences in tertiary industry development and their d composition.Figure 4a illustrates the relative differences in industrial development with each region, Figure 4b shows the relative differences in industrial development betwe the three major regions, and Figure 4c displays the contribution sources to the overall d ference in tertiary industry development.

Intra-Group Gini Coefficient
Inter-Group Gini Coefficient  From the perspective of intra-regional differences, the tertiary industry exhibits significant differences within the Guangzhou-Foshan-Zhaoqing economic circle and the Shenzhen-Dongguan-Huizhou economic circle, which have remained relatively stable.In contrast, the intra-regional differences within the Zhuhai-Zhongshan-Jiangmen economic circle are relatively small and show a fluctuating upward trend.This indicates that the Guangzhou-Foshan-Zhaoqing economic circle centered around Guangzhou and the Shenzhen-Dongguan-Huizhou economic circle centered around Shenzhen have consistently exhibited significant development disparities.As growth poles in regional development, their ability to radiate and drive surrounding cities is inadequate.On the other hand, the Zhuhai-Zhongshan-Jiangmen economic circle lacks a strong growth pole.

S-D-H G-F-Z Z-Z-J S-D-H & G-F-Z S-D-H & Z-Z-J G-F-Z & Z-
From the perspective of inter-regional differences, there are significant differences between the Guangzhou-Foshan-Zhaoqing economic circle and the Zhuhai-Zhongshan-Jiangmen economic circle, as well as between the Shenzhen-Dongguan-Huizhou economic circle and the Zhuhai-Zhongshan-Jiangmen economic circle.However, the inter-regional differences between the Guangzhou-Foshan-Zhaoqing economic circle and the Shenzhen-Dongguan-Huizhou economic circle are relatively small, and the development differences among all regions are relatively stable.This indicates that the tertiary industry development level in the Zhuhai-Zhongshan-Jiangmen economic circle significantly lags behind that of the Guangzhou-Foshan-Zhaoqing economic circle and the Shenzhen-Dongguan-Huizhou economic circle.
From the decomposition of the overall difference, the inter-regional difference is the main source of the overall difference 3.1.2.Spatial Differentiation Based on Kernel Density The regional differences in development level and their sources in the three major regions are analyzed in detail using the Dagum Gini coefficient, but it only reflects relative differences.Therefore, Kernel density estimation is employed to further analyze the absolute differences and dynamic evolution of the development levels of the three types of tertiary industries in various regions of PRD urban agglomeration.
The kernel density of POI data of the service industry in the Pearl River Delta was analyzed and classified by the Jenks method.The kernel density value was divided into five grades: low, lower, medium, higher, and high.The higher the density value, the higher the agglomeration degree of the service industry, and vice versa.The analysis results are shown in Figure 5. Overall, the high-density areas and agglomeration areas are mainly located in the administrative and economic centers of the central urban areas or new areas where infrastructure construction is relatively complete, the population is more concentrated, and they are greatly influenced by the radiation of the central urban areas.Further exploration of the evolving characteristics of the spatial structure of various types of tertiary industries reveals that the public tertiary industry consistently maintains a monocentric spatial structure with Guangzhou as the core.The coverage range of each density level expands continuously, but there is relatively little change in the development levels of the main cities.In addition, the productive tertiary industry always maintains a bicentric spatial structure with Guangzhou and Shenzhen as the cores.The Pearl River estuary's eastern corridor is gradually improving, while the western bank lacks sufficient development levels.The trend of urban integration between Guangzhou and Foshan is evident.Additionally, the consumptive tertiary industry gradually transitions from a monocentric structure with one main center to a polycentric structure.The Shenzhen-Dongguan-Huizhou economic circle develops rapidly and forms multiple secondary centers.The urban integration between Guangzhou and Foshan continues to advance.

Spatial Differentiation Based on Spatial Autocorrelation
Local spatial autocorrelation analysis is performed on various tertiary industries in the PRD urban agglomeration (Figure 6), and the following characteristics are obtained: (1) Spatial agglomeration distribution and evolution of tertiary industry are generally consistent.H-H areas are mainly distributed in core cities of PRD urban agglomeration, such as Guangzhou, Shenzhen, Dongguan, and Foshan; L-L areas mainly surround the core cities, mainly in Zhaoqing, Jiangmen, and Huizhou; and L-H areas are more dispersed and mainly occur in peripheral areas of cities where H-H areas are distributed, with a small range.In addition, there is no distribution of H-L areas in the agglomeration of overall and various types of tertiary industry, indicating that the tertiary industry structure of PRD urban agglomeration is stable as a whole.(2) Positive correlation is the main spatial correlation pattern of tertiary industry, and the H-L cluster is more and more extensive.There is a certain similarity in the overall distribution of tertiary industry and the distribution of H-H areas and L-L areas for various types.Among these, the distribution of H-H areas and L-L areas is wider and more spatially heterogeneous for the consumptive tertiary industry, indicating that developed cities in the consumptive tertiary industry have a stronger driving effect on surrounding cities.

Spatial Differentiation Based on Spatial Autocorrelation
Local spatial autocorrelation analysis is performed on various tertiary industries in the PRD urban agglomeration (Figure 6), and the following characteristics are obtained: (1) Spatial agglomeration distribution and evolution of tertiary industry are generally consistent.H-H areas are mainly distributed in core cities of PRD urban agglomeration, such as Guangzhou, Shenzhen, Dongguan, and Foshan; L-L areas mainly surround the core cities, mainly in Zhaoqing, Jiangmen, and Huizhou; and L-H areas are more dispersed and mainly occur in peripheral areas of cities where H-H areas are distributed, with a small range.In addition, there is no distribution of H-L areas in the agglomeration of overall and various types of tertiary industry, indicating that the tertiary industry structure of PRD urban agglomeration is stable as a whole.(2) Positive correlation is the main spatial correlation pattern of tertiary industry, and the H-L cluster is more and more extensive.There is a certain similarity in the overall distribution of tertiary industry and the distribution of H-H areas and L-L areas for various types.Among these, the distribution of H-H areas and L-L areas is wider and more spatially heterogeneous for the consumptive tertiary industry, indicating that developed cities in the consumptive tertiary industry have a stronger driving effect on surrounding cities. (3) There is a local micro-difference in the spatial agglomeration distribution of different types.The agglomeration degree of productive tertiary industry in Huaiji County and Guangning County in the north of Zhaoqing City is significantly lower than that of consumptive and public tertiary industry, which is related to the city's industrial foundation and location.Another region with differences is the southern part of Foshan City, where the H-H areas of all types of tertiary industry are scattered in a point-like manner, but the H-H areas of public tertiary industry have a larger agglomeration range compared to the other two types, which is related to the principle of equal distribution of public tertiary industry under the current high level of urbanization.
Land 2024, 13, x FOR PEER REVIEW 16 of 23 (3) There is a local micro-difference in the spatial agglomeration distribution of different types.The agglomeration degree of productive tertiary industry in Huaiji County and Guangning County in the north of Zhaoqing City is significantly lower than that of consumptive and public tertiary industry, which is related to the city's industrial foundation and location.Another region with differences is the southern part of Foshan City, where the H-H areas of all types of tertiary industry are scattered in a point-like manner, but the H-H areas of public tertiary industry have a larger agglomeration range compared to the other two types, which is related to the principle of equal distribution of public tertiary industry under the current high level of urbanization. (

Factor Detection Results
The purpose of factor detection is to calculate the q values of the influencing factors.The q value ranges from 0 to 1, and the closer it is to 1, the stronger the explanatory power of the independent variable X is for the attribute, and vice versa.The results are shown in Table 6.Overall, the q value of each type of tertiary industry is consistent with that of the whole.X2, X4, and X5 generally have high explanatory power in determining the spatial  The purpose of factor detection is to calculate the q values of the influencing factors.The q value ranges from 0 to 1, and the closer it is to 1, the stronger the explanatory power of the independent variable X is for the attribute, and vice versa.The results are shown in Table 6.Overall, the q value of each type of tertiary industry is consistent with that of the whole.X 2 , X 4 , and X 5 generally have high explanatory power in determining the spatial distribution of the tertiary industry, indicating that the development of the tertiary industry requires a favorable economic environment, a market for services and consumption, and convenient transportation conditions.From the perspective of different types of tertiary industries, X 4 and X 5 have a significant explanatory power for the public tertiary industry, indicating that public tertiary industries, which often have a public policy attribute, are usually built, and equipped with public services and infrastructure based on population indicators.X 3 and X 6 have significant explanatory power for the productive tertiary industry, indicating that the productive tertiary industry tends to have characteristics such as large land occupation and are labor-intensive.X 1 and X 7 have significant explanatory power for the consumptive tertiary industry, indicating that the placement of the consumptive tertiary industry often occurs within urban areas and is greatly influenced by the government's financial expenditure.

Interaction Detection Results
The results of the interaction detection for the overall tertiary industry and each type of tertiary industry are shown in Table 7.The study finds that: (1) The results of the interaction detection between the overall tertiary industry and individuals show a non-linear enhancement or dual-factor enhancement, indicating that the influence of interactions on the spatial distribution of tertiary industry is greater than that of single factors.(2) X 4 and X 6 rank first in the interaction factors for all types of tertiary industry, further confirming that the tertiary industry depends on the consumer market and promotes employment.(3) After interaction, the explanatory power of factors with small q-values is greatly increased, such as X 3 ∩X 5 , X 6 ∩ X 7 .(4) The results of the interaction detection for the tertiary industry show that neither the overall nor the individual factors exhibit non-linear weakening, single-factor nonlinear weakening, or independence, indicating that the spatial distribution of tertiary industry in the PRD urban agglomeration is not caused by a single factor, but is rather the result of the comprehensive effects of different influencing factors.

Discussion
The optimization of the tertiary industry structure is an important prerequisite for achieving high-quality economic development [57].The development process of modern economies has shown that all economies that have entered the ranks of high-income countries have seen the tertiary industry's GDP ratio surpass the sum of the first and second industries while maintaining a relatively stable level [58].Such economies show a coordinated development pattern in the productive tertiary industry, consumptive tertiary industry, and public tertiary industry.To achieve high-quality economic development, it is not enough to simply increase the overall proportion of the tertiary industry.It is also necessary to reasonably allocate and optimize the structure of three types of productive tertiary industry, consumptive tertiary industry, and public tertiary industry according to their respective characteristics.
The regional development differences in the tertiary industry within the PRD urban agglomeration suggest that the higher-level economic circles within the urban agglomeration often experience significant internal development disparities, with limited spillover effects on the surrounding areas.This further amplifies the regional development gaps.On the other hand, although the lower-level economic circles have smaller internal disparities, they are trapped in a long-term low-level development state due to the lack of spillover effects from growth poles.Therefore, in the future, efforts should be made to consolidate and strengthen the roles of Guangzhou and Shenzhen as growth poles, promote the high-end development of the tertiary industry, and maximize their spillover effects and coverage range to facilitate regional integration.The lower-level economic circles should leverage their comparative advantages and accelerate the cultivation of new growth poles.
From the analysis of spatial structure and spatial agglomeration, it can be seen that the structure of tertiary industry as a whole is relatively stable, and there exists a certain consistency in the distribution of agglomeration and decentralization.However, there also exists a certain incoherence in the structure of the industry, such as the higher level of agglomeration of the conservative tertiary industry; the more balanced distribution of the public tertiary industry; and the relatively decentralized layout of the productive tertiary industry.Therefore, in the future, it will be necessary to continue promoting the high-level development of the consumptive tertiary industry, facilitate the development of the productive tertiary industry in the cities on the western bank of the Pearl River estuary, and strengthen government regulation over the public tertiary industry to ensure the equitable allocation of public service facilities.
The analysis of influencing factors reveals that the interaction between population and employment dominates the spatial structure of the tertiary industry in the PRD urban agglomeration.Against the background of China's gradually diminishing demographic dividend, the future development of the tertiary industry in the PRD urban agglomeration should no longer rely on population size as a stimulus for industry growth.Instead, more attention should be put on optimizing industrial structure and promoting regional coordinated development under government guidance.
Based on the above analysis, it is concluded that in the future, the high-quality development of the tertiary industry in the PRD city clusters should optimize the internal structure of the tertiary industry, clarify the development priorities, and strengthen government guidance.Specifically, in terms of the development structure, the construction of a radiation-diffusion service network with the Pearl River Delta (PRD), especially Guangzhou and Shenzhen, as the radiation center and the surrounding cities as sub-centers, extending to the city and county levels as well as to the townships, should be accelerated, and regional services with specialized characteristics should be developed in various places in accordance with different circumstances, so as to form a pattern of development of the tertiary industry that combines centralized radiation with specialized division of labor.In terms of development priorities, the PRD region should also adhere to the strategy of high-end development, accelerate the development of a modern tertiary industry system that is complementary to the advanced manufacturing industry, and actively lead the construction of the Greater Bay Area of Guangdong, Hong Kong, and Macao in a wider context.Specifically, for the two development cores of the city cluster, Guangzhou and Shenzhen should both enhance the development of productive service industries.Guangzhou should develop high-end chains serving the regional economy with its strengths in shipping logistics and headquarters economy, while Shenzhen should give full play to its strengths in securities, banking, insurance, and cross-border finance, and vigorously develop the information economy.In terms of top-level design, the government's planning and policy guidance should also be strengthened.Given that the level of tertiary industry development in the PRD city clusters varies greatly in many dimensions, the development of the regional tertiary industry should be coordinated, and differentiated tertiary industry development policies should be implemented in light of the characteristics of each region.It should also strengthen the comprehensive coordination of regional industrial planning, land use planning, and tertiary industry policies, and focus on the internal structure of the tertiary industry to provide policy and financial guidance for its development, so as to break the development constraints of the core-periphery structure of the city clusters [59].
Based on the classification of the tertiary industry, this study analyzes the spatial differentiation characteristics of the tertiary industry based on socio-economic data and POI data and explores the influencing factors of the tertiary industry using multi-source data, so as to propose strategic recommendations for adjusting and optimizing the structure of the tertiary industry.Current studies on the spatial structure of the tertiary industry in urban agglomerations mainly analyze the tertiary industry as a whole or a single type [60].However, the tertiary industry includes a variety of categories, and its sub-types also have significant differences, so it is difficult to accurately grasp the spatial layout of the tertiary industry and the spatial differentiation of each type by a generalized study on the whole industry.Meanwhile, studying a single type of tertiary industry will not provide an overall comprehensive grasp of the complex system of the tertiary industry.This study categorizes the tertiary industry into the productive tertiary industry, consumptive tertiary industry, and public tertiary industry, and analyzes the spatial distribution of the tertiary industry as a whole and by type, which helps to accurately grasp the spatial distribution and differences of the tertiary industry in urban agglomerations.In terms of research data and methods, current research on the industrial spatial structure of urban agglomerations using POI data is becoming increasingly abundant [61,62].However, in terms of further analysis of the influencing factors of their spatial layout, the existing studies are mostly qualitative analyses or quantitative analyses using socio-economic data to construct models.Conversely, using multi-source big data and spatialized socio-economic data as the data basis enables this study to reflect the actual situation of socio-economic activities more accurately in PRD urban agglomeration.Furthermore, by using the geodetector to analyze the factors influencing the spatial layout of the tertiary industry in urban agglomerations, the study can objectively reflect the influencing factors in different spatial dimensions while filling in the gap in existing research's inadequate consideration of spatial dimensions.In addition, existing research has placed more emphasis on analyzing the spatial differentiation and influencing factors of the tertiary industry, but has neglected to study how to optimize and enhance its development.This study analyzes the structural problems of the tertiary industry in PRD urban agglomeration by exploring its spatial differentiation and influencing factors, and then proposes corresponding strategic recommendations which are of guidance significance for the high-quality development of the tertiary industry in the PRD urban agglomeration.
The spatial study of the tertiary industry in urban agglomerations is not a new topic, and many studies have comprehensively analyzed the spatial layout and influencing factors of the tertiary industry in different urban agglomerations or different types of tertiary industry in China [63][64][65].Based on previous research, this study identifies the spatial differentiation and influencing factors of different types of tertiary industry in the PRD urban agglomeration and analyzes the important role of the internal structure of the tertiary industry in its overall development.In addition, this study spatializes socio-economic data and combines multi-source big data to quantitatively analyze the factors influencing the spatial differentiation of the tertiary industry at the kilometer grid scale, making the research results more objective and reliable.Finally, this study analyzes the characteristics and problems of the tertiary industry's internal structure in PRD urban agglomeration and proposes corresponding strategic recommendations.This study also extends the research content of the tertiary industry in city clusters, expands the practical significance of the spatial study of the tertiary industry, provides a reference for the development of the tertiary industry in the PRD city clusters, and can provide corresponding planning guidance for the development of the tertiary industry in China's future city clusters.
However, this study still has some limitations.On the one hand, POI data only estimate and simulate different attributes in geospatial space, which is still significantly different from the actual urban construction and population distribution.On the other hand, we lack the time dimension to analyze the factors influencing the tertiary industry in the PRD city clusters.In future studies, we will try to use longer time series data to determine the spatial differentiation of service industries and their influencing factors in the PRD city clusters and to provide a practical basis for the future development of the PRD city clusters.

Conclusions
The adjustment and optimization of the tertiary industry structure is an important prerequisite for achieving high-quality development of the national economy.This study classifies the tertiary industry into different types, spatializes socio-economic data, and uses multi-source big data for quantitative analysis of influencing factors to obtain the regional differences in development levels and influencing factors of the tertiary industry in the PRD urban agglomeration in general and each type.This study then proposes the adjustment and optimization recommendations for the internal structure of the tertiary industry in the PRD urban agglomeration.The results show that the overall regional differences in the level of tertiary industry development in the PRD city clusters have continued to rise, with inter-regional differences being the main source; the structure of the tertiary industry is generally stable, with uncoordinated layout in localized areas; and the interaction between population and employment factors is the main influence on the layout of the tertiary industry in the PRD city clusters.With the backdrop of the country's highquality development, given that the total amount and proportion of the tertiary industry in the Pearl River Delta consistently rank among the nation's leaders, its future development should no longer be focused on pursuing an increase in the total amount.Instead, the emphasis should be on optimizing and upgrading its structure.Therefore, it is necessary to re-examine the role and differences of different types of tertiary industry to promote their high-quality development.

Figure 1 .
Figure 1.Scope map of the study area.Figure 1. Scope map of the study area.

Figure 1 .
Figure 1.Scope map of the study area.Figure 1. Scope map of the study area.

Figure 2 .
Figure 2. Various data types of the Pearl River Delta urban agglomeration.

Figure 2 .
Figure 2. Various data types of the Pearl River Delta urban agglomeration.

Figure 4 .
Figure 4. Regional differences in the development level of service industry in the Pearl River Delta from 2012 to 2022.(a) Intra-regional differences, (b) inter-regional differences, (c) sources of differences.

Figure 5 .
Figure 5. Evolution of spatial structure of service industry in the Pearl River Delta from 2012 to 2022.

Figure 5 .
Figure 5. Evolution of spatial structure of service industry in the Pearl River Delta from 2012 to 2022.

Figure 6 .
Figure 6.LISA Diagram of spatial agglomeration evolution of the Pearl River Delta urban agglomeration as a whole and various types of service industries from 2012 to 2022.

Figure 6 .
Figure 6.LISA Diagram of spatial agglomeration evolution of the Pearl River Delta urban agglomeration as a whole and various types of service industries from 2012 to 2022.

Table 1 .
POI Classification Statistics of Various Types of Tertiary Industry.

Table 2 .
Basic information on various types of research data.

Table 2 .
Basic information on various types of research data.

Table 3 .
The detection indicator system for influencing factors.

Table 4 .
Overall difference in service industry development level in the Pearl River Delta from 2012 to 2022.

Table 6 .
Factor detection results for various types of tertiary industry.
Note: Each influencing factor passed the significance test of p-value < 0.05.

Table 7 .
Results of interaction detection for each type of tertiary industry.
Note: Green, light green, light red and red, emphasize the values from low to high for each type of industry.