Prevalence and Associated Risk Factor of COVID-19 and Impacts of Meteorological and Social Variables on Its Propagation in Punjab, Pakistan

The current study identifies the spatial distribution of COVID-19 cases and its association with meteorological and social variables in Punjab (densely populated province of Pakistan). To identify the COVID-19 propagation, the weekly growth, recovery, and deaths rate have also been calculated. The geographic information system (GIS) has used to determine COVID-19 impacts on gender (male/female), age groups, and causalities over an affected population (km−2) for the period of 11th March to 12th August, 2020 in each district of province. Our results show that 43 peak days (where daily positive cases were above 900) have been observed in Punjab during 27th May to 8th July, 2020. The high population density districts, i.e., Lahore and Islamabad, have been affected (five persons per square kilometers) due to COVID-19, whereas the maximum death tolls (> 50 persons per millions) have also been observed in these urban districts. The meteorological variables (temperature, humidity, heat index, and ultraviolet index) show negative significant relationship to basic reproduction number (R0), whereas daily COVID-19 cases are positively correlated to aerosols concentration at 95% confidence level. The government intervention (stringency index) shows a positive impact to reduce the COVID-19 cases over the province. Keeping in view the COVID-19 behavior and climatology of the region, it has been identified that the COVID-19 cases may likely to increase during the dry period (high concentration of aerosols) i.e., October–December, 2020 and post-spring season (April to June), 2021 in urban areas of Pakistan. This study provides an overview on districts vulnerability that would help the policy makers, health agencies to plan their activities to reduce the COVID-19 impacts.


Introduction
The outbreak of coronavirus disease  has severely affected many countries of the world. The first case of the virus was observed in early December in Wuhan, China. The World Health Organization (WHO) declared the disease a global pandemic on 11 th March, 2020, during which 66 territories outside China were reported 8568 confirmed cases of the virus, with highest numbers observed in Iran and Italy (Hsu et al. 2020). As of 11th March, 2021, the confirmed cases of COVID-19 around the globe were 118,638,190 including 2,632,074 deaths. So far, the United States has the highest number of COVID-19 cases with a total of 28, 862,124 and 542,191 deaths,while India (11,285,561) and Brazil (11,205,972) have the second and third highest COVID-19 cases, respectively (Worldometer 2020).
The spread of the virus was highest in the neighboring countries of Pakistan, i.e., China in the north and Iran in the west. Initially, the COVID-19 cases were the second highest (23,042) after China till 24th March, 2020 (Tuite et al. 2020;WHO 2020). In Pakistan, foreign transmission (pilgrims and tourists) was the highest (80%), especially 8000 pilgrims from Iran out of which 60% were infected with COVID-19 which afterwards reduced to 10% till 31st March, 2020 (Badshah et al. 2020). Due to the limited testing facilities, prompt isolation, etc. caused internal transmission, especially in the major cities of Pakistan as the majority of the foreign people returned to their homes (Wilder Smith and Freedman 2020). The first confirmed case of COVID-19 in Pakistan was reported in Karachi (Sindh) and Islamabad on 26th February, 2020 (Noreen et al. 2020a, b). According to the National Institute of Health (NIH), Pakistan's latest COVID-19 update (available on https:// covid. gov. pk/ stats/ pakis tan) on 10th March, 2021, the data from all provinces and territories reported a total of 583,916 confirmed cases, 13,013 deaths, 554,255 recovered, 9,246,827 total tests, and 1628 critical cases. Province-wise statistics of Pakistan have shown COVID-19 cases of 258,679 in Sindh,173,395 in Punjab,72,801 in Khyber Pakhtunkhwa,44,690 in Islamabad,19,076 in Balochistan, 10,319 in Azad Jammu Kashmir, and 4956 in Gilgit Baltistan. Among the provinces, the highest number of cases has been reported from Sindh followed by Punjab and Khyber Pakhtunkhwa, Balochistan, and Azad Jammu Kashmir, whereas Gilgit Baltistan was the least affected province (NIH 2021). Currently, Pakistan is the 14th mostaffected country with a fertility rate of 2.23 and a point prevalence of 8.3, while the peak was observed in the mid of June. Since then, a considerable decline in the number of new cases was noted due to government strategies for effective surveillance and climatic factors (NAPPR 2020).
Punjab is a densely populated province consisting of 36 administrative districts and contributes 52.95% (110,012,442) of the total population of Pakistan (Census report of Pakistan, 2017). Among 24,104 confirmed cases, 439 deaths with 126 critical cases and 6507 recoveries were reported, whereas the fertility and recovery rate remained 1.8% and 27%, respectively. The community transmission is reported to be 91% (Noreen et al. 2020a, b). Punjab's high prevalence and fertility rate can be attributed to a mixture of asymptomatic transmission and its initial untraced spread. The total stands at 2,192 cases with 23 deaths, representing an alarming hike in the number of daily cases since 18th May, 2020 (100 plus cases daily) with the highest cases (151) reported on 27th May, 2020. According to the census report (2017), the male population contributes almost 50.07-53.09% of the total population in all the districts of Punjab except Attock, Jhelum, Chakwal, and Khushab, where the female contributes 50.04 to 52.18% of the total population of these districts.
The world COVID-19 data analysis shows that men have been infected more by COVID-19 as compared to women (Nikpouraghdam et al. 2020). Moreover, a weaker immune system is highly prone to be affected by COVID-19, which makes the older demographic much more susceptible and more likely to be hospitalized due to COVID-19 related complications (CDC 2020). However, several studies from Pakistan depicted that the highest number of individuals which are affected by COVID-19 belong to the 20-29 and 30-39 age groups, while the lowest are in the above 80 (Abid et al. 2020;Noreen et al. 2020a, b;PIDE 2020;Ladiwala et al. 2021 etc.). These data seem paradoxical given that the elderly are more susceptible to the virus owing to a weakened immune system and poor health (Clark et al. 2020), but this disparity can be explained by looking closely at Pakistan's demographic and social structure. According to the United Nation's data, only 4% of Pakistan's population is above 65, while the average age is 22.5. This is in stark contrast to the highly affected countries such as Italy, where 23% of the population is above 65 and the average age is 47. Italy has the second highest oldest population in the world after Japan (OECD, 2021). Initially, the slow response, testing capacity, family culture, and old age population are the main reasons for high mortality rate in Italy and other European countries. An elderly and unhealthy individual is likely to suffer more fatally from the viral infection due to their weak immune system (CDCP 2020).
Numerous studies (Chan et al. 2011;Paz and Semenza 2016;Lee and Chowell 2017;Morin et al. 2018;Keilman 2019, etc.) have identified that the weather plays significant role to lessen the spread of infectious disease like influenza, pneumonia, and flue, etc. during winter and dengue fever, chicken pox, and typhoid in summer season. Theses all diseases are highly associated with the diurnal temperature Published in partnership with CECCR at King Abdulaziz University range. Recently, many researchers around the world conducted studies and identified the meteorological parameters (temperature, humidity) impact on COVID-19 spread (Qi et al. 2020;Livadiotis 2020;Shi et al. 2020;Iqbal et al. 2020). Similarly, Bashir et al. 2020 identified that the temperature and air quality suppress COVID-19 significantly. The UV light helps to minimize the spread of microbial disease and remain effective for viral disease like flu and influenza (Welch et al. 2018;Budowsky et al.1981;McDevitt et al. 2012;Ko et al.2000). A significant nonlinear association has been observed among daily COVID-19 cases, meteorological variables, and air quality index on different lag days over Dhaka (Islam et al. 2021). Similarly, the dust particles inhaled by humans can damage the respiratory system and immunity causing high vulnerability to COVID-19 infection (Chen et al. 2021).
Geographic information system (GIS) is very helpful in mapping of spatial monitoring, disease spread, disease cases, risk assessment etc. (Mollalo et al. 2018). GIS algorithms facilitate early detection and outbreak of disease over a spatial area (Meng, 2017). In early times, the health agencies were highly dependent on old traditional methods to control and monitor the infectious disease, which take long time to detect and identify the spread of disease. Many researchers have used GIS to identify the spatiotemporal disease mapping, risk assessment and surveillance, etc. (Franch-Pardo et al. 2020;Sarwar et al. 2020;Sugg et al., 2021 etc.). In this study, GIS helped to analyze the COVID-19 cases in genders, age groups as well as infected patients over a population density (km2) and mortality in each of the districts. The current study identifies the propagation, growth, recovery, and death rate of COVID-19 cases in the Punjab province. The study determines impacts of COVID-19 on gender, age group, and fatality over positive patient density (km −2 ) in each district of Punjab. Moreover, the meteorological and social variables are used to identify their impacts on COVID-19 propagation over the Punjab province. The relationship and impact of these variables on COVID-19 may help the scientific community, epidemiologist, policy makers, and health agencies to understand the behavior, response, transmission, spread, as well as propagation. The aptitude of detection of disease pattern in different districts of Punjab would help in the efficiently planning strategies for disease control and prevention as well as provide information for further epidemiological studies to analyze disease transmission.

Data and Methodology
The daily data of COVID-19 (tests, patients) for each district of Punjab have been collected from the National Institute of Health (NIH), Ministry of Health, Pakistan for the period of 11 th March to 12 th August, 2020. In addition, the daily observational data of temperature (°C) and relative humidity (%) are obtained from Pakistan Meteorological Department (PMD). We calculated heat index (HI) which is the combination of temperature and relative humidity over the province as determined by Steadman, 1979. Moreover, the daily ultraviolet Index (UVI) and Aerosols optical Depth (AOD) data, with a spatial resolution of 1° × 1°, have obtained from the Ozone Monitoring Instrument (OMI), onboard NASA's Aura spacecraft (Torres 2006), and available at the following web-link: https:// giova nni. gsfc. nasa. gov/ giova nni/. Some researchers have used MODIS satellite AOD data for sand and dust storm monitoring and try to develop its relationship to various meteorological parameters over Saudi Arabia (Butt et al. 2017;Butt and Mashat 2018;Labban and Butt 2021).
To identify the COVID-19 affected patients over gender and various age groups, we first converted the affected population of male and female over per million populations, i.e., the number of affected male/female over the total number of male/female then multiplied by 106 to convert into million. Similarly, the number of deaths per million as well as affected population density has also been calculated. Second, we classified the population into four age groups, i.e., teenage (0-19), young (20-40), aged (41-60), and senior citizen (> 60), and then calculated the percentage of affected patient for each group in entire districts of Punjab. GIS has been used to display the data over each of the districts.
The data of the stringency index (SI) are acquired from https:// ourwo rldin data. org/ graph er/ covid-strin gency-index. This index is created by the University of Oxford which determines the effectiveness of lockdown, and it ranges from 0 to 100. According to Hale et al. 2020, SI is based on various indicators like schools closure, workplace closures, cancel public events, restrictions on gatherings, close public transport, public information campaigns, stay at home, restrictions on internal movement, international travel controls, testing policy, contract tracing, face coverings and vaccination policy, etc. The higher values of SI indicate a high level of stringency over the region. The both indices, i.e., HI and UVI, are used to identify the impacts and their association with basic reproduction number (R0), whereas stringency index and aerosols concentration are correlated with daily positive COVID-19 cases over Punjab province. The geographical location and elevation of Punjab districts are shown in Fig. 1.
The positivity, growth, recovery, and mortality rate have been calculated. Positivity rate (P r ) is the ratio between total positive patients (P pt ) over the number of tests (T t ) (Eq. 1): Published in partnership with CECCR at King Abdulaziz University The growth rate (G r ) is the ratio between numbers of new cases 1 day ( ΔN d ) to the new cases the previous day ( ΔN d−1 ) or it is defined as the ratio between two successive changes (Ridenhour et al. 2018) (Eq. 2): Recovery rate (R r ) is the ratio between numbers of recovered patient ( R pt ) to the total positive patients (P pt ) (Eq. 3): Mortality rate (M r ) is the ratio between numbers of deaths ( D t ) to the total positive patients (P pt ) (Eq. 4): Basic reproductive number (R 0 ) is not a biological constant for a pathogen; however, it determines the expected number of cases directly generated by one case in a population where all individuals are susceptible to infection.
R 0 is the ratio between the fractions of individual susceptible per day (β) to the fraction of recoveries (γ) defined by Diekmann et al. 1990 (Eq. 5): The value of R 0 > 1 represents the spread of epidemic, whereas the R 0 < 1 shows the reduction in the spread of disease.
The correlation was calculated using Pearson methods, whereas t test has been used to determine the significance at 95% confidence level.

Results and Discussion
To determine the variation in COVID-19, the daily and weekly analysis was conducted over the Punjab province. The first COVID-19 case was reported on 15 th March 2020 in Punjab, while the tests started 4 days earlier (11th March, 2020). Initially, most of the cases were foreign transmission. The daily tests were below 3000 up to 30th April, 2020.
(5) R 0 = . The government provided free testing facility by setting up a dedicated COVID-19 cell in the hospital. Furthermore, the track and test approach was also used to identify the COVID-19 patients. Meanwhile, the daily number of new cases started to increase due to local transmission and the government increased the daily tests facilities in the province after April. The highest tests (11,739) were conducted on 17th June, whereas the highest positive patients (2,705) in day were reported on 12th June 2020 (Fig. 2a). Overall, about 785,764 COVID-19 tests were conducted and 94,865 patients reported positive, 86,389 recovered, and 2179 died during 11th March to 12th August, 2020 in the province. Forty-three peak days were experienced during 27th May to 8th July, 2020 over the whole province where the daily positive cases were above 900 (Fig. 2).
The total recovered patients from COVID-19 were 86,389 up to 12th August, 2020, due to which the number of active cases reduced to 8476. The maximum patient recovered during 27 th June to 13 th July, 2020, whereas the highest patient recovered in a day were 12,802 on 24th July, 2020 (Fig. 2b). The total deaths were 2133 during the whole data period. The daily death peak period was observed from 27th May to 8th July, 2020, with the highest deaths (86) in a day on 24th June 2020 (Fig. 2c).
The weekly positivity, growth, recovery, and mortality rate of COVID-19 have been calculated for the Punjab province since 11th March, 2020 (Fig. 3). The positivity rate (P r ) defines as the ratio between patients to test, whereas the growth rate (G r ) is the difference between two successive changes of COVID-19 cases. The peak time period of Pr remained between 10 and 18 th weeks (13th May and 14th July), whereas it was the highest (24.6) in 13th week (3-9th June), as shown in Fig. 3a. Similarly, G r was the highest (10.38) on the 2nd week (18-24th March) and it consistently dropped. Later on, G r was below one after the 14th week (17th June, 2020) except the 17th week, where it remained 1.07 (Fig. 3b). The recovery rate (R r ) was initially low (below 1) but increased after the 15th week on 17th June. The peak period of R r was between 16 and 20th week, while it was the highest (9.83) in the 20th week (22-28th July), as shown in Fig. 3c. The weekly mortality rate (M r ) remained 0.37-3.66 in the province. Figure 3d shows that the peak period remained between the 13 th and 20 th weeks with the highest Mr (3.66) in the 16 th week (24 th June).
The weekly positivity, growth, recovery, and mortality rate of COVID-19 have been calculated for the Punjab province since 11 th March, 2020 (Fig. 3). The positivity rate (P r ) defines as the ratio between patients to test, whereas the growth rate (G r ) is the difference between two successive changes of COVID-19 cases. The peak time period of P r remained between 10 and 18 th weeks (13th May and 14th July) whereas, it was the highest (24.6) in 13th week (3-9 th June), as shown in Fig. 3a. Similarly, G r was the highest (10.38) on the 2 nd week (18-24th March) and it consistently dropped. Later on, G r was below one after the 14th week (17th June, 2020) except 17th week, where it remained 1.07 (Fig. 3b). The recovery rate (R r ) was initially low (below 1), but increased after the 15th week on 17th June. The peak period of R r was between 16 and 20th week, while it was the highest (9.83) in the 20th week (22-28th July), as shown in Fig. 3c. The weekly mortality rate (M r ) remained 0.37-3.66 in the province. Figure 3d shows that the peak period remained between 13 and 20th weeks with the highest M r (3.66) in the 16th week (24th June).

Geo-spatial Analysis
Pakistan is experiencing a high number of COVID-19 cases, with the maximum share from Punjab, Pakistan (Saeed et al. 2021). The numbers of affected population (male and female) per million due to the COVID-19 have been determined in 36 districts of the province including Islamabad (Fig. 4). The results show that the highest male population (> 500) per million was affected in Lahore, Rawalpindi, and Multan due to COVID-19, whereas it was between 300 and 500 per million in Sialkot, Gujranwala, Hafizabad, Faisalabad, and DG Khan region. Similarly, the highest female affected population (> 500) per million was observed in Lahore, whereas it was between 300 and 500 per million in Rawalpindi and 200 and 300 in Multan and Faisalabad region.
It has been observed that the majority of the male population was affected due to COVID-19 in the province as the male family members have to go outside for work and their interaction with the community is more than female. However, more females were affected due to COVID-19 than males in Attock and Bakkar districts. The majority of the population is female in Attock, whereas most of the females work in agriculture fields and associated with agriculture activity in Bhakkar and have become more exposed to COVID-19. The districts of Lahore, Multan, and Rawalpindi are the business hub and the majority of the males traveled to these cities for the work. More  importantly, these districts have high population density and the medical and hospital facilities are much better.
The percentage of COVID-19-affected patients over various age groups (teenage, young, aged, and senior citizen) has been calculated in the entire districts of Punjab (Fig. 5). It has observed that 10-25% of teenagers were affected in the majority districts of Punjab. The highest (> 90%) of affected patients in the teenage group were observed in Bhakkar districts followed by Jhang and Khushab, where it was 50%. Similarly, almost 50% of the young people were affected in 22 districts and 75% in 14 districts, whereas no young age group patient was reported in Bhakkar district. The results

Fig. 5
Percentage of affected population: a Teenage, b Young, c Aged, and d Senior Citizen due to COVID-19 in various districts of Punjab province also show that 25% of the aged group has been affected in 13 districts and 50% in 23 districts including Islamabad. Moreover, 10% of senior citizens were affected in 23 and 25% in 13 districts of Punjab. The results also show that 25% of the aged group has been affected in 13 districts and 50% in 23 districts including Islamabad. Moreover, 10% of senior citizens were affected in 23 and 25% in 13 districts of Punjab. The COVID-19 has affected most of the females than males in Bhakkar (Fig-4). Normally, the female took their children (teenagers) along with them in the fields, so that they can look after them easily and their children may help them in daily routine activity carried out in the agriculture field.
May and June are the harvesting and sowing season of Rabi (Wheat) and Kharif (Cotton, rice and Maize) crops, respectively. Most of the people living in rural districts work in fields and ignored the standard operating procedure like using of masks and social distancing etc. due to unawareness and non-seriousness of the problem. In order to overcome this problem, the government should start awareness programs and strict compliance with its order should be implemented through local administration. The results showed that COVID-19-positive cases were higher during the months of May and June. More importantly, the provincial administration eased the lock down in the province and people living in and around the urban cities started to work in the city. Most of the people moved towards the major cities in the province to earn their livelihood. Due to which, the COVID-19 spike was experienced in these districts.
The maximum number of deaths per million has occurred in the districts of Lahore, Rawalpindi, and Islamabad where it was greater than 50 persons per million, whereas 50 persons per million died in the districts of Sialkot, Gujranwala, Sheikhupura, and Multan (Fig. 6). The death tolls are high in the urban districts with high population density. According to the world air quality index (2021), the air quality of major cities of Punjab, i.e., Lahore, Gujranwala, Faisalabad, Multan, Rawalpindi, and Sialkot, are very poor. However, Lahore ranks at the top among the rest of the cities where the air quality index remains worst throughout the year, whereas the pollen allergy is common during the spring season in twin cities, i.e., Rawalpindi and Islamabad. The respiratory disease is more common in the inhabitants of these areas, and they have become highly susceptible to COVID-19 impact which directly affects the lungs creating suffocation and causing death. The recent study shows that death rate Fig. 6 Number of deaths per million and populations affected (km −2 ) due to COVID-19 in districts of Punjab remains high in poor air quality region as people in that region are highly exposed to bad environment and become victims of chronic diseases (Wu et al. 2020).
Population affected per square kilometers (ratio between positive patients to the population density) of each district was calculated which showed COVID-19 transmission rate was high in densely populated regions (Fig. 6). However, the majority of the foreign transmissions are observed in DG Khan and Bahawalpur districts, where lots of pilgrims from the neighboring countries, especially Iran, are quarantined in quarantine centers. The highest affected population per square kilometers was observed in Lahore and Islamabad where it was greater than 5 persons per km 2 .

Impact of Weather and Government
Intervention on COVID-19 The monthly climatology of various meteorological variables over Punjab province is shown in Table 1. The results show that mean monthly air temperature remains the lowest (11.1 °C) in January and the highest in June (32.2 °C), whereas the relative humidity remains the highest (69.9) in December. However, the heat index is only observed during May-September, where the highest value is observed in July (36.1 °C). The UV index remains high during May-August, whereas the highest aerosols concentration (0.61) is observed during May, which is the direst (lowest humidity) month of the year. Based on the daily COVID-19 cases, the basic reproductive number was calculated for the Punjab province. It has found that basic reproduction number (R 0 ) ranged for 1.01-3.8 in the province with maximum at Lahore (1.01-4.2) as well as in Islamabad (1.01-4.1). The temperature has the highest range (38.4 to 40.9 °C) during 21st to 27th May and 14th to 20th June, 2020, whereas the relative humidity (RH) consistently increased (> 55%) after 1 st July, 2020 and it was the lowest (< 25%) between 17 and 28th May, 2020. Moreover, HI and UVI are found to be rapidly increased after 20th May, 2020 with a maximum 52.4 °C and 11.53 at Lahore respectively. The stringency index was the highest (96.3) during 26th March-14th April, 2020. The aerosol concentration was highest during 20th April-20th June, 2020 due to dry weather conditions and low humidity, which is the climate of the region. Most importantly, wheat harvesting remains at its peak which enhances the aerosols concentrations during this period. The low humidity increases the aerosols suspension in the atmosphere as compared to high humidity, where these aerosols are deposited rapidly in the ground.
Pearson correlation analysis has been done for temperature, RH, HI, and UVI with respect to R0, whereas aerosol concentration and stringency index with respect to daily positive COVID-19 cases over Punjab province (Fig. 7). The results show that temperature, RH, HI, and UVI is negatively correlated to R0. Moreover, the aerosols and stringency index were positively and negatively correlated to daily positive cases, respectively. The correlation is statistically significant at the 95% confidence level ( Table 2).
The high temperature is observed during May and June, which is known as the hottest months in Punjab, whereas the humidity remained lower than the rest of the months. The HI starts to develop from the mid of June till 20th August due to the increase in humidity. Moreover, the UV index is also observed high during summer (May-September), because the sun is high in the sky, and its rays arrive more directly and travel through much less zone than in winter where the sun is lower in the sky. These all meteorological parameters significantly impacted to lessen the COVID-19  Adnan et al. 2020 over the major cities of Pakistan. Moreover, the aerosol concentration was high during May and June due to which the COVID-19 cases increased rapidly during these two months. However, the COVID-19 cases declined due to the arrival of the monsoon system (July to September), which normally enhances the humidity and reduces the aerosols concentration in the atmosphere. The government timely and precisely intervention also favored the rising trend of COVID-19 over the Punjab province.
The results clearly show that COVID-19 daily positive cases and its basic reproduction number (R0) decreased   (Davis et al. 2016a, b). Similarly, the COVID-19 transmission significantly decreases due to increase in temperature and humidity as identified by many researchers (Abduljali and Abduljali, 2020;MA et al. 2020;Loeffelholz and Tang, 2020 etc.). Altamini and Ahmed (2019) identified that the MERS-CoV broke out in Riyadh, Saudi Arabia during summer due to low humidity. It has already been determined that higher humidity reduces the transmission and spread of respiratory droplets like influenza and settles them rapidly to the ground (Lowen et al. 2007).
The government made three interventions in terms of lockdown (partial, complete, and smart) during the COVID-19 onset to decline period. The COVID-19 cases were small in number during 16th to 22nd March, 2020 in which partial lockdown was implemented in which public/private schools, offices, parks, and shops were closed. The country was under a complete lockdown during 23rd March-14th April, 2020 in which free movement of the public was restricted and all the public/private schools, offices, parks, shops, airports, and public inter/outer city transport remained closed, while the groceries' stores and some government important offices with less staff were opened. Most importantly, the business hours duration was reduced from 8 h (9am to 5 pm) to 6 h (10am to 4 pm). Moreover, the track and trace program to identify the COVID-19 cases was launched. The smart lock was implemented after 15th April, 2020, and those areas were sealed where high positive COVID-19 was observed. The business centers, offices, markets, parks, schools, public transports, and shops were opened gradually. The daily business hours increased up to 9 am to 10 pm. The smart lockdown along with standard operating procedure (social distancing, use of mask, etc.) helped a lot to reduce the COVID-19 wave in the country as well as in Punjab.
The massive media campaign and awareness programs were also launched to bring comprehension to the general public. The mobile companies aired the messages (text and audio) regarding COVID-19 precautionary measures. More importantly, weather also played a significant role to subside the COVID-19 cases after the mid-June, 2020. Recent study shows heat index (temperature, humidity) and ultraviolet index played significant role to subside COVID-19 over major cities of Pakistan including Lahore (Adnan et al. 2020). The cities of Lahore and Islamabad are the hub of business, financial and government functionaries, and the majority of the population traveled there each day from the neighboring regions. Most importantly, the population density of these two districts is very high as compared to the rest of the districts, and these remained more vulnerable to COVID-19.

Conclusion
Punjab is the largest and densely populated province of Pakistan. Based on the above results, it is concluded that total 785,764 COVID-19 tests were conducted out of which 94,865 patients reported positive, 86,389 recovered, and 2179 died during 11th March-12th August, 2020 in the province. The COVID-19 peak period (daily positive cases > 900) remained during 27th May-8th July, 2020 over which 43 peak days were observed in the province. The weekly analysis show peak duration for the positivity rate was 10th to 19th week (13th to 14th July), whereas it was highest in 13th week (3-9th June). Growth rate was higher in the 2nd week (18-24th March), while maximum recovery rate was observed during the 16th to 20th week. Moreover, the mortality rate remained 0.37-3.66 with a peak in the 16th week (24 th -30th June). The districts of Lahore, Islamabad, Rawalpindi, and Multan were the highly vulnerable districts regarding COVID-19. Significant negative correlation (at 95% level) has observed for temperature (− 0.36), RH (− 0.33), HI (-0.52), and UVI (− 0.34) to R0, whereas daily COVID-19 cases were positively correlated to aerosols (0.30) and negative to stringency index (-0.48) over Punjab province. The highest increase in the daily positive cases was observed during May-June, 2020 when RH was the lowest and aerosol concentration was highest due to dry conditions. This study identified that the COVID-19 affected the urban areas more than rural areas, especially those where the population was congested and there was a high concentration of aerosols. Based on these results, it has been concluded that the COVID-19 may likely to increase during the dry period (high concentration of aerosols) of October-December, 2020 and post-spring season (April-June), 2021 in Pakistan. The industrial and economic growth has strongly affected around the world due to COVID-19 and special bail-out packages and relief operation were carried out to minimize its impacts to industrial, corporate, and commercial sectors as well as general public. It is highly recommended that the public should keep social distancing, use of face masks, and sanitizers, whereas government intervention should be made in high COVID-19 vulnerable areas in the form of smart lockdown that could reduce the COVID-19 transmission and help to fight the future outbreak.