key: cord-0867040-r1fe05ga
authors: Khan, Yousaf Ali
title: Risk of mortality due to COVID-19 and air pollution in Pakistan
date: 2021-08-07
journal: Environ Sci Pollut Res Int
DOI: 10.1007/s11356-021-15654-z
sha: b3995dbf6a600bea95049208dd023b126e6877e0
doc_id: 867040
cord_uid: r1fe05ga

The present research aims to investigate the impact of air pollution on the number of mortalities caused by COVID-19 per Pakistani province. To do so, for each independent area of Pakistan, the observed mortality due to COVID-19 has been standardized over the entire population using standard age groups ranging from 0 to 4, 5 to 9, 10 to 14,…, 65, and above years, supported by the 2017 state people census. The impact of air pollution and COVID-19 transience among Pakistani areas, Islamabad Capital Territory (ICT), and the Federally Administered Tribal Region (FATA) was analyzed by a multiple-linear regression model, while the broad collection of attributes was observed by the resources of local spatial autocorrelation indicators, including the spatial portion of COVID-19 association. The result indicates that the observed mortality rate is much higher than predicted in certain provinces, namely, the Khyber Pakhtunkhwa and Punjab provinces, and the prevalence of PM(10) was independently linked to mortality due to the corona virus. Additionally, the results of the local spatial autocorrelation indicators on the standardized mortality rate and PM(10) define a collection of very higher ideologies in the broad range of KPK and the southern part of Punjab province, respectively, with a definite degree of connection between the two distributions in the Khyber Pakhtunkhwa region. In brief, this research seems to find a justification for confirming the existence of a correlation between the possibility of COVID-19 mortality and air pollution, more precisely considering air pollutants (i.e., particulate (PM(10)) and land take-over. To this end, the need to mediate in favor of measures aimed at eliminating emissions in the environment will be reiterated by speeding up current proposals and policies aimed at all causes of atmospheric pollution: urbanization, water and manufacturing, home heating, and transportation.

Air pollution is one of the biggest public quality concerns in the country. According to the International Health Organization, it is attributed that there are 7 million deaths globally each year (WHO Timeline 2020). In Asia, an expected 852,000 people die hastily annually from experience with air contaminants (WHO 2017) . Carbon monoxide (CO), particulate matter (i.e., PM 10 and PM2.5), carbon dioxide (CO 2 ), and nitrogen-dependent contaminants are the major air contaminants taken into consideration (e.g. NO x ). These compounds are mainly extracted from anthropogenic practices (trafficking, burning, development and industrialization, deforestation and urbanization, population increase, etc.), livestock breeding, and cultivation that modify the composition of the environment (Ghorani-Azam et al. 2016; Rapporto Quali ISPRA 2016) . The disability induced by human health fluctuates and relies on the subjectivity of the person, the length of exposure, and the perception of contaminants. Various conditions may be related to revelation; moreover, they vary from judicious and temporary to irreversible figures. Furthermore, sustained experience with toxic intensities above the control restrictions has been shown to trigger a persistent, provocative, and hypermetric syndrome that can contribute to a greater preference for viruses and growth of indicative diseases (Shi et al. 2020; Ma et al. 2020) , as well as additional preexisting resistant shifts (Shi et al. 2020; Mazhar et al. 2020) . Livelihood in an atmosphere of elevated levels of pollutants may often render a human more fragile, and thus it is more likely for infectious agents to develop prolonged respiratory problems. Relatively higher then in 2005, because the WHO has haggard up strategies defining experience thresholds for the small and lengthy term. In view of this, global monitoring thresholds are not the same and differ from country to country (WHO 2016) . The new epidemic of the COVID-19 pandemic SARS-CoV-2 virus has posed several questions in the academic world as to why certain countries have seen larger dissemination of the virus and a higher mortality rate (Lai et al. 2020; Shi et al. 2020; Ma et al. 2020) . Pakistan is, in reality, one of the developed countries of Asia, the country most impacted by the corona virus pandemic. The pandemic has acquired a striking character in the world with 45,70,750 viruses and 10,192 death confirmations as of 31 December 2020, comparatively senior mortally than further Asian countries and several inequalities inside the same national territory as well (Frontera 2020; Mazhar et al. 2020) . Numerous suggestions have been placed ahead between the potential details for the principal increase of the virus surrounded by the area, including regarding recent findings showing different associations from the end in view of the geographical and public contact between China and Pakistan (Waris et al. 2020; Anjum et al. 2020) . With the trade links connecting Pakistan and China, and the combination of possibilities for the virus to be contacted and imported into our country, knowing why the virus has advanced so suddenly in Pakistan's territory is one of the problems currently boiling the international science discourse. In addition, many other features might be concerned with making easy extent of viruses into society. In reality, a traditionally established occurrence about the spread of extra viruses (i.e., measles) (Peng et al. 2020) indicates that amounts of air effluence and, in particular, particulate matter could serve as a vector within the territory for the spread of the virus. The presence of SARS-CoV-2 RNA in particulate matter, however, was recently demonstrated by Setti et al. (2020a) . Consequently, a mixture of air pollution causes, such as pollutants and a broad spectrum of environmental issues in general, may be deemed liable for the objectives within the respiratory territory and declining the at-risk community, thus enhancing their chances of being impacted by respiratory diseases such as COVID-19 (Frontera 2020; D'Alessandro et al. 2020) . Khyber Pakhtunkhwa, renowned for its extreme air pollution rate, is home to the most affected areas of Pakistan. This study aims to determine the causal relationship between air pollution and COVID-19 at the regional level from the beginning of these assumptions and a set of current studies showing an advanced case casualty incidence in the Khyber Pakhtunkhwa area (Anjum et al. 2020; Mazhar et al. 2020 ). More specifically, in relation to the number of deaths caused by corona virus per Pakistani province, the function of air pollution and a collection of ecological variables, preferred from current comments (Waris et al. 2020; Murgante et al. 2020; Dettori et al. 2021) , were studied.

Pakistan was partitioned from a managerial position of observation into five districts: Azad Jammu and Kashmir (AJK), Islamabad Capital Territory (ICT), and the Federally Controlled Tribal Area (FATA). Punjab province is divided into 2 independent regions with regional obligations. COVID-19 results, as declared by the Ministry of Public Health Services, Government of Pakistan, and as obtained by the Pakistan Statistics Division (PSD), take into account, as of 31 December 2020, the number of overall tainted persons and the numeral of deaths at the provincial level (Pakistan Statistics Division 2020). Data and computational environment knowledge is given in Sect. 2.3, while references to the description and data attainment applicable to environmental variables are defined in Table 1 .

The remainder of this research is structured as follows: Sect. 2 presents the methodology and data used in this research, Sect. 3 provides empirical results, while in Sect. 4 the discussion is given. Finally, Section 5 ends this research with useful policy implementations and future research work.

Based on the 2017 national population census, standard mortality from COVID-19 over the entire population age classes (14 groups) ranging from 0 to 4 years to 65 years and above was averaged for each zone of Pakistan. The indirect standardization process initially rendered it possible to predict national real mortality by age category, obtained by dividing the number of COVID-19 deaths by 31 December 2020 and verified by the Ministry of National Health (MNH) by the 14 age groups listed at http://www.covid.gov.pk. Thus, the number of estimated age groups is based on the 2017 demographic survey for the prespecified age groups. COVID-19 deaths in all provinces, ICT, AJK, and FATA areas were estimated according to the formula:

In each observed region, P i is the unique age group population, and M i is the resulting mortality rate for that region's specific age group. The standardized mortality ratio (SMR) was obtained by comparing, in each province, the number of reported COVID-19 deaths with the corresponding number of predicted deaths:

where D is the number of deaths observed for COVID-19; the number of deaths predicted is "a." Ninety-five percent confidence intervals (95 percent CI) were determined for all regions under consideration, following the recommendations of Vandenbroucke (1982) .

With the assistance of the Pearson correlation, the relationship between the environmental variables under consideration was evaluated and mathematically expressed as follows:

Local predictor spatial autocorrelation According to Anselin (1995) , local predictor spatial autocorrelation (LISA) can be defined as: A statistic that satisfies the following two conditions is called LISA:

1. LISA for each observation indicates the extent of significant spatial clustering of similar values around that observation.

The sum of all LISAs of observations is proportional to a global indicator of spatial correlation.

LISA for a variable x j observed at loction j as measure L i can be mathematically expressed as

where f is a function of parameters, and x k j are the observed values in the neighborhood k j of j. The values x used in the derivation of the statistic may be the actual observations or some standardization of these to avoid scale dependence of the local indicators. Similar to the exercise often taken for global indicators of spatial association.

The data used in this analysis were extracted from numerous sources. Data related to COVID-19, such as a number of confirmed cases, deaths, active cases, and area recoveries, was collected from the Ministry of National Health Services (MNHS 2020) available at https://www.covid.gov.pk. Realtime air quality monitoring project (AQICN 2020) available at https://www.aqicn.org/country/pakistan/, Ministry of Environment, Government of Pakistan available at http:// www.adaptation-undp.org, and air pollution data were collected from Pakistan Ministry of Climate Change available at http://www.mocc.gov.pk. Results and graphs mentioned in this article were carried out using R-studio computational environment.

With the aid of a multivariate linear regression model, the association between mortality due to COVID-19 and all variables causing air pollution across various regions of Pakistan was well explained. The functional relation between this study's explained variable (SMR) and explanatory variables can be expressed mathematically as follows.

Equation (5) can be rewritten as follows after adding a constant term (α 0 ) and an error term ℇ t wherein (β k ) (k=1,2,..,7) represent the coefficients of PM 10 , PM 2.5 , CP, PA, Soil, NO 2 , Trees, MC, UGA, and BT. 

The findings show that major air toxins such as PM 10 , NO 2 , and PM 2.5 and different factors unique to the municipal climate, such as trees, urban population, bike lanes, walking streets, soil, community and town green areas, motorcycles/ cars plus heavy transport, as an independent variable, are defined in detail in Table 1 . In a recent analysis (Waris et al. 2020; Murgante et al. 2020; Dettori et al. 2020) , these variables were previously analyzed and defined on the foundation of their straight and tortuous connection to the frequency of investigation of air contaminants, evaluated through the explained variable outcome (SMR). In specific, linear regression was approved by means of the interrupted averages from the official measurement systems of the air pollutants collected in this study. At 5% t-statistics, the significance number was determined. The related p values and 95 percent confidence intervals were determined using R-statistical software for each coefficient and provided in the computation portion of Table 4 .

Based on official national statistics, the standardized mortality rate (SMR) contrasted the COVID-19 mortality rate with what was predicted. The SMR, shown in Table 2 , is far higher than predicted in all regions, namely, the Khyber Pakhtunkhwa province, including the northern areas and the Federally Administered Tribal Territory (FATA). Table 2 displays the average annual values for PM 10 corresponding to each region. Table 2 contains data pertaining to SMR for all areas of Pakistan. Along with COVID-19 reported cases, the area under consideration is graphically depicted in Fig. 1 , while the SMR for all regions under examination is shown in Fig. 2 . 

The coefficients of association among the ecological variables analyzed are exposed in Table 3 . Study of regression between environmental variables (such as NO 2 , PM -10 , PM 2.5 , bike paths, trees, soil, pedestrian areas, urban population, urban green spaces, motorcycles/cars, and bus/trucks), and the real mortality rate is given in Table 4 . The existence of PM 10 was separately correlated with the case status (p = 0.0011, 95 percent CI: 0.069, 0.234). There was no important correlation with the other variables (p > 0.05) with case status. The interaction with the SMR is seen in Fig. 3 considering only the PM 10 variables. In addition, Fig. 4 displays the country-wide scenario, separated by province with respect to SMR, thus highlighting the regional distribution with respect to PM 10 . Five groups were considered for the SMR distribution over the entire region of the Pakistan table. The two groups (below average and average) have lower than predicted mortality rates correlated with COVID-19. The intermediate class has a value about unity and is marginally higher than unity, thus indicating an estimated or slightly higher mortality rate, even as the extra two outlier have values where the mortality associated with COVID-19 is considerably higher than expected. The variety stands for the WHO and Pakistan particulate matter limits supported on a manifold of 10 (i.e., 20 μg/mc is the WHO boundary, while 68.3 μg/m 3 is the Pakistan limit) in comparison to the PM 10 map.

This study has been carried out following the recommendation of Anselin (1995) spatial autocorrelation analysis. Regression analysis was compared for comparative purposes, in regard to spatial autocorrelation. The conclusions of the study of spatial autocorrelation are summarized and presented in Fig. 5 

The risk of mortality due to COVID-19 depends largely on the incidence of the disease and the rate of spread, and the main aspect is the human activity within the community. However, other variables could have a fun or a function in the promotion of this disease, such as hereditary issues (Frontera "*" and "**" indicate significance at 5 and 1% (Mehmood et al. 2018; Waris et al. 2020; Anjum et al. 2020; Dettori et al. 2020) . The subsistence of correlations amid air quality and COVID-19 humanity was first examined (Giangreco 2020; Setti et al. 2020a; Wu et al. 2020 ) with respect to air emissions from air pollutants such as ("PM 10 " and "PM 2.5 "). Nitrogen-dependent constituents and elements are produced from human activities. The core idea was the proliferation of air-related emissions would bring a strain on the health status of at-risk populations and establish requirements for the growth and complications of ventilator illness, which Dettori et al. (2021) in Po Valley, Italy, further implemented. In Pakistan, as of 31 December 2020, the current area-wise SMR in Table 2 indicates a higher mortality rate in Pakistan. An increase of 6.7 times more than predicted was considered to represent the collection of 8 regions in which KPK is one of the most polluted, with a value of 6.7. It can be shown that nearly all impacted regions of all regions are located in the KPK and are particularly characterized by a high average annual worth of PM -10 . In the density groups, which typically contain between 238 and 268 individuals per quadrangle kilometer, a connection is also evident. The COVID-19 humanity analysis shows a very simple image between the area of KPK and the other regions included in these studies, with values a great deal senior than anticipated in the southern region and mortality standards consistent through those predictable in extra districts of Pakistan, especially in southern Punjab. Therefore, on the basis of recently published studies, a set of environmental variables was chosen to test the relationship with humanity due to COVID-19 in all regions of Pakistan. Keeping in brain that the major basis of air pollution may be similar, following the advice of Dettori et al. (2021) , we performed a correlation study between the variables that cause air pollution in the field, and a few of these variables illustrated a very strong connection. In particular, the uppermost association coefficient between PM 10 and PM 2.5 was found to be equal to 0.8983. Similarly, the correlation between PM 10 and NO 2 is 0.6201. This is due to the fact that the particulate measurement schemes at this time in use revisit the average PM -10 attention, counting the PM 2.5 portion. The two variables are, thus, closely connected. However, as the multiple-linear regression shows, among those variables that have been found to have a close relationship to the figure of deaths dispersed for every Pakistani region, artificial by COVID-19, PM -10 , PM -2.5 , and NO 2 are among those variables that refute the results of Italy (Dettori et al. 2021) . Pakistan was ranked number two in the world's air quality nations, which is why all air emissions tended to be statistically important for Pakistan. Since the variables PM 10 , PM 2.5 , and NO 2 are statistically relevant, we consider PM 10 as an important risk stratification factor in this analysis because PM 10 did not decrease significantly during the lockdown, which is why PM 10 was considered an independent mortality predictor for COVID-19 and early epidemic recurrence pointer as indicated by Setti et al. (2020a) and Dettori et al. (2021) . In addition, certain spatial autocorrelation with respect to PM 10 and SMR can be observed, indicating a certain degree of similarities between the COVID-19 regions most affected in Pakistan and the highest reported PM 10 values. Many studies hypothesize a link between the prevalence of air contaminants and humanity. In addition to a current methodical Copat et al. 2020; Dettori et al. 2021) , which highlights the significant contribution of PM -10 , PM -2.5 , and NO 2 activates for the increase and lethality of COVID-19 in Italy, Xiao identified in his research that a slight increase in long-term exposure to PM -2.5 leads to a strong ink Similarly, in his report, Shi et al. (2020) identified the role of temperature in the dynamics of the outbreak of COVID-19 in China. In any event, these are observational studies, so it would be important for further evaluation to be carried out, especially from an etiological point of view (Setti et al. 2020b; Anjum et al. 2020 ). There is a superior necessity to carry out parallel research in emerging countries in Asia and Africa, where air pollution is at an extreme degree, and to equate the findings with a study established in developed countries to further improve the investigation. Furthermore, what has therefore appeared is that a deep resemblance is needed in the control of air pollution, in the scrupulous of the PM -10 , which did not diminish dramatically throughout the lockdown. In this information, the efficacy of the steps adopted at the local level to mitigate air emissions can be checked by monitoring. In addition, because of the sensitivity of environmental threats, the pandemic could provide an incentive to contrast group aggression (Carducci et al. 2019; Dettori et al. 2021; Anjum et al. 2020 ).

Regarding ambient emissions and both the dissemination of the virus in the development of eruptions and the rate of container casualty, this investigation tends to establish rudiments to support the presence of a correlation among air effluence, referred to as PM 10 particulate matter and COVID-19 death.

In 2020, some impressive circumstances could have worsened the environmental situation throughout Pakistan, especially in the area of KPK, making Pakistan a prominent country in South Asia. The spatial allocation of SMR-related COVID-19 deaths shows correlations in the spatial models tried, particularly with particles, as shown by multivariate regression. Our findings on air quality and the risk of corona virus mortality are in line with the findings of Dettori et al. (2021) , Ma et al. (2020), and Xiao Wu (2020) . Conceptually, our study results are the same as those previously done for the USA (Xiao Wu 2020), for China , and for Italy (Dettori et al. 2021 ) but differ in the point table as all air pollutants seem to be significant. Pakistan is an emerging Asian economy where the indices of air quality are high to the extent of risk. Pakistan is second among the dirty countries in the world, according to WHO (2020) . Because of this, all air pollutant variables such as NO 2 , PM 2.5 , and PM 10 tend to be statistically important and play a significant role in the likelihood of COVID-19 mortality. In Pakistan, there are solid reasons for considering PM 10 as an independent air quality measure. The only air pollution indicator on which the data can be available. Second, PM 10 has a good association with NO 2 and PM 2.5 to compare the findings with those obtained for advanced countries where air quality is at normal levels.

With respect to policy recommendations, we should reinforce the need to proceed in support of measures aimed at reducing atmospheric emissions by speeding up current plans and initiatives addressing all sources of ambient pollution: industry, domestic heating, and transport. While Pakistan is in the right direction with respect to reducing the country's air pollution, it will take time to normalize the extent of the serious danger of air pollution. The sitting government is highly worried about the climate. In this regard, in his first term of power in Khyber Pakhtunkhwa, the new Prime Minister launched the Billion Tree Planning Project (BTPP). While in his current term, the BTP project has been upgraded to 10 billion nationally. Imran Khan, Pakistan's prime minister, made the nation realize that the sitting government made legal regulations on air quality in the country without a sound atmosphere and barred all illegal brick manufacturing factories for this reason. Indeed, these actions taken by Pakistan's government would help the environment and reduce air pollution in the long term. Additionally, there is a greater need for investment in safe and green transport, and green building technologies should also be increased, beginning even with rapidly relevant steps, i.e., road cleaning, pollution-eating paints, façades and plants, paying special attention to housing sanitary facilities.

Exploration of this research work uncover many research questions. We are going to explain some of them here. First, there is a greater potential to replicate similar research work in Europe and Africa to investigate the risk of mortality and the spread of COVID-19 and compare the results with the one obtained from this study. Second, age and season variation are very influential variables playing an important role in the spread of COVID-19 and its updated version. One may consider it as point of future research work. Apart from this, there is a greater need to investigate the effectiveness of corona virus vaccine and its classification. Additionally, techniques for optimum allocation of corona virus vaccine in a short span of time are also the need of the present situation.

Abbreviations COVID-19, Corona virus; CP, Cycle paths; PA, Pedestrianized road surface (m 2 /inhabitant); UGS, Urban green spaces; Moto, Motorcycles; SMR, Specific mortality rate; BT, Heavy buses and trucks; FATA, Federal Administrated Tribal Area; NO 2 , Nitrogen dioxide (average yearly value (μg/mc)); PM 2.5 , Particulate matter ≤ 2.5 μm (average yearly value (μg/mc)); PM 10 , Particulate matter ≤ 10 μm (average yearly value (μg/mc)); WHO, World Health Organization; LISA, Local indicator spatial autocorrelation; BTP, Billion Tree Plantation Program

Author contribution NA.

Funding This research was financially supported by the Provincial Government of Khyber Pakhtunkhwa, Pakistan, under Grant. No. KP/ HE/2020/1171901.

The data used in this analysis were extracted from numerous sources. Data related to COVID-19, such as a number of confirmed cases, deaths, active cases, and area recoveries, was collected from the Ministry of National Health Services (MNHS 2020) available online at http://www.covid.gov.pk. Real-time air quality monitoring project (AQICN 2020) available at https://aqicn.org/country/pakistan/, Ministry of Environment, Government of Pakistan available at http://www. adaptation-undp.org., and air pollution data were collected from Pakistan Ministry of Climate Change available at http://www.mocc.gov. pk. All the results reported in this research were carried out in Rcomputational environment.

Ethics approval and consent to participate This article does not contain any studies with human applicants or animals accomplished by any of the authors.

Competing interests The authors declare no competing interests.

Consent for publication NA.

An emerged challenge of air pollution and ever-increasing particulate matter in Pakistan; a critical review

Local indicators of spatial association-LISA

Retrieved from air pollution in Pakistan real-time air quality monitoring project

Environment and health: risk perception and its determinants among Italian university students

The role of air pollution (PM and NO2) in COVID-19 spread and lethality: a systematic review

Air pollutants and risk of death due to COVID-19 in Italy

Environmental risks perception among citizens living near industrial plants: a cross-sectional study

COVID-19 and living space challenge. Well-being and public health recommendations for a healthy, safe, and sustainable housing

Regional air pollution persistence links to COVID-19 infection zoning

Case fatality rate analysis of Italian COVID-19 outbreak

Effects of air pollution on human health and practical measures for prevention in Iran

XIV Rapporto Qualit'a dell'ambiente urbano

Coronavirus disease 2019 (COVID-19): current status and future perspectives

Effects of temperature variation and humidity on the death of COVID-19 in Wuhan

The effects of regional climatic condition on the spread of COVID-19 at global scale

Variations of PM2.5, PM10 mass concentration and health assessment in Islamabad

Retrieved from COVID-19 health advisory platform by Ministry of National Health Services Regulations and Coordination

Geographical analyses of COVID-19's spreading contagion in the challenge of global health risks

The effects of air pollution and meteorological factors on measles cases in Lanzhou

Potential role of particulate matter in the spreading of COVID-19 in Northern Italy: first observational study based on initial epidemic diffusion

Searching for SARS-COV-2 on particulate matter : a possible early indicator of COVID-19 epidemic recurrence

Impact of temperature on the dynamics of the COVID-19 outbreak in China

A shortcut method for calculating the 95 percent confidence interval of the standardized mortality ratio

Geographical overlap between alpha-1 antitrypsin deficiency and COVID-19 infection in Italy: casual or causal?

COVID-19 outbreak: current scenario of Pakistan

Air quality guidelines for particulate matter, ozone, nitrogen dioxide and sulfur dioxide

Air pollution

Risk factors associated with acute respiratory distress syndrome and death in patients with coronavirus disease 2019 pneumonia in Wuhan, China

Retrieved from COVID-19 PM2.5 Exposure to air pollution and COVID-19 mortality in the United States: a nationwide

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations