How do greenspace landscapes affect PM2.5 exposure in Wuhan? Linking spatial-nonstationary, annual varying, and multiscale perspectives

ABSTRACT As an ambient atmospheric pollutant, fine particulate matter (PM2.5) has posed significant adverse impacts on public health around the world. To attenuate the population exposure risk to PM2.5 pollution, greenspace has been considered as a promising approach. Little is known, however, about the attenuating impacts of greenspace landscapes on PM2.5 exposure risks at various locations, scales, and exposure levels. This study employed hotspot analysis, weighted barycenter, and time-series clustering to investigate the spatiotemporal dynamics of PM2.5 exposure across Wuhan. In addition, the multi-scale geographically weighted regression (MGWR) was used to determine the relationships between greenspace landscape patterns and yearly PM2.5 exposure over four years (2000, 2005, 2010, and 2015). Results revealed that, between 2000 and 2016, the variations in PM2.5 exposure hotspot coverages within Wuhan showed an inverse U-shape trend. The K-DTW clustering differentiated the study area into seven spatial clusters with homogeneous temporal dynamics. In general, there were three stages of fluctuations in PM2.5 exposure in Wuhan: 2000–2005, 2006–2011, and 2012–2016. MGWR also disclosed associations between PM2.5 exposure and greenspace landscape parameters (AI, ED, SI, and PLAND). PLAND of green spaces can mitigate PM2.5 exposure at a broader scale (the average bandwidth was 1391), while AI, ED, and SI are generally associated with PM2.5 exposure reduction on local scales. In Wuhan, we also confirmed such relationships between four landscape metrics with varying levels of exposure risks. The results indicate that the attenuation effectiveness toward PM2.5 exposure risk by greenspace landscapes is not only site- and scale-dependent but also affected by exposure risk levels. The findings of this study can contribute to greenspace planning and management for mitigating PM2.5-attributable adverse health impacts.


Introduction
Ambient air pollution has been considered as the primary concern for public health all over the world.According to the World Health Organization (WHO), 91% of residents worldwide are exposed to levels of atmospheric pollutants exceeding the WHO air quality guidelines (de Sherbinin et al. 2014;World Health Organization 2004, 2016).Among all types of pollutants, PM2.5 (fine particulate matter with aerodynamic dimensions of 2.5 m or less) has posed the most serious threat to public health around the world, particularly in China (Kan, Chen, and Tong 2012;World Health Organization 2016;Zhang et al. 2019b).Human exposure to PM2.5 has been linked to a variety of negative health outcomes, including lung cancer (Guo et al. 2019a), cardiovascular and cerebrovascular diseases (Chuang et al. 2005;Kloog et al. 2011), missed abortion (Zhang et al. 2019a), and even premature death (Cohen et al. 2017;Liu et al. 2019a).Since China is the world's most populous country, PM2.5 pollution has already posed a great concern to public health in China (particularly in highly populated areas).For instance, in 2010, nearly 942 million people in China were exposed to PM2.5 pollution exceeding the WHO Interim Target-1 (35 μ/ m 3 ) (Lim et al. 2012).
Attenuating the risks of PM2.5 exposure is crucial for all countries across the globe, especially for China.Fossil fuel-driven socioeconomic expansion in China has inevitably promoted ambient PM2.5 pollution (Chen et al. 2018;Fan et al. 2020aFan et al. , 2020b;;Yang et al. 2018a).To address the increasingly severe atmospheric pollution and related public health risks, the State Council of China enacted the Air Pollution Prevention and Control Action Plan (Action Plan) and the bulk of stringent air quality control measures in 2013, primarily aimed at restricting industrial emissions and coal fuel consumption (Chen et al. 2013;Geng et al. 2019;Yang and Zhao 2022;Zhang et al. 2019a).It has been investigated that such stringent measures in 2013-2017 could benefit improved air quality and associated public health developments.While the ultimate approach for PM2.5 pollution mitigation is argued to eliminate the source emissions, CONTACT Huimin Liu hmliu@whu.edu.cnSupplemental data for this article can be accessed online at https://doi.org/10.1080/10095020.2022.2085187.
other proposals, such as nature-based solutions to attenuate PM2.5 exposures, are equally necessary for developing societies, particularly for China (Chen et al. 2016;Fan et al. 2017).
Optimizing greenspace landscape planning and management has been considered as one of the most advisable nature-based solutions to urban environmental issues, e.g.air pollution removal and urban heat island (UHI) mitigation (Chen et al. 2016(Chen et al. , 2019;;He and Zhu 2018;Li et al. 2022;Wu et al. 2018).As "natural filters", greenspace are capable of purifying the air by filtering deposition, reducing re-suspension, and absorbing gaseous precursors of PM2.5 (Nowak et al. 2013).In addition, greenspace has the ability to alternate the microclimatic conditions, which can subsequently affect the deposition and dispersion of atmospheric particulate matter (Jeanjean et al. 2017;Wu et al. 2018).For the above reasons, greenspace has been proven effective in reducing PM2.5 pollution and hence improving public health around the world (Fan and Liu 2018;Jeanjean et al. 2017;Jin et al. 2014;Nowak et al. 2013).
In addition to quantifying the PM2.5 attenuation capacity of greenspace, it is also important to examine how such capacities of greenspace landscapes vary with locality, scale, and exposure level.The effectiveness of PM2.5 pollution mitigation by greenspace has been proven as site-, time-, and scale-dependent utilizing both in-site PM2.5 measurement and satellitederived PM2.5 products, which are intermediated by meteorological forcing, socioeconomic activities, and surface elements (Chen et al. 2016;Escobedo and Nowak 2009;Wu et al. 2018).In Chinese cities, greenspaces with multiple compositions and configurations have been substantiated to be able to remove PM2.5 with varying effectiveness at various seasons, locations, and scales (Chen et al. 2016;Wu et al. 2018).On the one hand, whether the association between greenspace landscapes and PM2.5 exposure risks is correspondingly site-, exposure level-, and scaledependent remains to be examined.The adoption of spatially sparse in-site PM2.5 monitoring, on the other hand, reduces the potency of massive contemporary studies examining PM2.5 exposure-greenspace correlations for supporting urban greenspace optimization (Jeanjean et al. 2017;Jin et al. 2014).
Fortunately, in this new era of earth-observation (EO), the overwhelming flood of EO-derived gridbased PM2.5 and population products allows us to investigate PM2.5 exposure-greenspace associations at spatially continuous perspectives and finer resolution (Guo et al. 2020a(Guo et al. , 2020b;;Huang and Wang 2020;Li et al. 2013;Ma et al. 2016;Wu et al. 2015a;Yang et al. 2020a).In this study, with a geographic focus on Wuhan, the megacity of Central China, a comprehensive investigation of the spatiotemporal dynamics of PM2.5 exposure patterns across Wuhan in 2000-2016 has been conducted using an annual ground-level PM2.5 dataset and gridded LandScan™ population.Furthermore, the associations between PM2.5 exposure and greenspace landscape patterns have been established in four critical years (2000, 2005, 2010, and 2015).Specifically, this study intended to answer the following questions: (1) How did the patterns of human exposure to PM2.5 vary across Wuhan in the period of 2000-2016?(2) Is the link between PM2.5 exposure and greenspace in Wuhan site-, exposure level-, and scale-dependent or not?(3) How did the greenspace landscapes in Wuhan affect the spatiotemporal patterns of PM2.5 exposure?

The study area
Wuhan, a megacity located in central China, has been selected as the study case.Located in the middle and lower reaches of the Yangtze River, Wuhan has a humid subtropical monsoon climate with abundant rainfall, sunshine, and four distinctive seasons.In the past decades, remarkable socioeconomic and landscape urbanization has been witnessed in Wuhan (Zhan, Yue, and Xiao 2018).By the end of 2019, more than 11 billion people permanently lived in Wuhan with a Gross Domestic Product (GDP) of 1.62 trillion China Yuan, according to the report by the Wuhan Municipal People's Government (WMPG) (WMPG 2020).As the industrial and economic center of central China, Wuhan suffers from severe fine particle pollution along with its remarkable socioeconomic development (Guo et al. 2019b;Xu et al. 2017).There are currently 10 PM2.5 monitoring sites (Figure 1) operated by the China National Environmental Monitoring Center (CNMEC) across Wuhan.These 10 monitoring sites are listed in Table S.1.

PM2.5 concentration
The sparsely distributed ground-level monitoring stations cannot fully capture the spatial variability of PM2.5 concentrations, thus posing difficulty in the spatially explicit delineation of PM2.5 exposure risks (Lin et al. 2016;Song et al. 2019;van Donkelaar et al. 2016).Hence, an annual ground-level PM2.5 concentration dataset downloaded from the Socioeconomic Data and Application Center (SEDAC) was adopted in this study (van Donkelaar et al. 2016(van Donkelaar et al. , 2018)).Currently, users can access the gridded annual PM2.

Population distribution
The Oak Ridge National Laboratory (ORNL) LandScan™ population distribution product has been used as the demographic input in the PM2.5 population exposure estimation.as greenspace in this study for its complex and interactive role in the emission and absorption of PM2.5 and related precursor pollutants (Lu et al. 2018;Wu et al. 2015b).g

PM2.5 population exposure estimation
Population weighted PM2.5 concentration (PW_PM2.5)has been used to quantify PM2.5 population exposure risks across Wuhan at 1km � 1km resolution.In each 1km � 1km grid, the higher population-weighted PM2.5 concentration indicates that residents are faced with higher PM2.5 exposure risks and vice versa (Lin et al. 2016;Song et al. 2019).The population-weighted PM2.5 concentration can be calculated as: where i is a grid in the study area, n is the amount of the grid in the study area.POP i and PM2:5 i are population and PM2.5 concentration in the ith grid.The calculation of population-weighted PM2.5 concentration has been conducted on MATLAB 2018b software platform.

PM2.5 exposure weighted barycenter
To better delineate the evolution of PM2.5 exposure spatial patterns, the PM2.5 exposure weighted barycenter has been introduced in this study.Considering population-weighted PM2.5 concentrations, the PM2.5 exposure weighted barycenter is capable of revealing the spatiotemporal dynamics of PM2.5 exposure patterns (Yang et al. 2019a).The coordinates of PM2.5 exposure weighted barycenter are calculated as follows: where i and n share the same definitions as Equation (1), and (x i ,y i ) is the coordinate of ith grid, and (� x,� y) is the coordinate of PM2.5 exposure weighted barycenter.

Hotspot identification
The statistically significant PM2.5 exposure hotspots from 2000 to 2016 were identified using Getis-Ord Gi* in this study.To identify a specific geographical element as a hotspot, Getis-Ord Gi* requires that such an element not only demonstrates high value itself but also to be surrounded by agglomerated high-value elements (Ord and Getis 1995).The z-value in Getis-Ord Gi* is a measure of spatial agglomeration.The statistical significance of a hotspot is determined by its z-value.A z-value close to or equal to zero suggests that certain elements have no spatial agglomeration.The calculation of Getis-Ord Gi* can be checked in the references (Ord and Getis 1995).Getis-Ord Gi* is implemented in the spatial econometric software GeoDa version 1.14.0.

Time-series clustering
Although various research have attempted to investigate the temporal dynamics of PM2.5 at scales ranging from global to regional, the spatial variations of PM2.5 patterns have rarely been fully disclosed (Hu et al. 2013;Peng et al. 2016).Time-series clustering could be a viable approach to capture the temporal dynamics of PM2.5 and locate the specific regions with homogeneous time-varying patterns across Wuhan from 2000 to 2016 (Hu et al. 2013;Liu et al. 2018;Peng et al. 2016).In this study, K-Means clustering coupling DTW (K-DTW) has been utilized to capture and spatialize the temporal dynamics of PM2.5 exposure in Wuhan.K-Means has been acknowledged as the most efficient and robust tool for exploring the temporal dynamics of environmental challenges as a raw databased time-series clustering method (Aghabozorgi, Shirkhorshidi, and Wah 2015;Liu et al. 2019b;Warren Liao 2005).However, Euclidean distance in K-Means clustering performs poorly when processing time-series with small sample size and inconsistent phase (Aghabozorgi, Shirkhorshidi, and Wah 2015;Warren Liao 2005).As the extension of Euclidean distance, Dynamic Time Warping (DTW) has been proposed as an extension of Euclidean distance to uncover temporal commonalities concealed behind phase variation using a dynamic planning technique (Aghabozorgi, Shirkhorshidi, and Wah 2015;Liu et al. 2018;Warren Liao 2005).The mathematical calculation of K-Means is well documented (Aghabozorgi, Shirkhorshidi, and Wah 2015;Liu et al. 2019bLiu et al. , 2018;;Warren Liao 2005), so it will not be repeated here.The definition of DTW is as: is the local temporal distance, w n is the weighted function of DTW.In this study, the K-DTW time-series clustering algorithm has been implemented on MATLAB 2018b software platform.

Landscape metrics
Landscape metrics have been recognized as feasible and effective in quantifying LU/LC landscape patterns, i.e. composition and configuration (McGarigal 1995).In this study, 10 metrics reflecting the aggregation, density, edge, shape, and coverage of greenspace landscape patterns have been selected.The detailed information about these landscape metrics is reported in Table S.2.All the selected landscape metrics are calculated at 1km � 1km resolution using Fragstats 4.2 software.

Multi-scale geographically weighted regression
The multi-scale geographically weighted regression (MGWR) was used to model the relationships between greenspace landscapes and PM2.5 exposure.Before the creation of MGWR, the geographically weighted regression (GWR) has been widely used in modeling various geographical issues, e.g.house prices, air quality, and land surface temperature (Chen et al. 2018;Fotheringham, Brunsdon, and Charlton 2003;Liu et al. 2019c).In GWR, all geographical processes are assumed to operate at the same scale (i.e. the same bandwidth).Such an assumption in GWR has now been argued to be probably inappropriate under certain scenarios.There is no single "correct" scale where it can model and simulate all the unique geographic and environmental problems (McMaster and Sheppard 2004).By relaxing the assumption that processes operate at one spatial scale, MGWR has been verified as a more powerful model than GWR (Fotheringham, Yang, and Kang 2017;Fotheringham, Yue, and Li 2019).In this study, MGWR was used to examine the multi-scale and spatial non-stationary associations between PM2.5 exposure and greenspace landscape patterns across Wuhan.MGWR can be mathematically represented as: where u g ; v g À � is the coordinate of the location, g,y g and x ig are dependent and independent variables, and ε g is the residual.Furthermore, bw0 is the bandwidth of intercept, bwi is the bandwidth of ith independent variable.In this study, we selected the bi-square kernel function in MGWR.
where w ghi is the spatial weight of the specific observation at location g with neighboring observation located at g for hth independent variable, h is the distance between g and g.In this study, the MGWR model was implemented on the MGWR Version 2.0.1 software.

Spatiotemporal patterns of PM2.5 exposure across Wuhan from 2000 to 2016
As shown in Figure 2, from 2000 to 2016, the areas where residents have been exposed to severe PM2.5 pollution (more than twice the WHO guidelines (35 μg/m3)) have been expanding.People who suffered from such long-term PM2.5 exposure would face an additional risk of premature death up to 15% (Cohen et al. 2017;World Health Organization 2004, 2016).Previous studies have revealed that from 2010 to 2014, the substantial increase in PM2.5 exposure risks was linked to pollutant emissions and built-up area expansion during the period (Xu et al. 2016;Zhang et al. 2018).
K-Means coupling DTW (K-DTW) was used to spatially differentiate the study area into clusters with homogeneous temporal dynamics of PM2.(Han et al. 2018).Clusters 4 and 5 experienced significant high-level of PM2.5 exposure in period of 2004-2006 (41.97 μg/m 3 of cluster 4 and 60.89 μg/m 3 of cluster 5 on average), which is consistent with built-up land expansion and intensification in the corresponding regions (Zhan, Yue, and Xiao 2018).
Table 1 presents the detailed information of primary surface attributes covered by the seven timeseries clusters.Figure 4 shows that the time-series clusters of PM2.5 exposure patterns across Wuhan, as sorted according to the average populationweighted PM2.5 concentration.The general temporal variations of PM2.5 exposure over the land surface property were covered by each time-series cluster.The stages of landscape urbanization were not revealed in the corresponding time-series clusters, contrary to our earlier studies on land surface temperature (LST) in Wuhan (Liu et al. 2019b).Such inconsistency can be attributed to the speculation that, ambient air pollution, compared with LST, is  It is mainly recognized as primary business areas and high-rise high-density settlements.In Jianghan and Jiang'an District, the coverage of this cluster extends from Baofeng to Xunlimen and Huangpu Road.In Wuchang District, this cluster mainly covers Simenkou, Zhongnan, Jiedaokou, and Xudong sub-centers.
a product of more fickle and complex surface-societyatmosphere interactions.Although surface element changes have been identified as one of the external forcings of PM2.5 pollution, urban landscape transformations alone are insufficient to explain the spatiotemporal dynamics of PM2.5 pollution patterns (Chen et al. 2016;Yang et al. 2018b).
The weighted barycenters (Figure 5) and hotspots (Figure 6) of PM2.5 exposure were also used to present a better understanding of the spatial evolution of PM2.5 exposures in Wuhan.Generalized from timeseries clustering, PM2.5 exposure weighted barycenters, and PM2.5 exposure hotspots, the spatiotemporal evolution of PM2.5 exposure patterns in Wuhan can be distinguished into three periods.

Associations between PM2.5 exposure and green spaces landscapes
MGWR was used in this research to examine the associations between PM2.5 exposure and green space landscapes throughout four years (i.e. 2000, 2005, 2010, and 2015).On the one hand, collecting greenspace maps every year from 2000 to 2016 is arduous and unnecessary.The calibration and validation of MGWR, on the other hand, is time-consuming.As a result, the associations established on four representative years in three stages were believed to be sufficient to explore PM2.5 exposure attenuation by green spaces across years, sites, and pollution levels.
Following the guidelines in our previous study, Random Forest (RF) regression was adopted to select the optimal green space landscape metrics from nominated metrics (including PLAND, ED, TE, TCA, PD, SI, PARA, CONTIG, AI, and SHDI) in four years (Yang et al. 2019b).The optimal metrics were successively determined in this study when the increased mean square error (%IncMSE) almost remained unchanged and the out-of-bag (OOB) error no longer decreased.The variable screening comes to an end in the sixth iteration, as shown in Table S5, and PLAND, ED, LSI, and AI were chosen as the best landscape metric combination.Multicollinearity and crosscorrelation were removed with the combination of PLAND (VIF = 1.73),ED (VIF = 1.42),SI (VIF = 1.15), and AI (VIF = 1.68), as measured by the variation inflation factor (VIF).It is noteworthy that the optimal variable sets in this study were collectively affected by the study region, period, and scale of metrics.Thus, the optimal landscape metrics in other sites and times may not keep the same.
In each year, the association between PM2.5 exposure and green space landscape metrics was established as follows: Ordinary Least Square (OLS) regression and GWR were also implemented for comparison, and the regression performance was summarized in Table S.6.The performance of regressions was evaluated by R-square (R 2 ), Akaike Information Criterion correction (AICc), and Sum of Squared Residuals (SSR) (Fotheringham, Yang, and Kang 2017;Liu et al. 2019c).MGWR's superior performance to GWR and OLS (Table S.6,) revealed that MGWR's multi-scale assumption facilitated the examination of associations between green spaces and PM2.5 exposure.
In general, the decrease in AICc and SSR from OLS to GWR and MGWR indicated that model fitting had improved.From 2000 to 2015, R2 increased in OLS (0.41 to 0.48), GWR (0.66 to 0.79), and MGWR (0.75 to 0.87) correspondingly, demonstrating that PM2.5 exposure attenuation was better explained by greenspace landscapes.According to Section 4.3, the year of 2005 witnessed the highest level of population exposure to PM2.5 pollution in Wuhan among all four years.In 2005, however, the associations between PM2.5 exposure and greenspace landscapes were weaker than the other three years.That is to say, greenspace landscapes were found to be less effective in decreasing risks of PM2.5 exposure when such risks were themselves relatively higher.Furthermore, in Wuhan, population exposure and its spatiotemporal dynamics were not entirely explained by greenspace landscape improvements (R 2 = 0.78 on average), showing that greenspace landscape improvements can be a feasible but not exclusive solution for PM2.5 exposure reduction.
Table 2 shows that the multi-scale effects in the associations between PM2.5 exposure risks and greenspace landscapes were underlined by the considerable difference between the optimal bandwidths discovered by GWR and MGWR.Almost all landscape metrics (except PLAND) operated at a local scale, according to GWR and MGWR's results, implying the significance of spatial heterogeneity in associations between PM2.5 exposure and greenspace landscape metrics.GWR and MGWR's optimal bandwidths in 2005 were narrower than in 2000, 2010, and 2015.Such results revealed that exposure attenuation by greenspace landscapes operated at narrower scales under higher exposure risks.In general, the operational scale of PLAND remained the broadest, followed by AI, ED, and SI.The spatial heterogeneity of the association between PM2.5 exposure and greenspace landscapes found above was further highlighted in the spatial distribution of local R 2 of MGWR (Figure 7).Over the junction of Wuchang and Hongshan Districts, as well as the central areas of Huangpi, Jiangxia, and Hanan Districts, strong associations between PM2.5 exposure and green spaces were demonstrated in 2000.In 2005, the association between PM2.5 exposure and green areas was very weak across Wuhan, with the exception of Jiangxia District's downtown.Furthermore, the spatial patterns of local R2 in 2010 and 2015 were very similar.In regions around Wuhan Iron and Steel Cooperation in Qingshan District, Changfeng sub-district in Qiaokou District, and subcenters of Huangpi, Xinzhou, and Jiangxia District, green spaces had a greater impact on PM2.5 exposure.The regression coefficients of landscape metrics in 2000-2015 were visually represented in Figures 8-11.
From 2000 (7.32) to 2005 (6.18), the absolute value of the average AI coefficient declined by 1.14, subsequently increasing to 8. The spatial patterns of regression coefficient of SI demonstrated significant spatial heterogeneity across Wuhan (standard deviation ≥0.96), and the   The spatial patterns of the PLAND coefficient were quite smoothing in 2000-2015 (average standard deviation = 0.39), showing that PLAND's PM2.5 mitigation effect functioned at broader spatial scales than AI, ED, and SI, as shown in Figure 11.The PLAND coefficient was found to be higher in the downtown regions of Jianghan, Jiang'an, Qiaokou, Hanyang, and Wuchang District.Between 2000 and 2015, the absolute values of the mean regression coefficients for all four green space landscape metrics increased.Specifically, the absolute values of their average regression coefficients were the lowest in 2005, indicating that the mitigation ability of greenspace landscapes on PM2.5 exposure was suppressed when the risk of PM2.5 exposure was high.

Implications for greenspace landscape planning
Based on the multi-scale, temporal varying, and spatial non-stationary associations between PM2.5 exposure risks and greenspace landscapes investigated in Section 4.2, this study has several implications for greenspace development across Wuhan as follows: The strong associations between PLAND of green spaces and PM2.5 exposure imply that the increase in greenspace landscape proportion can significantly attenuate PM2.5 exposure risks.Additionally, the larger bandwidths of PLAND indicate that the attenuation effectiveness of improving the percent cover of green spaces was relatively consistent across the study area.Overall, increasing the coverage of green spaces across Wuhan is the most direct and practical approach to proactive precautions for PM2.5 exposure.
It is highly recommended for urban planners to create ecological corridors to bridge distinct greenspace patches.This study revealed that there were strong spatial relationships between AI of greenspace and PM2.5 exposure risks around the key PM2.5 population exposure hotspots, according to the MGWR regression results (Figure 8).As a result, greenspace patches that are spatially aggregated near primary PM2.5 exposure hotspots, such as urban functional centers and settlements, can help release people from severe PM2.5 exposure risks.
Greenspace landscapes with excessive edge density could only provide limited effect in attenuating PM2.5 exposure.The results of MGWR indicate that although the local attenuation effect of PM2.5 by ED of greenspace was witnessed to a certain extent around PM2.5 exposure hotspots, the attenuation effect by ED was not strong, and the scope of influence was quite limited.Thus, in greenspace planning and implementation, it is not recommended to attenuate exposure risks through improving the edge complexity of greenspace landscapes.
PM2.5 exposure risks can be reduced by increasing the diversity of greenspace configurations on the outskirts of exposure hotspots in downtown regions.Higher SI of greenspace can mitigate PM2.5 exposure risks in the periphery of PM2.5 exposure hotspots in downtown Wuhan at local scales, according to MGWR's spatial patterns of SI regression coefficient, while improving the shape richness of green spaces within hotspots has limited benefit for PM2.5 exposure risk reduction.
Despite the significant spatial heterogeneity of mitigation effects by greenspace on PM2.5 exposure, the local authorities are still strongly recommended to plan and construct greenspace (e.g.gardens, Aleetrees, and other green infrastructure) as many sites as possible within central downtown.Although the economics and synergistic ecological effects (e.g.urban heat island mitigation (Ouyang et al. 2020;Yang et al. 2020bYang et al. , 2017;;Yu et al. 2017)) of this approach are unclear, it has been widely demonstrated to significantly reduce anthropogenic air pollution within cities (Nowak et al. 2013).

Limitations
There are several limitations in this study that can be further discussed.Firstly, the spatiotemporal dynamics of PM2.5 exposure in 2000-2016 have been investigated annually at 1km � 1km resolution.However, in the context of heterogeneous urban landscapes, the 1km � 1km may not be the most ideal spatial resolution for the exploration of spatiotemporal dynamics and underlying driving mechanisms of PM2.5 exposure risks.Up to now, 1km � 1km PM2.5 estimations are still too coarse to perfectly depict the spatial patterns of PM2.5 exposure patterns and the underlying interactions with urban landscapes (Zhang et al. 2018).The resolution improvement of PM2.5 estimation can open the door for new insights into detailed PM2.5 spatial variations in heterogeneous urban landscapes (Hu et al. 2013;Zhang et al. 2018).Secondly, the annual interval prevented this study from capturing more detailed temporal variations of PM2.5 exposure risks, e.g.daily, monthly, and seasonal variations.The adverse impacts of PM2.5 pollution on human health are not only attributed to long-term exposure but also resulted from severe short-term PM2.5 exposure (Liu et al. 2019a).Thus, the implications for greenspace landscape improvement based on long-term PM2.5 exposure variations investigated in this study may not necessarily have the same attenuating effects on short-term PM2.5 exposure risks (Escobedo and Nowak 2009;Jeanjean, Monks, and Leigh 2016).Furthermore, the spatiotemporal dynamics and exposure attenuation by greenspace investigated in this study may not be applicable to other atmospheric pollutants.Nowadays, more atmospheric pollutant datasets (e.g.black carbon, PM10, NO x , SO x, and O 3 ) are accessible, suggesting that more efforts can be made to study the health impacts of other pollutants in future studies.

Conclusions
Without a doubt, understanding the spatiotemporal dynamics of PM2.5 exposure risks and how attenuating effects of green spaces vary with places, times, scales, and exposure risk levels is crucial in combating increasingly intense and frequent atmospheric pollution in urban environments.In this study, MGWR was performed to examine the spatiotemporal dynamics of PM2.5 population exposure and underlying associations between green space landscape patterns in Wuhan, one of the megacities in central China troubled by PM2.5 pollution.The specific findings can be summarized as follows: The expansion and intensification of hotspots were mostly observed across core downtown and suburban district downtowns when analyzing in combination with K-DTW clustering, exposure hotspots, and exposure weighted barycenter.The spatiotemporal variations of PM2.5 exposure across Wuhan may be divided into three periods based on the geographic evolution of exposure patterns and temporal dynamics: 2000-2005, 2006-2011, and 2012-2016. PM2. PM2.5 exposure-greenspace landscape associations were site-, exposure level-, and scale-dependent, according to the modeling results of MGWR.PLAND was associated with PM2.5 exposure risk at a broader spatial scale than AI, ED, and SI.In 2005, when the risk of PM2.5 exposure was at its peak, greenspace metrics have a suppressed mitigating effect on PM2.5 exposure.In addition, in this year, the impacts of all four landscape indices were performed on a relatively narrower spatial scale than the other three years.
This study contributes to the growing body of literature concerning atmospheric environment improvement through green spaces.It can facilitate the collaboration between urban air quality studies and urban green space landscape planning, designing, and management.Although the associations between PM2.5 exposure and green space landscapes were established in four representative years only, one should not devalue the significance of green spaces in attenuating long-term PM2.5 exposures.Future studies are called to examine the interactions between PM2.5 exposure and green spaces in cities under various meteorological and socioeconomic conditions.Additionally, future studies can further dedicate to exploring the most cost-effective and environmentally sound greenspace-based solutions to PM2.5 pollution, considering vegetation species combination and neighboring urban functional conditions.
5 concentration from 1998 to 2016 at the 0:01 � � 0:01 � resolution (approximately 1km � 1km in China).This study utilized a re-projected subset of the global dataset covering Wuhan at 1km � 1km resolution in 2000-2016.The validation accuracy of satellite-derived PM2.5 concentration data can be checked in the Supplementary Materials (Figure S.1).

Figure 1 .
Figure 1.The study area and the spatial distribution of the 10 air quality monitoring sites (①-⑩) within Wuhan.
5 exposure from 2000 to 2006.The most appropriate category number of PM2.5 exposure time-series of K-DTW clustering was determined using the explained variance (Figure S.2).Seven clusters (explained variance = 0.9019) have been proposed as the ideal number of K-DTW clusters based on the explained variance.
Figure 3 shows the K-DTW time-series clustering results (c).

Figure 3
Figure 3(c) shows the temporal variations of PM2.5 exposure patterns of seven time-series clusters, which show exceptional diversity and strong non-linearity (p0.001).The first time-series cluster, which mostly consisted of natural landscapes, showed an inverted U-shape trend that peaked in period of 2006-2008.In the second cluster, PM2.5 exposure risks reached an alltime high in 2007 and 2014 and stayed in a moderate high level from 2008 to 2016.Such temporal variation is in line with the renewal and development of urban landscapes in the corresponding regions (Zhan, Yue, and Xiao 2018).Specifically, significant built-up land expansions in Wujiashan, Zhuankou, and the Donghu Gaoxin sub-district covered by the second cluster were witnessed in 2005-2010 (Xu et al. 2016).Clusters 3, 6, and 7 demonstrated similar PM2.5 exposure risk temporal dynamics, which may be divided into three stages:2000-2005, 2006-2011, and 2012-2016.Since 2006, such clusters primarily covered downtown areas that had not experienced substantial landscape renewal.Exogenous PM2.5 transmission(Xu et al. 2017) and population rearrangement (average LandScan

Figure 3 .
Figure 3.The time-series clustering results by K-DTW was shown in Figure 8(c).Spatial distribution of (a) inter-annual average, (b) standard deviation, and (c) time-series clustering result of PM2.5 exposure patterns in Wuhan in 2000-2016.

( 2 )
The second time period was from 2006 to 2011.PM2.5 exposure risks increased significantly from 2006 to 2007 (28.31 μg/m3 to 38.87 μg/ m3), then dropped by 4.65 μg/m 3 from 2007 to 2011, particularly in time-series clusters 3, 6, and 7. From a spatial perspective, the urban outskirts (e.g. the Donghu Gaoxin zone in the Hongshan District in Figure 6) and downtown suburban districts have seen a great expansion of PM2.5 exposure hotspots (including Caidian, Huangpi, Xinzhou, and Jiangxia Districts).During the second stage, it is evident that the weighted barycenters have shifted back to the south (particularly in 2006-2007).(3) The third time period is from 2012 to 2016.The spatial patterns of PM2.5 exposure stayed almost unchanged in the third period, with slight shrinkage since 2014.PM2.5 exposure barycenters revealed a significant but ruleless aggregating pattern (Figure 5), which aligned with the evolution of spatial patterns of hotspots (Table

Figure 4 .
Figure 4.The average temporal variations of PM2.5 exposure risks in seven clusters and averaged across Wuhan from 2000 to 2016.
41 and 8.70 in 2010 and 2015.PM2.5 attenuation by agglomerated greenspace weakened in 2000-2005, then strengthened in 2010-2015, according to such variations in the absolute value of the average coefficient.Furthermore, between 2000 and 2015, the standard deviation of the regression coefficient increased (from 0.77 to 1.21), showing that spatial heterogeneity in PM2.5-greenspace associations was more significant in 2010 and 2015 than in 2000 and 2005.Within downtown Wuhan, AI had impacts on PM2.5 exposure in increasingly wider regions from 2000 to 2015.Furthermore, AI of green space landscapes generally had significant impacts on PM2.5 exposure around downtowns of suburban districts (e.g.Huangpi, Xinzhou, and Jiangxia District).The lower absolute value of the average regression coefficient indicates that ED had weaker effective impacts on PM2.5 exposure in comparison to AI.In terms of geographic patterns, the spatial patterns of ED coefficients remained relatively stable from 2000 to 2015, as evidenced by standard deviations of ED coefficients (0.61 of 2000, 0.59 of 2005, 0.60 of 2010, and 0.58 of 2015).Greenspace ED and PM2.5 exposure were shown to be more closely associated in the Zhongnan region of Wuchang, the Honggangcheng region of Qingshan, and the Jianghan, Jiang'an, and Qiakou Districts.

Figure 7 .
Figure 7. Spatial distribution of local determination coefficient (Local R 2 ) modeled by MGWR in the four years.
average regression coefficients of SI represented similar variations to those of AI, i.e. weakened from 2000 to 2005 and strengthened in 2010-2015.The average coefficient of SI increased from −6.32 to −5.71 in 2000-2005 and subsequently decreased by 1.59 in 2010-2015.High SI regression coefficients, unlike AI and ED, were not concentrated near PM2.5 exposure hotspots in downtown Wuhan.On the contrary, the strong associations between SI and PM2.5 exposure risks were mainly witnessed in the peripheral regions of PM2.5 exposure hotspots, e.g.regions near the Jinyi Lake in Dongxihu District and the Yanxi Lake in Qingshan District.

Figure 11 .
Figure 11.Regression coefficients of percentage cover of green spaces landscapes (PLAND) from 2000 to 2015.

Table 1 .
The primary surface properties of seven time-series clusters.Junctures between built-up land and neighboring natural landscapes, including the beach of the Yangtze River located in Wuchang and Qingshan District, regions between Sha Lake and the East Lake, as well as Changfeng sub-district in Qiakou District.4Low/middle-risehigh-density residential areas underwent significant renewal in 2004-2006, for instance, Wutaizha and Xujiapeng region in Wuchang, Yingwujie Qinduankou regions in Hangyang, as well as Dazhilu sub-district in Jianghan.
5 Middle/high-rise middle density residential areas without significant renewal after 2012, including the Jiyuqiao region in Wuchang, Yuehu region in Hanyang, and Beihu and Huaqiao sub-district in Jianghan and Jiang'an.6 Middle/high-rise high-density residential areas are located in urban centers of Wuhan, and there was no remarkable change since 2006.This cluster generally includes Luonan and Shouyi sub-district in Wuchang, Hanzhengjie sub-district in Jianghan, Gutian, and Hanshuiqiao subdistrict in Qiaokou.7

Table 2 .
Optimal bandwidths of GWR and MGWR from 2000 to 2015.