Next Article in Journal
Bioactive Compounds Produced in Leaves of Mulberry (Morus alba L.) Transplants under Modified Environments of Root and Aerial Zones
Next Article in Special Issue
Responses of Forest Net Primary Productivity to Climatic Factors in China during 1982–2015
Previous Article in Journal
Identification of Structural Features of Hydrocinnamic Acid Related to Its Allelopathic Activity against the Parasitic Weed Cuscuta campestris
Previous Article in Special Issue
Varying Responses of Vegetation Greenness to the Diurnal Warming across the Global
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Does the Spatial Pattern of Plants and Green Space Affect Air Pollutant Concentrations? Evidence from 37 Garden Cities in China

1
College of Landscape Architecture, Nanjing Forestry University, Nanjing 210037, China
2
Research Institute of Architecture, Southeast University, Nanjing 210096, China
3
Graduate School of Human-Environment Studies, Kyushu University, Fukuoka 819-0395, Japan
4
State Key Laboratory of Vegetation and Environmental Change, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China
*
Authors to whom correspondence should be addressed.
Plants 2022, 11(21), 2847; https://doi.org/10.3390/plants11212847
Submission received: 15 September 2022 / Revised: 9 October 2022 / Accepted: 24 October 2022 / Published: 26 October 2022

Abstract

:
Relevant studies have demonstrated that urban green spaces composed of various types of plants are able to alleviate the morbidity and mortality of respiratory diseases, by reducing air pollution levels. In order to explore the relationship between the spatial pattern of urban green spaces and air pollutant concentrations, this study takes 37 garden cities with subtropical monsoon climate in China as the research object and selects the urban air quality monitoring data and land use type data in 2019 to analyze the relationship between the spatial pattern and the air pollutant concentration through the landscape metrics model and spatial regression model. Moreover, the threshold effect of the impact of green space on air pollutant concentrations is estimated, as well. The results showed that the spatial pattern of urban green space was significantly correlated with the concentrations of PM2.5 (PM with aerodynamic diameters of 2.5 mmor less), NO2 (Nitrogen Dioxide), and SO2 (Sulfur dioxide) pollutants in the air, while the concentrations of PM10 (PM with aerodynamic diameters of 10 mmor less) pollutants were not significantly affected by the green space pattern. Among them, the patch shape index (LSI), patch density (PD) and patch proportion in landscape area (PLAND) of forest land can affect the concentration of PM2.5, NO2, and SO2, respectively. The PLAND, PD, and LSI of grassland and farmland can also have an additional impact on the concentration of SO2 pollutants. The study also found that there was a significant threshold effect within the impact mechanism of urban green space landscape pattern indicators (LSI, PD, PLAND) on the concentrations of PM2.5, NO2, and SO2 air pollutants. The results of this study not only clarified the impact mechanism of the spatial pattern of urban green space on air pollutant concentrations but also provided quantitative reference and scientific basis for the optimization and updating of urban green space to promote public health.

1. Introduction

The emergence and development of modern urban green spaces are closely related to the improvement of human living environments and public health [1]. As an important discipline for the development of human living environments, landscape architecture provides important support in the establishment of public safety systems in urban and rural spaces and promotes the health and wellness of residents [2]. With the accelerated global urbanization and rapid industrial development, including the frequent occurrence of public health incidents such as the new coronavirus pneumonia, the public health of residents is under unprecedented threat [3]. The severe challenges posed by the pandemic have prompted the entire society to deeply recognize the importance of public health and have triggered unprecedented attention to the urban built environment and the construction of healthy urban green spaces.
As a global environmental and public health problem, air pollution has severe adverse effects on human health. The pollutants can enter the body through the respiratory system and affect the lungs and heart, causing cardiovascular and respiratory diseases [4]. The World Health Organization (WHO) report demonstrated that respiratory diseases caused by air pollution in 2019 ranked fourth in the top 10 causes of death worldwide (https://www.who.int/zh/ (accessed on 15 July 2021)). Numerous studies have verified that short-term or long-term exposure to air pollutants, including PM2.5, PM10, NO2, and SO2, increased the risks of mortality and morbidity, thereby posing a severe threat to the public health of the residents [5]. As an important part of the urban ecosystem, in addition to providing residents with green spaces for recreation and entertainment, urban green spaces also improve air quality. The ability of different plant species, plant communities, and green spaces along roads to retain dust and mitigate air pollution has been widely demonstrated [6]. In addition, by optimizing the structure of urban green spaces, increasing the areas of green spaces and green coverage could effectively reduce the concentration of particulate matter and gas pollutants in the air, which effectively improves the urban environment, and ultimately play an important role in promoting public health [7].
Domestic and international research on the mitigation of air pollutants in green spaces at the micro-level has primarily focused on the ability of urban garden green spaces to reduce air pollutant concentrations, using individual plant microstructure processes such as sedimentation, retardation, adsorption, and absorption. For example, Latha and Highwood demonstrated that changes in structures, such as the roughness of plant leaf surfaces, affect the sedimentation pattern of dust particles [8]. Beckett et al. deduced that plants maintained a higher humidity in a certain range during transpiration, and dust sedimentation is more likely to occur when it increases in weight after absorbing moisture, while the ability of the leaves to adsorb dust increases with an increase in their humidity [9]. In addition, numerous studies have demonstrated the variability in the absorption of different gaseous pollutants by different landscape plants [10,11,12].
The meso-level was adopted to investigate the reduction effect of the size, shape, plant configuration, vertical structure, etc., of small- and medium-scale urban green spaces on different pollutants. In his study on the ecological mechanism of urban open space planning, Wang Shaozeng pointed out that the more practical the mix of green space levels, the better the filtration effect on the atmospheric particulate matter [13]. By analyzing the relationship between the three-dimensional green volume of green spaces along roads and PM2.5 concentration, Sheng found that high 3D green volume did not indicate low PM2.5 concentration, and such green spaces with the uniform vertical distribution of biomass and diverse vegetation were more effective in reducing PM2.5 concentration [14]. Fan analyzed the correlation between the daily PM10 and PM2.5 concentrations and particulate matter concentrations of seven typical land cover types and different scales of land cover patterns and determined that pavement-type and low-to-medium canopy density vegetation exerted a more significant effect on PM10 levels, while PM2.5 concentrations were more sensitive to the response of building-type and low-to-medium canopy density vegetation [15]. In several ways, the aforementioned studies verified that green spaces, as living media, are one of the most important vehicles for mitigating air pollution and play a significant role in the public health of residents and the urban environment.
Recently, studies on how to optimize the landscape pattern of urban green spaces to reduce air pollution concentrations at the macro-level have increasingly garnered considerable attention, and the current research results primarily focused on exploring the correlation between the concentrations of PM2.5, PM10, and other pollutants, including land use or the landscape patterns of land cover. Ye et al. [16] explored the relationship between PM2.5 growth and land use changes in China from 1998–2015 and inferred that PM2.5 concentrations were higher in the eastern plains and Taklamakan Desert in China, and higher PM2.5 concentrations existed on artificial land surfaces, croplands, and deserts, while forests, grasslands, and unused land usually contained lower PM2.5 concentrations. Simultaneously, the average annual increase in PM2.5 concentrations in a desert land and artificial land surfaces was higher than that of other land types. Yue et al. [17] studied the quantitative relationship between vegetation coverage and atmospheric particulate matter based on remote sensing inversion and deduced that vegetation coverages of ≤10% and >45% exhibited a significant effect on mitigating atmospheric particulate matter pollution. Lei et al. [18] explored the effect of green spaces on particulate matter pollution by studying the landscape patterns of urban green spaces at multiple spatial scales and determined that increasing the biodiversity of green spaces and increasing the number of large green spaces significantly reduced PM10 concentrations at almost all scales. Zhao et al. [19] conducted a land use regression (LUR) analysis on the green spaces of lakes and wetlands, including the surrounding 500-m built environment in Wuhan City, and their results demonstrated that the lakes, wetlands, and nearby greenery exerted a positive and significant effect on PM2.5 concentrations within a buffer zone of 300 m or closer. Via a multivariate linear regression modeling of the daily average PM10 concentration data with land use pattern information for the cities of Vienna and Dublin, McNabola et al. [20] determined that adding transboundary air pollution and traffic activity representations to the predictor variables significantly improved the accuracy of LUR-based methods. Zhang et al. [21] applied the LUR model to analyze the correlation between air pollution levels and childhood asthma hospitalization rates, by establishing a spatial distribution LUR model of the daily pollutant concentration data of PM10 and SO2 with associated influencing factors in Shenyang, China and deduced that the number of childhood asthma hospitalizations was highly correlated with PM10 and SO2 pollution levels. Lee and Koutrakis adopted satellite ozone monitoring instruments together with land use parameters to develop a mixed-effects model through which they estimated the daily NO2 concentrations in New England, and explored the source areas of emissions (e.g., high population or traffic areas) in the study area and elucidated the seasonal characteristics of NO2, based on NO2 spatial distribution patterns [22]. Lu et al. [23] analyzed the landscape pattern indices and PM2.5 data in China and demonstrated that differences exist in the significant impact indicators of different land use types on PM2.5 concentration and that the landscape pattern indices exhibit a significant effect on PM2.5 concentrations. De Jalón et al. [24] conducted a comparative analysis on the effect of the dry deposition of trees on air pollution reduction for different land use types in the Basque Country, and finally deduced that coniferous forests are the most effective in eliminating air pollution. These research results mainly elucidated the effect of urban land use changes on univariate air pollutants, and the pollutants were mainly PM2.5, PM10, and other fine particulate matter. A fewer number of studies have been conducted on NO2, SO2, and other gaseous pollutants, and these studies rarely involved any comprehensive air pollution research on the quantification and regulation strategies of the concentration of mixed pollutants of gases and particulate matter, including the spatial distribution of urban landscape patterns.
On this basis, this study takes 37 garden cities with subtropical monsoon climate in China as the research object, selects the urban air quality monitoring data and land use remote sensing data in 2019 to analyze the relationship between the landscape pattern index and the air pollutant concentration through the spatial regression model method. Also, a new threshold effect estimation method based on a polynomial model is designed to explore the threshold effect of green space landscape patterns on air quality [25,26].

2. Results

In this paper, regression modeling tools in GeoDa software are used to conduct SEM regression analysis, and the results obtained are shown in Table 1.

2.1. Regression Analysis of Landscape Pattern Indices and Air Pollutants

PM2.5 concentrations were significantly and negatively correlated with the LSI of forestlands; each unit increase in the landscape shape index of the forestlands was followed by a 0.68 unit decrease in the PM2.5 pollutant concentration. NO2 concentration was significantly and negatively correlated with the PD of forestlands. With each unit of increase in the PD of forestlands, the concentration of the NO2 pollutant was subsequently reduced by 258.409 units. However, Grasslands and farmlands had no significant correlation with PM2.5 and NO2 pollutant concentrations. SO2 concentrations were significantly and positively correlated with the PD of forestlands, PLAND, and LSI of grasslands, where the relationship between SO2 pollutant concentrations and the PD of forestlands and LSI of grasslands was highly significant (p < 0.01). SO2 concentration was significantly and negatively correlated with the PLAND of forestlands, PD of grasslands, and PLAND of farmlands, where the PD of grasslands and PLAND of farmlands were highly significantly correlated with the SO2 pollutant concentration (p < 0.01). When the PD of forestlands, PLAND, and LSI of grasslands increased by one unit, the SO2 concentration increased by 149.939, 0.752, and 0.429 units, respectively. When each unit of PLAND of forestlands, PD of grasslands, and PLAND of farmlands increased, the SO2 concentration was reduced by 0.073, 214.564, and 0.172 units, respectively. In contrast, SO2 concentrations were not significantly affected by the LSI of forestlands, and the PD and LSI of farmlands.
From the results of the spatial correlation analysis between the landscape pattern indices of green spaces and air pollutants in the 37 cities, it was determined that the landscape pattern of urban green spaces was significantly associated with the concentrations of PM2.5, NO2, and SO2 pollutants in the air, while the concentrations of PM10 pollutants were not significantly affected by the pattern of green spaces. Meanwhile, LSI, PD, and PLAND of forestlands influenced the concentration of the three types of air pollutants, respectively. In addition, the PLAND, PD, and LSI of grasslands and PLAND of farmlands exert an additional influence on the concentration of SO2 pollutants.

2.2. Threshold Effect of the Impact of Landscape Pattern Indices on Air Quality

To further validate the relationship between the landscape pattern indices and the concentrations of the four pollutants, this study conducted regression analysis and linear fitting on the four pollutant concentrations and their landscape pattern indices, with EXCEL revealing significant effects on them. After comparing the R2 of the function’s trend line, it was deduced that the overall effect of the polynomial fit was better than several other types of functions in the experiment; hence, the polynomial was chosen for linear fitting in this study.

2.2.1. Threshold Effect of the Impact of Landscape Pattern Indices on PM2.5

Table 1 indicated that the LSI of forest land has a significant negative effect on the concentration of PM2.5 (p < 0.1). Therefore, this study analyzes the influence trend of PM2.5 concentration change to determine the threshold range of forest land LSI. It is found that the concentration of PM2.5 will first decrease and then increase with the increase of forest land LSI. According to the fitting calculation of the polynomial function, the coordinate value of the polynomial vertex of the forest landscape shape index is [18.02, 37.94]. Therefore, when the landscape shape index of the forest is 18.02, the minimum concentration of PM2.5 reaches 37.94 μg/m3. Therefore, this study found that the LSI value of 18.02 is the threshold of forest land (Figure 1).

2.2.2. Threshold Effect of the Impact of Landscape Pattern Indices on NO2

From Table 1, it was clear that the PD of forestland exerted a significant effect on NO2 concentration (p < 0.1); therefore, the study determined the range of both thresholds by analyzing the trend of the effect of NO2 concentration changes. From the equation curve in Figure 2, it can be observed that the NO2 concentration exhibited a decreasing trend, and then started increasing with the increase in the PD of the forestland. From the function calculation, the value of the polynomial vertex coordinates of the PD of forestland was [0.072, 29.161]; hence, when the forestland patch density was 0.072, the NO2 concentration reached the minimum of 29.161 μg/m3. Therefore, 0.072 was deduced to be the threshold value for the PD of forestland in this study (Figure 2).

2.2.3. Threshold Effect of the Impact of Landscape Pattern Indices on SO2

PD and PLAND of the forest, PD, and LSI of grassland, and PLAND of farmland will affect the concentration of SO2 pollutants, significantly. Through linear analysis of SO2 concentration change influence trend, the threshold results of 6 indicators are determined as follows.
For forest land, the concentration of SO2 will increase first and then decrease in the forest. Through polynomial fitting, it is found that when the PLAND of forest value is 50%, the SO2 concentration reaches the maximum value of 8.7 μg/m3. When the PD value of forest land is 0.038, the SO2 concentration reaches the maximum value of 8.72 μg/m3. In consideration of air quality and public health of residents, the optimal range of PD and PLAND of forest land is greater than or less than 50% and 0.038.
For grassland, the concentration of SO2 will increase first and then decrease with the increase of PD and PLAND of grassland. When the PLAND value of grassland is 3.33%, the SO2 concentration reaches the maximum value of 8.72 μg/m3. When the PD value of grassland is 0.121, the SO2 concentration reaches the maximum value of 9.545 μg/m3. Also, the SO2 concentration will decrease first and then increase with the increase of LSI of grassland, when the grassland LSI value is 14.13, the SO2 concentration reaches the minimum value of 7.84 μg/m3.
For farmland, the concentration of SO2 will increase first and then decrease with the increase of PLAND of farmland. When the PLAND value of farmland is 32.56, the PLAND value of 8.78 μg/m3. Therefore, this study found that 32.56 is the threshold of the PLAND value of farmland (Figure 3).

3. Discussions

3.1. Impact Mechanism of Landscape Pattern Index on Air Pollutant Concentration

3.1.1. Impact Mechanism of Landscape Pattern Index on PM2.5

The results of spatial regression analysis showed that the LSI of forest land had a significant negative effect on the concentration of PM2.5, which was consistent with the results obtained by previous studies and indicated that increasing the contact area between the edges of the green space patches and surrounding urban areas at large spatial scales significantly reduces PM2.5 concentrations [18]. PM2.5 primarily originates from industrial emissions, traffic emissions, and the burning of biomass, and is mainly present on roads, factories, and surrounding farmland in built-up areas [27]. When the LSI of forestland was elevated, the complexity of its patch edge also increased, and the contact area between the forestland patch edge and urban construction land was subsequently elevated. As the contact area between vegetation and PM2.5 particulate matter increased, the effect of dust reduction via vegetation leaf surface villi retardation, stem adsorption, and stomatal absorption by plants was enhanced, while the concentration of PM2.5 pollutants reduced [28,29].

3.1.2. Impact Mechanism of Landscape Pattern Index on NO2

The results of spatial regression analysis showed that the PD of forest land had a significant negative effect on the concentration of NO2, which indicated that increasing the contact area of forestland with NO2 pollutants was beneficial to the reduction of NO2 pollutant concentration.
The main sources of NO2 in the atmosphere are industrial and motor vehicle exhaust emissions [30]. Therefore, NO2 air pollution mainly exists in urban industrial areas with large traffic emissions. Studies have demonstrated that plants can absorb NO2 through their leaves [31]. The forestland patches outside the urban constructed land have a large area; however, NO2 was mainly generated in the built-up area and blocked by high-density buildings; hence, the contact area between urban forestland and NO2 pollutant gas was small, and increasing the forestland area exhibited no significant effect on the weakening effect of the overall urban NO2 pollutant concentration. Increasing the density of forestland patches in the built-up area enhanced the number of public green spaces in the city, and measures were taken to insert greenery in the edges, especially the greening of industrial areas, and street greenery was appropriately sufficient in moving plant leaves to come into contact with more polluting gases while playing a more significant role in the improvement of urban NO2 pollutant gases.

3.1.3. Impact Mechanism of Landscape Pattern Index on SO2

The results of spatial regression analysis showed that the PLAND of forest land, PD of grassland, and PLAND of farmland had a significant negative effect on the concentration of SO2. And the PD of forest land, PLAND of grassland, and LSI of grassland had a significant positive effect on the concentration of SO2.
The production of SO2 in urban air pollution mainly originates from coal combustion and industrial emissions, and studies have demonstrated that SO2 concentrations were lower in areas with rich vegetation under different spaces in cities [32]; hence, the increase in the area share of forestland patches effectively enhanced the efficiency of SO2 absorption by urban plants and reduced the SO2 air content in cities. Simultaneously, as a gas, SO2 was more susceptible to wind speed, temperature, and humidity; hence, the single vertical structure of grass greenery was more conducive to the circulation of urban winds and the transportation of SO2 pollutants from the inner city to the outer city, to achieve a lower average SO2 concentration [33]. The diffusion of SO2 is mainly influenced by wind speed and direction, and in areas with high building density due to the blockage of tall buildings resulting in low internal wind speed, SO2 concentration was not easily diffused. Therefore, with the increase in the proportion of farmland area, the degree of SO2 diffusion was accelerated and the concentration was reduced [34]. Meanwhile, cities with relatively large areas of farmland generally have industries positioned as agricultural cities, such as Henan Province, which is a largely agricultural province in China, including Xinyang and Nanyang, and these areas have low average SO2 concentrations.
As the density of forestland patches increased, SO2 concentration increased, and related studies have demonstrated that the larger the average area of green space patches in urban landscapes, the lower the fragmentation index, and the greater the role of green spaces in air pollution purification. In addition, when the degree of fragmentation of forestland was more severe, it could not function as an effective urban green heart, and the fragility of the landscape structure reduced the absorption capacity of SO2 [35]. As the ratio of the area occupied by grass increased, SO2 concentration increased. It was speculated that the reason for this result might be that the decrease in the number of plant leaves (the main organ of SO2 absorption by plants) in grass patches resulted in the subsequent decrease in the efficiency of SO2 absorption by green spaces [36]. Therefore, it was inferred that in addition to improving the air purification capacity, increasing the density of grass patches with small areas and low edge complexity also reduces economic costs.

3.2. Threshold Mechanism of Landscape Pattern Index on Air Pollutant Concentration

3.2.1. Threshold Effect of PM2.5

The LSI of forest land is related to the edge complexity and patch density of green space patches. The increase in LSI value can enhance the energy flow and exchange between green patches and surrounding patches and create more interaction opportunities for source and sink landscapes, so as to absorb and settle more pollution particles and reduce the concentration of PM2.5 [37]. However, as the LSI of forest land continues to increase, the patch density then increases the degree of green space fragmentation increases, and the abatement between forest land edges and air pollutants cannot offset the increasing amount of air pollution due to urban green space fragmentation, and PM2.5 concentrations then continue to rise to show a trend of first decreasing and then increasing [38].

3.2.2. Threshold Effect of NO2

As the PD of forest land keeps increasing, NO2 concentration shows a trend of first decreasing and then increasing. When the density of forested land patches increases continuously from 0, the intra-urban green space gradually transitions from single-core large green space to multi-core green space. With the connection of streets, rivers, and other green channels, the multi-core urban green space plays a greater ecosystem service function, more contact area of air pollutants, and higher efficiency of material exchange and energy circulation, which is conducive to the reduction of NO2 pollutant concentration [39]. However, with the increasing density of forest land patches, the urban green space system transitions from multi-core green space to green space patch fragmentation, and the concentration of NO2 pollutants then increases [40].

3.2.3. Threshold Effect of SO2

Except for the LSI of grassland, with the increasing values of PLAND and PD of forest land, PLAND and PD of grassland, and PLAND of agricultural land, the trend between the landscape pattern index and SO2 concentrations of different green space types was first increasing and then decreasing.
Green space is a sink landscape for mitigating urban SO2 pollution, which can adsorb and absorb and deter SO2 through plant leaves and branches, while urban construction land is a source landscape for SO2 pollutants [41]. The decrease in the area share of forest land patches is often associated with cities with high levels of development. According to relevant studies showing more, developed cities with relatively well-developed environmental measures and more rational urban master plans, as well as industrial structures with the upgrading and transformation of heavy industries to light industries and high-tech industries, tend to have lower pollution levels. The opposite is often true in medium-sized developing cities. Such cities are still at the stage of economic development, with more types of industries and insufficient attention to the environment, resulting in the reduction of SO2 air pollution by vegetation in such cities being much lower than the SO2 emissions, and therefore the pollution concentration is higher. With the increasing area of forest land emissions and reductions gradually reaching a balance or even vegetation reduction exceeding the local pollution emissions air pollution levels are reduced again [42,43].

3.3. Implications for Urban Planning and Management Policies

In the context of optimizing and renewing the landscape patterns of urban green spaces, several studies have demonstrated that regulating land use patterns by carefully planning the morphology and layout of green space networks effectively improves air quality and enhances the public health of residents. Based on the aforementioned findings, we made the following recommendations for the improvement of air pollution in cities with subtropical monsoon climates.
Forestland, grassland, and farmland can inhibit the concentration of PM2.5, NO2, and SO2. The grassland area in the urban built-up area should be properly controlled, the forest coverage should be gradually increased, the restoration and reconstruction of damaged forest land should be accelerated, the integrity and stability of the ecosystem should be improved, the urban green space landscape pattern should be reasonably controlled through scientific basis, and the production and living ecological space layout should be coordinated.
For cities with NO2 and SO2 as the main pollutants, the pollutants can be reduced systematically by reasonably arranging the land use pattern of forest land, grassland, farmland, and other green areas, and cooperating with each other. When the PD value of forest land is about 0.072, the overall value of urban NO2 pollutant concentration reaches the optimum. When the land and PD values of forest land are away from 50% and 0.038, the land and PD values of grassland are away from 3.33% and 0.121, and the LSI of grassland reaches 14.13, the urban SO2 pollutant concentration reduction effect is the best. Therefore, reasonable urban green space planning should be carried out to optimize the density of green space patches and evenly distribute the types of green space patches, to achieve a close relationship between various land uses and better alleviate air pollution. Through the way of inserting green in the gap, the coverage of urban green space can be improved economically and efficiently, and the contact area between green space and air pollutants can be increased to better improve the air quality level.
PM2.5 mainly comes from industrial waste gas emission, traffic emission, and biomass combustion, and mainly exists in roads and factories in built-up areas. China’s subtropical monsoon climate zone is in the stage of rapid urban development, and it is difficult to reduce traffic emissions and industrial emissions. However, PM2.5 pollution can be improved by optimizing the landscape pattern of urban green space. When the LSI of forest land reaches 18.02, the overall value of PM2.5 pollution concentration is the best. It is suggested to reasonably arrange the urban green space area. For the streets and factories with large traffic and industrial emissions, we will focus on improving the greening level of the streets, appropriately increasing the contact area between the street green space and PM2.5, improving the vertical structure of the urban green belt, block, and absorb PM2.5 pollutants to the greatest extent, and reduce the transmission route. In addition, the treatment efficiency of polluted gas in the factory shall be strictly controlled to reduce air pollution.
The comprehensive management of unused land and inefficient land should be promoted to form a reasonable and efficient urban green landscape pattern. At the same time, it is suggested to promote the development of residents’ lifestyle and consumption concepts towards green, healthy, and low-carbon. The scientific achievements in the prevention and control of PM2.5, NO2, and SO2 pollution should be actively disseminated to the residents so that the relevant departments can implement the relevant pollution control measures more smoothly, and the residents can support the air pollution control.
In addition, the results of the spatial regression model in this study also show that air pollution has obvious spatial effects, so it is not feasible to carry out internal air pollution control for a single city. From the perspective of urban agglomerations, we should coordinate and plan the air pollution control policies among cities, jointly improve the regional air quality level and reduce the threat of air pollution to residents’ public health. At the same time, our results show that the careful planning of urban green space landscape patterns has brought some positive benefits to air pollution, but compared with traffic emissions, industrial emissions, and biomass combustion, the role of pollutant gas emissions is still limited. Therefore, improving the rationality of urban green space landscape patterns is on the one hand. On the other hand, we should also attach great importance to energy efficiency and traffic management, further promote the reform of industrial structure, and eliminate backward industries with high pollution and high consumption. Finally, it aims to increase the proportion of high-tech industries, promote China’s transformation from incremental expansion to stock revitalization, improve development quality and resource utilization efficiency, and reduce air pollutant emissions from the source.

3.4. Research Innovations and Limitations

As early as the end of the last century Wickham et al. [44] proposed an integrated assessment of the environmental condition of a large region, by combining data on land cover, population, roads, rivers, air pollution, and topography. In 2009, Rafiee et al. [45] applied a combination of remote sensing image classification, landscape indicator assessment, and vegetation indices, to explore changes in urban landscape patterns and provide an assessment of changes and trends in urban living environments. Recently, there has been an increasing interest in applying remote sensing images to investigate the role of changes in urban land use patterns on air pollution concentrations. In 2015, Wu et al. [46] adopted landscape indicators such as PLAND, PD, ED, SHEI, and CONTAG to explore the effect of urban landscape patterns on PM2.5 pollution and determined that vegetation and water bodies were significant landscape components that reduced PM2.5 concentrations. In 2016, Xu et al. [47] explored the quantitative relationship between land use and air quality (SO2, NO2, and PM10) through binary correlation analysis, and the results indicated that land use had a significant effect on air quality. For each standard deviation increase in construction land, NO2 concentrations increased by 2%. For one standard deviation increase in water bodies, SO2 or PM10 concentrations decreased by 3–6%. Although this study quantified the role of land use type on air pollution levels, it did not explore the quantitative relationship between land use effects on PM2.5. In 2020, Li et al. [48] studied the non-linear effects of land use distribution on PM2.5, using the boosted regression tree method to capture the effects of land use scale on PM2.5 in different seasons. It was inferred that when the grassland and forestland areas were below 8% and 20%, respectively, the air quality improved significantly with the increase in grassland and forestland areas. When the distribution of construction land was greater than approximately 10%, PM2.5 pollution increased significantly with the increase in the construction land area. This study was the first to identify the threshold for the effect of land use type on PM2.5 concentrations. However, the role of the effect on the three air pollutants such as PM10, NO2, and SO2 is yet to be addressed. Most of the aforementioned studies adopted linear regression models to explore the role of land use types on air pollutants, which could not comprehensively incorporate the spatial effects between cities into the study, and also failed to demonstrate how the effects of landscape patterns on multiple air pollutants varied with changes in landscape pattern indices. On this premise, this study employed a spatial regression model approach to regress air pollutant concentrations and landscape pattern indices in 37 subtropical monsoon climate garden cities of China in 2019, to explore the core landscape pattern indices that exhibited significant effects on air pollution concentrations and the optimal core indicator thresholds for mitigating urban air pollution, to quantify the spatial effects of the landscape pattern indices on air pollutants, using urban spatial effects as a starting point, and then a general framework that differed from existing studies was proposed. In addition to focusing on quantifying the negative or positive effects of landscape pattern indices on air pollutant concentrations, to a certain extent, this study also reflected the appropriate threshold values of landscape pattern indices for reducing air pollution concentrations, which provides quantitative reference and technical support for urban planning in a more targeted manner.
However, in general, this study had some limitations. Firstly, the spatial resolution of satellite remote sensing images in this study is low, and there may be a certain misclassification of patch types, which may cause deviation to the impact indicators of air pollution. At the same time, because the research object is a large city with a subtropical monsoon climate, it is not universal for other climate belt cities. In addition, the image data used in our study provide a reference for the overall land use pattern of the whole city. The conclusion is that it provides a reference for the overall planning of urban green space landscape patterns from a macro perspective. The role of small-scale urban green space in reducing pollutants needs to be further discussed. At the same time, in order to ensure the accuracy of the urban green space landscape pattern indicators, the image data selected in this study are all summer image data, but the air pollution data is the annual average value of pollutants. Due to the seasonal changes, the green space coverage of a few cities in winter is reduced, so the impact analysis of this study on air pollutants is still insufficient. Therefore, based on the above shortcomings, in future research, we will expand the sample number of research objects, improve the accuracy of remote sensing data, and analyze the urban landscape pattern and air pollutant concentration in different seasons and climate zones respectively, so as to reduce variables, improve universality and research accuracy, and put forward more targeted urban green space landscape pattern planning strategies.

4. Materials and Methods

In this study, firstly, we take 37 garden cities with subtropical monsoon climate as the research unit and take the annual average concentration of PM2.5, PM10, NO2, and SO2 and the landscape pattern index of these 37 cities in 2019 as the dependent variable and independent variable, respectively, to carry out the spatial regression model. Based on the output of the model, the impact mechanism between landscape pattern index and PM2.5, PM10, NO2, and SO2 pollution was explored. Finally, a new threshold effect estimation method based on a polynomial model is designed to explore the threshold effect of green space landscape patterns on air quality (Figure 4).

4.1. Study Region

The subtropical monsoon climate zone is one of the four major climate zones in China, covering approximately a quarter of China’s land area. It supports almost 600 million people in China, is the most densely populated region in China, and mainly includes the Yangtze River Delta region, Pearl River Delta region, and middle reaches of the Yangtze River, including other areas of China’s economic center, where the level of economic development and urban landscaping are at a relatively high level in the country. With the advancement of urbanization, the original urban ecological network structure has been destroyed, owing to the drastic changes in urban land use patterns caused by the population concentration and rapid industrial development; in addition, the subtropical monsoon climate region has also become the region most threatened by air pollution in China, which exposes the public health of residents to tremendous pressure [49,50]. Considering the substantial influence of natural climate on urban air quality, to improve the accuracy of this study, 37 “national garden cities” with a population size of more than 4 million people and complete air quality data for 2019 were selected from the subtropical monsoon climate zone of China as the research objects of this study. According to the “National Garden City Standards” adopted by the Ministry of Housing and Urban-Rural Development of the People’s Republic of China, garden cities are cities with balanced distribution, practical structure, perfect functions, beautiful landscapes, fresh and comfortable human living, and ecological environments, safe, and pleasant. In addition, they were adopted as examples to investigate the effect of landscape patterns of urban green spaces on air pollution levels (see Figure 5, Appendix A).

4.2. Data Sources

Air quality index data were obtained from the “2019 National Air Quality Monthly Report” published by the Ministry of Ecology and Environment of the People’s Republic of China (http://www.mee.gov.cn/ (accessed on 13 September 2020)). In addition, the specific distributions of the air pollution index data (PM2.5, PM10, SO2, and NO2) for the 37 cities are presented in Figure 6, while the population data were obtained from the “2020 City Yearbook” for each city provided by the National Bureau of Statistics of China (http://www.stats.gov.cn/ (accessed on 1 September 2020)), and the number of permanent residents in 2019 was selected as the population index (Figure 6). The 37 urban land use classification maps in this study were obtained from the MCD12Q1.006 land use type product of the 2019 MODIS/Terra, provided by the official USGS online platform with a resolution of 500 m (https://www.usgs.gov/ (accessed on 20 September 2020)). The urban land use map was converted to a tiff format by ArcGIS10.3 and then imported into FRAGSTATS for calculations, to obtain the data for patch-level landscape pattern indices (Appendix B).

4.3. Landscape Pattern Indices

To some extent, the landscape pattern indices capture the intrinsic spatial structure of the environment and enhance the interpretation of spatial patterns and characteristics of the landscape and are now widely used to measure landscape patterns [51]. Based on the results obtained from previous related studies and validation experiments, we selected three landscape indicators to measure the urban landscape pattern of the study region, these indexes are patch proportion in landscape area (PLAND), patch density (PD), and landscape shape index (LSI). The reason why these three indicators are selected is that they are not only widely used to describe the fragmentation, irregularity, and complexity of urban landscape patterns, but also have a more direct macro control effect on the adjustment and optimization of urban green space landscape pattern in the later stage, which is conducive to clarifying and standardizing the follow-up policymaking. Patch proportion of landscape area (PLAND) is the proportion of different landscape patch types in the overall land area percentage measurement, which can help us judge the proportion of this patch type in the overall spatial pattern The patch density (PD) represents the number of patches of a certain type within 100 hectares and can reflect the spatial pattern of the landscape. Its value has a positive correlation with the fragmentation of patch types. The greater the density index, the higher the fragmentation of the patch type. The landscape shape index (LSI) is a robust index used to describe the complexity of urban morphology through the ratio of urban patch perimeter. The larger the LSI value, the more irregular the patch shape, the higher the landscape complexity, and the lower the stability. The above indicators are calculated using FRAGSTATS 4.2 software (Table 2).

4.4. Spatial Regression Modeling Methods

In this study, to investigate the quantitative influence of the landscape patterns of urban green spaces on air pollutants, the PLAND, PD, and LSI of three land use types (forestland, grassland, and farmland) were selected as the independent variables of air pollutant concentration that influence factors for spatial regression analysis.
First, in the selection of the model, after the regression analysis. The calculation of Moran’s I index for air pollutants using the GeoDa software revealed that the Moran’s I index was 0.57, and the model exhibited a significant spatial correlation. This indicated that air pollutants exhibit spatial interactions, and their effects could spread through adjacent regions, which confirmed that local and regional landscape pattern changes directly or indirectly affected their surrounding air pollutant concentrations. Therefore, the spatial regression model can better reveal the relationship between air pollutant concentration and landscape green space pattern.
Secondly, the regression results of the spatial autoregression model (SAR) and spatial error model (SEM) were compared in Table 3. The evaluation indicators for the goodness of fit of the spatial regression model included the coefficient of determination R2, the natural logarithm of the likelihood function (log-likelihood, logL), Akaike information criterion (AIC), and Schwarz Criterion (SC). Among them, the value range of R2 is (0,1), and the closer R2 is to 1, the better the regression fit of the model. In addition, the higher the logL value, and the smaller the AIC and SC values, the better the regression effect of the spatial regression model. For the selection of the optimal spatial regression model, the R2 and logL values of the SEM model were higher than those of the SAR model in the regression model information of the four pollutants, while the AIC and SC values of the SEM model were smaller for the remaining three pollutants, except for PM10 of the SEM model, which was 0.04 greater than the SAR model, thus implying that the regression results of the SEM model simulating air pollutants were better. thus, implying that the regression results of the SEM model simulating air pollutants were better. After comprehensive analysis and comparison, the SEM was selected for the spatial regression of the four atmospheric pollutants in this study. Therefore, this research established an SEM using the air pollutant concentrations of 37 Chinese cities as the dependent variables and the landscape pattern indices of the green spaces of each city as the independent variables [52,53].
Y = β1X1 + β2X2 + β3X3 + …… + β9X9 + ε
where Y denotes the dependent variable, i.e., pollutant concentration; Wy is the spatial weight matrix; ρ denotes the spatial regression coefficient of the spatial weight matrix WY, which is adapted to indicate the spatial interaction of air pollution; X1–X9 represents the impact factors of landscape pattern indexes (PLAND, PD, LSI) for three types of urban green space: forest land, grassland, and farmland, respectively (The independent variables); β1–β9 denote the regression coefficients of independent factors such as landscape pattern indices; ε denotes the random error. This was ultimately utilized to establish the SEM spatial regression model between pollutant concentrations and landscape patterns in 2019.

5. Conclusions

Urban green space is an effective tool to improve air quality. Quantifying the relationship between urban green space landscape patterns and air pollution concentration is of great significance to promote the public health level and sustainable development of residents in high-density population areas. This study takes 37 garden cities with subtropical monsoon climate in China as the research object, selects the urban air quality monitoring data and land use remote sensing data in 2019, carries out regression analysis on urban air pollutant concentration and landscape pattern index through spatial regression model method, and explores the relationship between landscape pattern index and air pollutant concentration, According to the above regression analysis results, the landscape pattern index threshold with significant correlation with air pollutant concentration was explored. The specific conclusions are as follows:
(1)
The landscape pattern of urban green space was significantly correlated with the concentrations of PM2.5, NO2, and SO2 pollutants in the air, while the concentrations of PM10 pollutants were not significantly affected by the green space pattern.
(2)
Among them, the patch shape index (LSI), patch density (PD), and patch proportion in landscape area (PLAND) of forest land can affect the concentration of PM2.5, NO2 and SO2, respectively. The PLAND, PD, and LSI of grassland and farmland can also have an additional impact on the concentration of SO2 pollutants.
(3)
The study also found that there was a significant threshold effect on the impact mechanism of urban green space landscape pattern indicators (LSI, PD, PLAND) on the concentrations of PM2.5, NO2, and SO2 air pollutants. When the PD value of forest land is about 0.072, the overall value of urban NO2 pollutant concentration reaches the optimum. When the land and PD values of forest land are away from 50% and 0.038, the land and PD values of grassland are away from 3.33% and 0.121, and the LSI of grassland reaches 14.13, the urban SO2 pollutant concentration reduction effect is the best.
The results of this study not only clarify the impact mechanism of the landscape pattern of urban green space on air quality but also propose a polynomial-based threshold effect estimation method, which provides quantitative reference and scientific basis for the optimization and updating of urban green space landscape patterns to promote public health.

Author Contributions

C.W.: Conceptualization, Methodology, Software; M.G.: Writing—Reviewing and Editing; J.J.: Data curation, Writing—Original draft preparation; Y.Y.: Visualization, Investigation; Y.R.: Software, Methodology; Y.W.: Validation; J.C.: Supervision. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Program of China (Funder: Jun Jin; Funding No.: 2020YFD1100405); the National Natural Science Foundation of China (Funder: Jiajie Cao; Funding No.: 32071832); the Humanity and Social Sciences Youth Foundation of the Ministry of Education of China (Funder: Chengkang Wang; Funding No.: 18YJCZH167) and the Natural Science Foundation of Jiangsu Province (Funder: Chengkang Wang; Funding No.: BK20190751).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The effect of landscape patterns of urban green spaces on air pollution levels.
Table A1. The effect of landscape patterns of urban green spaces on air pollution levels.
Greening IndicatorsAir Quality IndicatorsSocial Indicators
Serial NumberNameGreening RatePM2.5 (μg/m3)PM10 (μg/m3)NO2 (μg/m3)SO2 (μg/m3)Population Ten Thousand
1Guiyang38.2026.6647.0021.169.66497.14
2Taizhou38.2026.6648.6721.664.41614.00
3Nanning34.0830.3352.7529.339.08734.48
4Ningbo37.8728.5047.3335.507.91854.20
5Wenzhou34.4928.0052.9233.837.58930.00
6Jiangmen41.6226.6673.8332.087.00463.03
7Zhaoqing37.9931.7548.0033.509.33415.17
8Jinhua38.1532.4154.434.417.16562.40
9Yichun44.6535.4259.4123.7513.33558.26
10Yancheng39.1039.5066.0824.004.75720.89
11Jiaxing36.1525.3356.2532.756.66480.00
12Nantong40.0036.6655.0031.5810.50731.80
13Dongguan37.5031.9148.4236.589.66846.45
14Shaoxing37.1938.0061.5831.586.91505.70
15Foshan42.5029.7556.1641.089.08815.86
16Nanchang38.5135.569.3333.508.58560.05
17Guangzhou39.9130.0052.6644.836.831530.59
18Changde39.3147.6660.0022.668.00577.13
19Xiaoganunavailable43.3372.2521.257.08492.10
20Huanggangunavailable40.2573.0824.919.66633.30
21Jiujiang44.7545.9263.5029.9110.91492.03
22Yueyang39.4143.567.8326.758.75577.13
23Taizhou38.643.9267.0828.007.50463.61
24Changsha35.2547.0857.4133.087.08839.45
25Xinyang38.0048.2576.0823.916.33646.00
26Hangzhou37.2337.9166.4141.256.751036.00
27Zhuzhou34.9147.2565.2533.3310.75402.85
28Yichang36.1052.2572.8329.167.08413.79
29Wuxi39.8439.0068.7539.918.16659.15
30Hefei39.7543.5865.7538.086.16818.90
31Nannjing41.0039.7569.1641.759.83850.00
32Jingzhou33.1046.5083.0032.089.25557.01
33Chuzhou42.8648.2572.0835.089.66414.70
34Changzhou39.2046.8371.0041.0010.41473.60
35WUhan34.4745.2570.7544.168.831121.20
36Xiangyang33.3460.3384.5831.5810.50568.00
37Nanyang38.1059.6693.1628.836.331003.00

Appendix B

Table A2. The data for patch-level landscape pattern indices.
Table A2. The data for patch-level landscape pattern indices.
NumberCityForest-PLForest-PDForest-LSIGrass-PLGrass-PDGrass-LSIWater-PLWater-PDWater-LSIFarm-PLFarm-PDFarm-LSIConstruction-PLConstruction-PDConstruction-LSI
1Guiyang67.910.0214.360.150.005.620.450.003.8627.900.0318.403.570.007.29
2Taizhou78.940.0113.000.800.0111.812.540.0110.719.300.0423.858.360.0213.31
3Nanning57.400.0223.861.390.0224.821.110.0111.8837.160.0227.252.940.0115.06
4Ningbo66.700.0214.071.420.0214.653.210.0112.4511.360.03820.9517.220.0212.41
5Wenzhou85.190.018.591.000.0114.131.600.0110.374.620.0218.277.500.0110.79
6Jiangmen63.610.0219.432.010.0318.945.430.021422.740.0222.096.040.0110.31
7Zhaoqing92.560.007.931.500.0221.042.170.0112.421.980.0113.481.770.0111.32
8Jinhua83.640.0112.860.180.017.350.320.007.188.010.0321.797.830.0217.34
9Yichun72.370.0218.351.930.0224.090.830.0112.2423.260.0224.021.600.0111.61
10Yancheng5.120.0327.571.860.0220.315.840.0110.8083.390.008.803.530.0111.79
11Jiaxing16.130.0723.681.310.016.854.450.019.2457.230.0111.8920.870.0312.25
12Nantong10.850.0424.880.310.017.951.830.005.6478.030.019.158.920.0110.51
13Dongguan26.500.0411.615.430.0714.274.420.039.570.380.014.4562.620.016.71
14Shaoxing76.630.0112.291.380.0111.920.890.018.4810.430.0318.5210.660.0212.65
15Foshan44.000.0213.592.680.0414.498.070.0314.020.830.015.3842.850.029.32
16Nanchang21.780.0522.994.140.0522.7217.850.0211.9549.340.0116.506.790.017.06
17Guangzhou64.120.0215.593.010.0315.852.080.0213.697.110.0317.6523.470.0111.66
18Changde61.780.0222.290.920.0118.444.000.0117.9732.210.0124.541.080.008.46
19Xiaogan22.240.0422.674.320.0422.053.580.018.9767.430.0114.382.430.018.29
20Huanggang22.240.0422.671.560.0221.734.920.0116.6234.450.0228.232.350.0111.59
21Jiujiang75.830.0114.011.400.0220.239.900.0011.3811.130.0221.651.660.0110.75
22Yueyang56.730.0319.112.170.0219.6311.530.0113.2528.050.01723.501.480.007.71
23Taizhou21.620.0624.276.090.0521.137.020.0110.5554.430.01813.4010.760.01510.47
24Changsha82.300.0112.730.710.0111.180.630.007.5811.050.0323.635.290.018.13
25Xinyang28.200.0213.880.740.0116.471.540.0110.9468.010.019.651.500.019.47
26Hangzhou83.470.019.680.360.0112.213.810.0011.941.970.0116.2510.360.0111.57
27Zhuzhou90.530.008.780.390.019.680.660.008.506.510.0218.691.900.007.06
28Yichang87.350.017.360.570.0115.341.210.0012.809.810.0012.801.030.008.33
29Wuxi26.430.0521.932.470.0213.2619.740.015.3318.940.0316.5432.390.0210.37
30Hefei9.520.0321.981.120.0215.708.410.015.9174.660.0110.696.280.017.56
31Nannjing31.670.0527.871.790.0314.85310.580.018.3638.990.0321.3516.850.029.40
32Jingzhou17.560.0533.331.810.0321.8314.920.0116.3663.470.0116.362.210.019.52
33Chuzhou7.350.0319.661.770.0321.674.800.0113.7183.470.0011.462.590.0110.91
34Changzhou34.230.0523.642.370.0211.879.090.016.1129.420.0318.1224.880.028.55
35Wuhan35.310.0432.204.430.0524.2515.110.0216.6033.690.03224.6211.290.0110.32
36Xiangyang52.550.0212.001.450.0221.201.230.0112.9641.860.0114.722.900.0117.15
37Nanyang33.640.0112.492.130.0225.681.520.006.3459.120.0115.143.570.02528.08

References

  1. Gianfredi, V.; Buffoli, M.; Rebecchi, A.; Croci, R.; Oradini-Alacreu, A.; Stirparo, G.; Marino, A.; Odone, A.; Capolongo, S.; Signorelli, C. Association between Urban Greenspace and Health: A Systematic Review of Literature. Int. J. Environ. Res. Public Health 2021, 18, 5137. [Google Scholar] [CrossRef] [PubMed]
  2. Reyes-Riveros, R.; Altamirano, A.; De La Barrera, F.; Rozas-Vásquez, D.; Vieli, L.; Meli, P. Linking public urban green spaces and human well-being: A systematic review. Urban For. Urban Green. 2021, 61, 127105. [Google Scholar] [CrossRef]
  3. Barrett, E.; Apodaca, K. Response to “Alcohol and COVID-19: How Do We Respond to This Growing Public Health Crisis?”. J. Gen. Intern. Med. 2021, 36, 2476. [Google Scholar] [CrossRef] [PubMed]
  4. Jaafari, S.; Shabani, A.A.; Moeinaddini, M.; Danehkar, A.; Sakieh, Y. Applying landscape metrics and structural equation modeling to predict the effect of urban green space on air pollution and respiratory mortality in Tehran. Environ. Monit. Assess. 2020, 192, 412. [Google Scholar] [CrossRef] [PubMed]
  5. Kim, S.Y.; Kim, S.H.; Wee, J.H.; Min, C.; Han, S.M.; Kim, S.; Choi, H.G. Short and long term exposure to air pollution increases the risk of ischemic heart disease. Sci. Rep. 2021, 11, 51085. [Google Scholar]
  6. Luo, J.; Niu, Y.; Zhang, Y.; Zhang, M.; Tian, Y.; Zhou, X. Dynamic analysis of retention PM2.5 by plant leaves in rainfall weather conditions of six tree species. Energy Sources Part A Recovery Util. Environ. Eff. 2020, 42, 1014–1025. [Google Scholar] [CrossRef]
  7. Cai, L.; Zhuang, M.; Ren, Y. A landscape scale study in Southeast China investigating the effects of varied green space types on atmospheric PM2.5 in mid-winter. Urban For. Urban Green. 2020, 49, 126607. [Google Scholar] [CrossRef]
  8. Latha, K.M.; Highwood, E. Studies on particulate matter (PM10) and its precursors over urban environment of Reading, UK. J. Quant. Spectrosc. Radiat. Transf. 2006, 101, 367–379. [Google Scholar] [CrossRef]
  9. Beckett, K.P.; Freer-Smith, P.H.; Taylor, G. Urban woodlands: Their role in reducing the effects of particulate pollution. Environ. Pollut. 1998, 99, 347–360. [Google Scholar] [CrossRef]
  10. Sheng, Q.; Zhu, Z. Physiological Response of European Hornbeam Leaves to Nitrogen Dioxide Stress and Self-recovery. J. Am. Soc. Hortic. Sci. 2019, 144, 23–30. [Google Scholar] [CrossRef] [Green Version]
  11. Thomas, M.D.; Hill, G.R., Jr. Absorption of sulphur dioxide by alfalfa and its relation to leaf injury. Plant Physiol. 1935, 10, 291. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, H.; Sun, Y.-Q.; Zhang, M.; Liu, W.J.; Yang, L.; Zhou, X.W. Adsorption ability of air pollutants by indigenous tree species in ta-pieh mountains. Fresenius Environ. Bull. 2019, 28, 2908–2915. [Google Scholar]
  13. Wang, S.; Li, M. Study on the Principle of Urban Open Space Ecological Planning. Chin. Landsc. Arch. 2001, 5, 33–37. [Google Scholar]
  14. Sheng, Q.; Zhang, Y.; Zhu, Z.; Li, W.; Xu, J.; Tang, R. An experimental study to quantify road greenbelts and their association with PM2.5 concentration along city main roads in Nanjing, China. Sci. Total Environ. 2019, 667, 710–717. [Google Scholar] [CrossRef] [PubMed]
  15. Fan, S.; Li, X.; Dong, L. Field assessment of the effects of land-cover type and pattern on PM10 and PM2.5 concentrations in a microscale environment. Environ. Sci. Pollut. Res. 2019, 26, 2314–2327. [Google Scholar] [CrossRef] [PubMed]
  16. Ye, L.; Fang, L.; Tan, W.; Wang, Y.; Huang, Y. Exploring the effects of landscape structure on aerosol optical depth (AOD) patterns using GIS and HJ-1B images. Environ. Sci. Process. Impacts 2016, 18, 265–276. [Google Scholar] [CrossRef]
  17. Yue, F.; Dai, F.; Guo, X. Study on Correlation Between Air Pollutants and Vegetation Coverage in Wuhan Based on Remote Sensing Inversion. Landsc. Arch. 2019, 26, 76–81. [Google Scholar]
  18. Lei, Y.; Davies, G.M.; Jin, H.; Tian, G.; Kim, G. Scale-dependent effects of urban greenspace on particulate matter air pollution. Urban For. Urban Green. 2021, 61, 127089. [Google Scholar] [CrossRef]
  19. Zhao, L.; Li, T.; Przybysz, A.; Guan, Y.; Ji, P.; Ren, B.; Zhu, C. Effect of urban lake wetlands and neighboring urban greenery on air PM10 and PM2.5 mitigation. Build. Environ. 2021, 206, 108291. [Google Scholar] [CrossRef]
  20. Alam, M.S.; McNabola, A. Exploring the modeling of spatiotemporal variations in ambient air pollution within the land use regression framework: Estimation of PM10 concentrations on a daily basis. J. Air Waste Manag. Assoc. 2015, 65, 628–640. [Google Scholar] [CrossRef]
  21. Zhang, H.; Zhao, Y. Land use regression for spatial distribution of urban particulate matter (PM10) and sulfur dioxide (SO2) in a heavily polluted city in Northeast China. Environ. Monit. Assess. 2019, 191, 712. [Google Scholar] [CrossRef] [PubMed]
  22. Lee, H.J.; Koutrakis, P. Daily ambient NO2 concentration predictions using satellite ozone monitoring instrument NO2 data and land use regression. Environ. Sci. Technol. 2014, 48, 2305–2311. [Google Scholar] [PubMed]
  23. Lu, D.; Xu, J.; Yue, W.; Mao, W.; Yang, D.; Wang, J. Response of PM2.5 pollution to land use in China. J. Clean. Prod. 2019, 244, 118741. [Google Scholar] [CrossRef]
  24. de Jalón, S.G.; Burgess, P.; Yuste, J.C.; Moreno, G.; Graves, A.; Palma, J.; Crous-Duran, J.; Kay, S.; Chiabai, A. Dry deposition of air pollutants on trees at regional scale: A case study in the Basque Country. Agric. For. Meteorol. 2019, 278, 107648. [Google Scholar] [CrossRef]
  25. Anselin, L. Local Indicators of Spatial Association—LISA. Geogr. Anal. 1995, 27, 93–115. [Google Scholar] [CrossRef]
  26. Getis, A.; Ord, J.K. The Analysis of Spatial Association by Use of Distance Statistics. Geogr. Anal. 2010, 24, 189–206. [Google Scholar] [CrossRef]
  27. Alastuey, A.; Querol, X.; Plana, F.; Viana, M.; Ruiz, C.R.; De La Campa, A.S.; De La Rosa, J.; Mantilla, E.; Dos Santos, S.G. Identification and chemical characterization of industrial particulate matter sources in southwest Spain. J. Air Waste Manag. Assoc. 2006, 56, 993–1006. [Google Scholar] [CrossRef] [Green Version]
  28. Chen, L.; Liu, C.; Zhang, L.; Zou, R.; Zhang, Z. Variation in Tree Species Ability to Capture and Retain Airborne Fine Particulate Matter (PM2.5). Sci. Rep. 2017, 7, 3206. [Google Scholar] [CrossRef] [Green Version]
  29. Liang, D.; Ma, C.; Wang, Y.-Q.; Wang, Y.-J.; Chen-Xi, Z. Quantifying PM2.5 capture capability of greening trees based on leaf factors analyzing. Environ. Sci. Pollut. Res. 2016, 23, 21176–21186. [Google Scholar] [CrossRef] [Green Version]
  30. Gasmi, K.; Aljalal, A.; Al-basheer, W.; Abdulahi, M. Analysis of NOx, NO and NO2 ambient levels as a function of meteorological parameters in Dhahran, Saudi Arabia. WIT Trans. Ecol. Environ. 2017, 211, 77–86. [Google Scholar]
  31. Morikawa, H.; Erkin, Ö.C. Basic processes in phytoremediation and some applications to air pollution control. Chemosphere 2003, 52, 1553–1558. [Google Scholar] [CrossRef]
  32. Guo, Z.; Jiang, H.; Chen, J.; Cheng, M.; Wang, B.; Jiang, Z. The relationship between atmospheric SO2 column density and land use in Zhejiang, China. In Proceedings of the 2009 Joint Urban Remote Sensing Event, Shanghai, China, 20–22 May 2009; pp. 1–6. [Google Scholar] [CrossRef]
  33. Li, C.; Xu, Y.; Liu, M.; Hu, Y.; Huang, N.; Wu, W. Modeling the Impact of Urban Three-Dimensional Expansion on Atmospheric Environmental Conditions in an Old Industrial District: A Case Study in Shenyang, China. Pol. J. Environ. Stud. 2020, 29, 3171–3181. [Google Scholar] [CrossRef]
  34. Li, Y.; Shi, J. Characteristics of SO2 Dispersion in Urban of Taiyuan‚ Shanxi. J. Shanxi Univ. 2011, 34, 153–157. [Google Scholar]
  35. Shao, T.; Zhou, Z.; Wang, P.; Tang, W.; Liu, X.; Hu, X. Relationship between urban green-land landscape patterns and air pollution in the central district of Yichang city. Ying Yong Sheng Tai Xue Bao J. Appl. Ecol. 2004, 15, 691–696. [Google Scholar]
  36. Spedding, D.J. Uptake of Sulphur Dioxide by Barley Leaves at Low Sulphur Dioxide Concentrations. Nature 1969, 224, 1229–1231. [Google Scholar] [CrossRef]
  37. Li, K.; Li, C.; Liu, M.; Hu, Y.; Wang, H.; Wu, W. Multiscale analysis of the effects of urban green infrastructure landscape patterns on PM2.5 concentrations in an area of rapid urbanization. J. Clean. Prod. 2021, 325, 129324. [Google Scholar] [CrossRef]
  38. Yang, Y. The effect of urban form on PM2.5 concentration: Evidence from china’s 340 prefecture-level cities. Remote Sens. 2021, 14, 7. [Google Scholar]
  39. Rao, M.; George, L.A.; Rosenstiel, T.N.; Shandas, V.; Dinno, A. Assessing the relationship among urban trees, nitrogen dioxide, and respiratory health. Environ. Pollut. 2014, 194, 96–104. [Google Scholar] [CrossRef]
  40. Klingberg, J.; Broberg, M.; Strandberg, B.; Thorsson, P.; Pleijel, H. Influence of urban vegetation on air pollution and noise exposure—A case study in Gothenburg, Sweden. Sci. Total Environ. 2017, 599, 1728–1739. [Google Scholar] [CrossRef]
  41. Geng, Y.C.; Li, J. The value evaluation on the forest resource of forest industry zone of heilongjiang province in acid rain prevention. In Proceedings of the Asia-Pacific Power and Energy Engineering Conference, Wuhan, China, 25–28 March 2011; IEEE Computer Society: Washington, DC, USA, 2011. [Google Scholar]
  42. He, Z.; Shi, X.; Wang, X.; Xu, Y. Urbanisation and the geographic concentration of industrial SO2 emissions in China. Urban Stud. 2017, 54, 3579–3596. [Google Scholar] [CrossRef]
  43. Wickham, J.D.; Jones, K.B.; Riitters, K.H.; O’Neill, R.V.; Tankersley, R.D.; Smith, E.R.; Chaloud, D.J. An integrated environmental assessment of the us mid—Atlantic region. Environ. Manag. 1999, 24, 553–560. [Google Scholar] [CrossRef] [PubMed]
  44. Rafiee, R.; Mahiny, A.S.; Khorasani, N. Assessment of changes in urban green spaces of Mashad city using satellite data. Int. J. Appl. Earth Obs. Geoinf. 2009, 11, 431–438. [Google Scholar] [CrossRef]
  45. Wu, J.; Xie, W.; Li, W.; Li, J. Effects of Urban Landscape Pattern on PM2.5 Pollution—A Beijing Case Study. PLoS ONE 2015, 10, e0142449. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Xu, G.; Jiao, L.; Zhao, S.; Yuan, M.; Li, X.; Han, Y.; Zhang, B.; Dong, T. Examining the Impacts of Land Use on Air Quality from a Spatio-Temporal Perspective in Wuhan, China. Atmosphere 2016, 7, 62. [Google Scholar] [CrossRef] [Green Version]
  47. Li, C.; Zhang, K.; Dai, Z.; Ma, Z.; Liu, X. Investigation of the impact of land-use distribution on PM2.5 in Weifang: Seasonal variations. Int. J. Environ. Res. Public Health 2020, 17, 5135. [Google Scholar] [CrossRef]
  48. Garcia-Menendez, F.; Monier, E.; Selin, N.E. The role of natural variability in projections of climate change impacts on U.S. ozone pollution. Geophys. Res. Lett. 2017, 44, 2911–2921. [Google Scholar] [CrossRef] [Green Version]
  49. Lu, D.; Mao, W.; Yang, D.; Zhao, J.; Xu, J. Effects of land use and landscape pattern on PM2.5 in Yangtze River Delta, China. Atmos. Pollut. Res. 2018, 9, 705–713. [Google Scholar] [CrossRef]
  50. Xiao, H.; Huang, Z.; Zhang, J.; Zhang, H.; Chen, J.; Zhang, H.; Tong, L. Identifying the impacts of climate on the regional transport of haze pollution and inter-cities correspondence within the Yangtze River Delta. Environ. Pollut. 2017, 228, 26–34. [Google Scholar] [CrossRef]
  51. Parker, D.C.; Evans, T.; Meretsky, V. Measuring Emergent Properties of Agent-Based Landcover/Landuse Models Using Spatial Metrics; Society for Computational Economics: Dallas, TX, USA, 2001. [Google Scholar]
  52. Lee, W.J.; Park, C. Prediction of apartment prices per unit in Daegu-Gyeongbuk areas by spatial regression models. J. Korean Data Inf. Sci. Soc. 2015, 26, 561–568. [Google Scholar] [CrossRef] [Green Version]
  53. Tu, M.; Liu, Z.; He, C.; Fang, Z.; Lu, W. The relationships between urban landscape patterns and fine particulate pollution in China: A multiscale investigation using a geographically weighted regression model. J. Clean. Prod. 2019, 237, 117744. [Google Scholar] [CrossRef]
Figure 1. Threshold of significant landscape pattern indices for PM2.5 pollutant.
Figure 1. Threshold of significant landscape pattern indices for PM2.5 pollutant.
Plants 11 02847 g001
Figure 2. Threshold of significant landscape pattern indices for NO2 pollutant.
Figure 2. Threshold of significant landscape pattern indices for NO2 pollutant.
Plants 11 02847 g002
Figure 3. Threshold of significant landscape pattern indices for SO2 pollutant ((af) represents the threshold effect of SO2 and forest land PLAND, forest land PD, grassland PLAND, grassland PD, grassland LSI, and farmland PLAND landscape pattern index).
Figure 3. Threshold of significant landscape pattern indices for SO2 pollutant ((af) represents the threshold effect of SO2 and forest land PLAND, forest land PD, grassland PLAND, grassland PD, grassland LSI, and farmland PLAND landscape pattern index).
Plants 11 02847 g003
Figure 4. Flow chart Schemes follow the same formatting.
Figure 4. Flow chart Schemes follow the same formatting.
Plants 11 02847 g004
Figure 5. Distribution map of the selected 37 cities.
Figure 5. Distribution map of the selected 37 cities.
Plants 11 02847 g005
Figure 6. Distribution map of PM2.5, PM10, NO2, and SO2 concentration. (a) Distribution map of PM2.5 concentration in the study city. (b) Distribution map of PM10 concentration in the study city. (c) Distribution map of NO2 concentration in the study city. (d) Distribution map of SO2 concentration in the study city.
Figure 6. Distribution map of PM2.5, PM10, NO2, and SO2 concentration. (a) Distribution map of PM2.5 concentration in the study city. (b) Distribution map of PM10 concentration in the study city. (c) Distribution map of NO2 concentration in the study city. (d) Distribution map of SO2 concentration in the study city.
Plants 11 02847 g006
Table 1. Conclusion of thresholds for landscape pattern indices.
Table 1. Conclusion of thresholds for landscape pattern indices.
VariablepThresholdpThresholdpThresholdpThreshold
PM2.5PM10NO2SO2
Forest landPLAND0.569 0.1610.774 0.09 *50
PD0.41 0.5130.091 *0.07180.003 ***0.038
LSI0.059 *18.0180.3230.204 0.493
GrasslandPLAND0.847 0.7640.816 0.031 **3.33
PD0.571 0.7510.9902 0.000 ***0.121
LSI0.182 0.4990.818 0.000 ***14.13
Farm landPLAND0.247 0.5370.592 0.001 ***32.56
PD0.652 0.9210.624 0.925
LSI0.309 0.5920.774 0.139
*** p < 0.01, ** p < 0.05, * p < 0.1. Correlations significant at the level of 0.1 are marked in bold.
Table 2. Landscape pattern indices.
Table 2. Landscape pattern indices.
Metrics (Abbreviation)Calculation FormulaDescription
Percentage of landscape
types (PLAND)
PLAND = 1 A j = 1 n a i j PLAND quantifies the proportional abundance of each patch type in the landscape (percent)
Patch density (PD)PD = 1 A j = 1 M N j PD expresses number of patches on a per unit area for considered class (number per 100 hectares)
Landscape shape index
(LSI)
LSI = 1 A j = 1 n E LSI expresses the larger LSI value is, the more complex landscape shape is.
Table 3. Information comparison of SEM and SAR models.
Table 3. Information comparison of SEM and SAR models.
ModelR2LogLAICSC
PM2.5SEM0.747−116.730265.46291.662
SAR0.716−116.772267.543295.382
PM10SEM0.744−123.79279.581305.782
SAR0.73−122.789279.577307.416
NO2SEM0.549−108.694249.39275.591
SAR0.525−111.464256.928284.767
SO2SEM0.6339−60.004152.01178.211
SAR0.41−67.446168.894196.733
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, C.; Guo, M.; Jin, J.; Yang, Y.; Ren, Y.; Wang, Y.; Cao, J. Does the Spatial Pattern of Plants and Green Space Affect Air Pollutant Concentrations? Evidence from 37 Garden Cities in China. Plants 2022, 11, 2847. https://doi.org/10.3390/plants11212847

AMA Style

Wang C, Guo M, Jin J, Yang Y, Ren Y, Wang Y, Cao J. Does the Spatial Pattern of Plants and Green Space Affect Air Pollutant Concentrations? Evidence from 37 Garden Cities in China. Plants. 2022; 11(21):2847. https://doi.org/10.3390/plants11212847

Chicago/Turabian Style

Wang, Chengkang, Mengyue Guo, Jun Jin, Yifan Yang, Yujie Ren, Yang Wang, and Jiajie Cao. 2022. "Does the Spatial Pattern of Plants and Green Space Affect Air Pollutant Concentrations? Evidence from 37 Garden Cities in China" Plants 11, no. 21: 2847. https://doi.org/10.3390/plants11212847

APA Style

Wang, C., Guo, M., Jin, J., Yang, Y., Ren, Y., Wang, Y., & Cao, J. (2022). Does the Spatial Pattern of Plants and Green Space Affect Air Pollutant Concentrations? Evidence from 37 Garden Cities in China. Plants, 11(21), 2847. https://doi.org/10.3390/plants11212847

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop