Spatiotemporal Analysis of Air Pollution and Climate Change Effects on Urban Green Spaces in Bucharest Metropolis

Being an essential issue in global climate warming, the response of urban green spaces to air pollution and climate variability because of rapid urbanization has become an increasing concern at both the local and global levels. This study explored the response of urban vegetation to air pollution and climate variability in the Bucharest metropolis in Romania from a spatiotemporal perspective during 2000–2024, with a focus on the 2020–2024 period. Through the synergy of time series in situ air pollution and climate data, and derived vegetation biophysical variables from MODIS Terra/Aqua satellite data, this study applied statistical regression, correlation, and linear trend analysis to assess linear relationships between variables and their pairwise associations. Green spaces were measured with the MODIS normalized difference vegetation index (NDVI), leaf area index (LAI), photosynthetically active radiation (FPAR), evapotranspiration (ET), and net primary production (NPP), which capture the complex characteristics of urban vegetation systems (gardens, street trees, parks, and forests), periurban forests, and agricultural areas. For both the Bucharest center (6.5 km × 6.5 km) and metropolitan (40.5 km × 40.5 km) test areas, during the five-year investigated period, this study found negative correlations of the NDVI with ground-level concentrations of particulate matter in two size fractions, PM2.5 (city center r = −0.29; p < 0.01, and metropolitan r = −0.39; p < 0.01) and PM10 (city center r = −0.58; p < 0.01, and metropolitan r = −0.56; p < 0.01), as well as between the NDVI and gaseous air pollutants (nitrogen dioxide—NO2, sulfur dioxide—SO2, and carbon monoxide—CO. Also, negative correlations between NDVI and climate parameters, air relative humidity (RH), and land surface albedo (LSA) were observed. These results show the potential of urban green to improve air quality through air pollutant deposition, retention, and alteration of vegetation health, particularly during dry seasons and hot summers. For the same period of analysis, positive correlations between the NDVI and solar surface irradiance (SI) and planetary boundary layer height (PBL) were recorded. Because of the summer season’s (June–August) increase in ground-level ozone, significant negative correlations with the NDVI (r = −0.51, p < 0.01) were found for Bucharest city center and (r = −76; p < 0.01) for the metropolitan area, which may explain the degraded or devitalized vegetation under high ozone levels. Also, during hot summer seasons in the 2020–2024 period, this research reported negative correlations between air temperature at 2 m height (TA) and the NDVI for both the Bucharest city center (r = −0.84; p < 0.01) and metropolitan scale (r = −0.90; p < 0.01), as well as negative correlations between the land surface temperature (LST) and the NDVI for Bucharest (city center r = −0.29; p< 0.01) and the metropolitan area (r = −0.68, p < 0.01). During summer seasons, positive correlations between ET and climate parameters TA (r = 0.91; p < 0.01), SI (r = 0.91; p < 0.01), relative humidity RH (r = 0.65; p < 0.01), and NDVI (r = 0.83; p < 0.01) are associated with the cooling effects of urban vegetation, showing that a higher vegetation density is associated with lower air and land surface temperatures. The negative correlation between ET and LST (r = −0.92; p < 0.01) explains the imprint of evapotranspiration in the diurnal variations of LST in contrast with TA. The decreasing trend of NPP over 24 years highlighted the feedback response of vegetation to air pollution and climate warming. For future green cities, the results of this study contribute to the development of advanced strategies for urban vegetation protection and better mitigation of air quality under an increased frequency of extreme climate events.