key: cord-0750865-ita1krnz
authors: Broomandi, Parya; Crape, Byron; Jahanbakhshi, Ali; Janatian, Nasime; Nikfal, Amirhossein; Tamjidi, Mahsa; Kim, Jong R.; Middleton, Nick; Karaca, Ferhat
title: Assessment of the association between dust storms and COVID-19 infection rate in southwest Iran
date: 2022-01-21
journal: Environ Sci Pollut Res Int
DOI: 10.1007/s11356-021-18195-7
sha: 2eedf26ed3ca23f51a42718b85ca1e78720688ca
doc_id: 750865
cord_uid: ita1krnz

This study assesses a plausible correlation between a dust intrusion episode and a daily increase in COVID-19 cases. A surge in COVID-19 cases was observed a few days after a Middle East Dust (MED) event that peaked on 25(th) April 2020 in southwest Iran. To investigate potential causal factors for the spike in number of cases, cross-correlations between daily combined aerosol optical depths (AODs) and confirmed cases were computed for Khuzestan, Iran. Additionally, atmospheric stability data time series were assessed by covering before, during, and after dust intrusion, producing four statistically clustered distinct city groups. Groups 1 and 2 had different peak lag times of 10 and 4-5 days, respectively. Since there were statistically significant associations between AOD levels and confirmed cases in both groups, dust incursion may have increased population susceptibility to COVID-19 disease. Group 3 was utilized as a control group with neither a significant level of dust incursion during the episodic period nor any significant associations. Group 4 cities, which experienced high dust incursion levels, showed no significant correlation with confirmed case count increases. Random Forest Analysis assessed the influence of wind speed and AOD, showing relative importance of 0.31 and 0.23 on the daily increase percent of confirmed cases, respectively. This study may serve as a reference for better understanding and predicting factors affecting COVID-19 transmission and diffusion routes, focusing on the role of MED intrusions. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11356-021-18195-7.

Natural pathways are the primary sources of airborne particles on Earth, including marine aerosols and dust from arid areas Solomon et al. 2007 ). Levels of atmospheric particulate matter (PM) exceeding air quality standards and health guidelines are transported from deserts to inhabited areas under certain weather conditions (Marsham et al. 2013; Shahsavani et al. 2020) . For example, Shamal winds and cyclones cover the Middle East with dust clouds, leading to the deterioration of atmospheric conditions during the summer (Alizadeh-Choobari et al. 2016) . PM levels during such Middle East Dust (MED) events can be very high in western Iran: in the city of Ahvaz, a level of >9000 μg/m 3 was recorded on 27 th January 2017, and 5338 μg/m 3 on 3 rd June 2010; and at Sanandaj a PM 10 level of 5616 μg/m 3 was recorded on 5th July 2009 (Rashki et al. 2021; Salmabadi et al. 2020; Shahsavani et al. 2012) . The primary sources of dust particles for these events were located in Northern Iraq and along the Iraq-Syria border. In this arid climate, the vulnerability level of wind erosion is increased due to reduced vegetation coverage and major soil disturbance linked to human activities such as conflicts and land-use dynamics (Modarres et al. 2017) . Meteorological elements such as reductions in precipitation, low soil moisture level, high temperature, wind speed, and relative humidity (RH) intensify these MED events (Klingmüller et al. 2016) .

Airborne particles are linked to long-term and shortterm impacts on human health, such as cardiopulmonary and respiratory diseases (Middleton 2020) . The American Cancer Society recorded 8-18% more cardiovascularrelated deaths due to an increase of 10.0 μg/m 3 in PM 2.5 for 552,000 participants over 16 years in from US metropolitan areas (Fromme et al. 2008; Goudarzi et al. 2018; Niu et al. 2010) . It is worth noting that health conditions caused by natural PM (mineral dust) may differ from anthropogenic PM WHO 2007) . Some epidemiologists observed higher mortality rates on higher PM level dusty days than non-dusty days (Neophytou et al. 2013; Sandra et al. 2011) .

In the Iranian cities of Ahvaz and Kermanshah, some studies have shown increases in respiratory hospitalizations and cardiorespiratory mortality during dust storms (Delangizan and Jafari Motlagh 2013; Geravandi et al. 2017) . Ahvaz city in Khuzestan, which experiences major MED intrusions, experienced a 3.3% rise in daily deaths associated with a 10.0 μg/m 3 increment of PM 10 concentration during dusty days as compared to a 1.0% increase on nondusty days ). The authors found that any PM exposure during an MED intrusion may adversely impact human health in arid and downwind affected areas (receptors) .

During atmospheric transport, pollutants may be integrated into dust clouds and reach areas (receptors) downwind. Local gaseous pollutants (e.g., SO 2 ) can also condense on dust particles (Rodriguez et al. 2011) . During dust episodes, the accumulation, transport, and condensation of pollutants could make dust more dangerous to human health. Dust storms can transport bioaerosols on a large scale and affect ecosystems and populations downwind. Dust and sandstorms are capable of introducing many foreign microorganisms into the global system. Microorganisms found in desert dust are typically very resilient and highly resistant to desiccation, temperature extremes, high salinity, and exposure to ultraviolet radiation (Behzad et al. 2018; Broomandi and Rashidi 2018; Gonzalez-Toril et al. 2020; Goudarzi et al. 2014; Kai et al. 2017; Neisi et al. 2019; Nourmoradi et al. 2015; Reche et al. 2018; Weil et al. 2017) .

Particulate matter can provide a substrate, allowing a virus to maintain a long-term (hours or days) presence in the ambient atmosphere. Environmental parameters are the key players in viral inactivation. For example, high temperature and solar radiation may accelerate inactivation rates of viruses, while high humidity can promote diffusion rate Després et al. 2012; Wigginton and Boehm 2020) . Recent studies noted an association between ambient PM concentrations and presence of viruses among exposed populations (Chen et al. 2010) .

In China, a positive correlation (p < 0.001) was found between childrens' infection rates of RSV (respiratory syncytial virus) with PM 2.5 and PM 10 concentration. RSV deeper penetration in the respiratory system is reportedly promoted by particle-based transport (Ye et al. 2016) . Another investigation found a positive correlation between viral disease and exposure to high PM 2.5 levels in China (Chen et al. 2017) . Zhu et al. (2020) evaluated the association of short-term exposure to higher levels of air pollution with increases in risk of infection. The authors of this paper showed a substantial positive correlation between four days of exposure (23 rd January 2020 to 29 th February 2020) to concentrations of PM 2.5 and PM 10 and number of new daily positive cases in 120 cities in China. They also found that a 10.0 μg/m 3 increase in PM 2.5 and PM 10 levels was positively correlated with a 1.76% (95% CI: 0.89 to 2.63) and 2.24% (95% CI: 1.02 to 3.46) rise in confirmed daily positive cases . A study based in the USA found that longterm exposure to PM 2.5 was positively correlated with increased mortality risk of COVID-19 (Wu et al. 2020) . They reported that a 1.0 μg/m 3 incremental increase of PM 2.5 positively correlated with an increase in the SARS-CoV-2 mortality rate of 8% (Wu et al. 2020) .

In Istanbul, Turkey, a study of the associations of meteorological parameters and air quality with daily COVID-19 case numbers found that air quality indices, temperature, and relative humidity were associated with the spread of disease in the population (Shahzad et al. 2021) . Temperature and NO 2 concentrations were the primary factors associated with a substantial spike in the infection and death rates of COVID-19 in Istanbul. Moreover, they found a significant association for PM 10 , PM 2.5 , O 3 , and relative humidity with COVID-19 transmission (Shahzad et al. 2021) .

A study conducted in the state of NJ, USA, explored correlations between meteorological parameters and new COVID-19 case counts, finding that temperature was negatively correlated, and relative humidity and air quality index were positively correlated with new case counts. PM 2.5 , population density, human development index, and health security index were also correlated with spread of the disease (Doğan et al. 2020) . Moreover, autoregressive distributed lag (ARDL) analysis showed that relative humidity, air quality, and infections had lagged effects with COVID-19 spread across New Jersey (Doğan et al. 2020) .

Viruses such as SARS-CoV-1 are likely to spread in the air within the formation of tiny liquid droplets and then could freely transport with their viral content over tens of meters from their origin in the air Morawska and Cao 2020; Rahmani et al. 2020) . Researchers have concluded that COVID-19 can spread person-to-person through the air (Bourouiba 2020; Domingo et al. 2020; Hadei et al. 2020; Jayaweera et al. 2020; Morawska and Cao 2020) . Recently published studies suggested two mechanisms of direct and indirect aerosol transmission routes of SARS-CoV-2. With direct dispersion human-exhaled aerosols are transmitted in close proximity, while with indirect dispersion ambient particles act as cargo to spread the virus. Additionally, the presence of particulate matter can stimulate the expression of Transmembrane Serine Protease 2 (TMPRSS2) and Angiotensin-Converting Enzyme 2 (ACE-2), causing an increase in the SARS-CoV-2 binding sites and the facilitation of infection efficiency (Cao et al. 2021; Li et al. 2020; Paital and Agrawal 2020; Yao et al. 2020a) .

To the best of our knowledge, there is no study examining the association between MED dust storms and new case counts of SARS-CoV-2. This study aims to assess plausible links between MED dust intrusion and COVID-19 infection counts by assessing dust storm measures and disease occurence. This research is founded on the hypothesis of this possible relations in two ways: (1) the population becomes more susceptible to viral exposure during epidemics due to the damaging effect of MED events on the respiratory system, and (2) atmospheric clouds of dust can act as a carrier and increase atmospheric transportation of the virus (SARS-CoV-2). During the pandemic, a MED dust storm occurred on April 25-26, 2020, in Khuzestan, Iran, with the number of new cases in several cities increased dramatically followingr the event. To confirm the given hypotheses, the relationship between AOD levels and daily percentage increase in positive cases of COVID-19 are statistically examined in detail in 20 cities in Khuzestan, Iran.

Khuzestan province covers 63,213 km 2 in southwestern Iran and is home to 4.7 million inhabitants. It is located between 31°N and 32°N latitudes and 48°E and 49.5°E longitudes (Fig. 1) . The topographic altitudes change from 0 to 3740 m. The weather ranges from humid to arid. While the southern regions experience a tropical climate, the northern regions have a colder climate. Summertime begins in April and lasts to September, and wintertime lasts from October to March. During summer, the annual mean maximum temperatures are close to 50 °C, with the minimum value occurring in March (close to 9°C). The annual precipitation levels range 995-1100 mm in the north and 150-256 mm in the south. Roughly 70% of annual rainfall is from February to April. The annual amount of evaporation is 2000-4000 mm. The direction of the prevailing wind speed is from west to east and from northwest to southeast (Zarasvandi et al. 2011) .

There are two main origins of dust storms in Khuzestan, including neighboring countries such as Iraq and Saudi Arabia (Mohammadpour et al. 2020; Salmabadi et al. 2020 ) and susceptible to land to wind erosion in the west of Khuzestan and desiccation of water bodies such as Karkheh river and Hourazim wetland (Malamiri et al. 2019; MalAmiri et al. 2022) . According to the long-term analyses of meteorological-dust data, the highest levels of dusty days for both Iran and the Middle East occurred in the cities of Khuzestan province, Iran (Salmabadi et al. 2020) , with a maximum (Klingmüller et al. 2016) , and the continued spike of dust-AODs in downwind areas such as Khuzestan (Hamzeh et al. 2021a) . Dusty days are more frequent in the late spring and summer, with dust storms having horizontal visibility below 100 m most common in June and July (Mohammadpour et al. 2020; Salmabadi et al. 2020 ). The number of dusty days also increased in February and March due to the frontal dust storms over the Iraqi plains, with the highest impact on the southwest of Iran (Gholamzade Ledari et al. 2020; Hamidi 2019; Hamzeh et al. 2021b ).

Regional dust storms and other phenomena such as stable weather conditions and inversions can increase ambient PM levels. In our study, dusty days due to regional dust intrusions were statistically separated from non-dusty days by analyzing air pollutant levels. The threshold value of 100 μg/m 3 for PM 10 concentration was applied to distinguish dusty days from regular non-dusty days, as utilized by other researchers (Escudero et al. 2007; Givehchi et al. 2013 ). Multi-pollutant behavior analyses were used to distinguish days influenced by either inversion/atmospheric conditions or regional dust storms. During dust intrusion, notable increases in the PM levels are expected with no noteworthy increment of gaseous pollutants' concentrations, mainly formed from primary combustion like Carbon Monoxide (CO). To separate episodic days, the duration of notably elevated concentration of PM 10 , i.e., the daily average concentration with at least one STD (standard deviation) above the mean value, was chosen. During the study, when the daily carbon monoxide concentrations were at least one STD below the mean, the day's measurements were included in the analysis, as applied by Givehchi et al. (2013) . Among the determined dusty days in Khuzestan province, the dust storm outbreaks causing elevated PM 10 concentrations (daily PM 10 concentration above 200 μg/m 3 ) were selected for includion in the study. A statistical procedure by Givehchi et al. (2013) was followed to separate the amount of PM 10 increases caused by an intensive dust intrusion from local anthropogenic sources. Firstly, the initial correlation coefficients (R 2 with p values < 0.0001) of hourly-based PM 10 concentrations in the studied time frame among all pairs of available air quality stations were calculated during dust episodes. The pairs with a remarkably lower correlation coefficient (<0.6) were excluded from the studied cities since it is expected that they would be primarily affected by local emission sources, even during the dust intrusion.

In the current study, the hourly PM10 data were obtained from the National Air Quality Information System (Iranian meteorological organization). Also, the daily data of number of new cases and deaths were obtained from the database of the official coronavirus site enabled by the Ministry of Health (Food and Drug Association).

Since cross-correlation helps to evaluate the similarity between two-time series, it was implemented to explore the relationship between combined AOD and the daily percentage increase of SARS-CoV-2 infection (Wang et al. 2019 ). In its application, one of the main limitations is the different distribution characteristics of data sets correlated. During the growth or spread period, the pandemic usually follows an increasing trend; however, the AOD data is more random with some short-term autocorrelations due to the nature of local emissions and dust incursion characteristics. To provide correct correlation results, it is necessary to eliminate the typical growth characteristics of the COVID cases. In this study, a transformation procedure was used to normalize the virus infection growth characteristics, which typically follows an exponential function, while log transformation or original data was sufficient to have a normal distribution for AOD data. The percentage of daily increase values COVID infection numbers were calculated regarding the preceding day only (e.g., =

, where COVID t is the number of new cases in a day, and COVID t−1 is the number of new cases in a preceding day).

The positive lags denote correlations among the AOD at time t with the increasing percentage of COVID-19 at time t+1, t+2. The lag with the highest correlation coefficient was identified when the combined AOD and daily percentage increase of COVID-19 infection match the best.

COVID-19 infection data were received from the Iranian Ministry of Health and Medical Education. Cross-correlations between the daily combined AOD and the increase in the percentage of COVID-19 infection numbers were calculated for 20 cities in the province starting from 20 th April 2020 to 5 th May 2020. Before 20 th April, the infection data were not reported for each city, so it was impossible to include the data of earlier pandemic periods. The studied cities are listed in Table S1 in the supplementary material.

To investigate the possible association between weather conditions and vulnerability to infection during dusty days, it was decided to assign each city into a cluster using t h e K-m e a n s c l u st e r i n g m et h o d . A d a t a s et D = {O 1 , O 2 , …, O n } has n objects (instances) each O i is an object as a p-dimensional explanatory variable in the dataset (e.g., Pasquil Stability Class (PSC), Planetary Boundary Layer Height (PBLH), Wind speed, AOD, Duration of MED, and Population Density). The studied domain needs division into a combination of k clusters, representing a vector CKM = {C 1 , C 2 , …, C k } with denoted centroids by μ = {m 1 , m 2 , …, m k }. The beginning step is the randomly assigning of k points as cluster centers. For each data point, O i is a distance to each cluster centroid m j and was calculated through one distance metric like Manhattan, Euclidean, Minkowski, or Chebyshev distance. The argminj dist (O i m j ) is used to figure the nearest cluster for the respective instance to be assigned. The new midpoint of clusters is calculated by

where n j , m j , and O i the number of objects in cluster j, the centroid of cluster j, and the dataset instances, respectively. Iteratively, this process proceeds to where no data point alters cluster membership (Grace et al. 2016; Tüysüzoğlu et al. 2019 ).

Random Forest Analysis (RFA) was used in the current study to evaluate the relative importance of parameters including (AOD, temperature, PBLH, surface pressure, relative humidity, and wind speed) affecting the daily increase percent of COVID-19 infection. RFA helped in the selection of the most critical predictors affecting the response based on a classification from zero (no importance) to one (highest importance) (Molinaro et al. 2011 ).

Deep blue AOD (Aerosol Optical Depth) data from the Aqua-MODIS 550 nm Collection 6 MODIS and Terra-MODIS 550 nm Collection 6 MODIS were employed to investigate the changes in daily PM levels in the atmospheric column over 20 cities. At a 10-km resolution, the standard MODIS Level 2 (L2) AOD products are distributed. MYD04_L2 data from the Aqua-MODIS 550 nm L2 Aerosol Products and MOD04_L2 data from Terra-MODIS 550 nm L2 Aerosol Products for the study duration were extracted from LAADS (MODIS L1 Level 1 Atmosphere Archive and Distribution System) (http:// ladsw eb. nascom. nasa. gov/).

For improving the AOD spatial coverage in our study, DB_ DT AOD from Terra and Aqua (MODIS-carrying satellites) were combined by averaging. The method of averaging the Aqua and Terra measurements successfully estimates the daily AOD average mentioned in other studies (Lee et al. 2011; Nabavi et al. 2018 ). When either Aqua or Terra is missing, AOD with the help of simple linear regression can be estimated (Eqs. (1) and (2)) among daily Aqua and Terra AOD values (Hu et al. 2013; Nabavi et al. 2018) . and τ are estimated and available DB_DT AOD, respectively.

ERA5 reanalysis data, produced by C3S at ECMWF, as the current atmospheric reanalysis and regard to the 2016 version of IFS (Integrated Forecasting System), was employed to investigate the probable impact of meteorological data (including daily mean surface pressure (KPa), daily mean wind speed (m/s), daily mean relative humidity (%), daily mean temperature (°C), daily mean planetary boundary layer height (m)) on the increment of SARS-CoV-2 infection in Khuzestan province, Iran.

Stability time series by covering before-during-after dust intrusion for the studied cities, including Pasquil stability class (PSC) and vertical mixing Planetary Boundary Layer Height (PBLH), were downloaded from https:// www. ready. noaa. gov/ READY amet. php and were studied for any significant commonalities or differences between cities based on the dust incursion times, respectively.

The Visible/IR images from Spinning Enhanced Visible and Infrared Imager (SEVIRI), carried by MSG (Meteosat (1)

, if Terra is missing

, if Aqua is missing

Second Generation) satellites, are employed in our study to investigate dust storms. For continuous dust monitoring, the EUMETSAT (European Organisation for the Exploitation of Meteorological Satellites) recommends RGB images (available free of charge every hour: http:// www. eumet sat. int/ Home/ Main/ Image_ Galle ry/ Real_Time_Imagery/index. htm). The infrared channel data from SEVIRI is based on RGB image compositions, where dusty pixels show up in pink colors, and is used to monitor the dust event evolution during both day and night in deserts (Martínez et al. 2009 ).

For analyzing the sources and trajectories of dust particles, the Online HYSPLIT version 4 (Hybrid Single Lagrangian Integrated Trajectory) model was employed with 1°×1° resolution meteorological data (Draxler and Hess 1997 ). An HYSPLIT ensemble trajectory calculates the multiple pathways from one position by all-feasible deviations in X, Y, and Z. For the reduction of the model uncertainty, ensemble trajectories are used to compute all possible pathways (Draxler and Hess 1997) .

Identification of a regional dust invasion event during the pandemic growth period In Khuzestan, several air quality stations record hourly PM 10 concentrations. The recorded data from Ahvaz, Shush, Shushtar, Ramshir, Hamidiyeh, Abadan, Hendijan, Andimeshk, Dezful, and Omidiyeh were considered. Others were excluded due to the limited data coverage and unreliability with several hourly/daily gaps during the studied period (1 st April 2020-30 th April 2020). The resulting numbers of dusty days were 2, 2, 1, 1, 5, 4, 1, 1, and 3 for Shush, Ahvaz, Shushtar, Ramshir, Hamidiyeh, Abadan, Hendijan, Andimeshk, Dezful, and Omidiyeh, respectively during April 2020. Among the determined dusty days, one dust episode outbreak on 25 th -26 th April 2020, causing the elevated levels of PM 10 (daily PM 10 concentration above 200 μg/m 3 ), occurred and was selected for our research. To distinguish the increases of PM 10 caused by dust intrusion from local anthropogenic sources, the initial correlation coefficient (R 2 with p values less than 0.0001) of hourly PM 10 concentration among all pairs of air quality stations during dust outbreak was calculated. Correlation values above 0.75 showed that any of the studied stations were influenced by local sources during dust intrusion, and the amount of PM 10 increases during dust outbreak was caused by the dust intrusion. Fig. 2 illustrates the time series of daily measured PM 10 concentration and combined daily extracted deep blue AOD over the mentioned cities above.

Globally, the AErosol RObotic NETwork project (AER-ONET) provides distributed observations of spectral aerosol optical depth (AOD), precipitable water in diverse aerosol regimes, and inversion products (https:// aeron et. gsfc. nasa. gov/). PM 10 data (9 cities out of 20 cities), the correlation among the combined AOD data, and ground-based PM 10 measurements were acceptable with the R values ranging from 0.60 to 0.78 (all p values < 0.0001) (details of regression results presented in Supplementary Material). Both series in Fig. 2 represent a peak in PM 10 and AOD values during dust intrusion on 26 th April 2020.

To monitor MED outbreaks over land, their migration, and the corresponding changes in PM 10 levels in Khuzestan province, HYSPLIT backward trajectory modeling and the Visible/IR images of SEVIRI with a temporal resolution of 60 min and spatial resolution of 3 km×3 km were implemented in our study (Figs. 3 and 4) . Transport pathways of dust particles were tracked through 6-h time intervals up to 24 h before dust episodes reaching the study locations using HYSPLIT backward trajectories on 25 th April 2020. Fig. 3 shows the main corridors of dust transport originating from arid and semi-arid areas in northern Saudi Arabia and central parts of Iraq to Iran's west and southwest regions.

The Visible/IR images of SEVIRI are employed to support HYSPLIT results and identify dust storms (Fig. 4 ) that were consistent with the spatial and temporal coverage. In the thermal infrared part of the electromagnetic spectrum (8.7 to 12.0 μm), atmospheric dust can create a cooling anomaly in 'clear-sky' conditions (Hennen et al. 2019) . Relative RGB (Red/Green/Blue) beam strengths are rendered from inter-channel Brightness Temperature Differences (BTD), configured with specific limits (cf. Table 1 in Hennen et al. 2019 ) to distinguish 'thermally insulating' atmospheric components (such as soot and clouds from biomass burning) from atmospheric dust (Brindley et al. 2012; Hennen et al. 2019 ). In the 'Dust RGB' product, meteorological clouds appear as red or brown, the dust appears magenta or pink, and bare surfaces appear as white or blue (Fig. 3) (Brindley et al. 2012; Hennen et al. 2019 ). According to Fig. 4 , dust storms originated in central parts of Iraq and northern Saudi Arabia, gradually moved to the west and southwest of Iran, including Khuzestan province. Along with the results of applied statistical analyses to the pollutant concentrations, both HYSPLIT backward trajectory modeling and the Visible/IR images of SEVIRI confirmed that a dust storm started on 25 th and peaked on 26 th April 2020, similar to a prior case that covered large areas of Iraq and northern regions of Saudi Arabia that afterwards reached Khuzestan (Ashrafi et al. 2014; Martínez et al. 2009 ).

Previous studies have shown two different dust storm corridors to Khuzestan province in Iran: (a) W-E direction, from central regions of Iraq to the southwest and west of Iran, and (b) NW-SE direction, from eastern Syria and northwestern Iraq (Baghbanan et al. 2020; Broomandi and Bakhtiar Pour 2017; Cao et al. 2015; Zarasvandi et al. 2011) . During warm periods from March to September, Shamal winds are responsible for dust transportation from eastern Syria, western regions of Iraq, and Jordan to Khuzestan province. The main dust events over Iran come from the semi-arid and arid region of the dust belt, mainly northern, eastern, and central parts of Syria and Iraq, and the northern part of Saudi Arabia (Aliabadi et al. 2015; Broomandi and Bakhtiar Pour 2017; Cao et al. 2015) .

The first confirmed COVID-19 case was reported in Iran in mid-February 2020. As a result, lockdown and precautionary measures were implemented from 21 st March to 21 st April in 2020 in the whole country. By the beginning of June 2020, the total number of infected, dead, and recovered people in Iran are 157,562, 7942, and 123,077, respectively, while in Khuzestan are 15,988, 595, and 14,591, respectively. It is worth mentioning that the Khuzestan population ratio to Iran is about 5.9%. However, the rates of total infected, death, and recovered people in this city are around 10%, 7.5%, and 12%.

According to Fig. 5 , the initial number of Khuzestan cases had been slightly increasing, while the numbers were rapidly decreasing in the whole country. Subsequently, however, case numbers massively increased in Khuzestan, while the numbers were decreasing in the country, to the extent that the increase rate in Khuzestan outstripped that of the country starting from 27 th April 2020. The infection ratio of new cases in the province to the country was 9.6 %, 13.2%, 15.0%, 19.0%, and 28.0% on 25 th April (the day of the end of the studied dust storm event), 27 th April (2 days after the dust intrusion), 29 th April (4 days after the dust intrusion), 30 th April (5 days after the dust intrusion), and 2 nd May (7 days after the dust intrusion), respectively. Right after a 7-day lag, the newly reported number of cases tripled in Fig. 2 The time series of (A) combined daily extracted deep blue AOD, and (B) measured daily PM 10 concentration and in Khuzestan, Iran during 1 st April 2020 to 30 th April 2020 this province, indicating a significant abnormal increase in the daily infection rate. Additionally, Fig. 5 shows a slight increase in the overall death rate in Iran, while there is a rapid increase in death rate in Khuzestan.

Meteorological variables such as temperature, planetary boundary layer height (PBLH), surface pressure, RH (relative humidity), wind speed, and direction could impact the formation and dispersion of dust storms (Miri et al. 2017; Pirsaheb et al. 2016 ). The Pearson correlation coefficients analyses show that daily combined AOD in studied cities of Khuzestan, Iran had a statistically significant positive correlation with daily average wind speed (R 2 =0.40, p value <0.05), temperature (R 2 =0.32, p value <0.05), PBLH (R 2 =0.13, p value <0.05(. However, no statistically significant correlation was observed among daily combined AOD and relative humidity and surface pressure. The elevated AOD levels can be due to the (a) any increase in air temperature; leading to a reduction in air humidity and increase in dust acceptance in air, (b) any decrease in relative humidity; preventing the precipitation phenomena and the act of washing and dust descending, (c) any increase in the wind speed; facilitating the transmission of dust particles from origins to downwind areas, (d) any reduction in surface pressure; lowering the chance of precipitation, air turbulence, and dust movement, (e) and any increase in PBLH; causing effective dispersion of dust particles in the vertical direction (Ashrafi et al. 2014; Guan et al. 2017; Namdari et al. 2018; Pirsaheb et al. 2016; Rashki et al. 2015; Zhu et al. 2018 ).

The relationship between virus spread and atmospheric dust is not straightforward, and possible latent relations with other independent factors, meteorological factors, and atmospheric stability characteristics should be considered. To find any significant commonalities or differences among the cities during the dust incursion times (1 st April 2020-30 th April 2020), atmospheric stability data time series (e.g., Pasquil Stability Class (PSC) and Planetary Boundary Layer Height (PBLH)) studied for before, during, and after dust intrusion periods. Table 1 illustrates the dust event's duration, the dust cloud arrival time, PSC, PBLH, Wind Speed, AOD, and Population Density. The highest PBLH value was observed in Andimeshk (2526 m) with a stability class of B (Moderately unstable conditions), while in Mahshahr, the lowest value of 216 m occurred with a stability The K-means clustering algorithm clustered the cities to characterize their features or attributes (Shafiee et al. 2016a, b; Shobha and Asha 2017) , and it measures different aspects of the cities' atmospheric stability characteristics. The clusters were extracted by considering attribute values of PSC, PBLH, AOD, wind speed, population density, and the duration of dust events. Table 2 summarizes the cluster analysis in Khuzestan, Iran. RFA analysis showed the critical role of wind speed and AOD with the relative importance of 0.31 and 0.23, respectively, influencing the daily increase percent of COVID-19 infection ( Table 3) .

The cities were divided into three clusters based on their similarities. Cluster 1 includes Hendijan, Mahshar, Ramshir, Omidiyeh, and Behbahan cities with PBLH values below 1000 m (excluding Behbahan) with the stability classes of C and D and AOD ranging of 0.50-0.58. In cluster 1, the wind speed was relatively higher than the other two clusters (above 8.0 m/s). Cluster 2 includes Ahvaz, Shushtar, Abadan, Karun, Bavi, Andimeshk, and Shush cities with PBLH values above 2000 m (excluding Abadan), stability classes of A and B, and AOD ranging of 0.60-0.86. In cluster 2, the wind speed varied between 5 and 9 m/s. Cluster 3 includes Rahmhormoz, Khoramshahr, Dezful, Izeh, Baghe Malek, Masjed Soleyman, Hoveyzeh, and Hamidiyeh with PBLH values below 1000 m, with stability class of C and D, and AOD ranging of 0.71-0.87. In cluster 3, the wind speed was relatively lower than others, varying between 1.0 and 8.0 m/s. The main differences in the observed relationships of the virus spread and the dust event for each city type are discussed in the following sections.

The SARS-CoV-2 virus has an incubation time ranging 2-14 days (Guan et al. 2020; Lauer et al. 2020) ; thus, any possible impacts of the infection spread following a dust event would be observed after several days. Cross-correlation is a method to evaluate two-time series if one is more correlated to a lagged measure of the other, helping to identify at which lag time the two-time series are most strongly correlated. In this study, several cross-correlations up to +10 days lag were applied to the studied cities. Due to a lack of April 2020 at 1:00 am, (B) 25 th April 2020 at 10:00 am, (C) 25 th April 2020 at 3:00 pm, (D) 25 th April 2020 at 7:00 pm, (E) 25 th April 2020 at 10:00 pm, and (F) 26 th April 2020 at 8:00 am data availability, analysis with a higher lag level, e.g., 14 days or more, was impossible. The event data sets consist of measurments for (1) the prior week before the event, (2) during the event (1-3 days), and (3) after the event, including all possible recorded days during the study period.

In Table 4 , cities were categorized into four groups, based on calculated cross-correlations among the daily combined AOD and increases of COVID-19 infection numbers (%) from 20/04/2020 to 09/05/2020. In group 1, cities of Masjed Soleyman, Khoramshahr, and Izeh with R values above 0.70 showed that dust incursions might have a significant impact on the spread of COVID cases in these cities since the peak correlations come around 10-day lag that is the typical incubation period (Bontempi 2020) . Group 2 includes Shushtar, Behbahan, Andimeshk, Abadan, Bagh-e Malek, Hamidiyeh, Ahvaz, Bavi, Dezful, and Ramshir, showing R values ranging from 0.40 to 0.68. In this group, the AOD levels are high (except Behbahan city with relatively high correlations). The peak correlations were observed around 4-5 days lag times, which might be related to the weakening of the infected persons' respiratory and immune systems (Bontempi 2020) . Group 3 includes Omidiyeh, Mahshahr, Hendijan, and Karun, with R values ranging from 0.31 to 0.43. Most of the AOD levels in this group were lower, so they were utilized as our control group cities. These cities, having lower correlations for dust intrusions with new case counts, expectedly show no significant association with any lag. Finally, in group 4, the cities of Rahmhormoz, Hoveyzeh, and Shush, with R values ranging from 0.23 to 0.34, exceptions to the hypothesis since they had a higher level of dust incursion, but no significant correlation with no identifiable confirmed COVID case count increases.

In reference to Fig. 1 and Table 2 , neighboring cities such as Bavi, Karun, and Shushtar had similar lag values. Correspondingly, the cities clustered together (Table 2) based on atmospheric stability and dust storm characteristics (e.g. Ahvaz, Shushtar, Andimeshk, and Abadan in Cluster 2 and Ramshir, Hendijan, Mahshahr, and Omidiyeh Hendijan in Cluster 1) fell into the same groups in The reason behind this different grouping (based on crosscorrelation among combined AOD and daily increase percent of COVID-19 infection) might be due to poor medical surveillance such as insufficient testing in Shush, Rahmhormoz, and Hoveyzeh (grouped in G.4) . Despite the similar atmospheric stability and dust storm characteristics, they were exceptions to the hypothesis, not showing a significant correlation.

Another notable exception was observed in Cluster 3. Khorramshahr (grouped in G.1) and Dezful & Bagh-e Malek & Hamidiyeh (grouped in G.2) were clustered together because they shared similar cross-correlation values. All cities in cluster 3 share similar atmospheric stability and dust storm properties, while they have different elevations. Khorramshahr is located near the sea and at a lower elevation, about 3 m above sea level. It is speculated that it was more vulnerable to the dust storm intrusion than Dezful, Bagh-e Malek, and Hamidiyeh, located farther inland, at 921 m, 734 m, and 52 m above sea level, respectively, even though their AOD levels were slightly higher than Khoramshahr.

There is sparse research concerning COVID-19 epidemic phenomena and their possible association with PM levels and diffusion. Recently, studies showed a meaningful linkage among Saharan dust intrusion and observed COVID-19 outbreak in early 2020 on the Canary Islands, Spain, found that higher levels of PM 2.5 caused by Saharan dust intrusions and air temperature inversions can boost severe COVID-19 outbreaks (Rohrer et al. 2020) . Another study conducted in Spain, which examined the impact of Saharan dust intrusions on the incidence and severity of COVID-19, found that dust intrusion had an additional influence on COVID-19 incidences and hospital admission rates (Linares et al. 2021) . Desert dust events can increase the levels of ambient particulate matter (PM 2.5 ) as well as being a vector for fungal illnesses, which can intensify the mortality and morbidity related to SARS-CoV-2 (Verweij et al. 2020) . The overburdening of hospitals and health services may also be linked to PM 2.5 peak levels and particular meteorological conditions favorable for spreading and enhancing the virulence of COVID-19 (Rohrer et al. 2020) . Additionally, the observed pattern of COVID-19 cases concentrated in the 30-50N latitude area suggests that dust carrying the virus may have been spread by a circum-global northern sub-tropical jet-stream in the high-altitude troposphere over northern parts of China in early Spring, 2020 (Wickramasinghe et al. 2020) .

The authors speculated that dust particles responsible for carrying the COVID-19 virus were transported from considerable virus sources generated in Wuhan, China, toward southern parts of the USA, consequently across the Atlantic to Portugal and further states to the east (Wickramasinghe et al. 2020) . Their study suggested that the primary deposition of the dust particles carrying the virus depends on the interaction between the jet-stream and regional meteorological systems, leading to COVID-19 outbreaks in different countries along this latitude belt (Wickramasinghe et al. 2020) . The case of Brazil, which exceptionally is outside the 30-50N latitude belt, was suggested to being caused by the Azores cyclonic system affecting part of the jet-stream west of Portugal into the south-westerly trade winds with these winds entering Brazil in Spring 2020 (Wickramasinghe et al. 2020) . It is recommended not only to monitor the prevalence of the virus in the future but also to consider the occurrence of meteorological conditions that can cause unexpected and uncontrollable SARS-CoV-2 outbreaks.

Some other researchers suggested that transmission of infectious diseases can occur rapidly with increased air pollution levels (Bontempi 2020; Chen et al. 2010; Peng et al. 2020; Wu et al. 2020; Yang et al. 2020; Ye et al. 2016; Zhu et al. 2020) . Yang et al. (2020) discussed the susceptibility of people to viral disease due to exposure to high levels of air pollution with peaks of particulate matter (PM) concentration. Being exposed to elevated levels of PM reduces the effectiveness of the immune system and allows microorganisms to become more invasive . Current studies show that atmospheric particles could up-regulate the SARS-CoV-2 receptor and related protease (Cao et al. 2021) . Any exposure to aerosol pollution can be related to different cardiovascular and respiratory diseases (Pun et al. 2017) via various mechanisms such as the up-regulation of Angiotensin-Converting Enzyme 2 (ACE-2) and Transmembrane Serine Protease 2 (TMPRSS2) (Lin et al. 2018; Paital and Agrawal 2020) . The protective mechanism of ACE-2 upregulation helps the respiratory system maintain the Renin Angiotensin System (RAS) and decreases the inflammatory reaction (Ye and Liu 2019) . It is not only abundantly expressed in the lungs, but in the glandular cells of duodenal, gastric, and rectal epithelia of the infected people by COVID-19 (Cao et al. 2021; Paital and Agrawal 2020; Wu et al. 2020) . The ACE-2 acts as the primary receptor protein of SARS-CoV-2.

Additionally, the high affinity between the virus's synaptic glycoprotein and ACE-2 in host cell targets was reported (Vankadari and Wilce 2020) . Also, TMPRSS2 acts as a protease, cleaving viral spike protein and combining it with host cells to accelerate the infection (Kaur et al. 2021) . When particulate matter, especially PM 2.5 , attacks the human body, it can increase SARS-CoV-2 susceptibility for the human body by promoting the AEC-2 expression and its co-factor, TMPRSS2 . A recent in vivo study in mice showed a significant increase in the level of the AEC-2 expression after being exposed to PM 2.5 (Cao et al. 2021; Lin et al. 2018) . Statistical analysis suggests a significant association between PM concentrations, ACE-2 expression, and severity of COVID-19 infection (Cao et al. 2021; Paital and Agrawal 2020) . Therefore, it is recommended to conduct further studies paying particular attention to the level of air pollution in areas highly hit by COVID-19 and introduce suitable measures by decision-makers to decrease the level of air pollution. Bontempi (2020) reported notable registered positive cases in Italy (mainly in Brescia and Bergamo areas, Lombardy) after a PM 10 episode exceeding the set safety limit of 50 μg/m 3 for several days at the beginning of March 2020. However, a direct linkage between PM 10 peak and COVID-19 diffusion was not reported (Bontempi 2020) . Other related studies, summarized in Table 5 , supported the hypothesis that any increase in the spread of COVID-19 and related mortality could be due to poor air quality. Specifically, a positive association was observed among PM, focusing on PM2.5 and COVID-19 mortality or morbidity. Based on their studies, long-term exposure to a relatively high level of particulate matter may be responsible for the spread and pathogenicity of SARS-CoV-2.

This paper suggests that the captured impacts in our study can be related primarily to co-morbidities, inflammation, pre-existing Particulate Matter-related cellular damage, and up-regulation of ACE-2 & TMPRSS2 in the host cells (Cao et al. 2021; Ciencewicki and Jaspers 2007; Kim et al. 2020; Paital and Agrawal 2020; Pope 3rd et al. 2016; Tsai et al. 2019) . Previous studies also reported that air pollution exposure could negatively affect early responses of the immune system to the infection, which leads to later increases in worse prognosis and inflammation (Ciencewicki and Jaspers 2007; Conticini et al. 2020; Lambert et al. 2003 ), a possible alternative explanation for our findings.

A great concern is the presence of pathogenic microorganisms in the air and their transportation by ambient PM. It is essential to identify the potential impacts of airborne virus exposure caused by PM. During a MED event, these findings in Khuzestan, in combination with recently conducted studies of viral interactions with particulate matter, raise suspicion about the probable airborne transmission of COVID-19. As a result, any association between human health and air pollution, especially particulate matter, is vital to formulate positive strategies by policymakers to decrease atmospheric aerosols and potentially reduce the spread of future epidemic viruses and illnesses.

Unfortunately, during a pandemic, it is impossible to design a study and collect ideal quality temporal and spatial data to minimize sources of bias. As a result, our research has some limitations, and further investigations are required. We only focused on the correlations among AOD levels (representing inhalable particle levels in the studied cities) and the daily percentage increase of COVID-19 infection before, during, and after periods of a dust incursion event; this study did not focus on the association of air pollution on the COVID-19 pandemic. The data did not include medical history (no information on background health issues such as respiratory and cardiovascular conditions and diabetes), smoking status, or the age-specific or gender-specific confirmed COVID-19 cases, so no subgroup analyses on these factors could be conducted. We also could not include other factors in the analyses, such as indoor confinement duration, healthcare system capacity, and case identification procedures and practices. Also, our findings are not necessarily generalizable globally to other regions with dust intrusion events, given that we only had data in one region of Iran. Data on potential confounding factors such as the susceptibility of the population, surveillance data on respiratory infections, patterns of social relationships, public health measures, use of masks and sanitizers, hygiene, social distancing, virus resistance, mobility, urban density, socioeconomic variables were also not available for our study. Adjusting for these potential confounders can strengthen the findings for future studies.

Moreover, due to the lack of data, the current study was not able to construct generalized linear models with control variables such as autoregressive components, seasonality, and trend. Therefore, further studies are required with a sufficient length of data introducing the different waves of the pandemic to investigate the impacts of the combination of air pollution, climate variability, and other factors extrinsic to COVID-19. Since dust intrusion is associated with a decreased mixing layer height, we recommend in further studies to investigate the impact of other pollutants and also constitutes of particulate matter (PM 10 and PM 2.5 ), which can also be linked to the severity and incidence of COVID-19 (Frontera et al. 2020; Pandolfi et al. 2014; Salvador et al. 2019; Yao et al. 2020b; Zoran et al. 2020) . It is worth mentioning that reduced mixing layer height not only increases Karun-Abadan-Shushtar-Bavi-Shush-Ahvaz-Andimeshk 3

Rahmhormouz-Korramshahr-Dezful-Izeh-Bagh-e-Malek-Masjed Soleiman-Hoveizeh-Hamidiyeh the level of primary pollutants but also can make the dispersion of other pollutants more difficult . Health and environmental data on the ground in real time is almost always incomplete, if not entirely unavailable. Yet, to respond promptly to the demand of pandemics and to conduct short-term projections of disease incidence and death that can inform the planning for required beds, ventilators, and medicines, we require that data. The innovative utilization of satellite remote sensing for making projections following dust storms may provide that data and save lives. Our study is the first such application as a potential tool to inform the planning for the current pandemic and future epidemics in a timely manner.

Sparse earlier research has demonstrated a positive correlation between MED intrusion and health impacts in downwind areas. The subsequent global growth of the number of confirmed COVID-19 cases, the potential causal association of diffusion of COVID-19 with air pollution and dust intrusion was proposed. For the first time, this paper investigated the correlation between dust intrusion and daily increase in COVID-19 infection counts in Khuzestan, Iran, following the Middle East Dust (MED) intrusion of April 2020, when PM 10 pollution occurred for several days. In Khuzestan, only after a 7-day lag, the recorded number of new cases tripled, showing a significant abnormal increase in the province's daily disease rate. A direct correlation between the elevated daily combined AOD and an increase in the percentage of COVID-19 infection numbers was evident in the region, starting from 20/04/2020 to 09/05/2020. A Pearson correlation coefficient analysis shows that combined daily AOD in the studied cities of Khuzestan, Iran, was significantly positively associated with diurnal average wind speed, temperature, and the Planetary Boundary Layer. Also, RF analysis showed the critical role of wind speed and AOD with the relative correlation coefficient values of 0.31 and 0.23, respectively, influencing the daily increased percent of COVID-19 disease counts. This first study may serve as a reference to better understanding and predicting the factors affecting COVID-19 transmission and diffusion routes, focusing on the MED intrusions' role. Moreover, the findings can serve to support the adoption of public health measures to protect vulnerable population groups such as people with previous pathologies and the elderly by advising them to reduce their exposure to outdoor air pollution during dust intrusions.

The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s11356-021-18195-7. Table 5 The summarized conducted studies showing the association among particulate matter (PM 2.5 and PM 10 ) and COVID-19 infection and/or mortality in different regions around the world 

Evaluation and monitoring dust storm by remote sensing (case study: west and southwest of Iran)

Temporal variations in the frequency and concentration of dust events over Iran based on surface observations

Short-term effect of dust storms on the risk of mortality due to respiratory, cardiovascular and all-causes in Kuwait

Dust storm simulation over Iran using HYSPLIT

Impact of climate and ambient air pollution on the epidemic growth during COVID-19 outbreak in Japan

Temporal long-term variations in the occurrence of dust storm days in Iran

The effects of particulate matter sources on daily mortality: a case-crossover study of

Global ramifications of dust and sandstorm microbiota

Environmental pollution and COVID-19 outbreak: insights from Germany

First data analysis about possible COVID-19 virus airborne diffusion due to air particulate matter (PM): the case of Lombardy (Italy)

Turbulent gas clouds and respiratory pathogen emissions: potential implications for reducing transmission of COVID-19

A critical evaluation of the ability of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) thermal infrared red-green-blue rendering to identify dust events: theoretical analysis

Dust source identification using physical-chemical characterization and numerical modeling in Masjed Soleyman

The effect of dust storm on the microbial quality of ambient air in Ahvaz city and dust source identification using numerical modeling

Identification of dust storm source areas in West Asia using multiple environmental datasets

Multiple relationships between aerosol and COVID-19: a framework for global studies

Ambient influenza and avian influenza virus during dust storm days and background days

Is short-term exposure to ambient fine particles associated with measles incidence in China? A multi-city study

Air pollution and respiratory viral infection

Can atmospheric pollution be considered a co-factor in extremely high level of SARS-CoV-2 lethality in Northern Italy?

Dust phenomenon affects on cardiovascular and respiratory hospitalizations and mortality A case study in Kermanshah

Primary biological aerosol particles in the atmosphere: a review

Investigating the effects of meteorological parameters on COVID-19: case study of New Jersey, United States

Influence of airborne transmission of SARS-CoV-2 on COVID-19 pandemic. A review

Description of the HYSPLIT_4 modelling system

A methodology for the quantification of the net African dust load in air quality monitoring networks

Role of the chronic air pollution levels in the Covid-19 outbreak risk in Italy

Chemical and morphological properties of particulate matter (PM10, PM2.5) in school classrooms and outdoor air

Severe air pollution links to higher mortality in COVID-19 patients: the "double-hit" hypothesis

A comparative study of hospital admissions for respiratory diseases during normal and dusty days in Iran

Evaluation of the 13

Contribution of the Middle Eastern dust source areas to PM10 levels in urban receptors: case study of Tehran

Impacts of Saharan dust intrusions on bacterial communities of the low troposphere

Particulate matter and bacteria characteristics of the Middle East Dust (MED) storms over Ahvaz

Health risk assessment on human exposed to heavy metals in the ambient air PM10 in Ahvaz, southwest Iran

Air pollution analysis using enhanced K-Means clustering algorithm for real time sensor data

Dust storms in northern China: long-term spatiotemporal characteristics and climate controls

Clinical characteristics of 2019 novel coronavirus infection in China

A letter about the airborne transmission of SARS-CoV-2 based on the current evidence

Atmospheric investigation of frontal dust storms in southwest Asia

Atmospheric dynamics and numerical simulations of six frontal dust storms in the Middle East region

Spatial and temporal variability in dust storms in the Middle East

An assessment of SEVIRI imagery at various temporal resolutions and the effect on accurate dust emission mapping

Estimating ground-level PM2.5 concentrations in the southeastern US using geographically weighted regression

Transmission of COVID-19 virus by droplets and aerosols: a critical review on the unresolved dichotomy

The association between COVID-19 deaths and short-term ambient air pollution/meteorological condition exposure: a retrospective study from Wuhan, China

Characterization of atmospheric bioaerosols along the transport pathway of Asian dust during the Dust-Bioaerosol 2016 Campaign

Tmprss2 specific miRNAs as promising regulators for SARS-CoV-2 entry checkpoint

Diesel particulate matter 2.5 induces epithelial-to-mesenchymal transition and upregulation of SARS-CoV-2 receptor during human pluripotent stem cell-derived alveolar organoid development

Aerosol optical depth trend over the Middle East

Long-term exposure to air-pollution and COVID-19 mortality in England: a hierarchical spatial analysis

Effect of preexposure to ultrafine carbon black on respiratory syncytial virus infection in mice

The incubation period of coronavirus disease 2019 (COVID-19) from publicly reported confirmed cases: estimation and application

A novel calibration approach of MODIS AOD data to predict

Aerosol and boundary-layer interactions and impact on air quality

Acute lung injury in patients with COVID-19 infection

Instillation of particulate matter 2.5 induced acute lung injury and attenuated the injury recovery in ACE2 knockout mice

Impact of environmental factors and Sahara dust intrusions on incidence and severity of COVID-19 disease in Spain. Effect in the first and second pandemic waves

Identification of desert dust sources in the West Khuzestan Province using of brightness temperature algorithms

Kaskaoutis DG (2022) Mineralogical, geochemical, and textural characteristics of soil and airborne samples during dust storms in Khuzestan, southwest Iran

Meteorology and dust in the central Sahara: observations from Fennec supersite-1 during the

Use of SEVIRI images and derived products in a WMO Sand and dust Storm Warning System

Effect of environmental pollutants PM-2.5, carbon monoxide, and ozone on the incidence and mortality of SARS-COV-2 infection in ten wildfire affected counties in California

Health in dust belt cities and beyond-an essay by Nick Middleton

Human health impact assessment of exposure to particulate matter: an AirQ software modeling

Types of accelerated testing and modeling concepts

Spatiotemporal regionalization of atmospheric dust based on multivariate analysis of MACC model over Iran

Power of data mining methods to detect genetic associations and interactions

Airborne transmission of SARS-CoV-2: the world should face the reality

Assessing PM2.5 concentrations in Tehran, Iran, from space using MAIAC, deep blue, and dark target AOD and machine learning algorithms

Impacts of climate and synoptic fluctuations on dust storm activity over the Middle East

The impact of dusty days on fungi spores: hot vs. cold seasons of Ahvaz, Iran

Particulate matter concentrations during desert dust outbreaks and daily mortality in Nicosia, Cyprus

Evaluation of airborne particulate matter and metals data in personal, indoor and outdoor environments using ED-XRF and ICP-MS and co-located duplicate samples

The Effect of Dust Storm on the Microbial Quality of Ambient Air in Sanandaj: A City Located in the West of Iran

Air pollution by NO2 and PM2.5 explains COVID-19 infection severity by overexpression of angiotensinconverting enzyme 2 in respiratory cells: a review

Effect of atmospheric mixing layer depth variations on urban air quality and daily mortality during Saharan dust outbreaks

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

Evaluation of the potential relationship between Particulate Matter (PM) pollution and COVID-19 infection spread in Italy

Evaluating the effect of meteorological parameters (humidity,temperature,wind speed and pressure) on the dust phenomenon-case study

Exposure to fine particulate air pollution is associated with endothelial injury and systemic inflammation

Long-term PM2.5 exposures and respiratory, cancer and cardiovascular mortality in American older adults

Sampling and detection of corona viruses in air: a mini review

Dust-storm dynamics over Sistan region, Iran: seasonality, transport characteristics and affected areas

Dust storms in Iran-distribution, causes, frequencies and impacts

Deposition rates of viruses and bacteria above the atmospheric boundary layer

Transport of desert dust mixed with North African industrial pollutants in the subtropical Saharan Air Layer

Peaks of fine particulate matter may modulate the spreading and virulence of COVID-19

Transport routes and potential source regions of the Middle Eastern dust over Ahvaz during 2005-2017

Synergistic effect of the occurrence of African dust outbreaks on atmospheric pollutant levels in the Madrid metropolitan area

Saharan dust and associations between particulate matter and daily mortality in

Short-term effects of air pollution on daily single-and co-morbidity cardiorespiratory outpatient visits

The potential role of particulate matter in the spreading of COVID-19 in Northern Italy: First Evidence-based Research Hypotheses

Numerical analysis of the beam position monitor pickup for the Iranian Light Source Facility

Analysis of de-noising methods to improve the precision of the ILSF BPM electronic readout system

The evaluation of PM10, PM2.5, and PM1 concentrations during the Middle Eastern Dust (MED) events in Ahvaz

Short-term effects of particulate matter during desert and non-desert dust days on mortality in Iran

Does environmental quality and weather induce COVID-19: case study of Istanbul, Turkey. Environ Forensic 1-12

Monitoring weather based meteorological data: clustering approach for analysis

Climate change 2007: the physical science basis

Effects of short-and long-term exposures to particulate matter on inflammatory marker levels in the general population

Majority voting based multitask clustering of air quality monitoring network in Turkey

Emerging WuHan (COVID-19) coronavirus: glycan shield and structure prediction of spike glycoprotein and its interaction with human CD26

European Society for Clinical Microbiology and Infectious Diseases Fungal Infection Study Group; ESCMID Study Group for Infections in Critically Ill Patients (2020) Diagnosing COVID-19-associated pulmonary aspergillosis

Temporal cross-correlations between ambient air pollutants and seasonality of tuberculosis: a time-series analysis

Legal immigrants: invasion of alien microbial communities during winter occurring desert dust storms

Health relevance of particulate matter from various sources

Environmental engineers and scientists have important roles to play in stemming outbreaks and pandemics caused by enveloped viruses

Exposure to air pollution and COVID-19 mortality in the United States: a nationwide cross-sectional study

Acute effects of air pollution on the incidence of hand, foot, and mouth disease in Wuhan, China

A systematic review of lopinavir therapy for SARS coronavirus and MERS coronavirus-a possible reference for coronavirus disease-19 treatment option

Association of particulate matter pollution and case fatality rate of COVID-19 in 49 Chinese cities

ACE2 exhibits protective effects against LPSinduced acute lung injury in mice by inhibiting the LPS-TLR4 pathway

Haze is a risk factor contributing to the rapid spread of respiratory syncytial virus in children

Spatio-temporal occurrences and mineralogical-geochemical characteristics of airborne dusts in Khuzestan Province (southwestern Iran)

The characteristics of abnormal wintertime pollution events in the Jing-Jin-Ji region and its relationships with meteorological factors

Association between short-term exposure to air pollution and COVID-19 infection: evidence from China

Assessing the relationship between surface levels of PM2.5 and PM10 particulate matter impact on COVID-19 in

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

The authors are most grateful to the Department