Introduction

Many health, social, and economic problems have emerged with the outbreak of coronavirus disease 2019 (COVID-19) in Wuhan, China. In addition to causing pneumonia, the virus damages the heart, liver, and kidneys, as well as the immune system as a whole, as a result of which COVID-19 patients die due to multiple organ disorders (Huang et al. 2020). To date, no proven and effective treatment method has been identified against the virus (Cortegiani et al. 2020). The pandemic has spread rapidly from a single city to an entire country in as little as 30 days (Wu and McGoogan 2020). Since the virus is transmitted rapidly from person to person (Chan et al. 2020), numerous countries have called on their citizens to stay at home and apply social distancing rules to prevent the pandemic from spreading. Various lockdown measures have been implemented to flatten the pandemic curve, such as shutting down industries, halting vehicular traffic, increasing social distance, and stopping non-essential business activities (Bherwani et al. 2020). These measures, in turn, affect economic production and consumption activities. In the COVID-19 era, industrial activities have slowed down, vehicle use has decreased, the demand for imported goods has decreased, and many countries have suspended air travel—both international and domestic. Economic activities and environmental pollution, especially air pollution, are closely related. Particle matter 2.5 (PM2.5), one of the air pollution indicators, causes cardiovascular disease and lung cancer (Khan et al. 2017) and increases psychological distress, and for these reasons, the air pollution problem needs to be solved urgently (Xu and Liu 2020). In this respect, COVID-19 can be a solution for reducing PM2.5 emissions. Transport activities are significantly affected by COVID-19 lockdowns, resulting in less energy consumption and oil demand. As a result of lockdowns, it can be easily said that air quality has improved (Gautam 2020). However, COVID-19 also has some negative effects on the environment.

Zambrano-Monserrate et al. (2020) noted that despite the negative effects of COVID-19 on the environment due to increased waste and reduced recycling activities, the pandemic has also had a number of positive effects—for example, clean beaches, reduced environmental noise, and a reduction in nitrogen dioxide and particulate matter. The Copernicus Atmosphere Monitoring Service (CAMS) examines PM2.5 emissions in the atmosphere using satellite images from various countries. According to CAMS (2020), China’s average PM2.5 emissions in February 2020 were between 20 and 30% lower than the averages for the same month in 2017, 2018, and 2019. Similarly, some researchers have noted that the pandemic has reduced air pollution in various countries. For instance, Tobías et al. (2020) observed that the pandemic resulted in a reduction in emissions of nitrogen dioxide (NO2), black carbon, and particulate matter with a diameter of less than 10 (PM10) during the 2-week lockdown in Barcelona, Spain. Kerimray et al. (2020) concluded that the COVID-19 pandemic reduced PM2.5, CO (carbon monoxide) and NO2 emissions during the 27-day lockdown in Almaty, Kazakhstan. Moreover, Dantas et al. (2020) reported that the pandemic reduced CO and NO2 emissions in Rio de Janeiro, Brazil, from March 12, 2020 to April 16, 2020. Sharma et al. (2020) found that the pandemic decreased PM10, PM2.5, CO, and NO2 emissions in 22 Indian cities during the lockdown period. However, understanding the effect of the COVID-19 quarantine processes on environmental pollution requires more than the use of satellite data (Bao and Zhang 2020). Therefore, Asna-ashary et al. (2020) empirically investigated the pollution–COVID-19 nexus and reported a negative relationship between PM2.5 emissions and positive shocks in COVID-19 cases in 31 Iranian provinces.

The pandemic has also reduced air pollution in the USA, the world’s largest economy, whose production and consumption activities cause a high rate of air pollution. The first patient identified with COVID-19 in the USA was seen in Washington State on January 20, 2020, and as the pandemic spread rapidly, the USA became the country with the highest number of both cases and deaths. As of April 13, 2020, at least one COVID-19-related death had occurred in each of the 50 states of the USA. As of May 4, 2020, 1,212,000 cases and 69,921 deaths had been reported in the USA. On the same date, the global number of cases was 3,639,000 and the global number of deaths 252,240 (European Center for Disease Prevention and Control 2020). The USA thus accounted for 33% and 27% of worldwide COVID-19 cases and deaths, respectively. The number of cases and deaths continues to increase in the USA and the rest of the world, and the spread of the COVID-19 virus in one country can adversely affect other countries. However, the virus may have more of a positive effect on the environment during lockdown in places with a high population. For both reasons, we empirically analyzed the impact of worldwide COVID-19 cases and deaths on PM2.5 emissions in eight US cities with populations over 1 million.

The rest of the paper is organized in the following manner: “Data and methodology” provides the data and methodology, while “Methodology” presents and discusses the empirical results obtained in the study. Finally, “Results and discussion” gives a summary of the findings and concludes the study.

Data and methodology

Data

In this study, we examined the effect of COVID-19 on environmental pollution in eight US cities (New York, Los Angeles, Chicago, Phoenix, Philadelphia, San Antonio, San Diego, and San Jose) for the period of January 15, 2020 to May 4, 2020. Since environmental pollution data for Houston and Dallas were not available, these cities were excluded from the analysis. Data relating to worldwide COVID-19 cases and deaths were obtained from the European Center for Disease Prevention and Control (2020), while the PM2.5 data (daily per cubic meter air, μg/m3) for the eight US cities were collected from the United States Environmental Protection Agency (2020). The data utilized in the study are converted into natural logarithm to obtain a more stable data variance. The descriptive statistics of the data are illustrated in Table 1.

Table 1 Descriptive statistics of the variables

In terms of PM2.5 emissions, information from Table 1 illustrates that Los Angeles has the highest mean and median values followed by Chicago and San Diego. On the contrary, Phoenix and San Jose have the lowest PM2.5 emissions. Moreover, Skewness statistics demonstrate that 7 out of 10 variables are skewed (except PM2.5 emissions in Los Angeles, Phoenix, and Philadelphia), and Kurtosis statistics illustrate that all variables are leptokurtic. The Jarque-Bera test statistics indicate that PM2.5 emissions in six US cities are normally distributed based on a 5% significance level (except Chicago and San Diego). Regarding the COVID-19 variables, both the number of the cases and deaths do not follow the normal distribution.

After investigating the characteristics of the variables, we applied Fourier Lagrange multiplier (LM) unit root and asymmetric Fourier causality tests.

Methodology

Fourier Lagrange multiplier unit root test

Enders and Lee (2012) developed the LM-based Fourier unit root test on the basis of Gallant’s (1981) Fourier approximation. This approximation captures smooth structural shift using a small amount of low frequency information. The first step to implement the Fourier LM unit root test is shown in Eq. (1):

$$ \kern1.25em {\mathit{\Delta x}}_t={\beta}_0+{\beta}_1\Delta \sin \left(\frac{2\uppi kt}{T}\right)+{\beta}_2\Delta \cos \left(\frac{2\uppi kt}{T}\right)+{z}_t $$
(1)

In the first-differenced regression, represents the difference operator, β0 indicates the constant term, k denotes a particular frequency, and β1 and β2 illustrate the amplitude and displacement of the frequency approximation. With the estimated coefficients β0, β1, and β2, the detrended series is formed as in Eq. (2):

$$ \kern1.25em {\overset{\sim }{S}}_{\mathrm{t}}={x}_t-\overset{\sim }{\psi }-{\overset{\sim }{\beta}}_0\mathrm{t}-{\overset{\sim }{\beta}}_1\sin \left(\frac{2\uppi kt}{T}\right)-{\overset{\sim }{\beta}}_2\cos \left(\frac{2\uppi kt}{T}\right),\kern0.5em \mathrm{t}=2,\dots .,\mathrm{T} $$
(2)

where \( \overset{\sim }{\psi }={x}_1-{\overset{\sim }{\beta}}_0-{\overset{\sim }{\beta}}_1\sin \left(\frac{2\uppi kt}{T}\right)-{\overset{\sim }{\beta}}_2\cos \left(\frac{2\uppi kt}{T}\right) \), and x1 is the first observation of xt. At the last stage, the Fourier LM unit root test was performed using the detrended series:

$$ \kern1.25em {\Delta x}_t=\upvarphi {\overset{\sim }{S}}_{t-1}+{\alpha}_0+{\alpha}_1\Delta \sin \left(\frac{2\uppi kt}{T}\right)+{\upalpha}_2\Delta \cos \left(\frac{2\uppi kt}{T}\right)+\sum \limits_{i=1}^k{\vartheta}_i\Delta {\overset{\sim }{S}}_{t-\mathrm{i}}+{v}_t $$
(3)

In Eq. (3), the null hypothesis of unit root (H0: φ = 0) is tested using the t-statistic. The test statistic (τLM) depends only on the frequency k. Therefore, the critical values tabulated by Enders and Lee (2012) are a function of k. In addition, the authors used F-statistics to test the significance of the Fourier component as follows:

$$ \kern0.75em {\mathrm{F}}_{\upmu}(k)=\frac{\left({\mathrm{SSR}}_0-{\mathrm{SSR}}_1\right)/\mathrm{q}}{{\mathrm{SSR}}_1(k)/\left(T-k\right)} $$
(4)

where q indicates the number of regressors, SSR0 denotes the sum of squared residuals from the regression without Fourier approximation, while SSR1 represents SSR from the regression containing the trigonometric terms. When the F-statistic is greater than the critical value, it is convenient to use the Fourier LM unit root test; otherwise, more reliable and powerful results can be obtained using conventional unit root tests without a Fourier term.

Asymmetric Fourier causality test

Researchers began to investigate causal relations between macroeconomic variables by using the Granger (1969) causality test. However, the Granger and many other causality tests in the literature, such as that by Toda and Yamamoto (TY; Toda and Yamamoto 1995), neglect structural breaks that may occur in the series. To compensate for this negligence, Enders and Jones (2016) and Nazlioglu et al. (2016) proposed the Fourier Granger and Fourier TY causality tests, respectively. These tests are performed by adding Fourier functions to the equation, just like the Fourier LM unit root test. The authors referred to above stated that the null hypothesis could be rejected more accurately using this approach. Nazlioglu et al. (2016) relaxed the assumption that the constant term does not change over time. The model used for the Fourier TY causality test is shown in Eq. (5):

$$ \kern2.5em {y}_t={\alpha}_0+{\gamma}_1\sin \left(\frac{2\uppi kt}{T}\right)+{\gamma}_2\cos \left(\frac{2\uppi kt}{T}\right)+{\beta}_1{y}_{t-1}+\dots +{\beta}_{p+d\max }{\mathrm{y}}_{t-\left(p+d\max \right)}+{u}_t\kern7em $$
(5)

In the equation, yt represents the vector containing the variables of COVID-19 cases and deaths, and PM2.5 emissions, β is the coefficients matrix, t is the trend, T denotes the number of observations, γ1 and γ2 are the coefficients of the Fourier approximation that smooth structural shifts are captured, and dmax is the maximum integration degree of the series that can be determined by a unit root test. In our study, the optimal lag length p and the Fourier frequency k are determined by the Akaike information criterion (AIC). In the single-frequency Fourier TY causality test, the null hypothesis of no causality is tested as H0 : β1 = …βp = 0.

The reactions of people, firms, and decision units to positive and negative shocks are different. Analyzing the effects of both shocks as a whole leads to the hidden causal relationships being ignored. In order to reveal hidden causal relationships, Hatemi-J (2012) suggested separating variables into positive and negative shocks and applying the causality test to these shocks. Therefore, following Yilanci et al. (2019), we performed the asymmetric Fourier causality test by considering the cumulative positive and negative shocks of the variables. This test is applied by adding the shocks to the yt in Eq. (5). All other procedures are the same as those of a single-frequency Fourier TY causality test. This test offers two main advantages. First, it allows the causal relationships between the positive and negative shocks of the variables to be examined separately. Second, in the asymmetric Fourier TY causality test, structural breaks with an unknown number, form, and date are taken into account in the analysis. Due to the advantages, the asymmetric Fourier TY causality test provides a more accurate rejection of the null hypothesis of no causality. While applying this approach, a variable can be transformed into positive and negative components, as in Eq. (6):

$$ {\mathrm{PM}}_{2.5_t}={\mathrm{PM}}_{2.5_{t-1}}+{\varepsilon}_{1t}={\mathrm{PM}}_{2.5_{1,0}}+\sum \limits_{i=1}^t{\varepsilon_1}_i^{+}+\sum \limits_{i=1}^t{\varepsilon_1}_i^{-} $$
(6)

where \( {\mathrm{PM}}_{2.5_{1,0}} \) indicates the initial value of the relevant variable and \( {\varepsilon_1}_i^{+} \) and \( {\varepsilon_1}_i^{-} \) represent positive and negative shocks, respectively. This process is carried out in the same way for each variable analyzed. The shocks are then included in the Fourier causality equation. In this study, we investigated whether there is a causal relation from positive shocks of COVID-19 to positive and negative shocks of PM2.5 emissions. Because there is no cure for COVID-19, and therefore no negative shock for deaths and cases, we can only examine the positive shocks of COVID-19.

Results and discussion

In the first phase of the analysis, we investigated the stochastic properties of the variables to determine the maximum order of integration (dmax). Column 4 of Table 2 demonstrates that F-statistics are significant for all variables. Therefore, we decided to use trigonometric terms in unit root analysis and applied the Fourier LM unit root test to obtain more robust findings. According to the τLM statistics presented in Table 2, the raw data on COVID-19 cases and deaths contain a unit root. These variables are stationary in their first difference. At the same time, the PM2.5 emissions of eight cities are stationary at level I(0). To save space, the results of negative and positive components are not presented in the table. Positive and negative shocks of COVID-19 cases and deaths are also non-stationary at level.

Table 2 Fourier LM unit root test results

After determining the order of integration of the variables as 1, we analyzed the effects of the worldwide COVID-19 cases on the PM2.5 emissions of eight cities. According to the results presented in Table 3, we determined that the number of COVID-19 cases cause PM2.5 emissions only in Chicago. However, when we used the asymmetric causality test, the findings changed significantly. The findings of the asymmetric Fourier causality test illustrate that an increase in the number of cases reduces PM2.5 emissions in Los Angeles, Chicago, Phoenix, Philadelphia, San Antonio, and San Jose. Since there is no decrease in the number of cases and deaths, these variables do not have negative components. In addition, there is no relationship between the positive components of COVID-19 cases and PM2.5 emissions. This is not surprising, since COVID-19 reduces use of fossil fuels such as oil and coal, which are primary air polluters.

Table 3 The results of asymmetric Fourier causality test for COVID-19 cases

The causal relationships between worldwide COVID-19 deaths and PM2.5 in US cities are displayed in Table 4. According to symmetric causality test results, COVID-19 deaths are the cause of PM2.5 emissions in New York and Los Angeles. The results of the asymmetric Fourier causality test demonstrate that an increase in the number of deaths reduces the release of PM2.5 emissions in New York, San Diego, and San Jose. Overall, an increase in the number of cases affects air pollution more than an increase in the number of deaths. Therefore, it can be said that an increase in COVID-19 cases caused people to take more precautions and thus slow down economic activities.

Table 4 The results of asymmetric Fourier causality test for COVID-19 deaths

To sum up, COVID-19 deaths and cases positively affect environmental quality by reducing economic and social activities in eight cities with high populations in the USA. In line with the findings of Pata (2019), who stated that the 2001 economic crisis reduced carbon emissions in Turkey, our results indicate that the COVID-19 crisis reduced PM2.5 emissions in the selected US cities. Perhaps the only positive aspect of economic crises or pandemics is that they reduce human pressure on the environment. Humans destroy nature for the sake of their economic and social interests, and when human activities cease, nature can return to its balance. For this reason, humankind must review existing production and consumption activities with a view towards ensuring a cleaner environment and a more sustainable future.

Conclusions

This study investigates the effect of COVID-19 deaths and cases on environmental pollution in the USA. The results of the asymmetric Fourier causality test demonstrate that COVID-19 reduces PM2.5 emissions in US cities. An increase in the number of cases of COVID-19 affects pressure on the environment more than an increase in the number of deaths. Another important finding of the study is that positive and negative shocks should be taken into consideration. When shocks are not studied separately, a unidirectional causality from COVID-19 to PM2.5 emissions is found for New York, Los Angeles, and Chicago. However, an increase in positive shocks of COVID-19 causes negative shocks of PM2.5 in the eight high-population cities studied. In other words, an increase in worldwide COVID-19 deaths and cases causes a reduction in PM2.5 emissions. The rapid spread of the virus in the USA, especially in New York, led the lockdown process to be introduced. This resulted in the industry and service sectors largely ceasing their production activities. The slowdown in economic activities led to a reduction in environmental pollution. This demonstrates that environmental pollution is a man-made phenomenon, and that people are harming the natural environment in which they live. For this reason, perhaps, COVID-19 will assume its place as a pandemic that increased human awareness about environmental issues. Moreover, the effects of COVID-19 lockdown on air pollution will create an important opportunity to evaluate future air quality policies.

The overall findings imply that PM2.5 emissions can decrease when people are prevented from harming the environment. PM2.5 emissions are generally due to transport activities and the use of fossil energy sources such as oil and coal. The US government can improve air quality by promoting substitution of fossil fuels with renewables, by enforcing strict execution of air quality control plans, and by implementing awareness-raising programs on environmental issues. In addition, the US government and the private sector can expand remote working opportunities brought by the COVID-19 in the coming period, thereby reducing air pollution. Following the COVID-19 pandemic, if the regulatory authorities in the US take the necessary measures, PM2.5 emissions in cities can be reduced, and thus, environmental quality can be improved.