What Factors Dominate the Change of PM2.5 in the World from 2000 to 2019? A Study from Multi-Source Data

As the threat to human life and health from fine particulate matter (PM2.5) increases globally, the life and health problems caused by environmental pollution are also of increasing concern. Understanding past trends in PM2.5 and exploring the drivers of PM2.5 are important tools for addressing the life-threatening health problems caused by PM2.5. In this study, we calculated the change in annual average global PM2.5 concentrations from 2000 to 2020 using the Theil–Sen median trend analysis method and reveal spatial and temporal trends in PM2.5 concentrations over twenty-one years. The qualitative and quantitative effects of different drivers on PM2.5 concentrations in 2020 were explored from natural and socioeconomic perspectives using a multi-scale geographically weighted regression model. The results show that there is significant spatial heterogeneity in trends in PM2.5 concentration, with significant decreases in PM2.5 concentrations mainly in developed regions, such as the United States, Canada, Japan and the European Union countries, and conversely, significant increases in PM2.5 in developing regions, such as Africa, the Middle East and India. In addition, in regions with more advanced science and technology and urban management, PM2.5 concentrations are more evenly influenced by various factors, with a more negative influence. In contrast, regions at the rapid development stage usually continue their economic development at the cost of the environment, and under a high intensity of human activity. Increased temperature is known as the most important factor for the increase in PM2.5 concentration, while an increase in NDVI can play an important role in the reduction in PM2.5 concentration. This suggests that countries can achieve good air quality goals by setting a reasonable development path.


Introduction
In the context of global change issues, air quality has received increasing attention from researchers around the world. Fine particulate matter (PM 2.5 ) is particulate matter with an aerodynamic equivalent diameter less than or equal to 2.5 microns in ambient air, which has the characteristics of a long residence time in the atmosphere and long transmission distance, and has a great impact on human health and atmospheric environmental quality [1,2]. Many epidemiological studies have shown that many human health problems related to respiratory and lung diseases are associated with the concentration of PM [3,4]. In terms of the Global Burden of Disease (GBD), PM 2.5 contributes to millions of premature deaths around the world [5,6]. Despite many countries and international organizations making many efforts to solve it, the threat of PM 2.5 to the global environment and human health remains [7].
Numerous studies have been conducted regarding near-ground retrieval of PM 2.5 concentration, and many achievements have been made. Currently, as the demand for PM 2.5 research increases, the requirements for the high resolution and accuracy of near-surface MODIS on Terra and Aqua can capture information on vegetation, and the resulting global vegetation index product is provided by the 16-Day MODIS-Terra vegetation index from C6 (MOD13A2). The algorithm for this product chooses the best available pixel value from all the acquisitions from the 16 days. MOD13A2 has a spatial resolution of 0.01 × 0.01 • and a time horizon of 2000 to the present. MOD13A2 includes two primary vegetation index layers, NDVI and EVI. In this study, we calculated the global annual NDVI using the Google Earth Engine (GEE) (https://developers.google.com/earth-engine/datasets/. Accessed on 16 January 2023). Specifically, the data for 2000 is an incomplete dataset for the year; this dataset starts in February, thus, the global average NDVI in 2000 was calculated from February to December.

Meteorological Datasets
Changes in PM 2.5 are sensitive to meteorological conditions, and meteorological data play a significant role in changes in PM 2.5 concentrations [25]. The datasets of monthly averaged data on single levels from 1959 to the present produced by ERA5 were used in this study (https://cds.climate.copernicus.eu/. Accessed on 16 January 2023). Reanalysis combines model data with observations from across the world using the related laws of physics to complete a globally consistent dataset. These reanalysis datasets provide more than 260 meteorological variables, where the atmospheric variable spatial resolution is 0.25 × 0.25 • , and ocean wave variable spatial resolution is 0.5 × 0.5 • . Finally, total precipitation, 10 m u-component of wind, 10 m v-component of wind, boundary layer height and temperature at a height of 2 m were chosen for inclusion in the research.

Nighttime Light Dataset
Nighttime light (NTL) satellite data can visualize the intensity of human activities at night. Because of its strong correlation with urban structure and socioeconomic characteristics, an increasing number of researchers are using NTL remote sensing data to study the economic development of cities [26,27]. The NTL dataset is commonly used by the Defence Meteorological Satellite Program Operational Linescan System (DMSP-OLS) and Suomi National Polar-Orbiting Partnership-Visible Infrared Imaging Radiometer Suite (NPP-VIIRS). The two NTL datasets differ in time span and quality assessment criteria, making it difficult to obtain perfect data for long time-series analysis. There are also some studies that extend NTL data for long time-series by integrating two datasets [28]. In this study, we use the 2000-2020 annual average NTL dataset (https://dataverse.harvard.edu/dataset.xhtml?persistentId=doi: 10.7910/DVN/YGIVCD. Accessed on 16 January 2023) integrated by the autoencoder (AE) model using convolutional neural networks. The dataset has a spatial resolution of 0.01 × 0.01 • and shows good agreement at both the pixel scale and city scale, with remote production validation with R 2 values of 0.87 and 0.95, respectively [29].

Global Statistical Data
Population density, urban ratio and GDP per capita data can reflect socioeconomic attributes and regional development levels. The World Bank provides a crucial tool to support key management decisions and key statistical information regarding the Bank's operational activities in the context of the increasing demand for high-quality statistics. The World Bank database contains a large number of indicators covering statistics on various aspects of urban development, education, environment, economy, etc. Most of the data are sourced from its member countries, and therefore, the quality of the data is determined by the individual country systems. The population density, urban ratio and GDP per capita data used in this study are obtained from the World Bank database (https://data.worldbank.org.cn/. Accessed on 16 January 2023).

Theil-Sen Median Trend Analysis and the Mann-Kendall Test
Theil-Sen median trend analysis is a robust trend calculation method which uses nonparametric statistics. The method is a linear regression through the median of the slopes of all lines at paired points. It is computationally efficient and insensitive to measurement errors and outlier data. In addition, it is often used for trend analysis of long time-series data [30], so it can objectively reflect spatiotemporal changes in PM 2.5 concentrations.
where K x represents the slope of the fitted equation and x j and x i represent the value of the PM 2.5 concentrations in years j and i, respectively. When K x > 0, the PM 2.5 concentration shows an increasing trend. In contrast, when K x < 0, it shows a decreasing trend. The Mann-Kendall (MK) test is a nonparametric statistical test used to determine the significance of trends in time series data [31,32]. The statistic Z of the trend test is calculated as follows: S and Var(S) in equal 2 were calculated as follows: Var(S) = n(n − 1)(2n + 5) 18 (5) where S is a statistical variable, sgn x j − x i is the sign function, and Var(S) is the variance of S. Statistical Z is a bilateral trend test that can be combined with the normal distribution table to determine the significance level. In the study, statistic Z and the significance level correspond to Table 1. According to the first law of geography, everything is spatially correlated with each other; specifically, the closer things are to each other, the greater the spatial correlation. Therefore, when regression analysis is performed on spatial data, its spatial location relationship should be considered. In order to explore the spatial nonstationarity of spatial data, the GWR model was proposed [33]. Distinct from global regression, the results calculated by the GWR model can be used to analyze the different relationships between the dependent and explanatory variables of the study object at different spatial scopes. Its essence is to obtain spatially continuous varying regression coefficients by weighting distances. The classical GWR considers spatial heterogeneity, but it uses a constant bandwidth in whole study area which can cause poor fitting results for some variables in local range.
The MGWR relaxes the assumption that all processes to be modeled are on the same spatial scale and overcomes the limitations of GWR very well. It captures the most appropriate bandwidth for each variable, so that the results of the regression obtained are more reliable [34]. The MGWR can be described as follows: where Y i is a dependent variable representing the PM 2.5 concentration value at the spatial represents the intercept distance of the regression model, β bwk (u i , v i ) and X ik represent the coefficient and value of the explanatory variable k, respectively, and ε i is the random error. The core of the MGWR model is the spatial weight function, and the selection of this function is crucial to the correct estimation of the regression parameters. In this paper, a Gaussian function is used as the weight function, which can be described as follows: where d ij represents the distance between i and j, and b represents the bandwidth used to describe the relationship between distance and weight. According to the MGWR model, the relationship between PM 2.5 and explanatory variables can be established to better explain the spatial distribution characteristics of PM 2.5 .

Global Spatiotemporal Patterns of PM 2.5 in 2000-2019
The global spatiotemporal trends in PM 2.5 concentration over the past twenty years were calculated, and the significance was tested (Figure 1). In most parts of the world, PM 2.5 concentration showed a decreasing trend. The United States, eastern Canada, Japan, western Russia, northeast Argentina, the Tibetan Plateau region of China and a great number of countries in Europe show an extremely significant decreasing trend (p < 0.01), of which most are developed countries. Regions of western and eastern Africa (e.g., Namibia, Zambia, Nigeria, etc.), Turkmenistan, China's coastal region and the Korean Peninsula are either generally significant or slightly significant (p < 0.05 and p < 0.1). In contrast, some countries or regions show increasing PM 2.5 concentrations. Most regions of the Middle East area and India showed extremely significant increases (p < 0.01). Countries or regions such as central Canada, eastern Russia, Australia, Chile, and southern Argentina show different degrees of increase in the concentration of PM 2.5 (p < 0.05 and p < 0.1). It could be seen that the PM 2.5 concentration in developed countries has been declining more significantly than in other regions over the past two decades.
insula are either generally significant or slightly significant (p < 0.05 and p < 0.1). In trast, some countries or regions show increasing PM2.5 concentrations. Most regions o Middle East area and India showed extremely significant increases (p < 0.01). Countri regions such as central Canada, eastern Russia, Australia, Chile, and southern Argen show different degrees of increase in the concentration of PM2.5 (p < 0.05 and p < 0. could be seen that the PM2.5 concentration in developed countries has been declining m significantly than in other regions over the past two decades. During the twenty years, overall global PM2.5 concentration showed a nonmono trend ( Figure 2). It showed an increased from 2000 to 2003 then highly fluctuated ove next 12 years; during this time, the global economy was in a rapid development p However, in the past five years, it has begun to show a remarkable decreasing trend decreasing trend may be because the world has come to realize that environmental p tion should also be taken seriously as the economy develops. The decrease in PM2.5 centration is significantly influenced by stringent policies adopted [35][36][37]. The fitting of the time series shows that the trend line has decreased slowly, with a drop rate of − μg/m 3 /yr over the past 20 years. The highest annual average PM2.5 concentration appe in 2008, with concentrations above 20.2 μg/m 3 . The WHO issued air quality guidelin 2021, which set the annual average PM2.5 concentration at 5 μg/m 3 , and noted that w the concentration is more than 35 μg/m 3 , human health will be seriously threatened During the twenty years, overall global PM 2.5 concentration showed a nonmonotonic trend ( Figure 2). It showed an increased from 2000 to 2003 then highly fluctuated over the next 12 years; during this time, the global economy was in a rapid development phase. However, in the past five years, it has begun to show a remarkable decreasing trend; this decreasing trend may be because the world has come to realize that environmental pollution should also be taken seriously as the economy develops. The decrease in PM 2.5 concentration is significantly influenced by stringent policies adopted [35][36][37]. The fitting line of the time series shows that the trend line has decreased slowly, with a drop rate of −0.055 µg/m 3 /yr over the past 20 years. The highest annual average PM 2.5 concentration appeared in 2008, with concentrations above 20.2 µg/m 3 . The WHO issued air quality guidelines in 2021, which set the annual average PM 2.5 concentration at 5 µg/m 3 , and noted that when the concentration is more than 35 µg/m 3 , human health will be seriously threatened. At present, it seems that PM 2.5 concentrations are moving closer to the WHO standard, but there is still a long way to go.

Spatial-Temporal Patterns of Natural Driving Factors
NDVI is often used to reflect the state of terrestrial vegetation. There are a large number of studies showing that green spaces such as urban forests and parks effectively re-

Spatial-Temporal Patterns of Natural Driving Factors
NDVI is often used to reflect the state of terrestrial vegetation. There are a large number of studies showing that green spaces such as urban forests and parks effectively remove particle pollutants [38,39]. From 2000 to 2019 (Figure 3), the overall global NDVI showed an increasing trend. Before 2008, there was a slow decline, with an annual average decrease in NDVI of 0.01 over the 8 years, after which the trend of NDVI increased. The folding line of the global annual average value of NDVI was similar to a U-type; NDVI and PM 2.5 showed opposite trends over the 20 years. Correspondingly, the highest value of PM 2.5 and the lowest value of NDVI occurred in the same year. Regarding spatial distribution, most regions showed a trend towards becoming greener ( Figure 4). The NDVI shows an extremely significant increase (p < 0.01) in central and eastern China, northwestern India, Turkey, Greece, southern Brazil, and southern Chile. The United States, Mexico and many countries in central Africa are also experiencing significant increases (p < 0.05). In contrast, many regions in the Middle East, South Africa, Nigeria, Canada and Australia observed various significant decreasing trends. It is worth noting that the relationship between NDVI and PM 2.5 concentrations in some regions is opposite to that of the overall relationship. However, NDVI does not distinguish between agricultural land and other green spaces; at the same time, deforestation and agricultural expansion are occurring in many regions [40], such as India, where and the increasing trend is mainly contributed to by cropland [41]. For Canada, Guinea, Nigeria and Benin, although the NDVI showed a downwards trend, the PM 2.5 concentrations in these regions also showed a downwards trend. This suggests that the variation in PM 2.5 concentrations in each region may not be well explained from the perspective of trend changes in NDVI alone, and that such results may be dominated by changes in other factors.     Wet deposition is one of the most efficient ways to remove pollutants from the mosphere [42], and precipitation is the most common type of wet deposition metho During this twenty-year period, more than half of the areas did not show a significa increase or decrease in precipitation. However, in other areas, there was a large spatia stratified heterogeneity in the trends of precipitation across the globe ( Figure 5). Preci tation in the northern and eastern regions of the U.S. and Alaska showed increasing tren of varying significance, and the corresponding regional PM2.5 concentrations were on t decline, suggesting an important local correlation between precipitation and PM2.5 co centration. A similar situation occurred in Canada, Finland, Vietnam, and parts of weste Africa. However, many areas, such as India, Saudi Arabia and Indonesia, showed trend of increasing precipitation, but also an increasing trend in PM2.5 concentration; t increase in precipitation did not seem to have a visual impact on the decrease in PM Wet deposition is one of the most efficient ways to remove pollutants from the atmosphere [42], and precipitation is the most common type of wet deposition method. During this twenty-year period, more than half of the areas did not show a significant increase or decrease in precipitation. However, in other areas, there was a large spatially stratified heterogeneity in the trends of precipitation across the globe ( Figure 5). Precipitation in the northern and eastern regions of the U.S. and Alaska showed increasing trends of varying significance, and the corresponding regional PM 2.5 concentrations were on the decline, suggesting an important local correlation between precipitation and PM 2.5 concentration. A similar situation occurred in Canada, Finland, Vietnam, and parts of western Africa. However, many areas, such as India, Saudi Arabia and Indonesia, showed a trend of increasing precipitation, but also an increasing trend in PM 2.5 concentration; the increase in precipitation did not seem to have a visual impact on the decrease in PM 2.5 concentration from Figures 1 and 5. The main areas of reduced precipitation were Australia, eastern China, eastern Russia, northwestern Kazakhstan, central Africa, South Africa, and eastern and southern South America. The changes in PM 2.5 concentration in these areas with reduced precipitation also show different trends. The time-series line chart shows a slowly fluctuating increase in global annual mean precipitation with a growth rate of 0.0025 mm/yr ( Figure 6). There may be an important link between the increase in precipitation and global warming, but the increasing trend in precipitation may allow wet deposition to play a better role.  Figure 6). There may be an important link between the increase precipitation and global warming, but the increasing trend in precipitation may allow w deposition to play a better role.    Global warming is becoming an increasingly prominent issue, and has significant implications for the Earth's natural environment. The generation of the greenhouse effect has led to changes in the characteristics of meteorological factors, resulting in changes in natural conditions such as regional precipitation and temperature, and it also has had an Global warming is becoming an increasingly prominent issue, and has significant implications for the Earth's natural environment. The generation of the greenhouse effect has led to changes in the characteristics of meteorological factors, resulting in changes in natural conditions such as regional precipitation and temperature, and it also has had an impact on the growth of terrestrial vegetation [43,44]. Near-surface temperature is an important component of climate conditions, and there is an inextricable relationship between PM 2.5 concentration and temperature [45]. From 2000 to 2019, almost all regions showed a significant increase in temperature on a global scale, and even if there are very few regions, such as parts of Canada, the United States and Pakistan, with a decreasing trend, the overall temperature increase cannot be reversed (Figures 7 and 8). In the 20 years, the temperature rose 0.88 K. The effect of this on PM 2.5 is complicated because an increase in temperature has different effects on the different components of the air [46]. Therefore there is variability in the relationship between the warming effect and the change in PM 2.5 concentration in different regions. impact on the growth of terrestrial vegetation [43,44]. Near-surface temperature is a portant component of climate conditions, and there is an inextricable relationship bet PM2.5 concentration and temperature [45]. From 2000 to 2019, almost all regions sh a significant increase in temperature on a global scale, and even if there are very fe gions, such as parts of Canada, the United States and Pakistan, with a decreasing t the overall temperature increase cannot be reversed (Figures 7 and 8). In the 20 year temperature rose 0.88 K. The effect of this on PM2.5 is complicated because an incre temperature has different effects on the different components of the air [46]. Ther there is variability in the relationship between the warming effect and the change in concentration in different regions.

Spatial-Temporal Patterns of Socioeconomic Driving Factors
Over the past two decades, the world economy has experienced a period of rapid development. Global GDP per capita showed a steady upwards trend, with a decline only

Spatial-Temporal Patterns of Socioeconomic Driving Factors
Over the past two decades, the world economy has experienced a period of rapid development. Global GDP per capita showed a steady upwards trend, with a decline only in 2001, 2008 and 2015 ( Figure 9). In 2001, over-saturated markets due to rapid technological development caused an economic crisis; in 2008, the United States triggered the subprime mortgage crisis, resulting in a hard hit to the world economy; in 2015, world industrial production experienced low growth, trade continued to slump, and financial market turmoil intensified. These three separate years of global events ultimately led to the economic downturn. Economic development and environmental protection are two inseparable topics. If the protection of the environment is neglected during economic development, the unrestricted emission of various pollution sources may lead to the increasing concentration of PM 2.5 , resulting in significant risks to human health.  Figure 9). In 2001, over-saturated markets due to rapid technological development caused an economic crisis; in 2008, the United States triggered the subprime mortgage crisis, resulting in a hard hit to the world economy; in 2015, world industrial production experienced low growth, trade continued to slump, and financial market turmoil intensified. These three separate years of global events ultimately led to the economic downturn. Economic development and environmental protection are two inseparable topics. If the protection of the environment is neglected during economic development, the unrestricted emission of various pollution sources may lead to the increasing concentration of PM2.5, resulting in significant risks to human health. NTL can be a good measure of the socioeconomic development of a region. Most developed countries, such as the United States, South Korea, and the European Union region, have shown an extremely significant increasing trend of NTL during the past 20 NTL can be a good measure of the socioeconomic development of a region. Most developed countries, such as the United States, South Korea, and the European Union region, have shown an extremely significant increasing trend of NTL during the past 20 years, but the trend of PM 2.5 concentration in these regions has not increased with the increase in human activity intensity. Instead, these areas show a decreasing trend, which is most likely because developed countries have invested more money in environmental management and protection and have achieved certain management results. In addition, different countries and regions have different development routes and priorities, resulting in economic development while failing to consider environmental pollution control. For this reason there is an increasing trend of PM 2.5 concentration, mainly occurring in India, the Middle East, Chile, Colombia, and other developing regions. Moreover, in Canada, Japan and other regions, NTL shows a decreasing trend, and the PM 2.5 concentration in these regions also shows the same trend; perhaps the reduction in the intensity of human activity has also effectively controlled the reduction in anthropogenic emissions ( Figure 10). Changes in NTL and GDP are closely related to population density and urban The urban ratio is defined as the ratio of the urban population in the region to the population in the region. As the urban ratio continues to increase, there is a large of rural population into urban areas, which eventually leads to a trend of increasing population density and decreasing population density in the periphery. The chan population density has an important impact on the intensity of human activities region, and it will also have a prominent effect on the PM2.5 concentration. Throug ures 11 and 12, the trends of population density and urban ratio show a high deg consistency over the 20-year period on the global scale. Within the first decade, t crease in population density and urban ratio may have contributed significantly increase in PM2.5 concentration, whereas within the second decade, PM2.5 concent did not change at the same time. The decline in PM2.5 concentration indicates that t fect of increasing population density and urban expansion on the change in PM2.5 co tration is gradually weakening. Changes in NTL and GDP are closely related to population density and urban ratio. The urban ratio is defined as the ratio of the urban population in the region to the total population in the region. As the urban ratio continues to increase, there is a large influx of rural population into urban areas, which eventually leads to a trend of increasing urban population density and decreasing population density in the periphery. The change in population density has an important impact on the intensity of human activities in the region, and it will also have a prominent effect on the PM 2.5 concentration. Through Figures 11 and 12, the trends of population density and urban ratio show a high degree of consistency over the 20-year period on the global scale. Within the first decade, the increase in population density and urban ratio may have contributed significantly to the increase in PM 2.5 concentration, whereas within the second decade, PM 2.5 concentration did not change at the same time. The decline in PM 2.5 concentration indicates that the effect of increasing population density and urban expansion on the change in PM 2.5 concentration is gradually weakening. ures 11 and 12, the trends of population density and urban ratio show a high degree of consistency over the 20-year period on the global scale. Within the first decade, the increase in population density and urban ratio may have contributed significantly to the increase in PM2.5 concentration, whereas within the second decade, PM2.5 concentration did not change at the same time. The decline in PM2.5 concentration indicates that the effect of increasing population density and urban expansion on the change in PM2.5 concentration is gradually weakening.

Comparison with Traditional Linear Models
In this study, a MGWR model was applied to better understand the influence of natural and socio-economic factors on the spatial and temporal distribution of PM2.5 at national scales worldwide. Apart from the MGWR model, we compared different traditional linear models, such as ordinary least squares (OLS) and classic GWR. Table 2 shows the results of the three models. The traditional linear regression model ignores the spatial heterogeneity and the accuracy of the calculation was the worst. Variance inflation factor (VIF) was also calculated and it was found that the VIF values of the selected factors were all less than three. It is considered that there was no covariance between each factor, therefore GWR and MGWR can be used for regression. The GWR model takes into account the spatial heterogeneity of the study area, and the accuracy of the regression was also improved, however, because GWR defaults to all explanatory variables having the same spatial scale in the calculation process, this leads to less credible regression results in local areas. So, our study finally adopted the MGWR model, which also significantly improved the accuracy of the regression results, in addition to the interpretability of the model.

Comparison with Traditional Linear Models
In this study, a MGWR model was applied to better understand the influence of natural and socio-economic factors on the spatial and temporal distribution of PM 2.5 at national scales worldwide. Apart from the MGWR model, we compared different traditional linear models, such as ordinary least squares (OLS) and classic GWR. Table 2 shows the results of the three models. The traditional linear regression model ignores the spatial heterogeneity and the accuracy of the calculation was the worst. Variance inflation factor (VIF) was also calculated and it was found that the VIF values of the selected factors were all less than three. It is considered that there was no covariance between each factor, therefore GWR and MGWR can be used for regression. The GWR model takes into account the spatial heterogeneity of the study area, and the accuracy of the regression was also improved, however, because GWR defaults to all explanatory variables having the same spatial scale in the calculation process, this leads to less credible regression results in local areas. So, our study finally adopted the MGWR model, which also significantly improved the accuracy of the regression results, in addition to the interpretability of the model.  It can be assumed that with the increase in urbanization rate, the spatial layout within each country is more rational, which makes the production and living of people and the centralized management of the environment more efficient and orderly [47,48]. According to the urbanization rate statistics provided by the World Bank, countries in North and South America have relatively high urbanization rates globally, and they are already in a relatively mature state of urban development, maintaining a stable level of impact on PM 2.5 compared to other factors. For countries in the Asian region, the negative effect of urbanization rate on PM 2.5 concentrations mostly diminished over time, probably as a result of the increased influence of other factors on PM 2.5 concentrations during the development process. In Europe, the opposite trend to Asia was shown, with urbanization rate playing an increasingly important role in the change of PM 2.5 concentrations. Africa is relatively less developed than other regions and has relatively uneven development between countries, leading to greater heterogeneity. Population density, which reflects the density of the population living in a certain area, has a very important impact on PM 2.5 concentration [49]. Generally speaking, the higher the population density, the higher the intensity of human activity. Similarly, dense urban and industrial areas are often accompanied by high population densities, which will also lead to increased emission of air pollutants. We calculated that population density had a negative effect on PM 2.5 concentrations in most countries in 2000, except for central and southern Africa and the Middle East. However, in 2010, the negative effect of this factor on PM 2.5 concentrations weakened, and reversed in 2019. The year-on-year increase in population density is illustrated in Figure 11, and the high population density further leads to more anthropogenic emissions, which can provide a significant influence on the increase in PM 2.5 concentration [50]. In 2000, the effect of GDP per capita on PM 2.5 concentration had significant geographical differences, with negative effects in all countries in the European and American region and positive effects in countries in the Asian and African regions. Such positive effects were more pronounced in developing countries, while some economically developed countries in the United States, Canada, and Europe showed negative effects. However, in 2010, the effect of GDP per capita on PM 2.5 in all countries of the world showed a negative effect, with regression coefficients ranging from −0.025 to 0, which also indicates that its effect is relatively small. Years before 2010, as the global average annual PM 2.5 concentration was at a high level, various problems caused by the increase in environmental pollution has made countries around the world also began to pay attention to the monitoring and management of PM 2.5 , and air pollution prevention and control bills introduced by various countries one after another may be an important reason for this shift [51,52]. In 2019, the negative impact on PM 2.5 concentrations increased in most countries around the world as investment in environmental protection increased, but a positive impact has developed again in the Middle East and in eastern Africa, which may also contribute to the 20-year trend of increasing PM 2.5 concentrations in these regions.

Conclusions and Limitation
The Theil-Sen median trend and MGWR approach for analyzing the spatiotemporal distribution of PM2.5 concentrations and the contributions of driving factors have improved our understanding of PM2.5. Currently, PM2.5 concentrations show a decreasing trend in the more economically developed countries, while many relatively backwards developing countries show an increasing trend, especially in India and the Middle East. From the perspective of natural and socioeconomic factors, we can see that developed and developing countries show opposite results under the same driving factors and developed and developing countries are at distinct stages of urban development. Developed countries have largely completed the urbanization process, while developing countries are experiencing rapid urbanization. This has a lot to do with the higher economic, technological and management levels of developed countries, which have experienced the lesson of "pollute first, treat later" and have invested more in environmental protection, and transferred many heavy polluting enterprises to developing countries, resulting in two opposite trends. The trends of global average NDVI and global PM2.5 concentrations show "Utype" and inverse "U-type" distributions, respectively, suggesting that the increase in green areas has the power to explain reduced PM2.5 concentrations. In conclusion, the increase in population density, socioeconomic development and urban expansion in the context of global warming and increasing climate extremes tend to lead to an increase in PM2.5 concentration. However, through effective management and by setting a reasonable development path, for example, in the process of urban planning, reasonably arranging the spatial layout of different functional industries to facilitate centralized management of industries that emit air pollutants, and then increasing investment in high-tech development to accelerate the goal of conversion to clean energy, the goal of synergistic development of economic development and environmental protection can be achieved in the end.
In this study, a limited number of drivers were selected to analyze their relationship with PM2.5 concentration. Emission sources that can have a direct impact on PM2.5 concentrations were not considered, making the discussion of the detailed causes of PM2.5 concentration changes somewhat limited. In addition, although the MGWR model is better than the OLS model and the GWR model in terms of fitting accuracy, there are still some limitations. The MGWR model is also a type of linear model, and a nonlinear expression may yield better results under complex spatial distributions. Moreover, the explanation of a single factor is not completely reliable; the process of PM2.5 concentration change The influence of natural factors on PM 2.5 concentration cannot be ignored. In general, moisture in the air is adsorbed and collected by suspended PM 2.5 , which eventually falls to the ground in the form of precipitation, while the occurrence of precipitation also leads to an increase in moisture in the air, which has a good effect on the removal of PM 2.5 from the atmosphere [53]. In 2000, our results also show spatial distribution characteristics that are consistent with the above conditions. Precipitation has a negative effect on PM 2.5 concentrations in regions with abundant mean annual precipitation, and this effect is particularly evident in regions with sufficient annual precipitation such as the Americas, Australia, southern Asia, and western Africa. In contrast, regions with low precipitation, such as southern and eastern Africa and northern Asia, show different results. Although seasonal precipitation is high in some of these regions, this study is based on annual averages at the national scale, weakening its spatial heterogeneity and certain characteristics. In addition, the intensity and frequency of precipitation can affect PM 2.5 concentration variability [54]. In 2010 and 2019, the regression coefficients for precipitation were all between −0.05 and 0, indicating that the effect of precipitation on PM 2.5 concentration appears to be weakened relative to other factors. The process of temperature influence on PM 2.5 concentration is complex and a positive correlation between them is usually considered, which is consistent with our calculations at the global scale and shows a trend towards a stronger positive effect in many countries, which may have some connection with the global warming trend. However, some studies have also shown that higher temperatures enhance air convection, which in turn dilutes and diffuses PM 2.5 and reduces local PM 2.5 concentrations [55,56]. His diffusion effect may not be significant at larger spatial scales. During the past two decades, the positive effect of temperature on PM 2.5 concentration has become increasingly significant due to the effects of global warming. Among the many factors, the regression coefficient of temperature is in the higher range, which means that the effect of temperature on PM 2.5 concentration is crucial. Vegetation cover in all countries of the world showed different degrees of negative effects on PM 2.5 concentration, which is consistent with the results of most studies, and it can be inferred that good vegetation cover can effectively reduce atmospheric PM 2.5 concentration and is a very effective factor in suppressing PM 2.5 concentration [57]. At these three time points, the variation in the regression coefficients of NDVI is relatively small worldwide, but shows a clear gradient, with the largest negative effects in Asia and Australia, followed by Africa, while in the Americas and Europe such negative effects are relatively weak. Boundary layer height determines the volume available for diffusion and transport of pollutants in the atmosphere [58]; when the BLH is large, the PM 2.5 concentration mixed in it will be relatively low. The results of our study show that BLH has a negative effect on PM 2.5 globally, and the spatial distribution is stable over the past 20 years, with the value of the regression coefficient remaining between −0.35 and −0.1. The degree of negative influence also increases gradually from west to east. Wind is an important condition in the process of diffusion and transport of pollutants in the atmosphere [59]. Wind speed and wind direction together affect the diffusion rate and direction of PM 2.5 concentration. In 2000 and 2010, the positive effect of wind speed on PM 2.5 concentrations was the largest in Africa and smallest in Asia, and by 2019, the effect of wind speed shifted to negative globally, probably due to the acceleration between ocean currents and sea breezes in the context of global warming, which further enhanced the diffusion of PM 2.5 in the atmosphere and made PM 2.5 concentration decrease.
Both socio-economic factors and natural factors lead to change in global PM 2.5 concentration, and due to inconsistent development stages and speed of countries around the world, there is a strong spatial heterogeneity. According to the calculation results for developed countries, Europe and the Americas have more advanced science and technology and more reasonable urban structural layout, so the regression coefficients of these factors are relatively balanced. Conversely, economically underdeveloped countries still have to continue their economic development at the cost of the environment in order to seek better development, and so PM 2.5 is more significantly affected with a positive influence. Throughout all three years, the positive effect of temperature on PM 2.5 concentration was the largest of all variables. Although temperature is one of the natural factors, combustion and emissions from human activities are decisive causes of the greenhouse effect. The negative effect of NDVI on PM 2.5 concentration is the largest in all three years. NDVI is influenced not only by the natural growth of vegetation, but also by human activities. The increase in artificial green areas can effectively reduce PM 2.5 concentration. This also proves that the optimization of urban structure layout, the development of high technology, and the increase in vegetation all have positive effects on the reduction in PM 2.5 concentration. The decrease in global average annual PM 2.5 concentration in the past decade with the joint efforts of countries around the world confirms that the formulation of reasonable air pollution prevention and control acts can effectively monitor and manage PM 2.5 .

Conclusions and Limitation
The Theil-Sen median trend and MGWR approach for analyzing the spatiotemporal distribution of PM 2.5 concentrations and the contributions of driving factors have improved our understanding of PM 2.5 . Currently, PM 2.5 concentrations show a decreasing trend in the more economically developed countries, while many relatively backwards developing countries show an increasing trend, especially in India and the Middle East. From the perspective of natural and socioeconomic factors, we can see that developed and developing countries show opposite results under the same driving factors and developed and developing countries are at distinct stages of urban development. Developed countries have largely completed the urbanization process, while developing countries are experiencing rapid urbanization. This has a lot to do with the higher economic, technological and management levels of developed countries, which have experienced the lesson of "pollute first, treat later" and have invested more in environmental protection, and transferred many heavy polluting enterprises to developing countries, resulting in two opposite trends. The trends of global average NDVI and global PM 2.5 concentrations show "U-type" and inverse "U-type" distributions, respectively, suggesting that the increase in green areas has the power to explain reduced PM 2.5 concentrations. In conclusion, the increase in population density, socioeconomic development and urban expansion in the context of global warming and increasing climate extremes tend to lead to an increase in PM 2.5 concentration. However, through effective management and by setting a reasonable development path, for example, in the process of urban planning, reasonably arranging the spatial layout of different functional industries to facilitate centralized management of industries that emit air pollutants, and then increasing investment in high-tech development to accelerate the goal of conversion to clean energy, the goal of synergistic development of economic development and environmental protection can be achieved in the end.
In this study, a limited number of drivers were selected to analyze their relationship with PM 2.5 concentration. Emission sources that can have a direct impact on PM 2.5 concentrations were not considered, making the discussion of the detailed causes of PM 2.5 concentration changes somewhat limited. In addition, although the MGWR model is better than the OLS model and the GWR model in terms of fitting accuracy, there are still some limitations. The MGWR model is also a type of linear model, and a nonlinear expression may yield better results under complex spatial distributions. Moreover, the explanation of a single factor is not completely reliable; the process of PM 2.5 concentration change should be a very complex process, and the interaction between these influencing factors is not clear.
In future studies, we may consider more comprehensive influencing factors and try more models to explain the relationship between influencing factors and PM 2.5 concentrations, and undertake further exploration of the interaction between factors.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Conflicts of Interest:
The authors declare no conflict of interest.