Structural Modeling on the Determinants of Effectiveness of SOPs Containing COVID-19 in Mass Gatherings

The study is aimed to analyze the determinants of the effectiveness of SOPs in mass gatherings for containing COVID-19. The overall design of the study involves a literature review, data collection by field survey, structural modeling, and analysis. The study is built on the experts' opinion of a focus group (representing people who recently participated in and are responsible for mass gatherings). The study uses the discussion of the literature review to identify the determinants, interpretive structural modeling (ISM) for developing and analyzing a structural model, and Matrice d'Impacts Croises Multiplication Appliquée a un Classement (MICMAC) for corroboration of results of the ISM/classification of determinants. From the literature review, a list of determinants is generated and verified by a panel of experts. The results of the ISM revealed that the determinants “legal environment of the country,” “practicability of SOPs,” “perceived benefit of adapting SOPs,” and “possibilities of avoiding gathering” occupied the top of the model, therefore, they are less critical determinants, whereas “nature of gathering” occupied the bottom of the model, and is thus the most critical determinant. The remaining determinants form the middle of the model, and are therefore moderately severe. The results of MICMAC show that the determinant “perceived benefit of adapting SOPs” is dependent, “nature of gathering” is independent, and all others are linked. The results of MICMAC implicitly substantiate the findings of the ISM. The overall results of the study show that “nature of gathering” is the key determinant. This research does not require a priori theory since it is a theory-building study that uses an inductive approach. It is based on real data and it is useful for local authorities, organizers, participants (attendees/visitors) of mass gatherings, health officials/regulators, researchers, and the community at large. This study has fundamental importance for planning and preparing for such events while ensuring the minimum risk of COVID-19 transmission.


INTRODUCTION
COVID-19 has wrought havoc around the globe, with 200,840,180 confirmed infected cases and 4,265,903 confirmed deaths worldwide (World Health Organization, 2021). The pandemic has pushed mankind to reshape the global order (Akon and Rahman, 2020;Rehman and Ahmad, 2020). The outbreak of the COVID-19 pandemic in 2019 (now considered as coronaviridae family) has redirected the thinking of policymakers. It also has forced governments and communities to revisit governance, social, and business practices. Discouragement of gathering and social mixing has emerged as one of the measures to contain the transmission of novel disease. Evaluation of this precaution, qua reality, has gained fundamental importance because in certain cases, social distancing and its SOPs are practicable, but in certain cases it is difficult if not impossible. Particularly in mass gatherings, the practicability of the SOPs of social distancing is highly questionable. Mass gatherings are events that are capable of amplifying the transmission of COVID-19 and potentially disrupt the country's response capacity to the COVID-19 pandemic. Mass gathering includes cultural events, religious gatherings, sports events, entertainment concerts, political procession/assemblies, congregations, and crowds. The gatherings may be indoor or outdoor, since as far as the SOPs are concerned they are questionable in both cases. It has become highly important to analyze and uncover the structures of the relationships of the determinants of the effectiveness of SOPs of social distancing to contain COVID-19, particularly in mass gatherings. Right from the outbreak of this disease, researchers started attempting to make some sense of the system. For example, Carlin et al. (2021) evaluated the impact of the COVID-19 pandemic on mass gatherings for hockey and basketball games in North America, Hughes et al. (2020) conducted a study on the framework of sports in Australia, Mat et al. (2020) evaluated and found COVID-19 contamination in Malaysia because of mass gatherings, Nyasulu et al. (2021) argued that the COVID-19 outbreak emerged in Malawi-Southeastern Africa as a result of sociopolitical mass gatherings, Teixeira (2020) highlighted a possible transmission of COVID-19 due to mass gatherings in Brazil, Yezli and Khan (2020) highlighted possible health risks from the transmission of COVID-19 due to mass gatherings in Jeddah-KSA. A lot of research has also been conducted on religious mass gatherings and religious tourism with the perspective of COVID-19, including but not limited to the studies of Ahmed and Memish (2020), Atique and Itumalla (2020), Ebrahim and Memish (2020), Hoang et al. (2020), Mat et al. (2020), Mubarak (2020), Mubarak and Zin (2020), Quadri (2020), Rodriguez-Morales et al. (2020), VanderWeele (2020), Yezli and Khan (2020), Hashim et al. (2021), Hsu et al. (2021), Sokhna et al. (2021), and Yezli and Khan (2021). Although some of these studies address the issue partially or indirectly, there is still a clear gap in the literature regarding real determinants of the implementation of SOPs for COVID-19 in mass gatherings. This study is designed to fill this gap by way of analyzing the determinants of the effectiveness of SOPs in mass gatherings for containing COVID-19. The objectives of the study are: (i) to identify the determinants of the effectiveness of the SOPs of containing COVID-19, (ii) to develop and analyze a structural model of the determinants, (iii) to corroborate the results of the structural model, and (iv) to classify the determinants. For achieving these objectives, an array of methodologies was considered (Shaukat et al., 2021) and interpretive structural modeling (ISM) augmented with Matrice d'Impacts Croises Multiplication Appliquée a un Classement (MICMAC) was chosen as the methodology. It has many advantages over its rival methodologies. It is an exploratory approach, appropriate for making some sense of the early stages of systems. It takes a smaller amount of data and converts fragile and unclear mental models into graphical models to simplify complex issues. It outperforms the statistical methodologies that take large statistical data but contribute little toward the literature. The rest of this study is divided into the Literature Survey, the Methodology, the Modeling, Analysis, Results, and Discussion, and the Conclusion.

LITERATURE SURVEY
A survey of the literature is important: to set out the context of the study, to account for the current literature, and to acknowledge contemporary work. A thorough survey of the literature was conducted. We explored the online databases of Wiley-Blackwell, JStor, Emerald, Springerlink, Sage, Taylor and Francis, Elsevier (Science Direct), etc. Keywords used for the search include: mass gathering, COVID-19, social distancing, SOPs of COVID-19, coronavirus disease, novel disease of coronavirus, determinants of SOPs of COVID-19, social distancing, etc. An influx of published/unpublished literature is found to be available online. A large number of articles were also found as "accepted for publication." Hundreds of relevant papers were reviewed, and only a small amount of literature was found directly concerning the issue. It is not out of place to report some of the studies that can provide foundations for understanding the current phenomenon. McCloskey et al. (2020) proclaimed that the banning of mass gathering events could result in plummeting further global spread of COVID-19. Memish et al. (2020) concluded that the COVID-19 epidemic demands internationally combined, lucid, and cooperative measures by communities, individuals, institutions, commercial bodies, and regulators to suspend gatherings and mitigate the transmission. Nyasulu et al. (2021) has argued that mass anti-government demonstrations and cross-border immigration are causing increased COVID-19 disease and mortality, particularly among people aged 50-59 years. Piovani et al. (2021) studied the longitudinal data of 37 countries to analyze the effect of banning mass gatherings on COVID-19 cumulative mortality. The results affirmed that a 1-day delay in exercising the banning of mass gatherings is associated with an increase in mortality of 6.97%. Chakraborty and Maity (2020) asserted that COVID-19 has had a devastating impact on society which could be controlled by restricting mass gatherings. Through our review of the literature, an initial list of determinants (Table 1) of the effectiveness of the SOPs for mass gatherings for containing

10
Possibilities of avoiding gathering How much is it practically possible for a person to avoid gathering? Ahmed and Memish (2020) 11 Seriousness of SOPs implementers How much is government or implementers (whatsoever) serious in implementing the SOPs? Hsu et al. (2021) COVID-19 was prepared and placed before a panel of experts for ratification. For the verification of the initial list of determinants, a panel of experts consisting of 10 Ph.D. doctors was constituted. The experts are individually approached and asked to (i) evaluate the determinants, (ii) add determinants if they were not already on the list, (iii) merge two or more determinants if appropriate, (iv) delete determinants that seem irrelevant, (v) modify determinants if necessary, and (vi) vote for the inclusion of determinants in the study. The experts added one of the determinants (available physical space at gathering place). Although the determinants have been gathered from the literature, it is still important to get our findings ratified by a panel of experts because the literature collected is diverse, meager, and assorted. Anecdotal evidence of the determinants demands that we verify qua relevance, appropriateness, and sufficiency. In this way, the list of determinants is finalized for inclusion in the study ( Table 2).
In the case of determinant 4 ("existence of penal clauses in SOPs, " marked gray in Table 2) a tie occurred in the voting of the panel of experts. The authors, therefore, themselves evaluated the determinant and decided to include it in the study since it has fundamental importance and logically qualifies as part of the phenomenon under study.

METHODOLOGY
This work is an exploratory inductive study which embraces a literature review, data collection survey, structural modeling, and analysis. The methodology of the study focuses on the discussion of the literature review for the identification of the determinants, ISM for developing and analyzing a structural model, and MICMAC for the corroboration of the results of the ISM/classification of determinants. Interpretive structural modeling is very useful a tool for qualitative analysis. It is appropriate for analyzing complex interdependent relationships in rapidly changing and entangled situations. It is capable of articulating complexities into visible, well-defined models with graphical presentations (Warfield, 1974;Sushil, 2012;Li et al., 2019). Since the COVID-19 pandemic is a new and complex phenomenon that is being explored, the method is not only justified but also the most appropriate. It outperforms statistical approaches in the cases like the one at hand. The population of the study consists of the people participating or responsible for mass gatherings. The sampling design of the study is based on a non-probability type of sampling that is a focus group (panel of experts) from within the population. The sample has been drawn on the basis of pre-determined criteria. The size of the sample consists of 10 experts on a panel. There are various techniques for finalizing the variables to be included in the investigation. Common techniques used include a literature review or expert opinion (Majumdar and Sinha, 2019), case study (Li et al., 2019), Delphi method (Bhosale and Kant, 2016), exploratory factor analysis (Li and Yang, 2014), presumption of the authors or meta-analysis (Lohaus and Habermann, 2019), idea engineering workshop and brainstorming session (Kumar et al., 2013), and interview content analysis (Xiao, 2018). This study used a literature review because this method has advantages over its rivals. It provides a foundation of knowledge about the phenomenon under study, prevents duplication, gives credit to other researchers, highlights questions left from other research, provides justification of inclusion or exclusion of the variables in the study, and indicates the potential relationships among variables sets at the outset of the study. The data have been collected from the experts using a matrix type questionnaire that contains precise instructions for respondents about completing the questionnaire. The respondents were instructed as follows: (i) fill only white cells (ij part only) and do not fill black and gray cells ( Table A1 in Appendix), (ii) contextual relationship = leads to, (iii) enter V when the row leads to the column, (iv) enter A when the column leads to the row, (v) enter O when there is no relation between the row and the column, and (vi) enter X when Possibilities of implementation of penal clauses Possibilities of implementation of penal clauses V A

10
Possibilities of avoiding gathering A

11
Seriousness of SOPs implementers row and column lead to each other. The study uses the classical form of ISM that progresses stepwise and develops a hierarchical structural model by way of exploiting the binary matrices. It also uses MICMAC for the classification of the determinants and the analysis of their driving power and dependence power on a Cartesian plane.

Panel of Experts
To generating the primary data, a panel of experts was constituted according to the norms of ISM studies. Panels are constituted where the data does not exist, data collection is expensive, available data is limited, and/or the data is unreliable. This is the case in our study, therefore, we find it appropriate to form a panel of experts as a source of eliciting the primary data. This panel of respondents is different from the aforementioned panel used for the verification of the factors. There are two types of respondent panel i.e., homogeneous and heterogeneous. The size of the panel depends on the type of panel. In the current study, we intended to collect data from people responsible for and/or participating in different types of gathering, the context demands a heterogeneous panel. Therefore, the study uses a heterogeneous panel. The usual size of a heterogeneous panel varies from 8 to 12 experts. In this study, the panel of experts consisted of 10 experts (Jena et al., 2017). When recruiting experts for a panel, theoretical knowledge of the phenomenon under study, expert knowledge, and practical experience are the key criteria. In this case, the experts were recruited on the basis of their knowledge and recent experience of some familial, political, religious, funereal, or accidental mass gathering. Apart from this criterion, all the respondents are well-educated (minimum university graduates) and held responsible positions within organization, family, and society. The experts were identified, a rapport was developed with them, then they were invited to participate in this study and asked to fill a questionnaire. It took us more than 2 months to complete this process. There are different methods to elicit the primary data from the minds of the people. These methods include Delphi, brainstorming, nominal group technique, repertorygrid interview technique, matrix type questionnaire, laddering interview, one-to-one, face-to-face in-depth interview, triadic sorting, and elect alternatives (VAXO) for every pair of relations (Niazi et al., 2019;Qazi et al., 2020). This study used the one-to-one, face-to-face interview approach coupled with elect alternatives (VAXO) for every pair of relations. The data were collected individually from each expert in the ij part of the matrix   using VAXO symbols. The data was aggregated in MS Excel sheet applying the "count if " function and "majority rule" in the case of every pair. This panel was involved at two stages: stage one at data collection and stage two at model verification. After developing the model, the panel of experts was approached again to review the structural model logically, theoretically, and conceptually.

ISM Model Building
From the aggregated data, a structural self-interaction matrix (SSIM) is developed ( Table 3). The SSIM is converted into an initial reachability matrix ( Table 4) by applying the classical rules for converting an SSIM into binary codes.
Every 0 on the initial reachability matrix is evaluated from the view point of transitive relations and on the basis of evidence of transitivity, some of the 0 s are replaced with 1 * , and the reachability matrix is revised to form the final reachability matrix (fully transitive reachability matrix) Table 5.
The final reachability matrix is partitioned through the iteration method, exploiting the elementary concepts of set theory (Warfield, 1973(Warfield, , 1974. The iterations have been performed following the classical rules of iteration (Warfield, 1973(Warfield, , 1974 and presented below in Tables 6-9.
Hierarchies that were determined at different levels as a result of the iterations were extracted in a conical matrix and a digraph was prepared but, being an optional step in the ISM procedure, not reported here for brevity. The ISM process is presented in a condensed form as Table 10.
The model captured on the diagonal of the conical matrix was converted into a graphical representation as Figure 1.
Close observation of the model revealed that the determinants coded as 6, 7, 8, and 10 fall on Level I. The determinants coded as 2, 5, 9, and 11 fall at Level II. The determinants coded as 3 and 4 fall at Level III and the determinant coded as 1 falls at Level IV. The model was generated from the conical matrix in the classical format of bottom-top. The bottom of the model is considered to be the most critical and the top to be the comparatively least critical. In this way, the factors that occupy the bottom of the model are critical determinants of the adoption

MICMAC Analysis
In order to corroborate the findings of the ISM method, MICMAC analysis (Godet, 1986) was performed, as shown in Figure 2. Following the scale centric approach, we divided the Cartesian plane into four quadrants i.e., independent, autonomous, dependent, and linkage.
Close observation of the driving-dependence diagram (Figure 2) revealed that the determinant coded as 8 is categorized in the dependent quadrant, the determinant coded as 1 is categorized in the independent quadrant, and determinants coded as 2, 3, 4, 5, 6, 7, 9, 10, and 11 are categorized in the linkage quadrant, while no determinant falls in autonomous.

Results
The objective of the study is to analyze the determinants of the effectiveness of SOPs in mass gatherings for containing COVID-19. During the course of the study relevant determinants were identified, regarding which the data was collected and ISM and MICMAC were applied. The results of the literature review revealed that there are a total of 11 important determinants of the effective implementation of SOPs for containing COVID-19 ( Table 1). The results of the ISM show that determinants legal environment of the country (6), practicability of SOPs (7), perceived benefit of adapting SOPs (8), and possibilities of avoiding gathering (10) fall at Level I (top level). The determinants that occupy the top level of the model are considered less critical. The factors at the top of the model are dependent or driven by other factors. The determinants type of gathering community (2), campaign of awareness regarding SOPs (5), possibilities of implementation of penal clauses (9), and seriousness of SOPs implementers (11) fall at Level II (middle level), and the determinants available physical space at gathering place (3) and existence of penal clauses in SOPs (4) fall at Level III (middle level). The determinants that occupy the middle of the model are moderately critical and are mediators or moderators. The determinant nature of gathering (1) falls at Level IV (bottom level). The determinants that occupy the bottom of the model are considered as the most critical factors, and they are independent and drivers and they affect the upper levels critically. The bottom of the model deserves the immediate attention of stakeholders, particularly policymakers. In the current study, "nature of gathering" is the most critical factor, deserving immediate attention of health policymakers. Results of MICMAC are as follows: Dependent: the factors (in this case the determinants) that fall in the dependent quadrant have a high dependence power but a low driving power, which    They deserve to be eliminated from the system but since they have a few powerful links they are included in the study. In ISM models, most of the time they stand disconnected from the others. No determinant falls in the autonomous quadrant, so all the determinants in the study are highly connected to the system and are very much part of the system under study. Linkage: the factors that fall in the linkage quadrant are considered ambivalent, agile, unsettled, or unbalanced. If most of the factors are classified into the quadrant namely 'Linkage' in MICMAC analysis, it necessarily means that the system under study is not settled and it is trying to make some sense. The determinants "type of gathering community (2), " "available physical space at gathering place (3), " "existence of penal clauses in SOPs (4), " "campaign of awareness regarding SOPs (5), " "legal environment of the country (6), " "practicability of SOPs (7), " "possibilities of implementation of penal clauses (9), " "possibilities of avoiding gathering (10), " and "seriousness of SOPs implementers (11)" fall in the linkage quadrant. The overall results of the study show that "nature of gathering (1)" is the key determinant. The results are summarized in Table 11. The results of MICMAC analysis affirmed the results of the ISM. According to the ISM, "nature of gathering (1)" (gray highlighted in Table 11) occupies the bottom of the model and is therefore important/independent, and this result is also confirmed by MICMAC analysis as this factor is categorized in the independent quadrant, having a high driving power and low dependence power. Other results are also almost the same in both methodologies.

Discussion
Since the main objective of the study is to analyze the determinants of the effectiveness of the SOPs for the COVID-19 pandemic in mass gathering, it is a fundamentally important and research worthy topic in today's world. Therefore, it is worthwhile to discuss the results of the study in contrast with the current literature. We have identified representative determinants, form them into a hierarchy, imposed a structure on interactive relations among them, and categorized them into different quadrants based on their driving-dependence. In the literature, there are few studies that also have addressed this issue, partially or remotely. A comparative analysis of the studies closer to the current study is represented in Table 12. Yezli and Khan (2020) conducted a theoretical study in the Kingdom of Saudi Arabia and emphasized the use of the "Jeddah Tool" as a standardized guide for handling hazards. It is a suggestive study without any empirical evidence. El Alaoui (2020) conducted a study in Japan focusing on the possibilities of estimation of infections in mass gatherings. It uses probability theory and illustrates the approximation of infections as a result of mass gatherings. Memish et al. (2020) advocated for suspending mass gatherings in view of the possible danger of rapid transmission of COVID-19. It is a theoretical study that builds logical arguments. Although these studies directly examine the phenomenon of mass gathering, the current study addresses altogether the different aspects of mass gathering. This study is different from them in context, in methodology (design, population, sampling, analysis, etc.), and in results.

CONCLUSION
This study has a fundamental importance for planning and preparing for mass gathering events while ensuring the minimum risk of COVID-19 transmission. The study aims to analyze the determinants of the effectiveness of SOPs in mass gatherings for containing COVID-19. The problem of unclear determinants of the implementation of the SOPs for preventing coronavirus, especially in mass gatherings, is studied in this research. The overall design of the study involves a literature review, data collection survey, structural modeling, and analysis. As a methodology, discussion of the literature review is used to identify the determinants, ISM is used to develop and analyze a structural model, and MICMAC is used for corroborating the results of the ISM/classification of determinants. The results of the literature review reveal that there are a total of 11 important determinants of effective implementation of SOPs of COVID-19 containment ( Table 1). The results of the ISM show that the determinants "legal environment of the country (6), " "practicability of SOPs (7), " "perceived benefit of adapting SOPs (8), " and "possibilities of avoiding gathering (10)" fall at Level I, the determinants "type of gathering community (2), " "campaign of awareness regarding SOPs (5), " "possibilities of implementation of penal clauses (9), " and "seriousness of SOPs implementers (11)" fall at Level II, tge determinants "available physical space at gathering place (3), " and "existence of penal clauses in SOPs (4)" fall at Level III and the determinant "nature of gathering (1)" falls at Level IV. The results of MICMAC show that the determinant "perceived benefit of adapting SOPs (8)" falls in the dependent quadrant, the determinant "nature of gathering (1)" falls in the independent quadrant, no determinants fall in the autonomous quadrant, and the determinants "type of gathering community (2), " "available physical space at gathering place (3), " "existence of penal clauses in SOPs (4), " "campaign of awareness regarding SOPs (5), " "legal environment of the country (6), " "practicability of SOPs (7), " "possibilities of implementation of penal clauses (9), " "possibilities of avoiding gathering (10), " and "seriousness of SOPs implementers (11)" fall in the linkage quadrant. The overall results of the study show that "nature of gathering (1") is the key determinant.

Contribution of the Study
The study contributed a list of determinants, an ISM model, a MICMAC diagram, and simplified information about the inter-relationship of the determinants toward the contemporary literature.

Practical and Theoretical Implications
This study is useful for local authorities to regulate mass gatherings in local communities since it develops their understanding of the inter-determinant relationships. It is useful for organizers and participants (attendees/visitors) of mass gatherings because it provides a lot of information for understanding the implications of following/not-following the SOPs of COVID-19 containment. It is helpful to health officials/regulators since it identifies key determinants, provides information about the relationships of these determinants, and prioritizes them on the basis of their importance of handling. It is also helpful for researchers by offering research frameworks since it categorizes the determinants into independent, dependent, autonomous, and mediators/moderators. It also proposes many relationships that can be hypothesized and tested statistically by future researchers. The study is useful for the community at large because it provides a lot of new information and understanding of the phenomenon under study.

Limitations of the Study and Directions for Future Studies
There are certain limitations of the study as well. First, it uses a qualitative approach to address the issue; therefore, future studies should use quantitative approaches to address the issue. Second, the data has been collected from Pakistan (a Muslim country); therefore, the results of the study should accordingly be generalized and interpreted. In this context, future studies may be conducted in different countries. Third, a limited list of determinants was investigated since the literature was not very rich in this regard; therefore, future studies should enrich the list and replicate the study. Fourth, the data has been collected from a limited number of experts (focus group considered to be representative), it is recommended that in the future the data should be collected from a statistical population in order to augment or enhance the results of the study.

Policy Recommendations
Keeping in view the results of the study ("nature of gathering" being the key determinant and most of the other determinants being agile, unbalanced, or unsettled) it is recommended for policymakers to devise gathering specific SOPs and carefully evaluate the effects of unbalanced determinants in order to achieve the desired results of precautionary measures of COVID-19 pandemic.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.