Abstract

During the outbreak of an epidemic, it becomes significantly essential to monitor the effects of containment measures and forecast the outbreak, including the epidemic peak. Many countries have either implemented strict lockdown to counter the spread of coronavirus disease or taken necessary preventive measures across the world to reduce the outbreak of this epidemic war. Several epidemic models have been presented across the world to examine the effects of public health-related strategies on mitigating the spread of current infectious disease, yet no reputable model has been presented for Pakistan as well as other South-Asian developing countries as per the authors’ knowledge. In this research, an actual coronavirus prediction in Pakistan is presented, which may guide the decision-makers as to how this pandemic has spread across the country and how it can be controlled. Furthermore, in the absence of targeted medicines, the analysis helps to develop a precise plan for the eradication of the outbreak by adopting the calculated steps at the right time. The mathematical phenomenological models have been adopted in this study to predict, project, and simulate the overall affected cases reflected due to the recent outbreak in Pakistan. These models predict the expected growth, and the estimated results are almost well matched with the real cases. Through the calibration of parameters and analyzing the current situation, forecast for the appearance of new cases in Pakistan is reported till the end of this year. The constant level of number of patients and time to reach specific levels are also reported through the simulations. The drastic conditions are also discussed which may occur if all the preventive restraints are removed. This research quantitatively describes the significant characteristics of the spread of corona cases. It acknowledges and provides an understanding of a short-term and long-term transmission of coronavirus outbreak in the country as three evolutionary phases. Therefore, this research provides a pathway to cope with the emerging threat of a severe outbreak in developing and nondeveloping countries.

1. Introduction

The number of identified cases increased gradually in the province of Hubei, Wuhan city, China, since the novel 2019 coronavirus was discovered in late December 2019. The Chinese government has implemented major spread prevention measures, such as the quarantine of the city of Wuhan and shutting off roads for all Hubei cities by the end of January 2020. Amid this measure, the outbreak spread globally and the World Health Organization has declared it an epidemic [1, 2]. Consequently, a rapid rise in outbreak cases became a global emergency and a major health issue. Nowadays, the epidemic spread is a collective global public health concern, which constitutes serious political, economic, and social challenges. Efforts are being made to restrict and minimize the spread of this virus worldwide. The characteristics of corona are high infectivity, persistent incubation period, identification challenge, and complexity in how the disease is spreading. Unity is, therefore, recognized and considered as an essential necessity to fight against coronavirus during this period.

Various research studies have been conducted on the epidemiological study of modeling and forecasting [3] to tackle the spread of epidemic disease in order to provide not only the analytical predictions according to their size but also the ending phase of spread. Kuniya [4] has predicted the peak of the current epidemic disease through the application of the compartmental model “Susceptible-Exposed-Infected-Recovered (SEIR)” to provide the feasible insight for the Summer Olympics of 2020 to be held in Japan. Al-qaness et al. [5] applied Flower Pollination Algorithm (FPA) modified model together with the Salp Swarm Algorithm (SSA) in order to estimate the exact cases of corona in China for ten days. Koczkodaj et al. [6] predicted that the cases of corona outside China will rise by 1 million by a specific date on the basis of WHO heuristic situation reports. Wu et al. [7] have measured the outbreak of Wuhan and calculated the magnitude of the epidemic’s international and domestic threats to public health based on the social and antidrug spreading interventions. Younes [8] used the Lotka-Volterra model coupled with Kalman Filter Algorithm model. It was predicted that the epidemic could be controlled in a few weeks by adhering to social distancing and banning all types of traveling; otherwise, more cases will rise in upcoming days. Wang et al. [9] applied the Generalized Logistic Growth Model (GLGM), Richards Growth Models (RGM), and Subepidemic Wave Model. They conducted a modeling study in Zhejiang and Guangdong provinces of China and predicted that the ratio of cases is slowing down in both provinces.

Pinter et al. [10] used Susceptible-Infected-Resistant (SIR) model, Adoptive Network-Based Fuzzy Interference System (ANFIS), and Multi-Layered Perception-Imperialist Competitive Algorithm (MLP-ICA) model to conclude that fewer number of cases were reported in Hungary and projected that overall cases and mortality rate will significantly decrease in Hungary at the end of May 2020. Reno et al. [11] applied SEIR model and concluded that, in Italy, preventive measures would strictly be imposed for a long term in Lombardy due to an increase in the number of hospitalized patients but in Emilia Romagna fewer cases were reported than Lombardy, Italy. Sanchez-Caballero et al. [12] applied the Verhulst equation-based model and concluded that China had started complete lockdown when the number of cases was 571 while Italy and France imposed partial lockdown when 12462 and 5678 cases were reported, respectively. In addition, the increased number of cases in Italy or France compared to China is just due to not adhering to the early preventive measures and delaying in lockdown. To date, scientists, engineers, and social scientists have been making rapid use of analytical data that are beneficial for decision-making in public and health organizations and to predict local and global challenges. Some of the previous research studies investigated epidemic spread, based on the projection and predicting models, as mentioned in Table 1. Different nations have made different preventive measures to respond to various outbreak patterns. Such interventions measures and their effect on spread should be reported if an effective containment strategy is to be enforced. Worldwide, epidemic has a negative impact on the education system, public life, economy, and the social life.

The WHO [16] officially declared Global Health Emergency for outbreak on January 30, 2020. Within two months of that declaration, the first major pandemic over a century has been driven by coronavirus [2]. The number of valid cases of COVID-19 and mortality continues to grow rapidly, with serious impacts on local and global economies without a definite end in sight. The governments of many countries have adopted lockdowns and other measures to limit the transmission of epidemics in several respects. Over the past few months, several models have been developed at the national and provincial levels to examine the possible effects of coronavirus transmission mitigation strategies. Pakistan is an independent state having a population of 220 million [17]. The first case of epidemic was reported on 26 February [1]. Most of the cases were transported from foreign travelers and not triggered by local transmission. Pakistan has therefore made intensive efforts to eradicate and quarantine the pandemic outbreak, maintain an efficient and robust health system, and reduce the time needed to reach an epitome of epidemic. The containment measures presented in Table 2 were taken from Pakistan’s official releases [1, 18].

The intervention measures are aimed at reducing the contact between people, thus limiting the outbreak. These interventions include decreasing the disease rate among the elderly and those with chronic health conditions, mitigating the risk of spread among care and medical workers, significantly reducing mortality rate, and maintaining the health-care system. The Pakistani government shall take suitable action as nonmedicinal initiatives to minimize the outbreak impact [2]. In addition, Polymerase Chain Reaction (PCR) tests by the Pakistan Health Ministry were conducted to diagnose all international visitors for transmission [32]. As a result of the containment measures, the number of diagnosed infected cases of corona gradually increased. Citizens coming back from Iran and the Kingdom of Saudi Arabia were mostly infected and later on, they became the cause of the spread of the transmission [18, 19]. Nevertheless, some of the cases that have been reported recently are locally transmitted. This resulted in a dramatic rise in the early days in the number of cases.

The Government of Pakistan introduced many protective actions to safeguard the community, including sanitary education, patient isolation, and quarantine facilities, closing of key entrance points on state borders, closure of educational institutions, business centers, and the start of clinical testing and epidemiological monitoring. Around each other, such efforts require measured steps for optimal allocation of precautionary measures, structure of manufacturing processes, health supplies, and gradual stabilization of economic growth across the region. In the present contribution, an analytical approach is developed based on phenomenological models (DEGM, DG2M, DGLM, and DGRM). The research effort has allowed an important forecast to control the country’s rapidly increasing pandemic in both short- and long-term zones. In addition, the possible impact on pandemic spread in Pakistan is analyzed, keeping in mind existing preventive measures. This research may classify the original contributions as (1) to present the mathematical phenomenology approach to simulate the outbreak spread in Pakistan, in order to estimate the scale of the epidemic, reproductive number, and final phase of the outbreak, (2) to project the evolutionary stages and provide an understanding of controlling its spread in the short-time and long-term domains, (3) to assess the current government’s quick response and their effectiveness in containing outbreak spread in Pakistan with various preventive measures, and (4) to provide a prediction of outbreak distribution amid the repatriation plan for the political authorities of Pakistan for effective intervention.

2. Materials and Methods

2.1. Phenomenological Mathematical Models

The phenomenological models (Exponential Growth Model [37], Generalized Growth Model [38], Generalized Logistic Growth Model [39], and Generalized Richards Growth Model) [40] are articulated by the following differential equations, respectively [4143]:

The mathematical models were thoroughly taken into account to describe past epidemics. The total estimated cases at time t are shown by C(t). Equation (4) then enhances the initial Richard Model of Growth [44] that explains the outbreak [4549]. In this case, r is the early growth rate; the overall size of the outbreak (spread) is described by K; the term describes the parameters which helps to extract the epidemic growth of the model; and the exponent’s measures are labeled by α which shows the variation from the single logistic curve structure. Equation (4) is the same as the original form of Richard formula ; that is,

In GL model, the value of α and is restricted to 1. Although data was built in terms of confirmed cases, modifications, and fatalities, the variations in term C can be taken into account over a unit time interval. The rate of change is expressed in C, that is, in C(t), differently as . The changes in positive patients are . So, in a discrete form, equation (4) is remodeled aswhere the overall cases are reflected by on the nth day. The DGRMs, which are designed for mapping the pandemic, link the accumulative data on day nth with (n + 1)th day. The stable points in equation (6) for a DGRM are , which means no patient was tested positive, and , which means the number of patients reaches the final size. The stable levels are similar to those reached by the GRM. Further, we write the DG2M [50] as

It brings an exponential increase for the value of , which is equal to 1, which is , linear growth for the parameter relative to n as , and for the subexponential growth.

3. Analysis and Results

In the absence of medicines to control and monitor the outbreak, strict measures need to be adopted at the right time; therefore, besides the data evaluation, a clear public awareness, life, and economic risk management are essentially required. In all the provinces of Pakistan, regular pandemic growth identifies an overall increase in reported corona recovered and active cases, including mortality rate (death) which is reflected through blue color (dots) in Figure 1, which is tested using DG2M, Discrete Logistics Growth Models (DGM), and DGRM mathematically. Such models have been used in the United States, Arab, Asian, and European countries to comprehend the epidemic [3739, 4143, 51].

3.1. Epidemic Spread and Evolutionary Phases in Pakistan
3.1.1. Phase 1: 26 February to 15 March

On 26 February 2020 in Pakistan, the first two confirmed cases have been registered. The pandemic expands significantly throughout the country due to having a lack of the governance and public awareness and absence in the complete lockdown strategy. Comparing the data (real time) with numerical model results confirms the exponential growth of the outbreak at an absolute rate of 0.85. Figures 1 and 2 reflect the DG2M results in the red-dashed line.

Figure 1 illustrates the growth of reported coronavirus patients in the country measured against the period of time (number of days, n). In the graph, blue dots reflect the real-time data. The DG2M curve for the value of (reflecting the red-dashed line) and r = 0.3 shows the total number of confirmed cases of COVID-19 in the first twenty days. In the first twenty days, overall reported cases follow DG2M curve for r = 0 : 3 and (red-dashed line). The DG2M curve for the value of (shown in blue color) and r = 0.74 represents the real-time data in phase 3.

3.1.2. Phase 2: 16 March to 31 March

In the next 15 days as shown in Figure 2, the entry for zairin from Iran, Tablighi Jamaat members, and the country citizens incoming through overseas on special flights induced the sudden rise in the number of confirmed cases. Lack of testing and quarantine facilities, returning of zairin to their home cities, and unexpected strain of people with previous international traveling background put significant pressure on economic resource limitations and health-care setup. Poor governance has led the country to combat the increasing epidemic cases a long way forward.

3.1.3. Phase 3: 03 April to 08 May

The governance tries to combat the outbreak across the country in the next 38 days of quarantining by combatting public movement, implementing social distancing, surveillance and controlling, shutting down educational institutions, playgrounds, and hotels, and restricting air and road travels. The exponential growth has been declined from 0.85 (value of phase 1) to the value of a new exponent 0.74 which represents the DG2M curve in blue line in Figure 2.

The confirmed cases in Pakistan have been represented by the log plot in Figure 2, in which along with days n is plotted on 26 February 2020. The blue dots reflect the real-time data. In the first twenty days, the number of confirmed cases is followed by DG2M curve with which is drawn in red-dashed line in the graph and the value of r = 0.3. The DG2M curve meets/followed the real-time data for denoted by blue line and (1st value) and r = 0.74 to the value of a new exponent 0.74 which represents the DG2M results in blue line in Figure 2.

The measures taken in phase 3 to remove the pandemic effect in the country resulted in the identification of 28736 cases, including 20291 active cases, 636 deaths, and 7809 cases rescued. The statistical analysis shows a red- and blue-dashed lines that correlate in real time to the second week of May. It reflects that, at stage 1, the number of reported epidemic cases today was the same without any lockdown and public awareness drive. In phase 2, the country allowed foreigners to have a previous traveling background to enter only after proper screening and passing through the quarantine period. These lockdown measures have damaged the country’s economy at stages 2 and 3. At the end of the day, they would have to pay the cost for entry of people with foreign history. It also revealed the poor situation in Pakistan in the early stages to cope with the pandemic.

3.1.4. Phase 4: 09 May to 15 May

The evaluation of coronavirus reported cases centered on DGRM and interpretation to 31 May indicates that the pandemic scale is expected to upsurge around sixty thousand cases over the next twenty-two days when the lockdown is implemented (as reflected in black-dashed narrow line in Figure 3). Furthermore, it indicates that we can significantly hold back the outbreak in the country up to August 2020 by accelerating lockdowns at this point. An enforced lockdown in this crucial stage will help to enable the country to handle the outbreak in the upcoming two months rather than relaxing lockdown. This limits the size of the epidemic to approximately 100,000 cases as shown in Figure 4. If the lockdown stops fully at phase 4, the scale of the outbreak will increase to 120,000 coronavirus confirmed cases in the next twenty-two days as graphed by the red line in Figure 3. The number of COVID-19 confirmed cases is to be somehow between sixty thousand and one hundred twenty thousand if the lockdown is partially released.

Figure 3 exemplifies the growth of COVID-19 confirmed patients followed by the DG2M curve (in black-dashed line) for K = 100000, α = 2, , and r = 0.74 plotted against the number of days, n. The blue dots reflect the real time and the red line shows the DG2M results for and r = 0.3.

3.1.5. Phase 5: 15 June to Continue

In an account of the holy Ramadan days and the required occupational needs of those days, the planned lockdown in phases 2 and 3 is relaxed in phase 4. In the next month, a detailed and calculated approach is required for pandemic management. The country's pandemic size shall not increase by a lockdown scenario from May 09 to June 31. Through perceiving the DGRM results, a strict lockdown in June would carry around three hundred and fifty thousand confirmed cases in September, with a rise of less than one percent in the number of new cases per day as shown in Figure 4 (black curved line). At the end of August, the outbreak peak is expected to beat by about 200,000 reported cases. A compromised lockdown would allow epidemic cases at phase 4 to be managed by the end of August, with the projected epidemic scale about 250,000 cases as reflected in Figure 4 by the big black-dashed line.

Figure 4 describes the growth of confirmed patients followed by the DG2M curve (by red line) for and r = 0.3 plotted against the number of days, n. The blue dots reflect the real-time data. The black-dashed line shows the results of DGRM curve for K = 100000, , r = 0.74, and α = 2. Furthermore, the black-big-dashed line and the black line reflect the DGRM results for K = 250000 and K = 350000.

4. Discussion

Various scenarios were analyzed to provide policymakers with insights into the implementation of the appropriate intervention to mitigate the spread of the epidemic in Pakistan. The study used mathematical phenomenological models to estimate the daily and weekly peaks, total infected cases, and the estimated duration of the outbreak. In addition, the effects of preventive variation in the number of susceptible and infected populations were also investigated. Phenomenological models evaluated the effectiveness of methods of control such as social separation and lockdown. The Government of Pakistan has taken various control measures, as outlined in Table 2. In response to the outbreak, the measures are gradually increased to more restricted ones. The first measure was the NAP and the other measures include sealing the Taftan and Afghanistan borders, closure of school and universities all over the country, applying half and full lockdown, postponing the PSL, active participation of the Pakistan CAA in the closure of all international and national flights, manufacturing of hand sanitizers and face masks as per WHO guidelines, extension of lockdown, and relief package from the Government of Pakistan for the needy people. According to the Pakistan Bureau of Statistics, as of 1 January 2019, the proportion of people between 5 and 24 years of age accounts for approximately 25% of the total population [52]. Even though the early analysis indicating the low percentage of this age, the issue is very important because it can spread the disease to other age groups. In fact, a new WHO study disclosed a steady growth in the proportion of cases aged 20–39 [1]. The evolutionary phases describe the real scenario of cases: in the first phase, the first two patients have been diagnosed and due to having a lack in applying lockdown and people awareness operations, the pandemic had a blowout in the country in a very drastic way; in phase 2, the study projected that new cases emerged due to the arrival of citizens from foreign countries, thus, putting colossal pressure on the budgetary sources and health-care proper setup and making it difficult for the country to combat an increased number of coronaviruses cases; in phase 3, the study predicted that exponential growth of cases had decreased from 0.085 (phase 1 value) to 0.074 (phase 3 value) due to the enforcement of quarantine strategy, combatting public movement, implementing social distancing, surveillance, shutting down educational institutions, and restricting train and road traveling. The elimination of the outbreak in Pakistan in phase 3 resulted in 28736 cases, including 20291 active, 7809 recovered, and 636 death cases. This indicates that, in phase 1, the rate of total cases has been increased, leading to a lack of public awareness and applying lockdown. The cases increased in phase 2 as the country has allowed the entry of people having foreign travel history without any adequate checkup and quarantine. The study predicts that, even after enforcing complete lockdown in phase 3, the number of corona cases would cross 60,000 at the end of May and it will further go on and reach 100,000 cases till the end of June. Either by lifting the lockdown or with smart lockdown, cases are likely to increase between 60,000 and 120,000. So instead of relaxing the lockdown, a strict lockdown is obligatory to combat the outbreak effectively. The study projects that, even after enforcing strict lockdown in the month of June (in phase 4), the confirmed cases will rise to almost 350,000 at the end of August or the start of the month of September with 1% increase in the newly identified cases per day. This number is likely to decrease after the end of December 2020. Although these estimates may vary with the passage of time, it will really help us to observe the most influential factors that cause an increase in the epidemic spread. On the basis of this analysis, competent authorities may design the most effective strategies in order to control the epidemic.

5. Conclusion and Recommendations

In this study, phenomenological mathematical models were used to assess the feasibility of the appearance of the recent outbreak in Pakistan as well as measuring the definitive number of patients according to the current situation. By comparing the model outcomes with confirmed cases, it has been established that the estimated values have realistic correspondence with the confirmed numbers. If the current pattern continues, Pakistan has to bear the burden of around 350,000 cases till the end of August or 1st week of September 2020 with an increase of 1% of newly identified cases per day. The situation may vary by the deviation in transmission rate, death rate, recovery rate, and further implementation of social distancing. The Government of Pakistan has relaxed the lockdown under some restrictions since 9 May, in order to support the business activity and reduce the economic crisis. Currently, the shopping malls, public markets, business centers, trains, and local transport are opening across the country. Since the opening of lockdown, the number of COVID-19 cases is likely to increase exponentially due to more chances of transmission among peoples. As a precautionary measure, the government has launched a public service campaign with the collaboration of media during the relaxation in lockdown periods, in order to provide preventive measures through video messages (1-2 minutes). Meanwhile, the emphasis is laid to strictly act upon prevention measures during the ongoing pandemic outbreak. The study recommends that health experts, motivational speakers, administrative authorities, media, and physiological experts should arrange online lectures and provide health education, to motivate the public to follow all the essential protective measures for their own and other people safety. The study recommends imposing a strict lockdown from June to July 2020 to control the contagious disease across the country. There is also a need to focus on the effective preparedness of the health sector and enhancing the technical expertise of all front line warriors. Besides, optimal allocation of preventive or precautionary measures, judicious utilization of available resources, and construction of all productive initiatives would get the country away from the global pandemic. On the other hand, mismanagement and delayed responses are expected to increase the spread of the pandemic and result in excessive economic loss in the coming months.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

This research was fully funded by the China Scholarship Council and National Natural Science Foundation of China (71971196).