1. Introduction
According to the aging classification standards established by the United Nations in 1956 [
1], when the proportion of individuals aged 65 and above exceeds 7% in a country or region, it signifies the transition into an aging society; surpassing 14% indicates a deeply aged society; and crossing the threshold of 20% denotes entry into a super-aged society. According to reports from the United Nations, the global population is entering an aging phase, with the number and proportion of elderly populations increasing in almost every country. Population aging may emerge as one of the most significant social trends of the 21st century. However, as the degree of aging deepens, the supply of elderly care services fails to meet the increasing demands of the elderly population. Issues such as the lack of comprehensive planning and the irrational spatial allocation of elderly care facilities are becoming increasingly prominent [
2,
3,
4,
5,
6]. Therefore, addressing population aging trends effectively and optimizing the layout of elderly care facilities to promote supply–demand balance holds paramount significance.
Spatial equilibrium has long been considered a cornerstone in the planning research of elderly care facilities [
7,
8,
9,
10]. The spatial equilibrium of elderly care facilities (SEECF) implies that the optimal state is one where the supply of elderly resources within a specific spatial range matches the demands of the elderly population. Any spatial imbalances resulting from supply–demand disparities can be rectified by restructuring spatial arrangements and optimizing existing spaces [
11]. Current research on the SEECF varies in terms of the types of elderly care facilities, research methods, and spatial scale. Firstly, many studies focus solely on service-oriented facilities catering to the daily living, care, and recreational needs of the elderly population. Neglecting the essential medical and nursing care facilities required for the elderly daily life [
12,
13,
14,
15,
16,
17]. Secondly, in terms of research methodologies, there is a tendency to downplay the actual needs of the elderly population in favor of emphasizing the supply of elderly care facilities. As seen in methods such as kernel density estimation [
18,
19,
20], service coverage rates [
21], the Huff model [
22,
23,
24,
25], and accessibility assessments [
26,
27,
28,
29,
30,
31,
32,
33]. Additionally, in terms of scale, precise population data are fundamental components of effective urban management and strategic planning [
34,
35,
36]. However, many existing studies utilize administrative divisions, such as districts [
16,
37,
38,
39], streets [
40,
41,
42,
43], or communities [
38,
44,
45,
46], to represent population centers. This leads to less accurate simulations of real-world scenarios. The deficiencies in different types of research have brought severe challenges to the spatiotemporal sensitivity measurement of regional SEECF, making it difficult to conduct a refined study of the SEECF. The refined analysis of the SEECF is the basis for providing a decision-making basis for the scientific planning and management of elderly care facilities [
47,
48]. It is necessary to propose a fine-scale spatial equilibrium model of elderly care facilities that takes into account the types of elderly care facilities and the relationship between supply and demand.
Population aging plays a significant role in the transition of China’s economic growth. That is from the previous emphasis on the gross domestic product (GDP) scale expansion to a shift towards high-quality economic development. There is a wealth of research on the relationship between population aging and economic development [
49,
50]. Conflicting conclusions exist concerning internal mechanisms and consequences due to varying research perspectives. Some studies [
51,
52,
53,
54] find a significant negative impact of population aging on economic growth. In contrast, others [
55,
56,
57] argue that population aging may contribute to GDP growth. Additionally, some research [
58] suggests that with population aging, the GDP growth rate initially increases before decreasing. Moreover, existing studies on the relationship between elderly population factors and economic growth rarely integrate different attributes of the elderly population into the same research framework. Most studies only explore quantitative correlations, with limited research on spatial associations. Furthermore, the relationships among population aging, SEECF, and economic development remain unclear. Therefore, based on the panel threshold model, this study considers that population factors and facility allocation equilibrium simultaneously affect economic development. Thus, studying their influence patterns to provide multi-dimensional references for facility planning and sustainable economic development. While there is extensive research on the relationship between GDP and urban land use categories (ULUC) [
59,
60,
61,
62,
63,
64,
65,
66,
67], much of it focuses solely on construction land type [
68,
69,
70,
71]. Hence, this study further explores the spatial associations between the SEECF, GDP, and ULUC.
At present, research on the spatial equilibrium model of elderly care facilities is limited to service-type elderly care facilities. In the process of model construction, the focus is on the supply capacity of elderly care facilities, and the actual needs of the elderly population are weakened. The spatial scale of the research is mostly based on administrative divisions. This large-scale regional setting makes the spatial equilibrium evaluation of elderly care facilities low in spatiotemporal sensitivity. Most of the research on the SEECF and economic development has only explored the relationship between numerical distribution and rarely involved the analysis of spatial association relationships. At the same time, the impact of ULUC on the SEECF has not been considered. Therefore, it is urgent to construct a fine-scale spatial equilibrium model of elderly care facilities that takes into account different types of elderly care facilities and supply–demand relationships and further quantitatively analyze its spatial association with economic development and ULUC.
To solve the above problems, this study first took into account different types of elderly care facilities, such as service-oriented and medical care-oriented, and considered the supply capacity of elderly care facilities and the actual needs of the elderly population to construct a spatial equilibrium model of elderly care facilities with high spatiotemporal sensitivity (SEM-HSTS). At the 100-m fine grid scale, considering the supply capacity of elderly care facilities and the actual needs of the elderly population, we proposed constructing two factors, i.e., the spatial accessibility rate of elderly care services (SARecs) and the spatiotemporal supply–demand ratio for elderly care services (STSDRecs). By employing the modified two-step floating catchment area (M2SFCA), we obtained the spatiotemporal availability of medical services (STAms) factor. Facing the needs of per capita resource fairness and facility accessibility efficiency, the coordination degree model was introduced to calculate the SEECF with high spatiotemporal sensitivity. On this basis, considering the high spatiotemporal sensitivity of the SEM-HSTS, this study further explored its comprehensive association with economic development. Considering that the elderly population factor and the SEECF will have an impact on GDP at the same time, we investigated the phased influence relationships among the population aging, SEECF, and GDP based on the threshold effect test. Through bivariate local spatial autocorrelation analysis and risk factor detection, we quantitatively analyzed the spatial associations among SEECF, GDP, and ULUC. This study can provide valuable insights into the precise planning and layout of elderly care facilities and furnish a scientific theoretical basis for the multi-dimensional management of such facilities within the context of joint economic development and urban spatial structure.
4. Results
In empirical research, we first computed the SEECF for 2010, 2015, and 2020 based on the model constructed in this study. We compared our approach with the classical accessibility assessment method, evaluating the sensitivity to spatiotemporal dynamics using three indicators, i.e., standard deviation, q-value of Geodetector, and the number of significant regions detected by local spatial autocorrelation. Furthermore, employing a panel threshold model, we identified the phased influence patterns of population aging, SEECF, and economic development. Lastly, utilizing factor detection of Geodetector along with bivariate spatial autocorrelation analysis, we investigated the spatial associations among SEECF, economic development, and ULUC.
4.1. SEM-HSTS
The spatial equilibrium of elderly care facilities (SEECF) is an index to quantify the balance of supply and demand distribution of elderly care facilities. It is a numeric variable with a value range of [0, 1]. The spatial equilibrium state of elderly care facilities (SESECF) is an index to evaluate the balance of supply and demand distribution of elderly care facilities, which is a text variable, and its correspondence with the SEECF is shown in
Table 6 below. The larger the value of the SEECF, the more imbalanced, and the smaller the value of the SEECF, the morebalanced. This indicates that there is a negative correlation between the SEECF and the SESECF.
To further quantitatively compare the spatial sensitivity of the two methods, this study evaluated three indicators, i.e., standard deviation, q-value of Geodetector, and the number of significant clusters in local spatial autocorrelation. These indicators represented three evaluation dimensions. Since the SEECF, obtained through computation, was dimensionless and ranged from [0, 1], normalization of the accessibility results was required for comparison. A higher standard deviation indicates greater dispersion, hence higher spatial sensitivity. The q-value of Geodetector represents spatial stratified heterogeneity, with a range of [0, 1]. Moreover, calculations were made using the Graphical detector tool in the QGIS software. A higher q-value indicates more pronounced spatial stratified heterogeneity, implying higher spatial sensitivity [
78]. The number of significant clusters in local spatial autocorrelation represents local spatial heterogeneity, with a greater quantity indicating more evident local spatial heterogeneity [
79], assessed using the Local Moran’s I method. Calculations were made using the Univariate Local Moran’s I tool in the GeoDa software. In addition to spatial sensitivity, this study also investigated the temporal sensitivity of SEECF and accessibility over different periods using the same indicators. The Differential Local Moran’s I method was employed for local spatial autocorrelation.
4.1.1. Distribution of SEECF and Comparison of Spatial Sensitivity between SEECF and Accessibility
The classic accessibility assessment employed the M2SFCA method [
72] for computation, contrasting with the spatiotemporal sensitivity of the SEM-HSTS proposed in this study. In ArcGIS, experimental results depicted the SEECF and accessibility distribution maps for the central urban area of Wuhan in 2010, 2015, and 2020 (
Figure 3).
Figure 3a–c illustrate the SEECF distribution, while
Figure 3d–f depict the accessibility distribution. The SEECF was categorized based on computed values, with thresholds delineating balanced from imbalanced conditions. According to international standards, the number of nursing home beds per thousand elderly individuals typically ranges from 40 to 50. The Chinese Ministry of Civil Affairs set a target of 40 beds per thousand elderly individuals in 2020, while the Wuhan city government aimed for five beds per hundred elderly individuals. Building upon the analysis using the M2SFCA method and based on relevant literature [
16], we further classified the accessibility of elderly care facilities into six levels based on the number of beds available per 100 elderly individuals, i.e., strong (>5), relatively strong (4–5), slightly strong (3–4), slightly weak (2–3), relatively weak (1–2), and weak (0–1). In traditional evaluation methods [
80], the degree of supply–demand matching was directly represented by the accessibility of elderly care facilities, corresponding to SEECF.
- (1)
Dissimilarity in the distribution of SEECF.
Figure 3a–c show that over time, the red areas representing spatial imbalance had increased. Overall, the SEECF exhibited a slight tendency towards imbalance. Calculations revealed that the proportion of imbalanced areas had risen from 61.89% to 84.24%. However, significant differences in spatial balance persisted among regions, with distribution remaining uneven. Specifically, in
Figure 3a, the core urban areas in 2010 were predominantly imbalanced, primarily due to a supply of elderly care facilities significantly lower than the demand from the population aging, with slightly better spatial balance observed in more peripheral areas. In
Figure 3b for 2015, the southwestern and some peripheral areas showed better balance compared to the central core areas, attributed to the continuous improvement in the quantity and accessibility of facilities alongside the development of elderly care facilities and road networks. However, the growth rate of the elderly population demand in these regions was slower than in the core areas. In
Figure 3c for 2020, the balance shifted towards the core areas, mainly due to policy inclinations.
- (2)
Comparison of spatial sensitivity between SEECF and accessibility.
While the overall findings align with existing research [
16,
40,
42,
80,
81] and accessibility experiments, there were also divergent results. Taking 2020 as an example, regions marked with green circles in
Figure 3c indicated spatial imbalance of SEECF, whereas strong accessibility in
Figure 3f implied spatial balance. Further analysis, incorporating elderly care facilities kernel density and elderly population distribution (
Figure 4), revealed these areas exhibited an abundance supply of beds but relatively low elderly population numbers, indicative of spatial imbalance resulting from resource inefficiency. Consequently, it preliminarily revealed that the spatial distribution of SEECF proposed in this study was more complex, with significant differences extending beyond the relatively vague phenomena observed at the village and street levels. This advantage stemmed not only from the analysis at the high spatial resolution but also from the careful consideration of the following three factors proposed in this study: SARecs, STSDRecs, and STAms, which accounted for the attenuation of service capacity of elderly care facilities at different scales and the actual needs of the elderly population, thus demonstrating high spatial sensitivity. Based on these findings, precise spatial planning and optimization of existing space for elderly care facilities can be supported with robust geographical information.
The evaluation of the spatial sensitivity of SEECF and accessibility (
Table 7) revealed that the standard deviation of SEECF was greater than that of accessibility each year, indicating higher dispersion compared to accessibility. Additionally, the q-value for SEECF was consistently higher than that for accessibility, signifying more pronounced spatial stratified heterogeneity. With a total grid count of 73,289, the number of significant clusters of SEECF detected by Local Moran’s I exceeded that of accessibility, indicating more pronounced local spatial heterogeneity compared to accessibility. Altogether, these findings demonstrated that the SEM-HSTS constructed in this study exhibited higher spatial sensitivity.
4.1.2. Comparison of Temporal Sensitivity between SEECF and Accessibility
The temporal sensitivity evaluation was conducted for the periods 2010 to 2020, 2010 to 2015, and 2015 to 2020 (
Table 8). The results indicated that for any given period, the standard deviation of SEECF exceeded that of accessibility, as did the q-value, and the number of significant clusters detected by Differential Local Moran’s I was greater than that of accessibility, implying greater temporal heterogeneity compared to accessibility. These findings collectively demonstrated the higher temporal sensitivity of SEM-HSTS. This advantage was primarily attributed to the consideration of changes in service capacity (bed count) and elderly population demand (elderly population) over time within the factors STSDRecs and STAms.
4.2. Analysis of the Associations between SEECF and Economic Development
4.2.1. Influence Relationships between SEECF and Economic Development
- (1)
Threshold effect test of population aging on GDP.
Experiments were performed in Stata software.
Table 9 presents the results of the threshold effect model test for the influence of population aging on GDP. Generally, a
p-value < 0.05 indicates the existence of a threshold. The results indicated a single threshold effect test with a
p-value of 0.31, suggesting the absence of a threshold effect in the impact of population aging on GDP. The F-value represents the noise level, where a higher value indicates a smaller proportion of noise to the signal. Critical values at 1%, 5%, and 10% indicate the LR statistic accepted at the corresponding levels [
49].
Although there was no threshold effect in the influence of population aging on GDP, implying no significant difference in the effect of population aging on GDP under different SEECF, the results (
Table 10) showed negative coefficients for the variables in the right column, indicating an inhibitory effect of population aging on GDP. Additionally, it could be observed that under the backdrop of population aging, the coefficient for population size was 0.8974, while for transportation infrastructure, it was 0.2012. This suggested that the influence of population size on GDP was more potent than transportation infrastructure’s.
- (2)
Threshold effect test of SEECF on GDP.
Table 11 presents the results of the threshold effect model test for the influence of the SEECF on GDP. The results indicated a single threshold effect test with a
p-value of 0.0133 and a double threshold effect test with a
p-value of 0.06, suggesting that the influence of SEECF on GDP existed only as a single threshold effect.
Table 12 displays the effect of the SEECF on GDP, with population aging as the threshold variable, where the coefficients are listed in the right column. When the degree of population aging was ≤0.0637 (indicating a pre-aging society), the SEECF positively affected GDP, with a coefficient of 0.0247. As the SEECF and the spatial equilibrium state of elderly care facilities (SESECF) are inversely related, SESECF exerted an inhibitory effect on GDP, with a more substantial inhibitory effect as it approached spatial balance. Conversely, when the degree of population aging was >0.0637 (approaching an aging society), the SEECF exerted an inhibitory effect on GDP, with a coefficient of −0.0079. In this scenario, the SESECF promoted GDP, with a more substantial promotional effect as it tends towards spatial balance. The panel threshold model utilized time-series data from 2010 to 2020, enabling the capture of detailed effects over a comprehensive period due to the high temporal sensitivity of SEM-HSTS. These conclusions not only revealed the relationships between the SEECF and GDP but also demonstrated the staged influence patterns under different stages of aging, providing insights for the planning of elderly care facilities considering economic development prerequisites and population aging background. Additionally, it was observed that under the backdrop of SEECF, the coefficient for population size was 0.9127, while for transportation infrastructure, it was 0.1973, indicating that the influence of population size on GDP consistently outweighed that of transportation infrastructure.
4.2.2. Spatial Associations between SEECF and Economic Development
Following the identification of the phased influence relationships between the SEECF and economic development, this study explored their spatial associations. Initially, based on the bivariate spatial autocorrelation analysis method, an analysis of the SEECF and GDP in the central urban area of Wuhan in 2020 was conducted. Experiments were performed in GeoDa software. The global spatial autocorrelation analysis yielded a Moran’s I value of −0.323 with a
p-value of 0.001. The distribution of points along the horizontal and vertical axes was relatively uniform, indicating a significant negative association between the SEECF and GDP in the central urban area of Wuhan in 2020 within a 99% confidence interval. The results also highlighted the importance of investigating the spatial associations between the two variables. Local Spatial Autocorrelation Analysis results in obtaining a BiLISA cluster map (
Figure 5). LISA (Local Indicators of Spatial Association) is a method for spatial data analysis to identify spatial cluster patterns. Based on the similarity of data in geographic space, it calculates the local spatial correlation index of each region so as to judge the local spatial autocorrelation.
Figure 5 revealed significant “high-high” clustering, indicating a correlation between high GDP and spatial imbalance in some core areas and transitional zones, accounting for 13.93% of the total. Conversely, “low-low” clustering, denoting a correlation between low GDP and spatial balance, was observed in peripheral areas, representing 4.91% of the total. Additionally, “high-low” clustering, indicating a correlation between high GDP and spatial balance, was prevalent in most areas of the core region, accounting for 21.79% of the total. Furthermore, “low-high” clustering, indicative of a correlation between low GDP and spatial imbalance, was observed in peripheral areas and some transitional zones, representing 17.57% of the total. These phenomena preliminarily reflected the complexity of the spatial correlations between the SEECF and GDP. The areas with high GDP and spatial balance were primarily concentrated in the core region, followed by areas with low GDP and spatial imbalance, showing an uneven distribution. This conclusion confirmed the previously mentioned influence relationships and demonstrated its specific spatial distribution. It was the high spatial sensitivity of SEM-HSTS that enabled the precise identification of the differential distribution associated with GDP. Based on the spatial correlations conclusions, GDP development can be promoted by improving the SESECF in regions.
After the previous analysis, it can be seen that there are phased hierarchical impact relationships and spatial correlation differences between economic development and SEECF. At present, many studies have explored and found that there is a significant spatial association between GDP and different ULUC, so ULUC are likely to affect the spatial allocation of elderly care facilities. In this study, Geodetector was used to study the relative effects of ULUC, GDP, and the interaction between the two on the SESECF. The temporal changes are depicted in
Figure 6. It could be observed that over time, and with changes in the SESECF, the influence of GDP on the SESECF was more significant than that of ULUC. Additionally, the combined interaction of ULUC and GDP had a more substantial influence on the SESECF than any single factor alone. This underscored the necessity of considering ULUC when exploring the spatial associations between the SEECF and economic development.
In addition, risk area detection, risk factor detection, and ecological detection all indicated significant differences in the influence of GDP versus ULUC on the SEECF, with the former exhibiting a more decisive influence. We further explored the spatial associations in 2020. The spatial distribution maps of SEECF, GDP, and ULUC in the central urban area of Wuhan in 2020 are depicted in
Figure 7.
Figure 7c showed that ULUC were diverse and spatially unevenly distributed, with variations in land area. Therefore, the high spatial sensitivity of the SEM-HSTS provided a superior and essential prerequisite for investigating the spatial associations under different ULUC, enabling precise and distinct conclusions. Statistical analysis of SEECF and GDP by ULUC yields
Table 13 and
Figure 8. Analyzing in conjunction with
Figure 7, it was evident that the type with the lowest average value of SEECF, indicating the most balanced, was transportation stations. Simultaneously, this type exhibited the highest average GDP, corresponding to the phenomenon of spatial balance in the core areas correlated with high GDP as mentioned earlier. Similar patterns were observed for business office type, commercial service type, and sport and cultural type. Conversely, the urban land type with the highest average value of SEECF, indicating the most imbalanced, was the airport facility type. Additionally, this type showed the lowest average GDP, corresponding to the phenomenon of spatial imbalance in the peripheral areas with low GDP correlation as mentioned earlier. Similar patterns were observed for medical type, administrative type, and park and green type. Industrial types near the Qingshan district exhibited spatial imbalance correlated with high GDP, while those near Hanyang district exhibited spatial balance correlated with low GDP. Residential type in the Jiang’an district exhibited spatial imbalance correlated with high GDP, whereas most residential type in other areas exhibited spatial balance correlated with low GDP. Educational type showed no significant correlation between SEECF and GDP. The maximum difference in average values of SEECF among all ULUC was approximately 0.15, with all being in a state of imbalance, indicating an insignificant difference. This also indirectly reflected that the influence of ULUC on SEECF was smaller than that of GDP.