Spatio-Temporal Heterogeneous Impacts of the Drivers of NO 2 Pollution in Chinese Cities: Based on Satellite Observation Data

: Rapid urbanization in China has led to an increasing problem of atmospheric nitrogen dioxide (NO 2 ) pollution, which negatively impacts urban ecology and public health. Nitrogen dioxide is an important atmospheric pollutant, and quantitative spatio-temporal analysis and inﬂuencing factor analysis of Chinese cities can help improve urban air pollution. In this study, the spatiotemporal analysis methods were used to explore the variations of NO 2 pollution in Chinese cities from 2005 to 2020. The ﬁndings are as follows. In more than half of Chinese cities, NO 2 levels remarkably decreased between 2005 and 2020. The effective NO 2 reduction strategies contributed to the signiﬁcant NO 2 reduction during the 13th Five-Year Plan (2016–2020). Moreover, we found that the pandemic of COVID-19 alleviated NO 2 pollution in China since it reduced the trafﬁc, industrial, and living activities. The NO 2 pollution in Chinese cities was found highly spatially clustered. The geographically and temporally weighted regression model was used to analyze the spatio-temporal heterogeneity of NO 2 pollution inﬂuencing factors in Chinese cities, including natural meteorological and socio-economic factors. The results showed that the GDPPC, population densities, and ambient air pressure were positively correlated with NO 2 pollution. In contrast, the ratio of the tertiary to the secondary industry, temperature, wind speed, and relative humidity negatively impacted the NO 2 pollution level. The ﬁndings of this research contribute to the improvement of urban air quality, stimulating the achievements of the sustainable development goals of Chinese cities.


Introduction
In recent years, the use of tremendous amounts of fossil energy in Chinese cities has resulted in nitrogen oxide pollutants into the atmosphere, severely worsening the air quality. Nitrogen oxides (NO x ) are made up of two compounds: nitrogen monoxide (NO) and nitrogen dioxide (NO 2 ) [1]. Notably, NO 2 has posed a huge threat to the air quality of Chinese cities. As a result, it has become a significant pollutant in urban air quality monitoring systems [2]. NO 2 has had the fastest concentration growth rate in China over the last two decades [3]. Numerous studies have demonstrated that long-term exposure to high NO 2 concentrations increases mortality from respiratory and cardiovascular diseases [4]. Additionally, NO 2 is a significant contributor to acid rain and photochemical smog [3]. Most significantly, NO 2 takes part in the formation of ozone and aerosol, affecting the local climate change [5].
The Chinese government's National Air Pollution Prevention and Control Joint Center has conducted a synergistic treatment of PM 2.5 and ozone [6]. In addition, in 2015, the United Nations established 17 Sustainable Development Goals (SDGs) as a plan of action to achieve global peace and prosperity by 2030 [7]. The mitigation of NO 2 pollution is closely related to the mapping of hazardous chemicals and pollutants in the air, water, and soil in SDG Target 3 "good health and well-being", mapping of air quality in SDG Target 11 "sustainable cities and communities", and the environmental variables for climate change models in SDG Target "climate action" [8,9]. Therefore, NO 2 has received increased attention from researchers and global governments. China, as the biggest developing country, also has suffered from NO 2 pollution for a long time. Hence, it is necessary to analyze the spatio-temporal distribution of NO 2 pollution in China and explore its influencing factors. NO 2 pollution is primarily measured based on the ground platform. However, we discovered that the ground observations over China were unevenly distributed within the short period of historical time series (only available after 2013). On the other hand, satellite remote sensing observation data have the advantages of broad spatial coverage, a long observation period, and spatial continuity [10,11]. Therefore, the tropospheric NO 2 vertical column densities (VCDs) data retrieved by the Ozone Monitoring Instrument (OMI) are widely applied to detect the long-term variations in NO 2 pollution over China and investigate the drivers from a spatio-temporal perspective [1,3,5].
During the last two decades, many researchers have used satellite observation technologies to study the various factors influencing NO 2 pollution over China [12]. Since industrial, transportation, and residential emissions are all anthropogenic sources of NO x emissions [13], fluctuations in NO 2 concentrations are highly correlated with human activities [14]. It has been found that civil vehicles, electricity consumption, total population, built-up areas, and coal use are closely correlated with NO 2 pollution levels [12]. Wang et al. [15] analyzed the spatial and temporal distribution of NO 2 columns over China using the simple linear regression model based on OMI satellite observations. The association between changes in NO 2 pollution and urbanization in China was also conducted [3,16]. Bucsela et al. [17] used satellite measurements to examine the impacts of income and urban spatial form on urban NO 2 levels. Since air pollution is a regional problem, the traditional linear regression method cannot solve the spatial autocorrelation problem of air pollution. However, spatial econometric models can effectively address the complex spatial interactions and spatial dependence factors in regression models. The link between NO 2 pollution and its natural and socio-economic factors in Chinese cities was quantified by incorporating spatial effects in an extended STIRPAT model (stochastic impacts by regression on population, affluence, and technology) [12]. Moreover, socio-economic or public health events that significantly influence human activities can affect local NO 2 pollution levels, such as regional economic recessions [18,19], the 2008 Beijing Olympic Games [20], the 2016 G20 Hangzhou Summit [21], and the COVID pandemic incidents [2,[22][23][24][25].
The studies mentioned above drew fruitful conclusions but neglected the variability of influencing factors of NO 2 pollution under different spatial and temporal conditions. Additionally, NO x emissions in eastern China grew fast between 2000 and 2011 and began to fall steadily [26,27]. Governments recently established a regional coordinated development strategy, in which economically underdeveloped regions in central and western China were encouraged to absorb the energy-intensive industry shifted from the eastern coastal regions [28]. Since the 12th Five-Year Plan (2011-2015), some cities in central China and western China remarkably increased their NO x emissions [3]. It has been evidenced that spatially varying socio-economic conditions and natural geographic factors contribute to varying levels of NO 2 pollution at the city level [12]. In other words, these existing studies disclosed the average impacts of the drivers of NO 2 pollution over China, and they paid little attention to the spatio-temporal variations of the drivers from one city to another, which could mask important information on NO 2 pollution prevention and control based on the characteristics of Chinese cities.
Hence, our study aims to fill this gap using a novel geographically and temporally weighted regression (GTWR) approach proposed by Huang et al. [29] to test the spatiotemporal heterogeneity of the influencing factors of NO 2 pollution in each prefecture-level city. Huang et al. [29] incorporated both spatial and temporal characteristics into a regression model on the basis of the classical geographically weighted regression (GWR) approach and proposed the effective construction method of the spatio-temporal distance. In recent years, the GTWR model has been widely used in many studies to explore the heterogeneous relationships between independent and dependent variables in time and space, such as the relationship between carbon emissions and urbanization [30], air pollutants and natural geographical conditions [31,32], and ecosystem services and human factors [33].
This study first analyzed the spatial and temporal variations of NO 2 pollution over China from 2005 to 2020 using satellite observations data. Then, the GTWR model was applied to explore the influencing factors on the NO 2 pollution at the city level over China. Finally, the relationship between NO 2 pollution and both metrological and socio-economic variables was investigated and quantified. The findings of this study may provide both potential solutions for China's urban air pollution prevention and control and scientific support for achieving sustainable development goals for China.

Data Description
The tropospheric NO 2 VCDs were retrieved by the Ozone Monitoring Instrument onboard the EOS-Aura satellite [34]. The satellite, launched in September 2004, is in a sunsynchronous orbit at 705 km altitude with a 99 min time period. It has the OMI pixel size of 13 × 24 km 2 at nadir in the global mode and 13 × 12 km 2 in the zoom mode, with a local passage time of approximately 13:40 [10]. We used the Royal Netherlands Meteorological Institute's (KNMI) monthly QA4ECV NO 2 long-term dataset (version 1.1) in this study (https://www.temis.nl/airpollution/no2col/no2regioomimonth_qa.php, accessed on 16 March 2022). The QA4ECV tropospheric NO 2 VCDs product has a spatial resolution of 0.125 • × 0.125 • . The updated NO 2 retrieval algorithms were referred to Boersma, Eskes, Dirksen, van der A, Veefkind, Stammes, Huijnen, Kleipool, Sneep, Claas, Leitão, Richter, Zhou, and Brunner [34] and Boersma et al. [35]. The monthly data have an uncertainty of approximately 10% but range from 15% to 30% in contaminated places [35]. We further removed the data with row anomalies and cloud radiance fractions of more than 50%. Then, the monthly gridded composite of tropospheric NO 2 VCDs was averaged into the annual NO 2 VCDs gridded dataset.
The independent variables in the GTWR model are introduced and processed as follows. The ground-based meteorological variables, including temperature (temp), wind speed (WS), ambient air pressure near the ground (Pres), and relative humidity (Humi), were obtained from the National Meteorological Information Center of China Meteorological Administration (http://data.cma.cn/, accessed on 16 March 2022). Given the daily station monitoring data, we applied an inverse distance-weighted (IDW) interpolated method to interpolate into grid data with the spatial resolution of 0.125 • × 0.125 • , keeping the consistency with the spatial resolution of the QA4ECV tropospheric NO 2 VCDs product. The daily dataset was then processed into an annual mean dataset. Meanwhile, the yearly socio-economic variables and indicators of Chinese cities from 2005 to 2019, including foreign direct investment (FDI), population density (PD), gross domestic product per capita (GDPPC), and the ratio of the tertiary to the secondary industry (TSRatio), were obtained from the China Statistical Yearbooks and the China City Statistical Yearbooks. In addition, the tropospheric NO 2 and meteorological parameters at the prefectural level were retrieved according to China's administrative boundary vector data (http://www.resdc.cn/, accessed on 16 March 2022) by using ArcGIS 10.8 software.
It should be noted that although the data for NO 2 VCDs of prefecture-level cities are extracted, due to the data unavailability of the explanatory variables for some cities, the sample size in the regression analysis is restricted to 271 prefecture-level cities. The descriptive statistics of the variables involved in the regression models (standard deviation (S.D.), mean, median, minimum (Min), and maximum (Max)) are summarized in Table 1.

Methodology
We note that NO 2 VCDs in adjacent cities tend to be similar. In other words, NO 2 pollution may exhibit spatial autocorrelation. Hence, global Moran's I is introduced to test if there is spatial dependence for NO 2 VCDs of Chinese cities. It is as follows: and Y are the NO 2 VCDs of city i, j, and are the average values of NO 2 VCDs of all samples. W ij denotes a spatial weights matrix, which describes the spatial arrangement of these cities. n is the number of all samples.
When the global Moran's I is significant and positive, it indicates that NO 2 pollution of one city is similar to that of its neighbors, namely a spatial clustering. When global Moran's I is significant and negative, it shows that NO 2 pollution is spatially dispersed. When global Moran's I equals zero, NO 2 pollution may be randomly distributed in space.
Anselin [36] proposed a local Moran's I, called the local indicator of spatial association (LISA), to test whether similar or dissimilar observations are clustered together in a local area. The local Moran's I value of each city can be calculated as follows: where w ij is also the spatial weights matrix; x i is the attribute of city I; and x is the mean value of the attributes. 2 is the variance of the attributes. In this study, to measure the importance of the explanatory factors affecting NO 2 VCDs, the independent variables can be standardized as follows: where X k is the standardized independent variable k. X k is the original independent variable k. X k_min and X k_max are the minimum and maximum values of the independent variable k, respectively.
We use the GTWR model to quantify the spatio-temporal heterogeneous impacts of different independent variables on NO 2 VCDs changes in Chinese cities. Compared with the cross-sectional data used in the traditional GWR model, the GTWR model incorporates temporal variations into the GWR model. In this study, the ordinary least squares (OLS), GWR, temporally weighted regression (TWR), and GTWR models are all conducted for is the unknown coefficient of influencing factors to be estimated, including meteorological conditions and socio-economic factors; ε i denotes a random error. Therefore, the coefficients β n (u i , v i , t i ) are estimated by the least-squares method. The estimator reads below: , and k is the number of the explanatory variables. The GTWR model is essentially determined by the ability of the kernel function in the weight matrix to solve for both the temporal non-smoothness and spatial non-smoothness. The weight matrix is computed by the Euclidean distance and Gaussian distance-decaybased functions, as described in Wu et al. [37]. The weight scheme is generated by a bandwidth parameter h [38], and the schematic diagram of the spatio-temporal distance d ST ij between cities i and j is shown in Figure 1 and represented as follows: where µ is the scale factor of the temporal and spatial distance; then the GTWR model is built and compared with the other models (OLS, GWR, TWR, and GTWR) by a series of goodness-of-fit statistics, i.e., corrected Akaike Information Criterion (AICc) [29]. We use the GTWR model to quantify the spatio-temporal heterogeneous impacts of different independent variables on NO2 VCDs changes in Chinese cities. Compared with the cross-sectional data used in the traditional GWR model, the GTWR model incorporates temporal variations into the GWR model. In this study, the ordinary least squares (OLS), GWR, temporally weighted regression (TWR), and GTWR models are all conducted for completeness and comparison on the ArcGIS 10.8 software. The GTWR model can be expressed as follows: where , , denotes city i at location , and year ; , , is an intercept, , , is the unknown coefficient of influencing factors to be estimated, including meteorological conditions and socio-economic factors; denotes a random error.
Therefore, the coefficients , , are estimated by the least-squares method. The estimator reads below: where , , , , … , , and k is the number of the explanatory variables. The GTWR model is essentially determined by the ability of the kernel function in the weight matrix to solve for both the temporal non-smoothness and spatial non-smoothness. The weight matrix is computed by the Euclidean distance and Gaussian distancedecay-based functions, as described in Wu et al. [37]. The weight scheme is generated by a bandwidth parameter ℎ [38], and the schematic diagram of the spatio-temporal distance between cities i and j is shown in Figure 1 and represented as follows: (6) where is the scale factor of the temporal and spatial distance; then the GTWR model is built and compared with the other models (OLS, GWR, TWR, and GTWR) by a series of goodness-of-fit statistics, i.e., corrected Akaike Information Criterion (AICc) [29].  Figure 2 presents the annual mean of NO2 VCDs at the prefectural city level from 2005 to 2020. The most severe NO2 pollution is mainly located in northern China during this sample period. Notably, we observe that some northern cities were highly polluted, for example, Jiaozuo (Henan Province), Handan (Hebei Province), Shijiazhuang (Hebei Province), and Xingtai (Hebei Province), in which NO2 VCDs are larger than 2000 × 10 13 molecule/cm 2 . Apart from the North China Plain, moderately polluting cities (larger than  Figure 2 presents the annual mean of NO 2 VCDs at the prefectural city level from 2005 to 2020. The most severe NO 2 pollution is mainly located in northern China during this sample period. Notably, we observe that some northern cities were highly polluted, for example, Jiaozuo (Henan Province), Handan (Hebei Province), Shijiazhuang (Hebei Province), and Xingtai (Hebei Province), in which NO 2 VCDs are larger than 2000 × 10 13 molecule/cm 2 . Apart from the North China Plain, moderately polluting cities (larger than 800 × 10 13 molecule/cm 2 ) were observed in the important urban clusters, such as the Yangtze River Delta, Pearl River Delta, and Sichuan-Chongqing urban agglomeration.

Analysis of Spatio-Temporal Variation of NO 2 Pollution
Remote Sens. 2022, 14, x FOR PEER REVIEW 6 of 16 800 × 10 13 molecule/cm 2 ) were observed in the important urban clusters, such as the Yangtze River Delta, Pearl River Delta, and Sichuan-Chongqing urban agglomeration.    Figure 3c, we noticed that the COVID-19 pandemic resulted in a decrease in NO 2 in 288 (86.2%) Chinese cities. In particular, NO 2 VCDs significantly fell by more than 30 percent in some cities of the North China Plain and Yangtze River Delta regions. This shows that the epidemic has also exerted a significant impact on the traffic, industrial, and living activities in the local cities, leading to the rapid reduction of the NO 2 levels. Furthermore, the cities in the Pearl River Delta region are the exceptions, with NO 2 levels declining from 2005 to 2020 as a result of the Guangdong and Hong Kong governments' joint emissions control efforts, which began in 2003 [16].   Besides, to further evaluate the impact of the COVID-19 pandemic on NO levels in prefecture-level cities, we compared the average values of February and 2020 with the average values of earlier two years, namely, 2018 and 2019, remov influence of the Chinese New Year on NO2 levels. As shown in Figure 3c, we notic the COVID-19 pandemic resulted in a decrease in NO2 in 288 (86.2%) Chinese c particular, NO2 VCDs significantly fell by more than 30 percent in some cities of th China Plain and Yangtze River Delta regions. This shows that the epidemic has erted a significant impact on the traffic, industrial, and living activities in the loca leading to the rapid reduction of the NO2 levels. Furthermore, the cities in the Pea Delta region are the exceptions, with NO2 levels declining from 2005 to 2020 as a r the Guangdong and Hong Kong governments' joint emissions control efforts, wh gan in 2003 [16].
A spatial autocorrelation among the NO2 pollution of these Chinese cities can ified. We then calculated the Moran's I value of annual mean NO2 levels at the ci for robustness check using ArcGIS 10.8 software. As shown in Figure 4, the resu cated that the global Moran's I values during the sample period were statistically cant and larger than zero (p-value < 0.05 at the 95% confidence level), indicating th pollution exhibits a significant positive spatial autocorrelation every year. A spatial autocorrelation among the NO 2 pollution of these Chinese cities can be verified. We then calculated the Moran's I value of annual mean NO 2 levels at the city level for robustness check using ArcGIS 10.8 software. As shown in Figure 4, the results indicated that the global Moran's I values during the sample period were statistically significant and larger than zero (p-value < 0.05 at the 95% confidence level), indicating that NO 2 pollution exhibits a significant positive spatial autocorrelation every year.   Figure 5. It was discovered that the high-high cluster regions in these years were located in the North China Plain and the Yangtze River Delta. The high-high cluster of NO2 pollution appeared in Pearl River Delta only in 2005. The low-low cluster regions of NO2 pollution were mainly located in western China and northeastern China. Addi-  Figure 5. It was discovered that the high-high cluster regions in these years were located in the North China Plain and the Yangtze River Delta. The high-high cluster of NO 2 pollution appeared in Pearl River Delta only in 2005. The low-low cluster regions of NO 2 pollution were mainly located in western China and northeastern China. Additionally, it is noted that the low-low cluster regions narrowed from 2005 to 2020.  Figure 5. It was discovered that the high-high cluster regions in these years were located in the North China Plain and the Yangtze River Delta. The high-high cluster of NO2 pollution appeared in Pearl River Delta only in 2005. The low-low cluster regions of NO2 pollution were mainly located in western China and northeastern China. Additionally, it is noted that the low-low cluster regions narrowed from 2005 to 2020.

Regression Results
We selected eight important socio-economic and natural independent variables for this study: namely, GDP per capita, population density, ratio of the tertiary to the secondary industry, foreign direct investment, temperature, wind, pressure, and humidity. The hypothesized spatio-temporal relationships between NO 2 pollution and the independent variables are demonstrated with an N*T estimated coefficient matrix. Otherwise, such variations may be averaged incorrectly in a global model (e.g., OLS). Hence, the GTWR is able to capture the spatio-temporal heterogeneous impacts of these explanatory variables locally on NO 2 pollution in each city. Table 2 presents the estimation results of four models, namely, OLS, TWR, GWR, and GTWR. We find that the GTWR model has the highest R 2 value (0.904) compared with that of GWR (0.879), TWR (0.804), and OLS (0.776), indicating that the GTWR model has the highest explanatory power. Additionally, the GTWR model has the lowest AICc value (170.898) and the residual sum of squares (218.057) among these four models, also showing that the GTWR model is the best fitted compared to the other models in terms of these goodness-of-fit statistics. In general, the GTWR models might be suitable for analyzing n the spatial and temporal variations of the influencing factors of NO 2 pollution in this study. Given the geographical disparities and temporal trends in NO 2 pollution, natural geographic conditions, and socio-economic factors prevalent in all Chinese prefecture-level cities, there are apparent spatial and temporal variances in the effects of each explanatory variable on NO 2 VCDs. As a result, the GTWR model was developed to examine the impacts of relevant factors on NO 2 pollution in each city at different times. Table 3 summarizes the descriptive statistics for the estimated coefficients of the GTWR model (i.e., mean, median, minimum (Min), and maximum (Max)). Overall, GDP per capita, population density, foreign direct investment, and air pressure are positively correlated with NO 2 pollution. In contrast, the ratio of the tertiary to the secondary industry, temperature, wind speed, and relative humidity are negatively correlated with NO 2 pollution. This means that the increases in GDP per capita, population density, and foreign direct investment exacerbate NO 2 pollution in most cities. A similar finding also applies to the air pressure variable. By contrast, increasing the ratio of the tertiary to the secondary industry or warmer, windier, or more humid climates could reduce NO 2 VCDs. Figures 6 and 7 depict the temporal and spatial variations in the coefficients of the explanatory factors in the estimated results of the GTWR model from 2005 to 2019. As illustrated in Figure 6, most GDP per capita coefficients are positive, with only a few negative coefficients. GDP per capita is generally used to indicate the economic development level. It is linked because it generates a large amount of NO 2 emissions due to the raw economic development pattern typically adopted in the early period. Meanwhile, the GDP per capita coefficients in these cities from southern and northwestern China are the highest, indicating that increasing GDP per capita in these cities can increase NO 2 pollution dramatically. On the contrary, the northeastern and North China Plain cities have the lowest GDP per capita coefficients. We can conclude that the GDP per capita coefficients trend gets smaller as time increases from a temporal perspective. This reflects that the local economy has begun to intensify, allowing for a drop in NO 2 pollution per unit of GDP per capita.
All coefficients of the population density variable are observed to be positive. This is because the urban population growth will increase energy consumption, resulting in increased NO 2 emissions. As can be seen, population density expansion plays a crucial role in improving local NO 2 pollution in northern cities of China, notably in the northeast. The cities with the lowest coefficients, however, are located in southern and eastern China.
From a temporal perspective, the effect of population density on NO 2 pollution became stronger in the early study period and then turned smaller in most cities across the country, demonstrating that people have become more environmentally friendly and green in energy consumption in recent years. Figures 6 and 7 depict the temporal and spatial variations in the coefficients of the explanatory factors in the estimated results of the GTWR model from 2005 to 2019. As illustrated in Figure 6, most GDP per capita coefficients are positive, with only a few negative coefficients. GDP per capita is generally used to indicate the economic development level. It is linked because it generates a large amount of NO2 emissions due to the raw economic development pattern typically adopted in the early period. Meanwhile, the GDP per capita coefficients in these cities from southern and northwestern China are the highest, indicating that increasing GDP per capita in these cities can increase NO2 pollution dramatically. On the contrary, the northeastern and North China Plain cities have the lowest GDP per capita coefficients. We can conclude that the GDP per capita coefficients trend gets smaller as time increases from a temporal perspective. This reflects that the local economy has begun to intensify, allowing for a drop in NO2 pollution per unit of GDP per capita.  All coefficients of the population density variable are observed to be positive. This is because the urban population growth will increase energy consumption, resulting in increased NO2 emissions. As can be seen, population density expansion plays a crucial role in improving local NO2 pollution in northern cities of China, notably in the northeast. The Cities in northeastern China and northwestern China had the most significant TSRatio coefficients. The TSRatio variable of Hebei, Henan, Shandong, Jiangsu, Yunnan, Guangxi, and Hainan cities is negatively correlated and NO 2 pollution. NO x emitted by the tertiary industry is mainly from service and transportation, and the industrial NO x emission comes from the secondary industry. The negative coefficients suggest that increasing the shares of the tertiary industry may contribute to reducing NO 2 pollution. From a temporal view, we notice that the coefficients in eastern China and northeastern China decreased from 2005 to 2020. However, in western China and central China, their coefficients were increasing. We can infer that the contribution of services and transportation to NO 2 pollution is more significant than the increase in pollution levels caused by the secondary industry between 2005 and 2020. On the other hand, this also reflects that most cities in China have significantly improved their capacity to reduce industrial NO x emissions in recent years.
Most positive coefficients of foreign direct investment are concentrated in eastern and central China. Since foreign direct investment can bring advanced technologies from industrialized countries, it can help improve urban air quality to a certain extent. However, the GTWR model results indicate that changes in foreign direct investment have a negligible influence on NO 2 pollution. Additionally, it can be observed from the time dimension perspective that the effect of foreign direct investment decreases with time for most cities.
As shown in Figure 7, meteorological variables, including temperature, ambient air pressure, wind speed, and relative humidity, are also strongly correlated with the tropospheric NO 2 VCDs in various cities. Generally, increasing temperature accelerates photochemical reactions and decreases the atmospheric lifetime of NO 2 ; increasing relative humidity reduces tropospheric NO 2 by increasing the rate of NO x conversion to secondary aerosols; wind speed affects the rate of pollution diffusion and dilution in the atmosphere; and increasing air pressure increases NO 2 levels by improving atmospheric stability. From Figure 7, we find that the GTWR model results present a positive association between air pressure and NO 2 . On the other hand, increasing temperature, wind speed, and humidity can help cities alleviate NO 2 pollution.
Additionally, we normalized each explanatory variable to examine the variable's contribution to NO 2 pollution. Table 4 summarizes the estimated coefficients for each explanatory variable after standardization. We compared the mean and median values of the variables' coefficients and discovered that population density had the highest positive coefficients, 0.609 (mean) and 0.611 (median). On the other hand, humidity has the most significant negative coefficients, −0.375 (mean) and −0.426 (median). As a result, we quantified the contribution of each variable to NO 2 pollution by ranking them according to the absolute magnitude of their coefficients: population density (0.609), humidity (−0.375), GDP per capita (0.210), air pressure (0.193), temperature (−0.189), the ratio of the tertiary to the secondary industry (−0.141), foreign investment (0.066), and wind speed (−0.058).

Conclusions
This study explored the spatial and temporal variations of NO 2 pollution in the prefecture-level city from 2005 to 2020, covering the 11th, 12th, and 13th Five-Year Plan periods. The most polluting cities were located in the North China Plain. NO 2 VCDs in more than half of cities in 2020 was lower than that in 2005. In addition, 86.2% of Chinese cities experienced a rapid reduction in the NO 2 VCDs in February and March of 2020 due to the COVID-19 pandemic. It indicates that both the NO 2 reduction strategies and the COVID-19 pandemic have led to the great NO 2 reduction in the 13th Five-Year Plan. The global Moran's I results suggested that NO 2 pollution in Chinese cities was highly spatially clustered. The local Moran's I results indicated that high-high NO 2 polluted cities were located in the North China Plain regions while the low-low cities were in western China and northeastern China.
Meanwhile, the results obtained from the GTWR model analysis allow us to bring the following policy insights. First, the high-density population in cities is not conducive to reducing NO 2 pollution, especially in the cities in the northwest and northeast. Therefore, appropriate adjustment of population size can help reduce air pollution. Meanwhile, people need to be encouraged to use more clean energy and raise environmental awareness of energy-saving and conservation. Secondly, although GDP per capita is positively correlated with NO 2 pollution, the coefficient of GDP per capita has been gradually decreasing in recent years, indicating that NO 2 pollution gradually starts to reduce as income increases. Finally, if the ratio of the tertiary to the secondary industry increases, it can effectively reduce NO 2 pollution, especially in the North China Plain and the Yangtze River Delta. It can achieve a good effect on alleviating air pollution.
However, this study still has limitations that can be solved in future work. For example, the satellite observed NO 2 data used in this study are the tropospheric NO 2 vertical column densities products, which have a certain non-linear relationship with the ground-level NO 2 concentrations. Future work needs to transform tropospheric NO 2 column concentrations to ground-level NO 2 concentrations with some specific technical means for further analyses and discussions. In addition, more accurate socio-economic data could be used to replace the data from traditional yearbooks.