key: cord-0920015-nrqqvppu
authors: Iqbal, Anam; Haq, Wajiha; Mahmood, Tahir; Raza, Syed Hassan
title: Effect of meteorological factors on the COVID-19 cases: a case study related to three major cities of the Kingdom of Saudi Arabia
date: 2021-11-12
journal: Environ Sci Pollut Res Int
DOI: 10.1007/s11356-021-17268-x
sha: 6496145a19bcfa107d10ab42a60cfe7ab30ced29
doc_id: 920015
cord_uid: nrqqvppu

The COVID-19 pandemic affected the world through its ability to cause widespread infection. The Middle East including the Kingdom of Saudi Arabia (KSA) has also been hit by the COVID-19 pandemic like the rest of the world. This study aims to examine the relationships between meteorological factors and COVID-19 case counts in three cities of the KSA. The distribution of the COVID-19 case counts was observed for all three cities followed by cross-correlation analysis which was carried out to estimate the lag effects of meteorological factors on COVID-19 case counts. Moreover, the Poisson model and negative binomial (NB) model with their zero-inflated versions (i.e., ZIP and ZINB) were fitted to estimate city-specific impacts of weather variables on confirmed case counts, and the best model is evaluated by comparative analysis for each city. We found significant associations between meteorological factors and COVID-19 case counts in three cities of KSA. We also perceived that the ZINB model was the best fitted for COVID-19 case counts. In this case study, temperature, humidity, and wind speed were the factors that affected COVID-19 case counts. The results can be used to make policies to overcome this pandemic situation in the future such as deploying more resources through testing and tracking in such areas where we observe significantly higher wind speed or higher humidity. Moreover, the selected models can be used for predicting the probability of COVID-19 incidence across various regions.

At the inception of COVID-19, health care providers in Wuhan, capital of Hubei, a province of China, reported numerous cases of pneumonia by a known virus to the World Health Organization (WHO) on December 31, 2019. In January 2020, a novel virus was identified, and the samples were taken from cases (WHO 2020b). In February 2020, the WHO named it as coronavirus disease (WHO 2020c ) and referred to as acute respiratory syndrome coronavirus-2 (SARS-CoV-2) (Public Health England 2020). The World Health Organization (WHO) has officially announced the COVID-19 as a pandemic on March 11, 2020, due to the high contagiousness and widespread (WHO 2020d) . The incubation period of the virus is argued to vary between 5 and 10 days Guan et al. 2020 ). The initial manifestation may include cough, fever, sore throat, myalgia, anorexia, and breathing difficulties (Guan et al. 2020; Jiang et al. 2020) . It is declared that the virus emanated in animals, and most of the cases came into view in the seafood market in Wuhan (HuangCL 2020; Rizwan et al. 2020; . The genetics analysis and some clinical discoveries can be found in studies Zhu et al. 2020) and references therein. It is not to deny that mRNA vaccine has been introduced in different parts of the world but yet the best knowledge about the behavior of the COVID-19 is to be investigated. Especially, it is yet to be known as our best knowledge that how meteorological factors can play an important role in the increase/decrease of infected cases.

Many studies on the COVID-19's influence have been concentrated on China (Gao et al. 2021; Lin et al. 2020) , USA (Bashir et al. 2020; Sharma et al. 2021) , and Europe (Menebo 2020; Yuan et al. 2021 ). There is a need to explore the virus's influence on the Middle East as well. This research focuses on the Kingdom of Saudi Arabia (KSA) that has been negatively affected. The Kingdom of Saudi Arabia (KSA) is the largest country in the Arabian Peninsula and thus had a higher risk in terms of infections and death from COVID-19 from the very beginning. After the identification of the first case on March 2, 2020, the number of cases increased rapidly, resulting in around 3,63,061 cases as of January 3, 2021. The Kingdom of Saudi Arabia (KSA) is located in the Middle East and has a 95% desert area. It has a warm and dry climate with high temperatures during the day, unexpected drops in temperature at night, and very low annual rainfall in most of the country. The number of confirmed cases of COVID-19 started to rise quickly in the last week of March. The Saudi government took protective measures to minimize the destructive effects of COVID-19, such as suspension of Umrah, travel ban, closure of mosques for prayers, and curfew. On June 23, 2020, the Ministry of Hajj and Umrah declared that Hajj 2020 would be held with a restricted number of pilgrims. The Ministry of Health took serious measures during the Hajj 2020; pilgrims were required to quarantine themselves at hotel rooms for a week before arrival to perform Hajj, and they were asked to be carried on a safe social distance of 2 m and were split into groups of 20 individuals per group, who were followed by a guide throughout their whole Hajj performance (WHO 2020a) . From March 2, 2020, to January 3, 2021, a total of around 3,63,061 COVID-19 cases have been reported in Saudi Arabia (Arab News 2020).

According to epidemiological studies, environmental circumstances may have influenced the spread of past outbreaks such as severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome disease (MERS). The SARS epidemic was gradually faded with the warming weather coming and was basically ended in July 2003 (Cao et al. 2016 Tan et al. 2005; Wang et al. 2006) . Some researchers suggested that climate change may have aided the emergence and spread of different viral diseases, including the SARS and MERS (Gale et al. 2010; Lofgren et al. 2007; Stott 2016) . The previous studies presented that meteorological factors play an important role in the eruption of infectious diseases like influenza (Li 2011; Lowen et al. 2007; Shaman and Kohn 2009) . Researchers in Korea discovered that low daily temperature and low/high relative humidity (RH) greatly enhanced the risk of influenza, but that diurnal temperature range was positively connected with influenza incidence (DTR) (Park et al. 2020) . As the coronavirus belongs to the same family of the influenza virus and has symptoms similar to colds, coughs, flu, and fever, it must be impacted by weather as well. Hence, it is essential to investigate the effects of many elements on COVID-19 instances, such as temperature, humidity, and wind speed. In this study, the authors aim to investigate the effect of weather-related variables on COVID-19 spread across the three most affected cities of KSA (i.e., Medina, Riyadh, and Jeddah) . The relationships between meteorological factors and COVID-19 case counts in three cities of the KSA are not known to the best of our knowledge.

The rest of the article is organized as follows. In the next section, we provide literature review. The "Material and methods" section contains the study design, data collection methods, and analytical techniques used in this study. The "Results and discussion" section presents the results obtained using the analytical methods implemented for the data analysis and some discussion. Finally, we provide conclusions in the "Conclusion" section.

COVID-19's most researched relationship is with meteorological conditions and air pollution, which influence coronavirus transmission and contribute to COVID-19's spread (Gao et al. 2021; Suhaimi et al. 2020) . Some studies proposed that the COVID-19 epidemic might be affected by temperature and its alterations (Cao et al. 2016; Tan et al. 2005; Wang et al. 2006) . Another research presents that the weather with low temperature, mild daytime temperature range, and low humidity favors the transmitting process of COVID-19 ). The regional climatic condition might influence the COVID-19 outbreak (Mazhar Iqbal et al. 2020) . It is suggested that weather is related to the spread of COVID-19, but the high temperature may not decrease the case counts without taking into account the execution of large-scale public health interventions ). The literature Wang et al. 2020) proposed that when examining a single-factor model in the higher temperature group, every 1 °C rise in minimum temperature is related to the decline of the cumulative number of COVID-19 by 0.86. In another study, it was noted that there was no association between COVID-19 outbreak and high temperature or ultraviolet (UV) radiation in Chinese cities . Oliveiros et al. (2020) reported that the temperature rise might not surely be related to the decline of COVID-19 case counts without taking into account containment measures.

Since the outbreak of COVID-19 in November 2019, several empirical studies have been conducted in China , the USA (Zangari et al. 2020) , Brazil (Auler et al. 2020) , and India (Jain and Sharma 2020) . Menebo (2020) analyzed the correlation between the COVID-19 pandemic and components of weather in the capital city of Norway. Fu et al. (2021) studied the association of COVID-19 transmission with temperature, absolute humidity, and diurnal temperature range in the four European countries. Islam et al. (2021) found a substantial positive correlation between the minimum and mean temperature, wind speed, relative humidity, and absolute humidity, but contact transmission had no direct relationship with the number of COVID-19 confirmed cases in Bangladesh. Sarkodie and Owusu (2020) developed conceptual tools to compare confirmed cases, deaths, and recovery cases using dew/frost point, temperature, disaggregate temperature, wind speed, relative humidity, precipitation, and surface pressure. In literature, most of the studies were designed to highlight the COVID-19 impact on KSA public and health practitioners (see studies Al-Ahdal et al. 2012; Alsubaie et al. 2019; El-Masry et al. 2021; Nasar et al. 2020; Obied et al. 2020 ; Qadah 2020 and references therein).

To investigate the relationships between COVID-19 and meteorological conditions, several statistical and scientific methods are used. For instance, Khursheed et al. (2021) applied generalized additive model (GAM) and examined a non-linear relationship between meteorological factors and COVID-19 transmission with penalized spline methods. Yuan et al. (2021) studied the impact of meteorological variables on daily COVID-19 cases employing log-linear generalized additive model. Pani et al. (2020) investigated the connection between COVID-19 and climatic indicators using Spearman and Kendall rank correlation tests. Zoran et al. (2020) conducted spatial analysis; quantile on quantile approach was used by Shahzad et al. (2020) . Iqbal et al. (2020) used the wavelet approach to discover that temperature has no effect on COVID-19 instances. The M-SEIR model is used by Shi et al. (2020) to explain why there is no significant relationship between humidity and COVID-19.

In literature, most of the studies concluded that the countries that are now suffering the most from the COVID-19 epidemic are those that are located in colder climates. In regions with a warm, dry climate, there are few research exploring the relationship between meteorological parameters and COVID-19 cases. The current work hypothesizes a null relationship between the COVID-19 (confirmed cases) and meteorological parameters (temperature, humidity, and wind speed) in the three most infected cities of KSA using count modeling. The quest to find the factors affecting the disease will help manipulate the significant factors to control the disease.

This section introduces the study design, data collection methods, and the analytical approaches used to analyze the effect of meteorological factors on COVID-19 case counts.

A conspicuous increase in COVID-19 case counts was observed in Medina, Jeddah, and Riyadh, due to which we chose these cities to study the relationship between weather and case counts. Medina is the second holiest city for Muslims in the western region of Saudi Arabia. It has a desert climate with long and humid summers and short and windy winters. It is mostly clear all over the year with an average temperature of 81.4°F. Riyadh is located on the eastern part of the Najd (central region) plateau with a hot desert climate with an average temperature of 77.7°F and low humidity throughout the year; its wettest month is April. Jeddah is situated in the western region of Saudi Arabia, the average high temperature is 84.2°F, and November is its most wet month. The COVID-19 confirmed case counts in three cities were extracted from King Abdullah Petroleum Studies and Research Center (KAPSARC) Data Portal [https:// datas ource. kapsa rc. org/] for the period of March 2 to July 15, 2020. The daily meteorological data containing average temperature, average relative humidity, and average wind speed (mph) were obtained from the website of weather underground [https:// www. wunde rgrou nd. com/] from March 2 to July 15, 2020.

A descriptive analysis was carried out to examine the characteristics of all variables in three cities, followed by a three-stage analysis. In the first stage, we fitted the Poisson, NB, ZIP, and ZINB distribution on the number of case counts for each city and estimated their parameters. Also, we performed diagnosis analysis based on the model selection criteria such as the Akaike information criterion (AIC) and Bayesian information criterion (BIC). In the second stage, cross-correlation analysis was conducted to identify the lag effects of meteorological factors that might be useful in predicting COVID-19 case counts. Then we adopted generalized linear models with Poisson and NB distribution. Moreover, we used zero-inflated versions of Poisson and NB models to estimate the city-specific case counts. The backward elimination method was applied to select the most significant meteorological factors in regression models. The structures of models used in this study are given below:

In count modeling, the most used technique is the Poisson regression model. It requires equidispersion that is the expected value of the response variable is equal to its variance. This technique is easy to use and permits zero counts and positive skewness (Cameron and Trivedi 2013) . Let (y i , x i ) be an observation and y i given x i has Poisson distribution. Thus, the model can be formulated as:

where i is the mean vector, X i represents covariates, X ′ i is the transpose of the vector X i , and is the vector of unknown parameters. According to the equidispersion property, the conditional mean and variance of the Poisson regression

i . The estimation method for the Poisson regression model is maximum likelihood estimation (MLE). So the log-likelihood function for the response variable is defined as:

Furthermore, the parameter vector can be estimated iteratively using a Newton-Raphson and a Fisher scoring technique (Agresti 2003) . The parameters are interpreted as: for one-unit change in an explanatory variable, the difference in the logs of the expected count is expected to change by

the corresponding coefficient, given the other explanatory variables in the model are constant.

The Poisson model assumption (i.e., mean equals the variance) is often violated due to the heterogeneity in the data. Commonly, the NB regression model is used to model the overdispersed count data. It can be taken as the generalization of Poisson regression because it has the same mean structure as the Poisson regression and has an extra parameter to model the overdispersion. Let (y i , x i ) be an observation and y i given x i has NB distribution. Thus, the model can be formulated as:

The conditional mean and variance of the NB model are

The index k represents the dispersion parameter. In NB regression model, the mean of the response variable is linked to a linear combination of explanatory variables by a log link function and given as:

where i is the mean vector, X i represents covariates, X ′ i is the transpose of the vector X i , and is the vector of unknown parameters. The parameters are estimated by the MLE method. So, the log-likelihood function for the response variable is given as:

Similar to the Poisson model, in the NB model, the parameter vector can also be estimated by using the Newton-Raphson and a Fisher scoring method. To fit the NB model, package "MASS" (Venables et al. 2002) is used in R. Given that the other explanatory factors in the model are constant, a one-unit change in an explanatory variable is anticipated to affect the difference in the logs of the expected count by the corresponding coefficient.

Sometimes, the cause of variation in the count data set might be an excess of zeros. In such a situation, the ordinary Poisson regression model is not suitable as it provides biased and inadequate estimates due to the violation of the equidispersion assumption. So, to get overcome the effect of excess zeros, an alternative to the ordinary Poisson regression model is proposed by Lambert (1992) . In the proposed ZIP model, all zeros are grouped into two parts such as (i) first part models a proportion (1 − )e − i of zeros coming from Poisson distribution and zero inflation by incorporating zeros with a proportion and (ii) the second part models the non-zero counts that follow the Poisson distribution. Let the response variable be independent, and then the ZIP model is defined as:

where Z i represents the vector of covariates describing the probability of excess zero i . The conditional mean and variance for the ZIP regres-

The ZIP model is formulated by adding the following two link functions:

In the ZIP model, i is modeled as:

where X ′ i is the vector of covariates, Z ′ i is the vector of covariates for defining the probability of excess zero i , and and are the vectors of unknown parameters. For estimation of parameters, the log-likelihood function for the response variable is defined as:

Let I i be an indicator variable which is defined as follows:

Then the log-likelihood function takes the form:

Furthermore, the expectation-maximization (EM) algorithm can be used to maximize the log-likelihood function iteratively (Lambert 1992) . To fit the ZIP model, package "pscl" (Jackman 2010) is employed in R. The interpretation of the parameter is same as in Poisson regression, but the parameter can be interpreted by the odds ratio (Lambert 1992 The zero-inflated negative binomial model

The data is mostly overdispersed due to heterogeneity in the data, and the ZINB distribution is used to model such data . The ZINB distribution is a continuous mixture of Poisson distribution, and it allows the Poisson mean to be gamma distributed to model the overdispersion. Similar to the ZIP model, in ZINB, all zeros are grouped into two parts such as (i) first part model a proportion (1 − ) 1 + i − of zeros coming from NB distribution and zero inflation by incorporating zeros with a proportion and (ii) the second part models the non-zero counts that follow ordinary NB distribution. The ZINB regression model is defined as follows:

where Z i is the vector of covariates describing the probability of excess zero i . The conditional mean and variance of t h e Z I N B m o d e l a r e

where is the dispersion parameter. The ZINB model is formulated by adding the following two link functions:

In the ZINB model, i is modeled as:

where ′ X i is the vector of covariates, ′ Z i is the vector of covariates for defining the probability of excess zero i , and and are the vectors of unknown parameters. The ZINB loglikelihood function is defined as:

The Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm is used to estimate the parameters of the ZINB model (Mahmood 2020) . To fit the ZINB model, package "pscl" (Jackman 2010) is employed in R. The interpretation of the parameter is the same as in Poisson regression, but the parameter can be interpreted by the odds ratio.

The autocorrelation is usually more robust in residuals from time series regression of infectious diseases (Cameron and Trivedi 2013); for that reason, in these models, the response variable (daily COVID-19 case counts) was introduced as an explanatory variable at 1-day lag. No obvious seasonal variation and long-term trend of case counts in three cities were perceived during this period (see Fig. 1 ).

The selection of a suitable model is an essential step for data modeling. The choice of a wrong model may lead to estimation errors and false conclusions. Therefore, in the

ln(Γ( )) + ln Γ 1 + y i − ln Γ + y i third stage of analysis, we have used a two-step procedure known as the LRT-Vuong method; the procedure is given in Mahmood (2020) 

The Vuong (Vuong 1989) test is used to choose between the Poisson model and the ZIP model or NB model and ZINB model. Let P Y i = y i |X i be the predicted probability of response variable from the standard model and P Y i = y i |X i , Z i be the predicted probability of response variable from the zero inflation model, and then the statistic m i be defined as:

Thus, Vuong's statistic for testing the nested hypothesis of standard model versus inflation model, i.e., H 0 ∶ E m i = 0 versus H 1 ∶ E m i ≠ 0 , is defined as: Under the null hypothesis, Vuong's statistic follows normal distribution. At 5% level of significance, when V > 1.96 , the standard model is chosen, when V < −1.96 , the zeroinflated model is chosen, and |V| < 1.96 , then both models are assumed equivalent.

For the assessment of overdispersion, we first applied both Poisson and NB models to the data of each city. Furthermore, the LRT test was utilized to examine whether the NB model gives the best fits or the Poisson model. The significant findings of the NB model give evidence of overdispersion in the data set. In the second step, we had used the Vuong test. In the case of overdispersion, the NB model was compared with the ZINB model; otherwise, the Poisson model was compared with the ZIP model. If the zero-inflated model yields a better fit as compared to the ordinary model, then it reveals zero inflation in the data. Hence, in the case of zero inflation, the zero-inflated model is effective for estimation and drawing conclusions. It is also to be noted that all analyses were conducted using R software with version 3.6.2.

The descriptive results for all variables are presented in Table 1 . A total of 95,115 cases were officially reported in three cities of KSA from March 2 to July 15, 2020. Among the three cities, the average daily case counts were high in Riyadh with a large variance. During our study period, the daily average temperatures for Medina, Riyadh, and Jeddah were 88.92°F, 86.98°F, and 85.44°F, respectively. Medina city had a smaller number of cases with the highest average temperature, and Jeddah city had the lowest average temperature but with large case counts as compared to Medina. The humidity and wind speed were relatively higher in Jeddah than in the other two cities.

The distribution of COVID-19 case counts in three cities can be seen in Fig. 2 . The distribution results for COVID-19 case counts in three cities are given in Table 2 . We found that AIC and BIC were smaller for ZIP as compared to Poisson distribution and smaller for ZINB as compared to NB distribution. Overall, ZINB distribution fitted better the case counts in all three cities.

In cross-correlation analysis (cf. Table 3 ), we found a significant positive association between daily cases and temperature on the same day and at 1-to 3-day lag for Medina. In Riyadh, daily cases were significantly and positively associated with the temperature on the same day and at 1-to 3-day lag. A significant negative association was noticed between daily case counts and humidity on the same day and at 1-to 3-day lag. Similarly, a significant negative association was also found between daily cases and wind speed on the same day and at 1-to 3-day lag, which is in consensus with the previous study by Rendana (2020) . We observed that, in Jeddah, there was a significant positive relationship between daily cases and temperature on the same day and at 1-to 3-day lag. So, our results confirmed that weather conditions influence COVID-19 case counts, and this finding is in line with previous studies on infectious diseases (Cao et al. 2016; Li 2011; Lowen et al. 2007; Tan et al. 2005) . All significant lagged variables were included in city-specific regression models.

The Poisson regression analysis results for Medina (cf. Table 4) revealed that the daily cases at 1-day lag had a significant positive impact on COVID-19 daily case counts. A 1°F increase in temperature at 2-day lag was associated with an increase in daily COVID-19 case counts. However, each 1°F increase in temperature at 3-day lag was related to the decline in the number of daily cases. From ZIP regression analysis, we found that daily cases at 1-day lag had a significant positive effect on daily case counts. Each 1°F increase in temperature at 2-day lag led to a rise in the number of cases. At the same time, each 1°F increase in temperature on the same day and at 3-day lag was related to declining in case counts. We also found that the probability of COVID-19 occurrence was negatively related to the temperature on the same day. From NB and ZINB regression analysis results (cf . Table 4 ), it was observed that one unit increase in daily cases at 1-day lag and 1°F increase in temperature at 1-and 2-day lag were likely to increase daily case counts. Also, the probability of COVID-19 occurrence was negatively related to the temperature on the same day. The studies (Liu et al. Wang et al. 2020; Yuan et al. 2021) reported that the spread of COVID-19 is related to temperature. In our study, as we can see from Table 4 , the temperature at 1-to 3-day lag had a greater impact on COVID-19 case counts than other variables in Medina. For example, Zoran et al. (2020) found that the maximum temperature might have a major effect on COVID-19 transmission. Behnood et al. (2020) indicated significant association between temperature and COVID-19 infection rates.

In Poisson and ZIP regression analysis for Riyadh (cf .  Table 5 ), we perceived that the number of daily cases at 1-day lag was linked with increased case counts. Each 1°F increase in temperature on the same day at 1-and 3-day lag was associated with the prevalence of COVID-19 cases, while each 1°F increase in temperature at 2-day lag was likely to decrease the cases. Each 1% rise in humidity on the same day, at 1-and 3-day lag, was related to the increase of daily cases. At the same time, a 1% increase in humidity at 2-day lag leads to a decline in cases. Furthermore, each 1 mph increase in wind speed on the same day at 1-and 3-day lag was connected with reduced case counts. But each 1 mph increase in wind speed at 2-day lag was likely to lead rise in the case of counts. It was also observed that the temperature at 1 and 3 days was related negatively to the probability of the occurrence of COVID-19. In NB and ZINB regression analysis (cf. Table 5), we found that number of cases at 1-day lag was associated with the rise in case counts. Each 1°F increase in temperature on the same day at 3-day lag was correlated to the rise of case counts. Also, each 1% increase in humidity and 1 mph increase in wind speed at 3-day lag were related to increased case counts. Also, the probability of the occurrence of COVID-19 case counts was negatively associated with the temperature at 1-and 3-day lag. Hence, we perceived that the temperature at 3-day lag and wind speed on the same day and 3-day lag had a greater impact on the COVID-19 epidemic than other variables in Riyadh. This finding is consistent with some previous studies (Auler et al. 2020; Guo et al. 2021; Menebo 2020; Şahin 2020) . For instance, Wei et al. (2020) reported that the temperature and wind speed significantly affected the confirmed cases of COVID-19 in the spring in Beijing. Gupta (2020) identified that wind speed promoted COVID-19 spread in Turkey.

The output of Poisson and ZIP regression analysis for Jeddah (cf. Table 6 ) showed that the rise in daily cases at 1-day lag was likely to lead to an increase in COVID-19 daily case counts. Each 1°F increase in temperature on the same day and temperature at 2-to 3-day lag was related to rising of COVID-19 daily case counts, while 1°F increase in temperature at 1-day lag was associated with the decline of cases. It seemed no explanatory variable had a significant impact on the probability of COVID-19 occurrence. The NB and ZINB regression analysis (cf. Table 6) indicated that daily cases at 1-day lag were likely to increase COVID-19 daily case counts. Moreover, we found that each 1°F increase in temperature on the same day and at 2-to 3-day lag was correlated with increased case counts. Also, it was observed that no explanatory variable had a significant effect on the probability of occurrence of COVID-19 case counts. These results showed that the temperature on the same day and at 1-day lag had a greater impact on COVID-19 case counts in Jeddah. Similar results were found in previous research, which showed the contribution of temperature to COVID-19 transmission (Chien and Chen 2020; Sarkodie and Owusu 2020; Yang et al. 2021) . For example, Bhardwaj and Agrawal (2020) presented that a modest rise in temperature will increase spatial distances and viral exposure risk due to a rise in the intensity of Brownian motion of air particles. Sharma et al. (2021) showed that the risk of coronavirus infection in people was affected by air pressure and temperature.

The LRT statistic results (c.f. Tables 4, 5, and 6) in three cities indicated that the P values were less than 0.01, which led to the rejection of the null hypothesis. So, the significant findings of the NB model gave evidence of overdispersion in the data set of each city. After finding the significant results of the overdispersion test, we used the Vuong test to test for zero inflation (NB vs ZINB). The findings (c.f. Tables 4, 5, and 6) revealed that the ZINB model was more appropriate as compared to the NB model at 1% level of significance for all three cities. Also, the other model selection criteria, such as AIC and BIC, were in favor of the ZINB model. This finding is consistent with early research findings; for example, Benita and Gasca-Sanchez (2021) 

As COVID-19 is a continuous and increasing phenomenon; thus, it was vital to see the effect of meteorological factors on the disease. We examined the impact of meteorological factors on COVID-19 cases after considering the excess of zero counts. Our data revealed a significant correlation between COVID-19 cases and climatic parameters for all three most contaminated cities under study. For temperature as a meteorological factor, we found a significant positive association between COVID-19 case counts and temperature on the same day and past days. We also found that humidity and wind speed were associated negatively at various lags with COVID-19 case counts. The epidemic of COVID-19 in Medina, Riyadh, and Jeddah indicated overdispersed daily cases and zero cases on some days. When the response variable is in a discrete nature, then the Poisson model is preferred. The ZIP and ZINB models are suitable for modeling a data set with excess zero counts (Cameron and Trivedi 2013; Lambert 1992) , and in the case of the overdispersion, NB model is considered appropriate (McCullagh and Nelder 1989; Nelder and Wedderburn 1972) . Hence, we built four models to fit COVID-19 daily case counts for each city. Besides, we compared these models and perceived that the ZINB model is the best-fitted model for COVID-19 daily cases in Medina, Riyadh, and Jeddah.

The Kingdom of Saudi Arabia has a warm, dry climate with high temperatures and low annual rainfall. We found that, in three cities of KSA, an increase in temperature was not necessarily associated with the decline of cases in contradiction to the hypothesis. Although our findings show an independent impact of meteorological factors on COVID-19 case counts, the public health protective measures play a vital role in controlling the epidemic of This research adds to the domain of pandemic management in both a practical and theoretical sense. The findings provide a better knowledge of the situation and may aid in adopting suitable measures to implement pandemic response policies and activities as needed in a fast and effective way. For instance, our results ascertained that temperature, humidity, and wind speed contribute significantly towards the incidence of COVID-19; therefore, it can be argued that government needs to impose strict regulations during high temperature, humidity, and low wind speed. The most important implication is that, in pandemics, where widespread and effective vaccination is not possible, government preparation should be longer term, taking into account weather effects, and establish public health and safety directions accordingly.

For instance, with the proper weather-based changes, mask and distancing guidelines should be followed. Also, a robust response to the pandemic can be created by combining meteorological and medical expertise. In conclusion, our findings shed light on the association between meteorological factors and COVID-19 case counts in three cities of KSA. Overall, the areas with high temperature, low wind speed, and high humidity favor the rise in COVID-19 case counts. This study could be useful for forecasting agencies in developing predictive tools for COVID-19 prevention in KSA. The selected models could be used to predict the probability of COVID-19 occurrence and the number of cases in three cities of KSA. Furthermore, experts may reconsider and alter policy approaches in order to lessen the total impact on the environment and economy while still adhering to sustainable development.

There are certain limitations to work. Other variables regulating disease transmissions, such as air pollution, water pollution, ultraviolet radiation, rainfall, atmospheric pressure, and geographical characteristics, were not examined. These influencing variables should be taken into account in models like the ones utilized in this study. For future studies, the present work can be extended to other parts of the world with extensive data set and varying weather situations. The optimal temperature to counter-effect the spread of disease is yet to be investigated to the best of our knowledge regarding COVID-19 vaccines.

Coronavirus respiratory illness in Saudi Arabia

Middle East respiratory syndrome coronavirus epidemic impact on healthcare workers' risk perceptions, work and personal lives

Evidence that high temperatures and intermediate relative humidity might favor the spread of COVID-19 in tropical climate: a case study for the most affected Brazilian cities

Correlation between climate indicators and COVID-19 pandemic

Determinants of the infection rate of the COVID-19 in the US using ANFIS and virus optimization algorithm (VOA)

The main factors influencing COVID-19 spread and deaths in Mexico: a comparison between phases I and II

Likelihood of survival of coronavirus in a respiratory droplet deposited on a solid surface

Analysis of spatiotemporal characteristics of pandemic SARS spread in mainland China

A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating personto-person transmission: a study of a family cluster

Meteorological impacts on the incidence of COVID-19 in the US

Novel coronavirus disease-related knowledge, attitudes, and practices among the residents of Al-Jouf region in Saudi Arabia

The association of voter turnout with county-level coronavirus disease 2019 occurrence early in the pandemic

Meteorological factors, governmental responses and COVID-19: evidence from four European countries

Assessing the impact of climate change on vector-borne viruses in the EU through the elicitation of expert opinion

Impact of the COVID-19 pandemic on air pollution in Chinese megacities from the perspective of traffic volume and meteorological factors

Clinical characteristics of coronavirus disease 2019 in China

Meteorological factors and COVID-19 incidence in 190 countries: an observational study

Effect of ambient temperature on COVID-19 infection rate

Missed takes towards a pandemic of COVID-19? A systematic literature review of Coronavirus related diseases in Pakistan

The pandemic paradox: domestic violence and happiness of women

Clinicalfea• tures of patients infected with 2019 novel coronavirus in Wuhan

The nexus between COVID-19, temperature and exchange rate in Wuhan city: new findings from partial and multiple wavelet coherence

Effect of meteorological factors on COVID-19 cases in Bangladesh

Classes and methods for R. Developed in the Political Science Computational Laboratory

Social and travel lockdown impact considering coronavirus disease (COVID-19) on air quality in megacities of India: present benefits, future challenges and way forward

Investigating the roles of meteorological factors in COVID-19 transmission in Northern Italy

Zero-inflated Poisson regression, with an application to defects in manufacturing

Environmentally mediated transmission models for influenza and the relationships with meteorological indices

A mechanism-based parameterisation scheme to investigate the association between transmission rate of COVID-19 and meteorological factors on plains in China

Impact of meteorological factors on the COVID-19 transmission: a multi-city study in China

Influenza seasonality: underlying causes and modeling theories

Influenza virus transmission is dependent on relative humidity and temperature

The role of absolute humidity on transmission rates of the COVID-19 outbreak

Generalized linear model based monitoring methods for high-yield processes

Jawed Iqbal M (2020) The effects of regional climatic condition on the spread of COVID-19 at global scale

Generalized linear models 2nd edition chapman and hall

Temperature and precipitation associate with Covid-19 new daily cases: a correlation study between weather and Covid-19 pandemic in Oslo

A review of the prevalence of COVID-19 in the Arab world

Generalized Linear Models

Arab News (2020) Saudi Arabia records major decline in virus tally since March

Containment of COVID-19: the unprecedented response of Saudi Arabia

Role of temperature and humidity in the modulation of the doubling time of COVID-19 cases

Association of COVID-19 pandemic with meteorological parameters over Singapore

Effects of temperature, humidity, and diurnal temperature range on influenza incidence in a temperate region

Public Health England Qadah T (2020) Knowledge and attitude among healthcare workers towards COVID-19: a cross sectional study from Jeddah city, Saudi Arabia

Impact of the wind conditions on COVID-19 pandemic: a new insight for direction of the spread of the virus

Current perspective on diagnosis, epidemiological assessment, prevention strategies, and potential therapeutic interventions for severe acute respiratory infections caused by 2019 novel coronavirus (SARS-CoV-2)

Impact of weather on COVID-19 pandemic in Turkey

Impact of meteorological factors on COVID-19 pandemic: evidence from top 20 countries with confirmed cases

Effects of climatological parameters on the outbreak spread of COVID-19 in highly affected regions of Spain

Absolute humidity modulates influenza survival, transmission, and seasonality

Meteorological factors, COVID-19 cases, and deaths in top 10 most affected countries: an econometric investigation

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

How climate change affects extreme weather events

Demystifying a possible relationship between COVID-19, air quality and meteorological factors: evidence from Kuala Lumpur

An initial investigation of the association between the SARS outbreak and weather: with the view of the environmental temperature and its variation

Modern applied statistics with S, 4th edn

Likelihood ratio tests for model selection and nonnested hypotheses

Spatial dynamics of an epidemic of severe acute respiratory syndrome in an urban area

Temperature significant change COVID-19 Transmission in 429 cities

Impacts of transportation and meteorological factors on the transmission of COVID-19

Coronavirus: WHO chief hails Saudi Arabia's COVID-19 measures during Hajj. Al Arabiya English WHO (2020b) Novel Coronavirus -China. World Health Organization WHO (2020c) WHO Director-General's remarks at the media briefing on 2019-nCoV. World Health Organization WHO (2020d) Coronavirus disease (COVID-19) outbreak situation

Influence of meteorological factors on the COVID-19 transmission with season and geographic location

No association of COVID-19 transmission with temperature or UV radiation in Chinese cities

Non-linear correlation between daily new cases of COVID-19 and meteorological factors in 127 countries

Air quality changes in New York City during the COVID-19 pandemic

A novel coronavirus from patients with pneumonia in 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