Descriptive statistics for all observed data are included on Table 1. According to the results, mean 19,753 people caught the disease daily and a total of 19,795 people died linked to COVID-19 during the time interval of the study. Within this period, the highest daily confirmed cases in Turkey reached 63,082 cases on 16 April 2021, with highest mortality for 394 people on 30 April 2021. Mean PM10, PM2.5, SO2, CO, NO2 and O3 were calculated as 38.31, 15.54, 11.80, 882.19, 29.47 and 52.38 µg/m3, respectively. Additionally, according to the Jarque-Bera test, apart from NO2 and O3, other datasets did not have normal distribution. For the period from 29 April-17 May when the strictest lockdown precautions were applied, a clear fall was experienced in air pollution indicators, as seen in Fig. 2. In this period, the PM10, PM2.5, SO2, CO and NO2 concentrations reduced by 6%, 15%, 14%, 27% and 11% rates, while the O3 amount increased by nearly 12% compared to the mean for the previous 17-day period. The results overlap with findings from other studies investigating air pollution in Turkey during the pandemic (Ali et al. 2021; Orak and Ozdemir 2021; Sari and Esen 2021).
Table 1
Variables | Mean | Median | Min. | Max. | Skewness | Kurtosis | Jarque-Bera |
Daily cases | 19753 | 9380 | 4418 | 63082 | 1.13 | 2.76 | 26.20* |
Daily deaths | 162.25 | 125.5 | 35 | 394 | 0.46 | 1.73 | 12.62* |
PM10 (µg/m3) | 38.31 | 37.25 | 23.35 | 66.21 | 0.87 | 4.10 | 21.45* |
PM2.5 (µg/m3) | 15.54 | 14.65 | 8.12 | 24.92 | 0.64 | 3.08 | 8.28** |
SO2 (µg/m3) | 11.80 | 11.63 | 8.54 | 17.01 | 0.68 | 3.60 | 11.27* |
CO (µg/m3) | 882.19 | 900.66 | 511.43 | 1425.24 | 0.11 | 1.90 | 6.40** |
NO2 (µg/m3) | 29.47 | 28.92 | 17.31 | 45.14 | 0.27 | 2.64 | 2.16 |
O3 (µg/m3) | 52.38 | 53.52 | 25.96 | 69.09 | -0.32 | 2.76 | 2.44 |
* and ** signifies indicates significance at 1% and 5%. |
Table 1
The results for Spearman and Kendall correlation tests used to investigate the relationship between air pollutants and COVID-19 in Turkey are presented in Table 2. According to both correlation tests, for confirmed cases and mortality, there were positive correlations with PM2.5, SO2, CO and NO2 and a negative correlation with O3. The variation in air pollution appeared to have a more significant correlation with those infected with the disease rather than reported deaths. Positive correlation coefficients were obtained between COVID-19 and air pollutants in the order SO2 > CO > PM2.5>NO2 > PM10. However, no significant correlation was encountered between PM10 and any dataset related to the pandemic (p > 0.1). Additionally, Spearman correlation coefficients were higher and had more significance compared to Kendal coefficients. Findings are similar to studies researching crowded cities like New York and those in California in the USA (Bashir et al. 2020c, b).
Table 2
Kendall and Spearman correlation coefficients.
| | PM10 | PM2.5 | SO2 | CO | NO2 | O3 |
Spearman's rho | New cases | 0.209** | 0.292* | 0.672* | 0.532* | 0.240* | -0.549* |
Daily deaths | 0.121 | 0.185** | 0.376* | 0.111 | 0.260* | -0.700* |
Kendall's tau | New cases | 0.136** | 0.179* | 0.479* | 0.381* | 0.154** | -0.364* |
Daily deaths | 0.077 | 0.109*** | 0.246* | 0.100*** | 0.170* | -0.493* |
*, ** and *** signifies shows significance levels (2-tailed) at 1%, 5% and 10%. |
Table 2
Figure 3 and Fig. 4 show the CCF results between daily COVID-19 cases and deaths with air pollution indicators. Here, positive lag [1;14] indicates the COVID pandemic is ahead, while negative lag [-1;-14] indicates air pollutants are ahead and lag [0] indicates both time series have no lag. For daily cases and PM10, PM2.5, SO2, CO, NO2 and O3, the highest values were obtained on the 11th day (0.316), day 0 (0.483), -7th day (0.751), -7th day (0.624), day 0 (0.331) and day 0 (-0.621), respectively. Additionally, as seen in the correlogram in Fig. 3, there was increasing CCF trend with PM2.5, SO2 and CO ahead and PM10, NO2 and O3 behind. When the CCF values between daily mortality and air pollutants in Fig. 4 are examined, all air pollutants were ahead of COVID-19 with an increasing profile dominant. The highest CCF values were calculated for PM10, PM2.5, SO2, CO, NO2 and O3 on the 2nd day (0.295), -14th day (0.42), -14th day (0.715), -14th day (0.596), -2nd day (0.331) and − 2nd day (-0.689). In other words, the variation in PM2.5, SO2, CO, NO2 and O3, especially, has more significant effects in the progressive periods of the pandemic. Research investigating the effect of meteorological factors revealed evidence that air pressure, wind speed and precipitation in India had a delaying effect on the spread of the pandemic, similar to air pollution (Kulkarni et al. 2021).
Figure 3
Figure 3. Cross-correlation factors of confirmed cases with air pollutants. All cross-correlograms were assessed for lags of − 14 to 14 days.
Figure 4
Figure 4. Cross-correlation factors of daily deaths with air pollutants. All cross-correlograms were assessed for lags of − 14 to 14 days.
The graphical presentation of the wavelet transform coherence analysis for COVID-19 with air pollutants is given in Fig. 5 and Fig. 6. The horizontal axis represents time, while the vertical axis represents periods on the scale of days. The cone of influence showing the region remaining at 5% significance level is separated by a fine black line. Areas with consistency within this region are surrounded by thick black contours and are filled with a blue-red color gradient indicating higher correlations toward red. The color correlation scale is given on the right side of each graph. The small arrows inside circles represent the direction of the correlation between COVID-19 and air pollutants.
Figure 5
Figure 5. Wavelet transform coherence (WTC) between new confirmed cases and air pollutants.
Figure 6
Figure 6. Wavelet transform coherence (WTC) between daily deaths and air pollutants.
Figure 5 shows the WTC graphs between confirmed daily case and PM10, PM2.5, SO2, CO, NO2 and O3. Four correlation circles appear for PM10. The first is in the 0–4 frequency band at the end of May beginning of June. The left-up arrows dominating this region indicated a lead non-phase correlation for PM10 according to daily cases (R2 > 0.8). The second is in the 4–8 frequency band at the end of June and the other two are in the 8–16 frequency interval in May and July, respectively. For PM2.5, there were three correlation circles in the 0–4, 4–8 and 8016 frequency domains. The left-up arrows in the islet at the end of May indicates a direct negative and high correlation between PM10 and daily cases. The right-up and left-down arrows in other areas show in-phase and out-of-phase correlations leading daily cases. For the SO2 graph, the 0–4 frequency band had two regions in May and July, respectively. The dominant left arrows in these regions indicate negative correlations. For CO there are three coherence regions with two in the 0–4 band and the other in the 4–8 frequency band. Here, the left and up arrows show CO leads and has non-phase correlation. The NO2 graph has six islets surrounded by thick contours. The left and up arrows in the 0–4 frequency band in May show NO2 had a negative correlation leading COVID-19 cases. The right arrows in the other regions indicate in-phase correlation and the down arrows represent NO2 leading, while the up arrow indicates daily cases leading. For O3, there are five large and small correlation islets. Here, there are correlation regions for the 0–4 frequency band in May, June and July, for the 4–8 frequency band in June and the 8–16 frequency band in April. The left arrows indicate non-phase correlations and left-down arrows show daily cases were leading.
Figure 6 shows the wavelet coherence graphs between daily mortality and air pollutants. Three high-correlation regions are notable for PM10. For the frequency domains of 0–4, 4–8 and 8–16 the right arrows for the regions during April and May indicate high positive correlations (R2 > 0.9). For the region from the end of April beginning of May in the 4–8 frequency band, the left-up arrows show out-of-phase extending from PM10 to daily deaths. The 0–4 and 4–8 frequency domains remaining within the cone of influence for PM2.5 contain five coherence regions. The dominant left arrows in these regions indicate negative correlations. Similarly, for SO2, there are several correlation islets in the 0–4 and 4–8 frequency bands and in different time intervals. Here the left-up and right-down arrows show the lead of daily deaths. For all time intervals on the CO graph, 0–4, 4–8 and 8–16 frequency bands have five coherence areas. The dominant left-down arrows for the region at the end of April encompassing the 4–8 and 8–16 frequency bands indicate daily deaths lead compared to CO. The small islet in May for the 8–16 frequency band shows CO leads and has positive correlation. There are islets in the 0–4 and 4–8 frequency interval for NO2. There is a mutual relationship for regions in both frequency areas in May. In July, an out-of-phase relationship is present. Additionally, for the frequency area above 16, there are down and left arrows showing positive correlation with NO2 leading for the region in May. Finally, the O3 graph appears to have many scattered islets in a broad time interval in the 0–4, 4–8 and 8–16 frequency intervals. The right arrows showing in-phase correlation are dominant in the 0–4 frequency band. The down arrows in the 4–8 and 8–16 frequency domains indicate O3 leads compared to deaths.
Results for wavelet coherence reveal that PM2.5 and NO2 are the main determinants for daily confirmed cases. PM2.5 has short-term effects and NO2 has short-medium-term effects on the spread of COVID-19. The other pollutants have effects in narrower time intervals. Between daily mortality and all air pollutants, links with high correlation were identified in different periods. While the effects of PM10 and O3 spread over a larger period, PM2.5, SO2, CO and NO2 have short- and medium-term effects on mortality.
Research investigating the short- and medium-term effects of atmospheric particulate matter on COVID-19 reported the highest mortality risk (relative risk: 1.23, 95% CI) was exposure to aerosols in the outdoors on the 14th day (Wang et al. 2020a). Martelletti and Martelletti (2020) reported a significant correlation between high concentrations of atmospheric aerosols (PM10 and PM2.5) with high transmission and mortality rates for northern regions of Italy most affected by COVID-19. Several studies revealed that particulate matter amounts play a key role in the spread of COVID-19 (Bashir et al. 2020a; Yao et al. 2020). Additionally, Zhao et al. (2019) reported that PM10 contained higher levels of avian influenza virus compared to PM2.5, but that PM2.5 spread more in air. According to a study encompassing 120 major cities in China, PM2.5 is more effective than PM10 on new COVID-19 cases (Zhu et al. 2020). The present study obtained findings showing PM2.5 had more effect compared to PM10.
Bashir et al. (2020b) found a significant correlation between PM10, PM2.5, SO2, NO2 and CO levels in California with the COVID-19 pandemic. A study applying the generalized additive model by Zhu et al. (2020) proposed that a 1 mg/m3 increase in CO caused a 15% increase in daily confirmed cases. Research including 10 provinces in Turkey found significant correlations between major air pollutants and the pandemic and the highest positive correlation was obtained for SO2 (Şahin 2020a). For people dying due to COVID-19 in Milan, Italy, Zoran et al. (2020) found a negative correlation with NO2 and a positive significant correlation with ground ozone levels. Additionally, they mentioned the significant negative effects on human immunity and respiratory systems of air pollutants, which are probable carriers of viral infections. According to Ogen (2020), long-term exposure to NO2 has severe effects on the mortality of COVID-19. Research investigating the effects of air pollutants on COVID-19 in Germany proposed that PM2.5, NO2 and O3 were the main determinants of the pandemic (Bashir et al. 2020a). In South America, PM10, NO2, CO and O3 were included among the most important pollutant parameters that need to be controlled in the struggle against the pandemic (Bilal et al. 2021). A study in Dhaka city revealed O3 was the essential pollutant determining the spread of the pandemic (Rahman et al. 2021).
According to empirical findings, high levels of air pollution can be said to contribute to the spread and mortal effects of COVID-19 by notable amounts. Many studies confirm that air quality has key importance in reducing the effects of the pandemic (Lolli et al. 2020; Xu et al. 2020; Zhang et al. 2020). When assessed from this aspect, it is considered greatly important that in addition to limiting human interactions, the authorities take sustainable environmental precautions about developing air quality.