Is Urbanization Good for the Health of Middle-Aged and Elderly People in China?—Based on CHARLS Data

The purpose of this paper is to test whether improved healthcare services can mitigate health hazards resulting from environmental pollution in the urbanization process. Specifically, using China Health and Retirement Longitudinal Study (CHARLS) data and official statistics, this paper constructs comprehensive urbanization indicators and healthcare service indicators by applying the fully arrayed polygonal graphical indication method. Then, we introduce healthcare service indicators, urbanization indicators, environmental pollution indicators, and the interaction term between environmental pollution and healthcare into an ordered-logistics regression model. Our results indicate that improvement in health conditions can decrease the health risks from multiplied emissions of industrial sulfur dioxide, industrial soot and dust, and industrial effluents, but it cannot counteract the adverse health effects of PM2.5. Furthermore, heterogeneity tests show that, when considering the multidimensionality of urbanization, the positive influence of healthcare is the greatest in residential surroundings urbanization and economic urbanization, which reduces the prevalence of chronic diseases by 18.4% and 14.9%, respectively. Among the diverse city types, mixed-economy cities have the most obvious positive effects, where healthcare has the greatest mitigating effect on the health damage caused by industrial sulfur dioxide and industrial soot and dust, decreasing the prevalence of chronic diseases among the middle-aged and elderly by 27.3% and 16.4%, respectively. When considering the regional impacts of urbanization, there is a large difference in the positive effects brought about by medical care, which is reflected mainly in eastern and western China. In eastern China, although healthcare does not offset the health damage of PM2.5, the increase in chronic diseases among the middle-aged and elderly is only 0.5%, while in western China, the increase rises to 22.4%.


Introduction
Since the 1990s, both the speed and scale of China's urbanization have accelerated [1]. During this period, the urbanization rate in China has increased rapidly from 30.48% to 60.60%, an increase of nearly 31%. It has been reported that the urbanization rate in China will reach 75% by 2030 (The China urbanization 2.0 report). According to that result, it can be seen that, in several decades, China will have achieved the urbanization progress that occurred in western countries over hundreds of years [2]. Such rapid urbanization growth poses large challenges to public health. Specifically, quick urbanization leads to serious problems of environmental pollution and unhealthy lifestyles [3]. Among these problems, the most widespread impact has been on air pollution [4], for example, increased levels of industrial sulfur dioxide [5] and particulate matter (PM 10 , PM 2.5 ) pollution [6]. PM 2.5 Sustainability 2021, 13, 4996 2 of 20 pollution has had a particular impact, according to a previous study; both long-term and short-term exposure to PM 2.5 increased the probability of chronic diseases for residents, such as respiratory disease and cardiovascular morbidity [7]. It is even reported that 48.6% of the Chinese population (nearly 100 million people) suffer from obstructive pulmonary disease (COPD), and 18.7% of COPD deaths are attributable to environmental PM 2.5 exposure (China Pulmonary Health Study). At the same time, the hazards of wastewater to human health cannot be ignored [8]. Current research indicates that water pollution is associated not only with acute waterborne diseases (cholera, diarrhea) [9] but also with cancer risk in severe cases [10]. The poor lifestyles and unhealthy diets resulting from urbanization also introduce many risks for residents' health [11], such as obesity [12], hypertension [13], and heart disease [14]. Despite the negative impacts of urbanization, it also leads to an improved standard of living, better healthcare resources, increased job opportunities, and higher education levels, which lead to greatly improved public health [6]. However, the most sensitive of the impacted groups are middle-aged and older adults. According to statistics, China's elderly population accounted for 18.1% of the total population in 2019, and the proportion is expected to rise to 20% by 2030. At that time, China will enter a heavily aging society. Furthermore, with an accelerating rate of increase in the proportion of the population that is aged [15], chronic diseases among the elderly are a particular problem in China [16]. According to statistics, nearly 80% of deaths in the elderly population are caused by chronic diseases, and the elderly are 3.2 times more likely to suffer from chronic diseases than other groups [17], with environmental pollution factors increasing the probability of chronic diseases in the elderly. Researchers have shown that in Wuhan, China, a 10 µg/m 3 increase in NO 2 is associated with a 1.6% increase in CVD mortality in the elderly [18]. Some scholars have even predicted the future trend of CVD in China; by 2030, the number of CVD events per year in China will increase by more than 50% due to population aging and population growth alone [19], of which environmental pollution is a factor that cannot be ignored. These statistics demonstrate the importance of studying the health effects of urbanization on middle-aged and elderly people in China. Based on this, are the positive or negative effects of urbanization on the health of older adults greater?
The influence that urbanization has on health is complex [20,21]. Some scholars argue that urbanization adversely affects the health of residents [20,22]. Previous research has considered the negative health effects of single factors such as changes in environmental pollution [3,23], unhealthy lifestyles, and socioeconomic status [22] caused by urbanization. Among these factors, researchers reported that environmental pollution has the greatest influence on health. Other scholars believe that in the process of urbanization rising living standards, improved healthcare resources, more job opportunities, and higher educational levels have dramatically enhanced public health [6,24]. In particular, it should be noted that the advancement of sanitary conditions brought about by urbanization has made it easier for people to enjoy the latest achievements in medical technology and to obtain better medical care, which has benefited city residents and promoted their health. In conclusion, it can be seen that there are disparate impacts of urbanization on health.
In studies of the implications of urbanization for health, most of the scholars have focused on single indexes such as the demographic urbanization rate [25], constructing the urbanization indicators of broad community characteristics [22], nighttime light data [26], and so on. Even when researchers have taken into account the complexity and multidimensionality of urbanization, they have only briefly adopted multiple indexes or have used data that were not representative. For example, Liu et al. [27] considered the effects of differing levels and rates of urbanization on health by using indexes of land-use transitions, the growth of the economy, population clustering, and health services, and Chen et al. [26] adopted nighttime databases to assess the effects of urbanization of different dimensions, development speed, and level of health in county-level regions. However, only some scholars consider the multidimensionality of urbanization and employ the entropy method [28] or the fully arrayed polygonal graphical indication method [29] to divide urbanization into the urbanization of the population, urbanization of the economy, urbanization of Sustainability 2021, 13, 4996 3 of 20 residential environments, and urbanization of residential conditions. Although the selected methods and indicators can objectively reflect the urbanization of all of China, most of them are only analyzed at the provincial level and not at the urban level. More importantly, these measurements of urbanization do not apply to the impact of urbanization on health. Therefore, this paper adopts the fully arrayed polygonal graphical indication method to reconstruct the comprehensive urbanization index for China's prefecture-level cities and incorporate it into a model of urbanization and health, as well as dividing it into four aspects: demographic urbanization, economy urbanization, residential surroundings urbanization, and residential condition urbanization.
Most scholars studying healthcare service indicators use single indexes, such as healthcare expenditures [30], the number of hospital beds [31,32], and the number of doctors in hospitals [33]. However, medical resources are not only a reflection of these singular indexes, but are also closely related to financial resources, material resources, and manpower for healthcare. Therefore, some researchers have employed the entropy method [34] to comprehensively evaluate healthcare resources, but these methods are mostly applied at the provincial scale in China; the data are not easily available at the municipal level. Thus, taking into account the availability of data and the applicability of the method, this paper employs the fully arrayed polygonal graphical indication method [29] to comprehensively evaluate the level of medical and health resources by selecting three indicators: the number of hospital beds per thousand people, the number of medical and health institutions per thousand people, and the number of hospitals per hundred square kilometers.
To summarize, this paper adopts the fully arrayed polygonal graphical indication method to calculate composite urbanization indexes and healthcare indicators. In this method, we construct the comprehensive urbanization indicators based on the four dimensions of demographic urbanization, economy urbanization, residential surroundings urbanization, and residential condition urbanization. Next, the urbanization indicators, the health service indicators, and the interaction terms of environmental pollution and health services are incorporated into an ordered logistic model to explore the advantages and disadvantages that change in healthcare conditions and environmental pollution caused by urbanization and see what effect they have on the health of middle-aged and elderly people. The heterogeneity of urbanization in multidimensionality, regional differences, and city types are considered as well. The primary contributions of this paper are as follows: First, we incorporate urbanization, environmental pollution, healthcare services, and health into the same research framework to analyze the dual effects and the magnitude of the impacts of environmental pollution and the improvement of healthcare service caused by urbanization. Then, the comprehensive urbanization and healthcare service indexes measured by the fully arrayed polygonal graphical indication method are integrated into an ordered logistics model of the impact of urbanization on health, and the interaction term between the environmental pollution and healthcare service is introduced to study the negative and positive impacts of urbanization on health and the extent of those effects. Finally, the comprehensive urbanization indicators calculated by the fully arrayed polygonal graphical indication method are used to divide the integrated urbanization into four dimensions: urbanization of the population, urbanization of the economy, urbanization of residential environment, and urbanization of residential conditions. This is then used to explore the heterogeneity of the influence of urbanization on health based on the different dimensions of urbanization, the regional variability of urbanization, and the variability of city types. The manuscript is arranged as follows: the second part contains the data and methods, the third part contains the results and discussion, and the last part contains the conclusion and policy implications. Most scholars apply the rate of demographic urbanization to calculate the urbanization level [25]. However, urbanization is a complex concept with a variety of dimen-sions and characteristics [20], including not only demographic aspects but also economic, cultural, social, and ecosystem aspects. Thus, the comprehensive index approach can better demonstrate the integrative characteristics of urbanization. This study reconstructs the consolidated indicator system for measuring urbanization based on the research of Wu et al. [28,29]. The index system can be classified into three layers: The first layer is the target level, that is, the comprehensive urbanization level. The second layer is the criterion level, which contains the indicators demographic urbanization, economy urbanization, residential surroundings urbanization, and residential condition urbanization. The third layer is the indices level, which covers 16 sub-indicators (as shown in Table 1). We measure the integrated urbanization layer using the fully aligned polygonal graphical indexing method. Table 1 presents detailed information on the comprehensive urbanization indicator system. As described in Section 2.1.1, the integrated urbanization approach employs the fully aligned polygonal graphical indexing method to deal with the composite urbanization in China.

Methods and
The fundamental principle of the method is that N indicators are set up and the higher limit value of these indexes is taken as the radius to construct a central N-side shape. In other words, the connecting line of each index value is used to construct an uneven central N-polygon. The vertex of the n-index is the full permutation of N indexes, which can form (n − 1)!/2 distinct non-regular central n-squares. The complex indices are specified as the ratio of the mean area of all these non-regular polygons to the area of the central polygon. The specific calculation process is as follows: (1) The hyperbolic standardization function is applied to standardize the index values: where a, b, and c represent the parameters of the hyperbolic function, respectively.
(2) The hyperbolic standardized function fulfills F(U) = 1, F(T) = 0, F(L) = −1, where U is the higher limit of the exponent X, L is the lower limit of X, and T is the critical value of the X. The threshold value can be expressed as the exponential mean.
From Equation (2), the standardized function, F(x), maps the index value located in [L, U] to [−1, 1]. The normalization procedure will lead to a fast-slow-fast nonlinear growth trend of the indicator values.
(3) For the ith index, the single index value, S i , is: The vertex of the nth-edge consists of the value at S i = 1, the hub point consists of the value at S i = −1, and the threshold value of the polygon index is the value at S i = 0. If the indicator value is above the threshold value, the value of each indicator is positive; otherwise, it is negative.
(4) The fully aligned polygonal graphical indexing method is indicated as follows:

The Model of the Relationship between Urbanization and Health
To explore the dual effects that the factors of improved medical care service and increased environmental pollution from urbanization have on chronic diseases in middleaged and elderly people, and considering that the dependent variable is a discrete order variable, this research applies the ordered-logistics model. The constructed baseline model is as follows: In this formula, H represents the disease status of the micro individual; Urban is the urbanization index, including comprehensive urbanization (CU), demographic urbanization (DU), economy urbanization (EU), residential surroundings urbanization (RSU), and residential condition urbanization (RCU); AP j represents environment pollutants (industrial sulfur dioxide, industrial soot and dust, industrial wastewater) and the yearly average concentration of PM 2.5 ; Mhs is the medical and health services index; X i is the vector of the individual-level control variables, including age, sex, education, material, etc.; δ i is a random perturbation term; and β 1 , β 2 , and β 3 represent the effects of urbanization, medical services, and environmental pollution, respectively, on the health status of the micro individual.
To further study the effect of environmental pollution and healthcare service on the health of middle-aged and elderly people, and consider which is more important, this paper introduces an interaction term of environmental pollution and medical health service to construct the following model: In the formula, AP j × Mhs is the interaction term of environmental pollution (industrial sulfur dioxide, industrial soot and dust, industrial wastewater) and medical care service.
To research the interaction term using the ordered-logistics model, our study calculates the partial effect models of environmental pollution and medical care service on the health status of middle-aged and elderly people to specify the magnitude of this effect. The specific bias effects model is as follows: health of middle-aged and elderly people, and consider which is more important, this paper introduces an interaction term of environmental pollution and medical health service to construct the following model: In the formula, × ℎ is the interaction term of environmental pollution (industrial sulfur dioxide, industrial soot and dust, industrial wastewater) and medical care service.
To research the interaction term using the ordered-logistics model, our study calculates the partial effect models of environmental pollution and medical care service on the health status of middle-aged and elderly people to specify the magnitude of this effect. The specific bias effects model is as follows: where, is the probability of the individual health rank and ℊ μ − X − ℊ μ − X are the partial effects of the transformation. In the positive and negative examination of ℊ μ − X − ℊ μ − X β , if the value is above zero, it indicates that the improved healthcare conditions cannot counter the adverse effects of environmental pollution on middle-aged and elderly person health, while if it is less than zero, it shows that the improved healthcare service can reduce the negative effects of environmental pollution.

Variables
As described in the above model and analysis, the core explanatory variables in the paper are urbanization, medical and health service, and environmental pollution. For the urbanization indicator, the fully arrayed polygonal graphical indication method is used to calculate comprehensive urbanization. Environmental pollution is defined based on per capita industrial wastewater emissions; per capita industrial sulfur dioxide emissions; per capita industrial soot and dust emissions; and the annual average concentration of PM2.5 [35] because it is not only closely related to economic activities such as industrial structure [36], energy consumption [37], and international trade [38] but is also a key factor in public health [39]. More importantly, in order to maintain the original features of the data and to reduce or eliminate heteroscedasticity in the data, the indicators are presented in their natural logarithms. The healthcare service level is measured using the fully arrayed polygonal graphical indication method for three indicators: the number of beds per thousand people, the number of healthcare institutions per thousand people, and the number of hospitals per hundred square kilometers. In addition, the individual characteristic var- health of middle-aged and elderly people, and consider which is more importa paper introduces an interaction term of environmental pollution and medical hea vice to construct the following model: In the formula, × ℎ is the interaction term of environmental pollution trial sulfur dioxide, industrial soot and dust, industrial wastewater) and medical c vice.
To research the interaction term using the ordered-logistics model, our stud lates the partial effect models of environmental pollution and medical care servic health status of middle-aged and elderly people to specify the magnitude of thi The specific bias effects model is as follows: where, is the probability of the individual health rank and ℊ μ − ℊ μ − X are the partial effects of the transformation. In the positive and n examination of ℊ μ − X − ℊ μ − X β , if the value is above zero, it in that the improved healthcare conditions cannot counter the adverse effects of e mental pollution on middle-aged and elderly person health, while if it is less than shows that the improved healthcare service can reduce the negative effects of e mental pollution.

Variables
As described in the above model and analysis, the core explanatory variable paper are urbanization, medical and health service, and environmental pollution. urbanization indicator, the fully arrayed polygonal graphical indication method to calculate comprehensive urbanization. Environmental pollution is defined based capita industrial wastewater emissions; per capita industrial sulfur dioxide emissi capita industrial soot and dust emissions; and the annual average concentration [35] because it is not only closely related to economic activities such as industrial s [36], energy consumption [37], and international trade [38] but is also a key factor i health [39]. More importantly, in order to maintain the original features of the dat reduce or eliminate heteroscedasticity in the data, the indicators are presented natural logarithms. The healthcare service level is measured using the fully arra lygonal graphical indication method for three indicators: the number of beds pe sand people, the number of healthcare institutions per thousand people, and the of hospitals per hundred square kilometers. In addition, the individual characteri µ j − X α β 1 + health of middle-aged and elderly people, and consider whic paper introduces an interaction term of environmental pollutio vice to construct the following model: In the formula, × ℎ is the interaction term of enviro trial sulfur dioxide, industrial soot and dust, industrial wastew vice.
To research the interaction term using the ordered-logistic lates the partial effect models of environmental pollution and m health status of middle-aged and elderly people to specify the The specific bias effects model is as follows: where, is the probability of the individual health ran ℊ μ − X are the partial effects of the transformation. In examination of ℊ μ − X − ℊ μ − X β , if the value that the improved healthcare conditions cannot counter the a mental pollution on middle-aged and elderly person health, wh shows that the improved healthcare service can reduce the ne mental pollution.

Variables
As described in the above model and analysis, the core ex paper are urbanization, medical and health service, and enviro urbanization indicator, the fully arrayed polygonal graphical to calculate comprehensive urbanization. Environmental pollut capita industrial wastewater emissions; per capita industrial su capita industrial soot and dust emissions; and the annual aver [35] because it is not only closely related to economic activities s [36], energy consumption [37], and international trade [38] but i health [39]. More importantly, in order to maintain the original reduce or eliminate heteroscedasticity in the data, the indicat natural logarithms. The healthcare service level is measured u lygonal graphical indication method for three indicators: the sand people, the number of healthcare institutions per thousan of hospitals per hundred square kilometers. In addition, the ind health of middle-aged and elderly people, an paper introduces an interaction term of enviro vice to construct the following model: In the formula, × ℎ is the interactio trial sulfur dioxide, industrial soot and dust, in vice.
To research the interaction term using the lates the partial effect models of environmenta health status of middle-aged and elderly peop The specific bias effects model is as follows: where, is the probability of the indivi ℊ μ − X are the partial effects of the t examination of ℊ μ − X − ℊ μ − X that the improved healthcare conditions cann mental pollution on middle-aged and elderly p shows that the improved healthcare service c mental pollution.

Variables
As described in the above model and ana paper are urbanization, medical and health ser urbanization indicator, the fully arrayed polyg to calculate comprehensive urbanization. Envir capita industrial wastewater emissions; per cap capita industrial soot and dust emissions; and [35] because it is not only closely related to econ [36], energy consumption [37], and internationa health [39]. More importantly, in order to main reduce or eliminate heteroscedasticity in the natural logarithms. The healthcare service lev lygonal graphical indication method for three sand people, the number of healthcare institut of hospitals per hundred square kilometers. In (7) ∂P i ∂Mhs = medical services, and environmental pollution, respectively, on the health status of the micro individual.
To further study the effect of environmental pollution and healthcare service on the health of middle-aged and elderly people, and consider which is more important, this paper introduces an interaction term of environmental pollution and medical health service to construct the following model: In the formula, × ℎ is the interaction term of environmental pollution (industrial sulfur dioxide, industrial soot and dust, industrial wastewater) and medical care service.
To research the interaction term using the ordered-logistics model, our study calculates the partial effect models of environmental pollution and medical care service on the health status of middle-aged and elderly people to specify the magnitude of this effect. The specific bias effects model is as follows: where, is the probability of the individual health rank and ℊ μ − X − ℊ μ − X are the partial effects of the transformation. In the positive and negative examination of ℊ μ − X − ℊ μ − X β , if the value is above zero, it indicates that the improved healthcare conditions cannot counter the adverse effects of environmental pollution on middle-aged and elderly person health, while if it is less than zero, it shows that the improved healthcare service can reduce the negative effects of environmental pollution.

Variables
As described in the above model and analysis, the core explanatory variables in the paper are urbanization, medical and health service, and environmental pollution. For the urbanization indicator, the fully arrayed polygonal graphical indication method is used to calculate comprehensive urbanization. Environmental pollution is defined based on per capita industrial wastewater emissions; per capita industrial sulfur dioxide emissions; per capita industrial soot and dust emissions; and the annual average concentration of PM2.5 [35] because it is not only closely related to economic activities such as industrial structure [36], energy consumption [37], and international trade [38] but is also a key factor in public health [39]. More importantly, in order to maintain the original features of the data and to reduce or eliminate heteroscedasticity in the data, the indicators are presented in their natural logarithms. The healthcare service level is measured using the fully arrayed polygonal graphical indication method for three indicators: the number of beds per thousand people, the number of healthcare institutions per thousand people, and the number of hospitals per hundred square kilometers. In addition, the individual characteristic var- medical services, and environmental pollution, respectively, on the health statu micro individual.
To further study the effect of environmental pollution and healthcare servic health of middle-aged and elderly people, and consider which is more import paper introduces an interaction term of environmental pollution and medical he vice to construct the following model: In the formula, × ℎ is the interaction term of environmental pollution trial sulfur dioxide, industrial soot and dust, industrial wastewater) and medical vice.
To research the interaction term using the ordered-logistics model, our stud lates the partial effect models of environmental pollution and medical care servic health status of middle-aged and elderly people to specify the magnitude of th The specific bias effects model is as follows: where, is the probability of the individual health rank and ℊ μ − ℊ μ − X are the partial effects of the transformation. In the positive and examination of ℊ μ − X − ℊ μ − X β , if the value is above zero, it i that the improved healthcare conditions cannot counter the adverse effects of mental pollution on middle-aged and elderly person health, while if it is less tha shows that the improved healthcare service can reduce the negative effects of mental pollution.

Variables
As described in the above model and analysis, the core explanatory variabl paper are urbanization, medical and health service, and environmental pollution urbanization indicator, the fully arrayed polygonal graphical indication method to calculate comprehensive urbanization. Environmental pollution is defined base capita industrial wastewater emissions; per capita industrial sulfur dioxide emiss capita industrial soot and dust emissions; and the annual average concentration [35] because it is not only closely related to economic activities such as industrial s [36], energy consumption [37], and international trade [38] but is also a key factor health [39]. More importantly, in order to maintain the original features of the da reduce or eliminate heteroscedasticity in the data, the indicators are presented natural logarithms. The healthcare service level is measured using the fully arra lygonal graphical indication method for three indicators: the number of beds p sand people, the number of healthcare institutions per thousand people, and the of hospitals per hundred square kilometers. In addition, the individual character µ j − X α β 2 + medical services, and environmental pollution, respectively, micro individual.
To further study the effect of environmental pollution an health of middle-aged and elderly people, and consider whi paper introduces an interaction term of environmental polluti vice to construct the following model: In the formula, × ℎ is the interaction term of envir trial sulfur dioxide, industrial soot and dust, industrial wastew vice.
To research the interaction term using the ordered-logisti lates the partial effect models of environmental pollution and health status of middle-aged and elderly people to specify th The specific bias effects model is as follows: where, is the probability of the individual health ra ℊ μ − X are the partial effects of the transformation. I examination of ℊ μ − X − ℊ μ − X β , if the value that the improved healthcare conditions cannot counter the a mental pollution on middle-aged and elderly person health, w shows that the improved healthcare service can reduce the n mental pollution.

Variables
As described in the above model and analysis, the core e paper are urbanization, medical and health service, and enviro urbanization indicator, the fully arrayed polygonal graphical to calculate comprehensive urbanization. Environmental pollu capita industrial wastewater emissions; per capita industrial su capita industrial soot and dust emissions; and the annual aver [35] because it is not only closely related to economic activities [36], energy consumption [37], and international trade [38] but health [39]. More importantly, in order to maintain the origina reduce or eliminate heteroscedasticity in the data, the indica natural logarithms. The healthcare service level is measured u lygonal graphical indication method for three indicators: the sand people, the number of healthcare institutions per thousan of hospitals per hundred square kilometers. In addition, the in medical services, and environmental pollutio micro individual.
To further study the effect of environmen health of middle-aged and elderly people, an paper introduces an interaction term of enviro vice to construct the following model: In the formula, × ℎ is the interacti trial sulfur dioxide, industrial soot and dust, in vice.
To research the interaction term using th lates the partial effect models of environmenta health status of middle-aged and elderly peo The specific bias effects model is as follows: where, is the probability of the indiv ℊ μ − X are the partial effects of the t examination of ℊ μ − X − ℊ μ − X that the improved healthcare conditions cann mental pollution on middle-aged and elderly shows that the improved healthcare service c mental pollution.

Variables
As described in the above model and ana paper are urbanization, medical and health se urbanization indicator, the fully arrayed poly to calculate comprehensive urbanization. Envi capita industrial wastewater emissions; per ca capita industrial soot and dust emissions; and [35] because it is not only closely related to eco [36], energy consumption [37], and internation health [39]. More importantly, in order to main reduce or eliminate heteroscedasticity in the natural logarithms. The healthcare service lev lygonal graphical indication method for three sand people, the number of healthcare institut of hospitals per hundred square kilometers. In where, P i is the probability of the individual health rank and of the individual-level control variables, including age a random perturbation term; and , , and re medical services, and environmental pollution, respe micro individual.
To further study the effect of environmental pol health of middle-aged and elderly people, and cons paper introduces an interaction term of environment vice to construct the following model: In the formula, × ℎ is the interaction term trial sulfur dioxide, industrial soot and dust, industria vice.
To research the interaction term using the order lates the partial effect models of environmental pollu health status of middle-aged and elderly people to s The specific bias effects model is as follows: where, is the probability of the individual h ℊ μ − X are the partial effects of the transfor examination of ℊ μ − X − ℊ μ − X β , if that the improved healthcare conditions cannot cou mental pollution on middle-aged and elderly person shows that the improved healthcare service can redu mental pollution.

Variables
As described in the above model and analysis, t paper are urbanization, medical and health service, an urbanization indicator, the fully arrayed polygonal g to calculate comprehensive urbanization. Environmen capita industrial wastewater emissions; per capita ind capita industrial soot and dust emissions; and the an [35] because it is not only closely related to economic a [36], energy consumption [37], and international trade health [39]. More importantly, in order to maintain th reduce or eliminate heteroscedasticity in the data, th natural logarithms. The healthcare service level is m lygonal graphical indication method for three indica sand people, the number of healthcare institutions pe of hospitals per hundred square kilometers. In additio of the individual-level control variab a random perturbation term; and medical services, and environment micro individual.
To further study the effect of e health of middle-aged and elderly paper introduces an interaction term vice to construct the following mod ℎ = + + In the formula, × ℎ is th trial sulfur dioxide, industrial soot a vice.
To research the interaction term lates the partial effect models of env health status of middle-aged and e The specific bias effects model is as where, is the probability of ℊ μ − X are the partial effects examination of ℊ μ − X − ℊ that the improved healthcare cond mental pollution on middle-aged an shows that the improved healthcar mental pollution.

Variables
As described in the above mod paper are urbanization, medical and urbanization indicator, the fully arr to calculate comprehensive urbaniza capita industrial wastewater emissio capita industrial soot and dust emis [35] because it is not only closely rela [36], energy consumption [37], and i health [39]. More importantly, in ord reduce or eliminate heteroscedastic natural logarithms. The healthcare lygonal graphical indication metho sand people, the number of healthc of hospitals per hundred square kilo µ j − X α β i are the partial effects of the β i transformation. In the positive and negative examination of trial sulfur dioxide, industrial soot and dust, industrial wastewater) and the yearly average concentration of PM2.5; ℎ is the medical and health services index; is the vector of the individual-level control variables, including age, sex, education, material, etc.; is a random perturbation term; and , , and represent the effects of urbanization, medical services, and environmental pollution, respectively, on the health status of the micro individual.
To further study the effect of environmental pollution and healthcare service on the health of middle-aged and elderly people, and consider which is more important, this paper introduces an interaction term of environmental pollution and medical health service to construct the following model: In the formula, × ℎ is the interaction term of environmental pollution (industrial sulfur dioxide, industrial soot and dust, industrial wastewater) and medical care service.
To research the interaction term using the ordered-logistics model, our study calculates the partial effect models of environmental pollution and medical care service on the health status of middle-aged and elderly people to specify the magnitude of this effect. The specific bias effects model is as follows: where, is the probability of the individual health rank and ℊ μ − X − ℊ μ − X are the partial effects of the transformation. In the positive and negative examination of ℊ μ − X − ℊ μ − X β , if the value is above zero, it indicates that the improved healthcare conditions cannot counter the adverse effects of environmental pollution on middle-aged and elderly person health, while if it is less than zero, it shows that the improved healthcare service can reduce the negative effects of environmental pollution.

Variables
As described in the above model and analysis, the core explanatory variables in the paper are urbanization, medical and health service, and environmental pollution. For the urbanization indicator, the fully arrayed polygonal graphical indication method is used to calculate comprehensive urbanization. Environmental pollution is defined based on per capita industrial wastewater emissions; per capita industrial sulfur dioxide emissions; per capita industrial soot and dust emissions; and the annual average concentration of PM2.5 [35] because it is not only closely related to economic activities such as industrial structure [36], energy consumption [37], and international trade [38] but is also a key factor in public health [39]. More importantly, in order to maintain the original features of the data and to reduce or eliminate heteroscedasticity in the data, the indicators are presented in their natural logarithms. The healthcare service level is measured using the fully arrayed polygonal graphical indication method for three indicators: the number of beds per thousand people, the number of healthcare institutions per thousand people, and the number of hospitals per hundred square kilometers. In addition, the individual characteristic var- trial sulfur dioxide, industrial soot and dust, industrial wastewater) and the yearly aver age concentration of PM2.5; ℎ is the medical and health services index; is the vector of the individual-level control variables, including age, sex, education, material, etc.; is a random perturbation term; and , , and represent the effects of urbanization medical services, and environmental pollution, respectively, on the health status of the micro individual.
To further study the effect of environmental pollution and healthcare service on the health of middle-aged and elderly people, and consider which is more important, this paper introduces an interaction term of environmental pollution and medical health ser vice to construct the following model: In the formula, × ℎ is the interaction term of environmental pollution (indus trial sulfur dioxide, industrial soot and dust, industrial wastewater) and medical care ser vice.
To research the interaction term using the ordered-logistics model, our study calcu lates the partial effect models of environmental pollution and medical care service on the health status of middle-aged and elderly people to specify the magnitude of this effect The specific bias effects model is as follows: where, is the probability of the individual health rank and ℊ μ − X − ℊ μ − X are the partial effects of the transformation. In the positive and negative examination of ℊ μ − X − ℊ μ − X β , if the value is above zero, it indicates that the improved healthcare conditions cannot counter the adverse effects of environ mental pollution on middle-aged and elderly person health, while if it is less than zero, i shows that the improved healthcare service can reduce the negative effects of environ mental pollution.

Variables
As described in the above model and analysis, the core explanatory variables in the paper are urbanization, medical and health service, and environmental pollution. For the urbanization indicator, the fully arrayed polygonal graphical indication method is used to calculate comprehensive urbanization. Environmental pollution is defined based on per capita industrial wastewater emissions; per capita industrial sulfur dioxide emissions; per capita industrial soot and dust emissions; and the annual average concentration of PM2. [35] because it is not only closely related to economic activities such as industrial structure [36], energy consumption [37], and international trade [38] but is also a key factor in public health [39]. More importantly, in order to maintain the original features of the data and to reduce or eliminate heteroscedasticity in the data, the indicators are presented in their natural logarithms. The healthcare service level is measured using the fully arrayed po lygonal graphical indication method for three indicators: the number of beds per thou sand people, the number of healthcare institutions per thousand people, and the number of hospitals per hundred square kilometers. In addition, the individual characteristic var µ j − X α β 4 , if the value is above zero, it indicates that the improved healthcare conditions cannot counter the adverse effects of environmental pollution on middle-aged and elderly person health, while if it is less than zero, it shows that the improved healthcare service can reduce the negative effects of environmental pollution.

Variables
As described in the above model and analysis, the core explanatory variables in the paper are urbanization, medical and health service, and environmental pollution. For the urbanization indicator, the fully arrayed polygonal graphical indication method is used to calculate comprehensive urbanization. Environmental pollution is defined based on per capita industrial wastewater emissions; per capita industrial sulfur dioxide emissions; per capita industrial soot and dust emissions; and the annual average concentration of PM 2.5 [35] because it is not only closely related to economic activities such as industrial structure [36], energy consumption [37], and international trade [38] but is also a key factor in public health [39]. More importantly, in order to maintain the original features of the data and to reduce or eliminate heteroscedasticity in the data, the indicators are presented in their natural logarithms. The healthcare service level is measured using the fully arrayed polygonal graphical indication method for three indicators: the number of beds per thousand people, the number of healthcare institutions per thousand people, and the number of hospitals per hundred square kilometers. In addition, the individual characteristic variables are introduced as control variables, including age, sex, per capita household income, smoking, alcohol use, material, residence, education, and health insurance (as shown in Table 2).

Data
The data used for the urbanization and environmental pollution indexes originates from The China City Statistical Yearbook, The China Urban Construction Statistical Yearbook, the National Bureau of Statistics of China, and the Provincial Bureau of Statistics. The health and individual characteristic variables are derived from data in the China Health and Retirement Longitudinal Study (CHARLS) [40], which began in 2011 and was repeated every two years and was only updated to 2015. The data primarily focused on China's middle-aged and elderly over 45 years of age. The study covered a wide range of topics, including information on individual backgrounds, households, health, healthcare and insurance, work, retirement, pensions, income, and basic community information, etc. More information about the CHARLS can be found in Zhao et al. [41]. The summarized data can help meet the needs of scientific research. Table 3 reports statistics about variables in the sample. The processing of the data in question is described in detail in the Supplementary Materials.

Variable Definition
Composite urbanization (CU) Calculated based on demographic urbanization economy urbanization, residential surroundings urbanization, and residential condition urbanization.

Results and Discussion
This section shows the results of the regression models in estimating the middle-aged and elderly population chronic disease and urbanization, healthcare service, and environmental pollution indexes. First, our paper reports the baseline models; that is, the regression results for health and comprehensive urbanization, healthcare service, and environmental pollution. Second, based on the multidimensionality of the urbanization (demographic urbanism, economy urbanism, residential surroundings urbanism, and residential condition urbanism), the regional difference, and the difference of the city types, the paper explores the heterogeneity on the influence of healthcare service and environmental population on middle-aged and elderly people's health. Finally, the health multidimensions and the mobility of the population are considered to test the resulting robustness. Table 4 reports the model results of Equations (5) and (6). As shown in Table 4, Models 1 to 4 present the effects of per capita industrial sulfur dioxide emissions, per capita industrial soot and dust emissions, per capita industrial wastewater emissions, and PM 2.5 concentrations on chronic diseases in middle-aged and older adults, and Models 5 to 8 introduce the interaction terms of per capita industrial sulfur dioxide emissions, per capita industrial soot and dust emissions, per capita industrial wastewater emissions, and PM 2.5 concentrations with healthcare service to investigate the dual effects of urbanization on health in middle-aged and elderly people. The results of Models 1 to 4 show that integrated urbanization significantly reduces the morbidity of chronic diseases in middle-aged and elderly adults at the statistics level of 1%. As expected, modeling results show industrial soot and dust having a significant negative effect on middle-aged and elderly people, but for industrial sulfur dioxide and wastewater, the results show positive effects for chronic disease, which is in disagreement with objective facts. Thus, the interaction terms were introduced in Models 5 to 8. The estimated results from these models indicate that the coefficients of the interaction terms for industrial sulfur dioxide, industrial soot and dust, and healthcare service are significantly negative, showing that the improvement of healthcare conditions can dramatically relieve the health risks to middle-aged and elderly people caused by increased industrial sulfur dioxide and industrial soot and dust. To depict the effects of healthcare on the mitigation of health hazards posed by environmental pollution, this paper drew marginal effects graphs, shown in Figure 1. The second row of the plots in Figure 1 shows different levels of mitigation effects of healthcare on the health risks produced by industrial sulfur dioxide and industrial soot and dust, with the former decreasing the prevalence of chronic diseases from 0.26 to 0.16 and the latter decreasing prevalence of chronic diseases from 0.31 to 0.16. More noticeably, the coefficient of interaction terms between PM 2.5 and healthcare service is markedly positive, indicating that the enhancement of the healthcare service fails to cancel the health problems from PM 2.5 , which increased the prevalence of chronic conditions among the elderly from 0.19 to 0.24 ( Figure 1, first row). It is possible that the smaller particle size of PM 2.5 makes it more likely to penetrate through a body's protective layer than other pollutants and harm public health [42]. However, the interaction term of industrial wastewater and healthcare has a positive effect on chronic diseases in middle-aged and elderly adults, but the effect is not significant. Notes: ***, **, and * represent the significance at 1%, 5%, and 10% levels, respectively; the parentheses represent rob standard errors; the income, industrial wastewater, soot and dust, sulfur dioxide, and PM2.5 values are logarithmic.

The Heterogeneity Analysis of Urbanization on Chronic Diseases in Middle-Aged an Elderly People
According to the baseline analysis described in Section 3.1, comprehensive ur tion can improve the health conditions of middle-aged and elderly people. Howev ing to the multidimensionality [26], spatial variability [43], and imbalance [44] in urbanization, there are also significant regional and urban differences in the im urbanization on health [45]. Section 3.2 reveals this influence to a greater extent by ing the heterogeneity of the dual effects of urbanization on health based on the m mensionality of urbanization, regional differences, and different city types.

The Heterogeneity Analysis of Urbanization on Chronic Diseases in Middle-Aged and Elderly People
According to the baseline analysis described in Section 3.1, comprehensive urbanization can improve the health conditions of middle-aged and elderly people. However, owing to the multidimensionality [26], spatial variability [43], and imbalance [44] in China's urbanization, there are also significant regional and urban differences in the impact of urbanization on health [45]. Section 3.2 reveals this influence to a greater extent by exploring the heterogeneity of the dual effects of urbanization on health based on the multidimensionality of urbanization, regional differences, and different city types.  Notes: ***, **, and * represent the significance at 1%, 5%, and 10% levels, respectively; the parentheses represent robust standard errors; the income, industrial wastewater, soot and dust, sulfur dioxide, and PM 2.5 values are logarithmic.

The Influence of the Development Trend of Urbanization on Chronic Diseases
Each city has a different direction in the process of urbanization. For example, some cities may focus more on economic development and population clustering in the primary stage of urbanization, while others concentrate more on the improvement of the living environment and living conditions in the mature stage of urbanization. These differences lead to variability in the impact of diverse dimensions of urbanization on health. To address variability, this paper adopted the fully arrayed polygonal graphical indication method to classify the comprehensive urbanization into four dimensions to express the features of the urbanization process: the urbanization of the population, the urbanization of economy, residential surroundings urbanization, and residential condition urbanization. To clearly depict the impact of urbanization on health in different dimensions, this paper not only carried out the corresponding regression (Table 5) but also reported the marginal effect of medical service and environmental pollutants on health (Table A1). The results indicate that demographic urbanization, economy urbanization, and residential surroundings urbanization all significantly reduce the incidence of chronic diseases among elderly people (as shown in Table 5); however, no matter which dimension of urbanization is focused on, improved healthcare fails to offset the health risks caused by increased PM 2.5 from the process of urbanization.
When economy urbanization is the focus dimension, the marginal effects of the interaction term between healthcare and PM 2.5 are the smallest (DU: 0.152; EU: 0.124; RSU: 0.144; RCU: 0.161) (Table A1); that is, when economy urbanization is the main orientation of the urbanization process, the increase in the prevalence of chronic diseases among the middle-aged and elderly adults is the lowest, at 12.4% (Table A1). In addition, in this circumstance, the marginal effects of the interaction between medical care and industrial soot and dust are the largest, which indicates that healthcare has the largest mitigation effect on the health damage caused by industrial soot and dust when economic urbanization is the focus dimension, reducing the incidence of chronic diseases among middle-aged and elderly people by 18.4% (Table A1). However, when the main development direction of urbanization is the urbanization of the residential surroundings, healthcare has the greatest reducing effects on the health dangers induced by industrial sulfur dioxide (Table 5), with the incidence of chronic diseases decreased by 14.9% among the middle-aged and elderly adults (as shown by the marginal effect value of 0.149 in Table A1). For industrial wastewater, this mitigation effect is most significant only when the main development direction is the urbanization of living conditions, and the significance level is 1% (Table A1).

The Impact of City Type on Chronic Diseases
The imbalance of China's urbanization levels leads to disequilibrium of the economic development level and the industrial structures, resulting in differing health impacts. This paper uses Nelson's classification method from the United States [46] to classify 114 cities into three different types, based on data from 2011, 2013, and 2015. A city for which the proportion of the secondary industry in GDP is above the average level plus one standard deviation is classified as an industrial city, one for which the proportion of the tertiary industry to GDP is higher than the average level plus one standard deviation is classified as a commercial city, and the rest are classified as a mixed-economy city. This paper uses three groups of regression models to estimate results for the three types of cities, as reported in Table 6. First, although improved healthcare fails to cancel the health risks caused by PM 2.5 , the coefficient and marginal effect of the interaction term between healthcare and PM 2.5 are smallest in a commercial city (commercial city: OR = 2.380, 95% CI: 0.969-5.848; mixed-economy city: OR = 2.250, 95% CI: 1.193-4.242). In other words, at 13.2%, a commercial city has the lowest increase in the morbidity of chronic diseases among middle-aged and elderly people (Table A2). Second, the mitigation effects of healthcare for the health problems induced by industrial sulfur dioxide and industrial soot and dust are greatest in a mixed-economy city (industrial soot and dust in an industrial city: OR = 0.174, 95% CI: 0.074-0.409; in a commercial city: OR = 0.786, 95% CI: 0.453-1.366; in a mixed-economy city: OR = 0.176, 95% CI: 0.123-0.251), which showed decreases in chronic diseases among the middle-aged and elderly people by 27.3% and 16.4%, respectively (Table A2). For the industrial wastewater index, the results reported in Tables 6 and A2 show that the mitigation effect of healthcare on the health hazards brought about by industrial wastewater is insignificant. A possible explanation for this is that the treatment rate of industrial wastewater is constantly increasing with the rapid urbanization process. It is reported that, by 2016, the treatment rate of industrial wastewater in China had reached approximately 95%, so the improvement of health care conditions has little influence on the health problems caused by industrial wastewater.  Notes: ***, **, and * represent the significance at 1%, 5%, and 10% levels, respectively; parentheses represent robust standard errors; the control variables are the same as in Table 1, but the results for the control variables are not reported due to space limitations.

The Regional Differences in the Impact of Urbanization on Chronic Diseases
Finally, this paper considers the spatial discrepancy of the urbanization level in China by separating the sample into three subsamples: Eastern, Central, and Western. In general, urbanization is highest in eastern China, followed by central and western China [47]. This classification method is frequently adopted and is useful for analyzing the regional differences in the impact of urbanization on health. The results shown in Table 7 indicate that, in three regions, consolidated urbanization significantly increased the morbidity of chronic diseases among middle-aged and elderly people, but the extent of influence is most obvious in western China (CU: OR = 23.288, 95% CI: 3.360-161.400). Additionally, although the medical care does not offset the health risks caused by PM 2.5 , the health risks are lowest in eastern China (Eastern: OR = 1.039, 95% CI: 0.423-2.552; Central: OR = 6.367, 95% CI: 1.540-26.320; Western: OR = 1.210, 95% CI: 0.626-2.341), which showed an increase of 0.5% in the incidence of chronic diseases among middle-aged and elderly people (Table A3). In contrast to PM 2.5 , healthcare is the most effective in health hazards caused by industrial soot and dust in eastern China (Eastern: OR = 0.463, 95% CI: 0.300-0.717; Central: OR = 0.316, 95% CI: 0.195-0.512; Western: OR = 0.197, 95% CI: 0.114-0.341), where the morbidity of chronic diseases among middle-aged and elderly adults decreased by 10.6% (Table A3) (Table A3). More interestingly, the enhancement of medical care most effectively reduced the health dangers caused by industrial wastewater in western China, where the prevalence of chronic diseases among the elderly population was reduced by 51.6% (Table A3).

The Robustness Test
To maintain the robustness of the results, it is important to consider the robustness of the results from various aspects. Therefore, this paper focuses on both the substitution dependent variable and the endogeneity; on the one hand, we account for the multidimensionality of health in middle-aged and elderly people, and on the other hand, personal avoidance behaviors are considered. The part of the robustness test mainly explores these two aspects. The first endogenous question is the residential health measure. It uses the objective health measure self-rated health as the dependent variable and considers the impact of the interaction term between medical care service and environmental pollution on this measure among elderly people. The results show in Table 8 estimated results that are generally consistent with the baseline model; that is, the improvement of the healthcare conditions can significantly reduce the negative effects of environmental pollution (sulfur dioxide, soot and dust) on self-rated health among middle-aged and elderly people. The only exception is that for, the environmental pollution index PM 2.5 , the impact of medical care and air pollution on the residential self-rated health is different from the baseline model. This could be due to sampling error or statistical error, or it is possible that respondents may incorrectly assess their health.
The second endogenous issue considered is individual avoidance behavior. Specifically, residents may move to areas with lower air pollution in order to avoid individual health damage. In order to consider population mobility, this paper further restricted the sample by excluding respondents who had lived outside of their permanent residence for more than 6 months. Compared to the baseline regression, the sign and significance of coefficients are roughly insignificant, but the magnitude of the interaction term coefficients is higher (Table 9), which suggests that the mitigation effects of healthcare service levels on the health risks induced by air pollution are likely to be underestimated without accounting for population mobility.   Notes: ***and * represent the significance at 1% and 10% levels, respectively; the parentheses are robust standard errors; the control variables are the same as those in Table 1, but only the results for the control variables are reported due to space limitations. Notes: *** represent the significance at 1% levels; the values in parentheses are robust standard errors; the control variables are the same as those in Table 1, but only the results for the control variables are reported due to space limitations.

Conclusions
This paper employs the fully arrayed polygonal graphical indication method to construct comprehensive urbanization indicators (based on the four dimensions of demographic urbanization, economy urbanization, residential surroundings urbanization, and residential conditions urbanization) and a medical treatment index. Then, we incorporate urbanization indicators, environmental pollution indicators, healthcare, and the interaction terms of environmental pollution and healthcare into the ordered-logistics model to explore the magnitude of the impacts of urbanization on health, including the negative impact of environmental pollution and the positive impact of improved medical care. We also study the heterogeneity effects based on multidimensional urbanization, regional differences, and different city types. Based on this analysis, we conclude the following: (1) The comprehensive urbanization index was constructed by applying the fully arrange polygon graphical index method from four dimensions: the urbanization of the population, the urbanization of economy, the urbanization of residential environment, and the urbanization of residential conditions. The results indicate that integrated urbanization can significantly decrease the rate of chronic diseases among the middleaged and elderly. However, although the enhancement of healthcare conditions can significantly reduce the prevalence of chronic diseases induced by industrial sulfur dioxide and industrial soot and dust, it cannot offset the health damage caused by PM 2.5 . (2) Each city has different development directions in the process of urbanization, which can be considered by the four dimensions of demographic urbanization, economy urbanization, residential surroundings urbanization, and residential conditions urbanization. When the main development direction is economic urbanization, although the medical care still does not fully counteract the health risks from PM 2.5 , the health damage is the lowest, with an increase in the morbidity of chronic diseases among middle-aged and elderly adults of only 12.4%. In addition, the mitigation effect of healthcare for health hazards caused by industrial soot and dust is the highest, with the largest decrease in the incidence of chronic diseases at 18.4% in the primary stage of urbanization. When the development direction is residential surroundings urbanization, the healthcare service has the greatest mitigating effect on the health problems caused by industrial sulfur and dioxide, reducing the incidence of chronic diseases in middle-aged and elderly people by 14.9%. With residential condition urbanization, healthcare can significantly decrease the incidence of chronic diseases caused by industrial wastewater, and the significance level is 1%. (3) Referring to Nelson's classification in the United States, the study cities are separated into three categories: industrial, commercial, and mixed-economy cities. The results suggest that although medical treatment cannot counteract the health risks induced by PM 2.5 , at 13.2%, the commercial cities have the lowest overall increase in the incidence of chronic diseases among middle-aged and elderly people. The greatest reduction in the incidence of chronic diseases in middle-aged and elderly people induced by industrial sulfur dioxide and industrial soot and dust is seen in mixed-economy cities, with incidences decreased by 27.3% and 16.4%, respectively. (4) In order to consider the regional imbalance of urban development in China, the sample is divided into three sections: eastern, central, and western China. The regression results suggest that in eastern China, although healthcare does not offset the health damage caused by PM 2.5 , the overall health risks are lowest, with the incidence of chronic diseases increase by 0.5%. The mitigation of health hazards caused by industrial sulfur and dioxide is the lowest in eastern China, with the prevalence of chronic diseases decreased by 7.4%. However, the reduction influence of medical treatment on health risks caused by industrial soot and dust is the largest in the central part of China, which shows an 18.6% decrease in the prevalence of chronic diseases. Healthcare can significantly decrease the incidence of chronic diseases induced by industrial wastewater.
In light of these conclusions, this paper argues that local governments should focus on controlling air pollutants (particularly PM 2.5 ) and industrial wastewater pollution. First, with regard to PM 2.5 , local governments can compensate for the health hazards of PM 2.5 by enhancing the economic development level and narrowing the development levels between regions. Furthermore, the government should take into account improving the living environment and economic level, balancing the development of secondary and tertiary industries, and strengthening the configuration of medical facilities to tackle the industrial sulfur dioxide and industrial soot and dust pollution. Finally, when treating industrial wastewater pollution, the local government can take some measures to improve the living environment and living conditions, balance the development of industries, and improve other measures to focus on the treatment of eastern and central China and cities dominated by the secondary and tertiary industries.
Although this study explores the dual impact that urbanization has on health through improved healthcare and increased environmental pollution by incorporating healthcare services, environmental pollution, urbanization, and health into the same framework, the data studied are cross-sectional and do not account for changes in dynamics. In future studies, other types of data may be considered to more precisely examine the relationship between the dynamics of urbanization and health.