Spatiotemporal heterogeneity and its determinants of COVID-19 transmission in typical labor export provinces of China

Determinants of health poverty vulnerability in rural areas of Western China in the post-poverty relief era: an analysis based on the Anderson behavioral model

BACKGROUND

Although China has eliminated absolute poverty, the effects of sickness still pose a threat to the prospect of returning to poverty in western rural areas. However, poverty governance extends beyond solving absolute poverty, and should enhance the family's ability to resist risks, proactively identify the existence of risks, and facilitate preventive measures to reduce the probability of falling into poverty again. This study aimed to assess the health poverty vulnerability of rural households in western China and decompose its determinants.

METHODS

Based on survey data from 2022, the three-stage feasible generalized least squares method was used to calculate the health poverty vulnerability index. Then, Anderson's health behavior theory model was extended to analyse various influencing factors using binary logistic regression, and the contribution of each influencing factor was decomposed using the Shapley index. Finally, Tobit regression and the censored least absolute deviations estimation (clad) method were used to test the model's robustness.

RESULTS

A total of 5455 families in the rural Ningxia region of western China were included in the study. The health poverty vulnerability index of the sample population in 2022 was 0.3000 ± 0.2223, and families with vulnerability ≥0.5 accounted for 16.9% of the sample population. From the Anderson behavioral model, the three models including propensity, enabling, and demand factors had the best fit, and the AIC and BIC values were the smallest. The Shapley decomposition showed that the dimensions of the propensity factor, number of residents, age and educational level of the household head, and dependency ratio were the most important factors influencing vulnerability to health poverty. Tobit regression and the clad method proved the reliability of the constructed model through a robustness test.

CONCLUSION

Rural areas still face the risk of becoming poor or falling into poverty owing to residents' health problems. Health poverty alleviation should gradually change from a focus on treatment to prevention, and formulate a set of accurate and efficient intervention policies from a forward-looking perspective to consolidate the results of health poverty alleviation and prevent widescale poverty return.

Spatiotemporal associations between hand, foot and mouth disease and meteorological factors over multiple climate zones

Prior studies of hand, foot, and mouth disease (HFMD) have often observed inconsistent results regarding meteorological factors. We propose the hypothesis that these meteorological associations vary in regions because of the heterogeneity of their geographical characteristics. We have tested this hypothesis by applying a geographical detector and Bayesian space-time hierarchy model to measure stratified spatiotemporal heterogeneity and local associations between meteorological factors and HFMD risk in five climate zones in China from January 2016 to December 2017. We found a significant spatial stratified heterogeneity in HFMD risk and climate zone explained 15% of the spatial stratified heterogeneity. Meanwhile, there was a significant temporal stratified heterogeneity of 14% as determined by meteorological factors. Average temperatures and relative humidity had a significant positive effect on HFMD in all climate zones, they were the most obvious in the southern temperate zone. In northern temperate, southern temperate, northern subtropics, middle subtropics and southern subtropics climate zone, a 1 °C rise in temperature was related to an increase of 3.99%, 13.76%, 4.38%, 3.99%, and 7.74% in HFMD, and a 1% increment in relative humidity was associated with a 1.51%, 5.40%, 2.21%, 3.44%, and 4.78% increase, respectively. These findings provide strong support for our hypotheses that HFMD incidence has a significant spatiotemporal stratified heterogeneity and different climate zones have distinct influences on the disease. These findings provide strong support for our hypotheses: HFMD incidence had significant spatiotemporal stratified heterogeneity and different climate zones had distinct influences on it. The study suggested that HFMD prevention and policy should be made according to meteorological variation in each climate zone.

Spatial analysis of overweight prevalence in China: exploring the association with air pollution

BACKGROUND

Overweight is a known risk factor for various chronic diseases and poses a significant threat to middle-aged and elderly adults. Previous studies have reported a strong association between overweight and air pollution. However, the spatial relationship between the two remains unclear due to the confounding effects of spatial heterogeneity.

METHODS

We gathered height and weight data from the 2015 China Health and Retirement Long-term Survey (CHARLS), comprising 16,171 middle-aged and elderly individuals. We also collected regional air pollution data. We then analyzed the spatial pattern of overweight prevalence using Moran's I and Getis-Ord Gi* statistics. To quantify the explanatory power of distinct air pollutants for spatial differences in overweight prevalence across Southern and Northern China, as well as across different age groups, we utilized Geodetector's q-statistic.

RESULTS

The average prevalence of overweight among middle-aged and elderly individuals in each city was 67.27% and 57.39%, respectively. In general, the q-statistic in southern China was higher than that in northern China. In the north, the prevalence was significantly higher at 54.86% compared to the prevalence of 38.75% in the south. SO2 exhibited a relatively higher q-statistic in middle-aged individuals in both the north and south, while for the elderly in the south, NO2 was the most crucial factor (q = 0.24, p < 0.01). Moreover, fine particulate matter (PM2.5 and PM10) also demonstrated an important effect on overweight. Furthermore, we found that the pairwise interaction between various risk factors improved the explanatory power of the prevalence of overweight, with different effects for different age groups and regions. In northern China, the strongest interaction was found between NO2 and SO2 (q = 0.55) for middle-aged individuals and PM2.5 and SO2 (q = 0.27) for the elderly. Conversely, in southern China, middle-aged individuals demonstrated the strongest interaction between SO2 and PM10 (q = 0.60), while the elderly showed the highest interaction between NO2 and O3 (q = 0.42).

CONCLUSION

Significant spatial heterogeneity was observed in the effects of air pollution on overweight. Specifically, air pollution in southern China was found to have a greater impact on overweight than that in northern China. And, the impact of air pollution on middle-aged individuals was more pronounced than on the elderly, with distinct pollutants demonstrating significant variation in their impact. Moreover, we found that SO2 had a greater impact on overweight prevalence among middle-aged individuals, while NO2 had a greater impact on the elderly. Additionally, we identified significant statistically interactions between O3 and other pollutants.

Integrative Effects between a Bubble and Seal Program and Workers' Compliance to Health Advice on Successful COVID-19 Transmission Control in a Factory in Southern Thailand

Applying health measures to prevent COVID-19 transmission caused disruption of businesses. A practical plan to balance public health and business sustainability during the pandemic was needed. Herein, we describe a "Bubble and Seal" (B&S) program implemented in a frozen seafood factory in southern Thailand. We enrolled 1539 workers who lived in the factory dormitories. First, the workers who had a high fatality risk were triaged by RT-PCR tests, quarantined and treated if they had COVID-19. Newly diagnosed or suspected COVID-19 workers underwent the same practices. The non-quarantined workers were regulated to work and live in their groups without contact across the groups. Workers' personal hygiene and preventive measures were strongly stressed. Between the 6th and 9th weeks of the program, the post-COVID-19 infection status (PCIS) of all participants was evaluated by mass COVID-19 antibody or RT-PCR tests. Finally, 91.8% of the workers showed positive PCIS, which was above the number required for program exit. Although no workers had received a vaccination, there was only one case of severe COVID-19 pneumonia, and no evidence of COVID-19 spreading to the surrounding communities. Implementation of the B&S program and workers' adherence to health advice was the key to this success.

A joint Bayesian spatiotemporal risk prediction model of COVID-19 incidence, IC admission, and death with application to Sweden

The three closely related COVID-19 outcomes of incidence, intensive care (IC) admission and death, are commonly modelled separately leading to biased estimation of the parameters and relatively poor forecasts. This paper presents a joint spatiotemporal model of the three outcomes based on weekly data that is used for risk prediction and identification of hotspots. The paper applies a pure spatiotemporal model consisting of structured and unstructured spatial and temporal effects and their interaction capturing the effects of the unobserved covariates. The pure spatiotemporal model limits the data requirements to the three outcomes and the population at risk per spatiotemporal unit. The empirical study for the 21 Swedish regions for the period 1 January 2020-4 May 2021 confirms that the joint model predictions outperform the separate model predictions. The fifteen-week-ahead spatiotemporal forecasts (5 May-11 August 2021) show a significant decline in the relative risk of COVID-19 incidence, IC admission, death and number of hotspots.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00168-022-01191-1.

Short-term effects of indoor and outdoor air pollution on the lung cancer morbidity in Henan Province, Central China

Lung cancer is one of the most common cancer types and a major cause of death. The relationship between lung cancer morbidity and exposure to air pollutants is of particular concern. However, the relationship and difference in lung cancer morbidity between indoor and outdoor air pollution effects remain unclear. In this paper, the aim was to comprehensively investigate the spatial relationships between the lung cancer morbidity and indoor-outdoor air pollution in Henan based on the standard deviation ellipse, spatial autocorrelation analysis and GeoDetector. The results indicated that (1) the spatial distribution of lung cancer morbidity was related to the geomorphology, while high-morbidity areas were concentrated in the plains and basins of Central, Eastern and Southern Henan. (2) Among the selected outdoor air pollutants, PM_2.5, NO₂, SO₂, O₃ and CO were significantly correlated with the lung cancer morbidity. The degree of indoor air pollution was measured by the use of heating energy, and the proportions of coal-heating households, households with coal/biomass stoves and households with heated kangs were highly decisive in regard to the lung cancer morbidity. (3) The interaction between two factors was more notable than a single factor in explaining the lung cancer morbidity. Moreover, the interaction type was mainly nonlinear enhancement, and the proportion of households with coal/biomass stoves imposed the strongest interaction effect on the other factors.

Methods Used in the Spatial and Spatiotemporal Analysis of COVID-19 Epidemiology: A Systematic Review

The spread of the COVID-19 pandemic was spatially heterogeneous around the world; the transmission of the disease is driven by complex spatial and temporal variations in socioenvironmental factors. Spatial tools are useful in supporting COVID-19 control programs. A substantive review of the merits of the methodological approaches used to understand the spatial epidemiology of the disease is hardly undertaken. In this study, we reviewed the methodological approaches used to identify the spatial and spatiotemporal variations of COVID-19 and the socioeconomic, demographic and climatic drivers of such variations. We conducted a systematic literature search of spatial studies of COVID-19 published in English from Embase, Scopus, Medline, and Web of Science databases from 1 January 2019 to 7 September 2021. Methodological quality assessments were also performed using the Joanna Briggs Institute (JBI) risk of bias tool. A total of 154 studies met the inclusion criteria that used frequentist (85%) and Bayesian (15%) modelling approaches to identify spatial clusters and the associated risk factors. Bayesian models in the studies incorporated various spatial, temporal and spatiotemporal effects into the modelling schemes. This review highlighted the need for more local-level advanced Bayesian spatiotemporal modelling through the multi-level framework for COVID-19 prevention and control strategies.

Identifying spatiotemporal patterns of COVID-19 transmissions and the drivers of the patterns in Toronto: a Bayesian hierarchical spatiotemporal modelling

Spatiotemporal patterns and trends of COVID-19 at a local spatial scale using Bayesian approaches are hardly observed in literature. Also, studies rarely use satellite-derived long time-series data on the environment to predict COVID-19 risk at a spatial scale. In this study, we modelled the COVID-19 pandemic risk using a Bayesian hierarchical spatiotemporal model that incorporates satellite-derived remote sensing data on land surface temperature (LST) from January 2020 to October 2021 (89 weeks) and several socioeconomic covariates of the 140 neighbourhoods in Toronto. The spatial patterns of risk were heterogeneous in space with multiple high-risk neighbourhoods in Western and Southern Toronto. Higher risk was observed during Spring 2021. The spatiotemporal risk patterns identified 60% of neighbourhoods had a stable, 37% had an increasing, and 2% had a decreasing trend over the study period. LST was positively, and higher education was negatively associated with the COVID-19 incidence. We believe the use of Bayesian spatial modelling and the remote sensing technologies in this study provided a strong versatility and strengthened our analysis in identifying the spatial risk of COVID-19. The findings would help in prevention planning, and the framework of this study may be replicated in other highly transmissible infectious diseases.

The Geographical Distribution and Influencing Factors of COVID-19 in China

The study of the spatial differentiation of COVID-19 in cities and its driving mechanism is helpful to reveal the spatial distribution pattern, transmission mechanism and diffusion model, and evolution mechanism of the epidemic and can lay the foundation for constructing the spatial dynamics model of the epidemic and provide theoretical basis for the policy design, spatial planning and implementation of epidemic prevention and control and social governance. Geodetector (Origin version, Beijing, China) is a great tool for analysis of spatial differentiation and its influencing factors, and it provides decision support for differentiated policy design and its implementation in executing the city-specific policies. Using factor detection and interaction analysis of Geodetector, 15 indicators of economic, social, ecological, and environmental dimensions were integrated, and 143 cities were selected for the empirical research in China. The research shows that, first of all, risks of both infection and death show positive spatial autocorrelation, but the geographical distribution of local spatial autocorrelation differs significantly between the two. Secondly, the inequalities in urban economic, social, and residential environments interact with COVID-19 spatial heterogeneity, with stronger explanatory power especially when multidimensional inequalities are superimposed. Thirdly, the spatial distribution and spread of COVID-19 are highly spatially heterogeneous and correlated due to the complex influence of multiple factors, with factors such as Area of Urban Construction Land, GDP, Industrial Smoke and Dust Emission, and Expenditure having the strongest influence, the factors such as Area of Green, Number of Hospital Beds and Parks, and Industrial NOx Emissions having unignorable influence, while the factors such as Number of Free Parks and Industrial Enterprises, Per-GDP, and Population Density play an indirect role mainly by means of interaction. Fourthly, the factor interaction effect from the infected person’s perspective mainly shows a nonlinear enhancement effect, that is, the joint influence of the two factors is greater than the sum of their direct influences; but from the perspective of the dead, it mainly shows a two-factor enhancement effect, that is, the joint influence of the two factors is greater than the maximum of their direct influences but less than their sum. Fifthly, some suggestions are put forward from the perspectives of building a healthy, resilient, safe, and smart city, providing valuable reference and decision basis for city governments to carry out differentiated policy design.

Heterogeneity of the COVID-19 Pandemic in the United States of America: A Geo-Epidemiological Perspective

The spread of the COVID-19 pandemic has shown great heterogeneity between regions of countries, e. g., in the United States of America (USA). With the growing of the worldwide COVID-19 pandemic, there is a need to better highlight the variability in the trajectory of this disease in different worldwide geographic areas. Indeed, the epidemic trends across areas can display completely different evolution at a given time. Geo-epidemiological analyses using data, that are publicly available, could be a major topic to help governments and public administrations to implement health policies. Geo-epidemiological analyses could provide a basis for the implementation of relevant public health policies. With the COVID-19 pandemic, geo-epidemiological analyses can be readily utilized by policy interventions and USA public health authorities to highlight geographic areas of particular concern and enhance the allocation of resources.

Tracking and Controlling the Spatiotemporal Spread of SARS-CoV-2 Lineage B.1.1.7 in COVID-19 Reopenings

Understanding why or how the emergence of SARS-CoV-2 variants has occurred and how to control them is crucial as regards the potential of global reopening. To explore and further understand the spatiotemporal dynamics of the B.1.1.7 spread in the 368 districts of Taiwan, a district-level geographic prediction model of the risk of COVID-19 symptom onset has been proposed. It has been found that, (a) the human mobility, epidemic alert measures, and vaccination rates all played an important role in the spatiotemporal heterogeneity of B.1.1.7 transmission; (b) for regions with high human mobility and low vaccination rates, the partial relaxation of entry quarantine measures for specific imported groups would, in fact, lead to a wide spread of B.1.1.7 with a consequent doubling of high-onset-risk areas and together with the overall onset risk, a further increase of more than 20% would occur; (c) compared with the closing of business places and public venues in all districts, both lockdown in those areas of high-onset-risk and the gathered control effects regarding other districts, the control of B.1.1.7 spread would be better enabled by an onset risk reduction of up to 91.36%. Additionally, an increase in the vaccination rate in each district by up to 5-10 times would further reduce the onset risk by 6.07%-62.22%.

Spatial heterogeneity of bacillary dysentery and the impact of temperature in the Beijing-Tianjin-Hebei region of China

Previous studies indicate that the incidence of bacillary dysentery is closely related to meteorological factors. However, the impact of temperature and the spatial heterogeneity of the disease in regions of unbalanced socioeconomic development remains unclear. Therefore, this research collected data for 29,639 daily bacillary dysentery cases in children under 5 years of age, as well as the meteorological variables from China's Beijing-Tianjin-Hebei region, to analyze the spatial pattern of bacillary dysentery and reveal its nonlinear association with temperature. The SatScan method was employed first, to detect the spatial heterogeneity of the disease risk, and then the distributed lag nonlinear model (DLNM) was used to analyze the relationships between the daily minimum, mean, and maximum temperatures and bacillary dysentery in the stratified heterogeneous regions. The results indicated that bacillary dysentery incidence presented statistically significant spatial heterogeneity. The area of highest risk was found to be Beijing and its neighboring regions, which have high population densities. There was also a positive association between bacillary dysentery and temperature. Hotter temperatures were accompanied by higher relative risks. In the most likely spatial cluster region, the excess risk (ER) values for a 1°C rise in minimum, mean, and maximum temperatures above the median were 4.65%, 11.30%, and 19.21%, respectively. The effect of temperature on bacillary dysentery peaked at a lag of 3 to 4 days. The findings of this study will aid risk assessments and early warning systems for bacillary dysentery.

Study of SARS-CoV-2 transmission in urban environment by questionnaire and modeling for sustainable risk control

Caused by SARS-CoV-2, COVID-19 has become a severe threaten to society and human health, its epidemic control emerges as long-term issue. A sustainable epidemic and environmental transmission risk control (SEERC) in urban area is urgently needed. This work aims to conduct a new investigation on the transmission risk of SARS-COV-2 as virus/hazardous material through various environmental medias, routes and regions in the entirely urban area for guiding the SEERC. Specifically, 5 routes in 28 regions (totally 140 scenarios) are considered. For a new perspective, the risk evaluation is conducted by the quantification of frontline medicals staffs' valuable experience in this work. 207 specialists responsible for the treatment of over 9000 infected patients are involved. The result showed that degree of risk was in the order of breath>contact-to-object>contact-to-human>intake>unknown. The modeling suggested source control as the prior measure for epidemic control. The combination of source control & mask wearing showed high efficiency in SEERC. The homeworking policy needed to cooperate with activity limitation to perform its efficiency. Subsequently, a new plan for SEERC was discussed. This work delivered significant information to researchers and decision makers for the further development of sustainable control for SARS-COV-2 spreading and COVID-19 epidemic.

Spatio-temporal variation in tuberculosis incidence and risk factors for the disease in a region of unbalanced socio-economic development

BACKGROUND

Previous research pointed to a close relationship between the incidence of tuberculosis (TB) in aging populations and socio-economic conditions, however there has been lack of studies focused on a region of unbalanced socio-economic development. The aim of this paper is to explore the spatio-temporal variation in TB incidence and examine risk determinants of the disease among aging populations in a typical region.

METHODS

Data on TB-registered cases between 2009 and 2014, in addition to social-economic factors, were collected for each district/county in Beijing, Tianjin and Hebei, a region characterized by an aging population and disparities in social-economic development. A Bayesian space-time hierarchy model (BSTHM) was used to reveal spatio-temporal variation in the incidence of TB among the elderly in this region between 2009 to 2014. GeoDetector was applied to measure the determinant power (q statistic) of risk factors for TB among the elderly.

RESULTS

The incidence of TB among the elderly exhibited geographical spatial heterogeneity, with a higher incidence in underdeveloped rural areas compared with that in urban areas. Hotspots of TB incidence risk among the elderly were mostly located in north-eastern and southern areas in the study region, far from metropolitan areas. Areas with low risk were distributed mainly in the Beijing-Tianjin metropolitan areas. Social-economic factors had a non-linear influence on elderly TB incidence, with the dominant factors among rural populations being income (q = 0.20) and medical conditions (q = 0.17). These factors had a non-linear interactive effect on the incidence of TB among the elderly, with medical conditions and the level of economic development having the strongest effect (q = 0.54).

CONCLUSIONS

The findings explain spatio-temporal variation in TB incidence and risk determinants of elderly TB in the presence of disparities in social-economic development. High-risk zones were located mainly in rural areas, far from metropolitan centres. Medical conditions and the economic development level were significantly associated with elderly TB incidence, and these factors had a non-linear interactive effect on elderly TB incidence. The findings can help to optimize the allocation of health resources and to control TB transmission in the aging population in this region.

Emergence of universality in the transmission dynamics of COVID-19

The complexities involved in modelling the transmission dynamics of COVID-19 has been a roadblock in achieving predictability in the spread and containment of the disease. In addition to understanding the modes of transmission, the effectiveness of the mitigation methods also needs to be built into any effective model for making such predictions. We show that such complexities can be circumvented by appealing to scaling principles which lead to the emergence of universality in the transmission dynamics of the disease. The ensuing data collapse renders the transmission dynamics largely independent of geopolitical variations, the effectiveness of various mitigation strategies, population demographics, etc. We propose a simple two-parameter model-the Blue Sky model-and show that one class of transmission dynamics can be explained by a solution that lives at the edge of a blue sky bifurcation. In addition, the data collapse leads to an enhanced degree of predictability in the disease spread for several geographical scales which can also be realized in a model-independent manner as we show using a deep neural network. The methodology adopted in this work can potentially be applied to the transmission of other infectious diseases and new universality classes may be found. The predictability in transmission dynamics and the simplicity of our methodology can help in building policies for exit strategies and mitigation methods during a pandemic.

Comparison of Public Health Containment Measures of COVID-19 in China and India

Objective

This study aimed to make a comparative analysis of the public health containment measures between China and India, explore the causes of the serious COVID-19 epidemic in India, and eventually to improve global infectious disease control.

Methods

We extracted publicly available data from official websites, summarized the containment measures implemented in China and India, and assessed their effectiveness.

Results

China has responded to the COVID-19 outbreak with strict public health containment measures, including lockdown of Wuhan city, active case tracing, and large-scale testing, ultimately preventing a large increase in daily new cases and maintaining a low mortality rate per million population (as of May 5, 2021, daily new cases were 11 and mortality rate per million population was 3.37). India, although imposing a national lockdown to control the pandemic, has not implemented strict testing, tracking, and quarantine measures due to the overburdened healthcare system. Combined with massive lockdown, it has accelerated human mobility and exacerbated the epidemic, resulting in a rapid increase in daily new cases and a high mortality rate per million population (as of May 5, 2021, daily new cases were 412,431 and mortality rate per million population was 166.79).

Conclusion

China and India implemented public health containment measures to contain the spread of the COVID-19 pandemic based on their national situations. Meanwhile, daily new cases and mortality of COVID-19 also were affected by environmental and socioeconomic. Countries make a comprehensive strategy not only in terms of the biological, pharmaceutical, health, and sanitation sectors but also based on sustainability science and environmental science.

Spatial pattern of the population casualty rate caused by super typhoon Lekima and quantification of the interactive effects of potential impact factors

BACKGROUND

Typhoons greatly threaten human life and property, especially in China. Therefore, it is important to make effective policy decisions to minimize losses associated with typhoons.

METHODS

In this study, the GeoDetector method was used to quantify the determinant powers of natural and socioeconomic factors, and their interactions, on the population casualty rate of super typhoon Lekima. The local indicator of spatial association (LISA) method was followed to explore the spatial pattern of the population casualty rate under the influence of the identified dominant factors.

RESULTS

Both natural and socioeconomic factors were found to have significantly impacted the population casualty rate due to super typhoon Lekima. Among the selected factors, maximum precipitation was dominant factor (q = 0.56), followed by maximum wind speed (q = 0.45). In addition, number of health technicians (q = 0.35) and number of health beds (q = 0.27) have a strong influence on the population casualty rate. Among the interactive effects of 12 influencing factors, the combined effects of maximum precipitation and ratio of brick-wood houses, the maximum precipitation and ratio of steel-concrete houses, maximum precipitation and number of health technicians were highest (q = 0.72). Furthermore, high-risk areas with very high casualty rates were concentrated in the southeastern part of Zhejiang and northern Shandong Provinces, while lower-risk areas were mainly distributed in northern Liaoning and eastern Jiangsu provinces.

CONCLUSIONS

These results contribute to the development of more specific policies aimed at safety and successful property protection according to the regional differences during typhoons.

Determinant Powers of Socioeconomic Factors and Their Interactive Impacts on Particulate Matter Pollution in North China

Severe air pollution has significantly impacted climate and human health worldwide. In this study, global and local Moran’s I was used to examine the spatial autocorrelation of PM_2.5 pollution in North China from 2000–2017, using data obtained from Atmospheric Composition Analysis Group of Dalhousie University. The determinant powers and their interactive effects of socioeconomic factors on this pollutant are then quantified using a non-linear model, GeoDetector. Our experiments show that between 2000 and 2017, PM_2.5 pollution globally increased and exhibited a significant positive global and local autocorrelation. The greatest factor affecting PM_2.5 pollution was population density. Population density, road density, and urbanization showed a tendency to first increase and then decrease, while the number of industries and industrial output revealed a tendency to increase continuously. From a long-term perspective, the interactive effects of road density and industrial output, road density, and the number of industries were amongst the highest. These findings can be used to develop the effective policy to reduce PM_2.5 pollution, such as, due to the significant spatial autocorrelation between regions, the government should pay attention to the importance of regional joint management of PM_2.5 pollution.