1 Introduction

Virology scientists believe that there are approximately 1.7 million undiscovered types of viruses that are hosted only by mammals and birds. These scientists have no precise information about their nature, how they are transmitted to humans and spread, or their accompanying symptoms (Carroll et al. 2018). Therefore, the probability of humans becoming infected with a previously unknown virus that could lead to a new pandemic is unknown. In December 2019, the world witnessed the emergence of a new virus that was found in the city of Wuhan, China. Since that time, the spread of coronavirus disease 2019 (COVID-19), which is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has increased exponentially. According to the John Hopkins University of Medicine Coronavirus Resource Center website, by the end of January 2021, there were approximately 108, 148,755 infected people and 2,267,768 deaths across 188 countries (JHU 2021). The common symptoms reported by infected persons include fever, shortness of breath and several respiratory disorders. Saliva droplets and physical contact are reported as the main modes of transmission. As a result of limitations in terms of epidemiological surveillance and partly, diagnostic capacity, confirmed cases are likely to be underestimated in most countries (Jayaweera et al. 2020).

Considering that the virus spreads very rapidly, it was deduced that at some point, healthcare systems around the world would collapse. COVID-19 forced governments to implement mitigation strategies from restricting non-essential activities, for example, by closing schools, shopping centres and amusement parks, to instituting severe curfews with strict travel and movement restrictions. While most lockdown protocols implemented by countries considered the need to maintain access to important services, differences in the extent to which other non-essential activities were closed or opened were found across countries. For instance, the Saudi Arabian government implemented a nationwide full curfew from April 6 until May 28 that included 24 h movement restrictions all over the country (Anil and Alagha 2020). Mexico implemented a mitigation strategy that comprised lockdowns that restricted all non-essential activities. Although the government encouraged the populace to voluntarily reduce its mobility and carry out only essential activities, compulsory restrictions were not imposed (Estévez-Soto 2020). Unlike the previously mentioned countries, Sweden did not impose strict lockdowns, and an important part of its economic sectors remained opened (Carbone and Montecucco 2020). It is clear that different levels of lockdowns yield different results. In this regard, there is plenty of literature that contrasts performance and outputs among different types of lockdowns, strategies and scenarios across countries (Dickens et al. 2020; Sardar et al. 2020; Shammi et al. 2020). Months after the emergence of SARS-CoV-2, the debate has moved to identify the best strategy for returning “new regularity” by minimizing risk and maximizing benefits. In this context, the main objective of this work is to propose a decision support system (DSS) that aids governments, policymakers, or decision makers for assessing the most relevant indicators while a reopening strategy is designed after COVID-19 lockdowns. By identifying a set of indicators that should be prioritized after Covid-19 lockdowns through replicable and evidence-based methodology, we make a significant contribution for achieving a deeper understanding of this pandemic. The main “take-away” from this work available to readers, practitioners and stakeholders are given by the flexibility and capability of our methodology, which can be adapted to different contexts and scenarios. The investigated indicators are grouped into four dimensions: economic growth, environmental preservation, well-being in society and individual health. Note that these indicators are utilized by worldwide organizations as International Labour Organization (ILO 2021), International Monetary Fund (IMF 2021), Food and Agricultural Organization (FAO 2021) for investigating, ranking, or comparing countries. In this form, indicators proposed by mentioned organizations represent a framework which facilitates the replicability of this methodology across countries.

The paper is structured on six sections. A literature review, which comprises definitions of indicator criteria and a revision of the existing work in the field, is provided in the next section. The proposed methodology is introduced in section three. The results are provided on section fourth. A discussion that contrasts our results with another research is presented in the fifth section. The last section presents the conclusions and future research directions.

2 Contextual framework (Problem description)

There is scientific evidence that lockdowns and curfews, which comprise strategies such as social distancing, working at home and school closures, are extremely useful to “flatten the curve” of the number of infected people and therefore play an important role in avoiding overdemand for health care services in a short period of time. Along 2020, mathematical models and data science projects proliferated around the world as support systems for local governments to make better informed decisions about lockdown objectives and scopes (Ibarra-Vega 2020; Sardar et al. 2020). When forecasting how a pandemic is evolving, the general approach in epidemiology is centred on three variables: infected, deaths and recovered, which are also are non-linear functions with parameters given by the mitigation policies. Comorbidities and sociodemographic factors can be included in the model as covariates (Medrek and Pastuszak 2020). As long as the number of infected people is lower than the healthcare capacity, the pandemic remains manageable. Considering this the main objective of lockdowns and curfews, parallel debates emerged in relation to the collateral effects from another perspectives: economic, environmental, public and individual health (Martin et al. 2020). Quarantines and lockdowns reduce pandemic intensity but make it last longer. Regarding collateral effects, studies have estimated that lockdowns shorten economic activity by at least 40% (Martin et al. 2020). The International Monetary Fund (IMF) expected that the global economy will contract by 4.7% by 2020 (IMF 2021). The International Labour Organization reported that unemployment reached its peak point by reaching levels of approximately 6.7% of the global force, which is approximately equal to 195 million full-time employees (ILO 2021). Similarly, COVID-19 lockdowns, mainly related to confinement and isolation, brought several mental health problems, such as stress, depression and suicide. Because COVID-19 is an emerging phenomenon, studies accurately reporting repercussions on individuals’ mental health are not fully available. In 2003, SARS-CoV-1 was associated with a 30% increase in suicide in those aged 65 and older. Approximately 50% of the influenza-recovered patients were diagnosed as stressed and anxious, and approximately 30% of the healthcare workforce at health services showed emotional distress (Lee et al. 2007). However, not all side effects of COVID-19 are negative or undesirable. Recent data released by NASA (National Aeronautics and Space Administration) and the ESA (European Space Agency) indicated that as a result of quarantines and mobility restrictions, air pollution in Wuhan, Italy and Spain was reduced by approximately 30% (Cicala et al. 2020). Another study demonstrated that CO2 levels showed important reductions (30.3–48.5%) in the metropolitan area of Rio de Janeiro Brazil as a result of the COVID-19 quarantine (Dantas et al. 2020).

While analysing the impact of lockdowns on human activities, it is important to follow a systematic approach which considers all outputs, either positive, negative or neutral. There is supportive evidence that systematic approaches yield to better-informed decision-making processes and more accurate results (Pileggi 2020). In this case, it is evident that outputs related to COVID-19 lockdowns might conflict or be negatively correlated. While the main purpose of quarantines is to minimize the number of deaths caused by the virus, this is achieved by incurring collateral costs. When the time for reopening economies and societies has arrived, decision makers (DMs) must select the most suitable indicators by analysing the negatively correlated variables or conflicting criteria. Consider the following illustrative example: let “A” and “B” represent the relevant indicators to be assessed to determine whether a lockdown should be lifted. From one perspective, “A” can be considered better than “B”, and from a different point of view, “A” is better than “B”. Conflictive scenarios occur when indicators can be grouped simultaneously in two or more opposite categories, named “benefits” or “costs”. In this context, multi-criteria decision-making (MCDM) is suitable for decision-making problems with conflicting, incomplete or contradictory information. Research that proposes different types of Decision Support Systems based on MCGDM methods in the context of COVID pandemic recently emerged. Ashraf and Abdullah (2020) propose a decision-making approach for investigating COVID pandemic under the spherical fuzzy environment. Their results include indicators for measuring the degree of global emergency and uncertainty due to the pandemic. An application of DEMATEL method with intuitionistic fuzzy sets for modelling relaxations in lockdowns, on Ocampo and Yamagishi (2020) is proposed. By the application of TOPIS, Shrestha et al (2021) investigate the consequences of COVID pandemic on globalization. Due to space limitations other works related to this field were not mentioned. Readers interested on more details are advised to get Garg and Kumar (2019).

Most MCDM problems require the participation of more than one DM in the process. Hence, many MCDM approaches are employed in a group decision-making (GDM) environment. Since the emergence of multi-criteria GDM (MCGDM) methods in the early 1970s, they have been applied in different disciplines and contexts. MCGDM methods have been demonstrated to be suitable for finding solutions to complicated problems with conflicting criteria (Park et al. 2008; Deshpande et al 2020; Rowley et al. 2020; Pileggi, 2020). During the past decade, methodologies such as the analytic hierarchy process (AHP), technique for order of preference by similarity to ideal solution (TOPSIS), evidential reasoning (ER) and Višekriterijumsko kompromisno rangiranje (VIKOR), among others, have been widely applied in different disciplines and contexts (Yang et al. 1994; Hsu and Hsu 2008; Park et al. 2008; Zhao et al. 2013; Büyüközkan and Göçer 2018; Toros and Gazibey 2018; Lei et al. 2020). AHP is a mathematical method that places decision objects into a simple and practical hierarchical structure (Saaty 1977). The hierarchical shape, which is characteristic of AHP, allows the decision objects to align with the overall goal. The hierarchy simplifies the comprehension of problems with multiple aspects that are often affected by complex relations and therefore enables DMs to assess relationships in a structured and systematic way (Abdullah and Najib 2016). By using crisp values, pairwise comparisons are performed, and ultimately, overall weights for each criterion are obtained. The first versions of AHP were based on crisp numbers, and later, its limitations for accurately describing real-world problems were noted (Abdullah and Najib 2016; Büyüközkan and Göçer 2018). Because real-world problems are mainly characterized by vagueness and uncertainty, intuitionistic fuzzy set (IFS) theory was incorporated into AHP. Then it was extended to the fuzzy environment by incorporating concepts from interval-value intuitionistic fuzzy (IVIF) set theory (Atanassov 1986). This method is popularly known as intuitionistic fuzzy AHP (IVIF-AHP) (Abdullah and Najib 2016). By differentiating membership, non-membership and hesitancy, the IVIF-AHP has proven to be a reliable method for representing DM preferences as intuitionistic fuzzy values rather than exact numbers (Merigó and Casanovas 2011).

According to Abdullah and Najib (2016), IVIF-AHP can use incomplete information to accurately standardize scores across DMs. To overcome this limitation, the fuzzy cross entropy of IVIF sets was introduced (Ye 2011). By considering the optimal weights of the indicators, DM preferences can be aggregated, and the product of each pair of interval-value intuitionistic fuzzy number (IVIFN) matrices can be obtained. Considering this framework, several authors proposed the utilization of IVIF-AHP, where the IVIF numbers were obtained from linguistic variables rather than natural numbers (Wang et al. 2015; Abdullah and Najib 2016). As illustrated in the next section, IVIF-AHP can be employed to formalize linguistic variables and therefore provides a more transparent decision-making process. Taking in to account the characteristics of the problem, DMs might either come from different areas of expertise or not. The entropy weight method is the most common method and is capable of determining the weights of the DMs in decision-making processes (Ye 2011).

Under this framework, it is relevant for policy makers, governments and DMs to utilize a DSS that assigns importance to indicators to design a reopening strategy once COVID-19 lockdowns have been lifted. The present work introduces a novel interval-value intuitionistic fuzzy sets-DSS (IVIFS-DSS) methodology based on MCGDM approaches to identify the most important indicators that should be considered when a reopening strategy is implemented. A total of 16 criteria, which are grouped into four “dimensions”, are evaluated with respect to their impact on a given reopening strategy (see Table 1).

Table 1 Definitions of the key variables included in this research

The reason for selecting these indicators and no others is based on a literature review carried out for these purposes. According to worldwide organizations such as the United Nations (UN), Organization for Economic Cooperation and Development (OECD), Bank for International Settlements (BIS) and Food and Agriculture Organization of the United Nations (FAO), among others, the cited indicators are widely utilized for ranking, comparing or investigating countries (FAO 2021b; OECD 2021a; OECD 2021b; UN 2021). Note that some of the mentioned organizations gather more than 100 countries. These indicators represent a solid framework which allows the replicability of our methodology across countries. It is beyond the scope of this word to provide a detailed operational definition of these indicators; however, summarized definitions are provided for these indicators in “Appendix 1”.

Our work aims to assist DMs and policy makers by identifying the best mix of indicators that should be considered when a reopening strategy is implemented. The novel contributions of this research are as follows:

  • This study identifies 16 indicators in four categories, the economy, environmental protection, societal well-being and individual health and ranks their importance for reopening strategies after COVID-19 lockdowns have been lifted.

  • This study employs a novel IVIFS-DSS methodology based on MCGDM approaches (entropy weighting and AHP) to help policy makers prioritize and select indicators that must be considered when a reopening strategy is carried out.

  • This study provides practical guidance that can help policymakers and planners make better decisions during the ongoing pandemic and future pandemics.

3 Materials and methods

This section gathers the theoretical foundations behind the proposed methodology. It is composed by three subsections. Formal definitions for IFS and IVIFS are firstly presented. They are followed by detailed explanations of linguistic variables and theirs respective IVIFN equivalences. At the end, the IVIF DSS methodology, which comprises eight steps, is presented. For sake of simplicity, note that section four follows the same eight-step sequence that is here presented at subsection three.

3.1 Definition of IVIFS

After ordinary fuzzy set theory was proposed by Zadeh (1965) to address the uncertainty of human judgement, an improved extension referred to as the IFS was introduced by Atanassov (1986). If \(\stackrel{\sim }{X}\) denotes an IFS bounded in domain Z, then it can be expressed in the following way.

$$\widetilde{ X} = \left\{ {\left\langle {z,\mu_{{\widetilde{{\tilde{X}}}}} \left( z \right), \nu_{{\tilde{X}}} \left( z \right)} \right\rangle /z \in Z} \right\}$$
(1)

In formula (1), \({\mu }_{\stackrel{\sim }{X}}\left(z\right)\) and \({\nu }_{\stackrel{\sim }{X}}\left(z\right)\) represent the membership degree and non-membership degree, respectively. Furthermore, \(z\in Z\) for any subset \(\stackrel{\sim }{X}\) of \(Z\), and \({\mu }_{\stackrel{\sim }{X}}\left(z\right)\) and \({\nu }_{\stackrel{\sim }{X}}\left(z\right)\) satisfy \(0\le {\mu }_{\stackrel{\sim }{X}}\left(z\right)+{\nu }_{\stackrel{\sim }{X}}\left(z\right)\le 1\). \({\pi }_{\stackrel{\sim }{X}}\left(z\right)\) represents the hesitancy degree, and therefore, \({\pi }_{\stackrel{\sim }{X}}\left(z\right)=1-{\mu }_{\stackrel{\sim }{X}}\left(z\right)-{\nu }_{\stackrel{\sim }{X}}\left(z\right)\). In this formula, any element \(z\in Z\) is assigned to three different categories, as illustrated in Fig. 1.

Fig. 1
figure 1

Membership, non-membership, and hesitancy relations

In some cases, it is not easy to exactly define the membership degrees for certain elements of the IFSs. Thus, Atanassov and Gargov (1989) proposed a generalized type of IFS, which is called IVIFS in the literature. Let A be an IVIFS, where \(z\in Z\); then, it can be calculated as follows:

$$A = \left\{ {z,\mu_{A} \left( z \right), \nu_{A} \left( z \right)/z \in Z} \right\}$$
(2)

where \({\mu }_{A}\left(z\right)\) and \({\nu }_{A}\left(z\right)\) take values between [0,1], under the condition that \(0\le \mathrm{sup}{\mu }_{A}\left(z\right)+\mathrm{sup}{\nu }_{A}\left(z\right)\le 1\) for any \(z \epsilon Z\). The closed intervals of \({\mu }_{A}\left(z\right)\) and \({\nu }_{A}\left(z\right)\) are denoted by the lower and upper bounds as \({\mu }_{A}^{L}\left(z\right)\), \({\mu }_{A}^{U}\left(z\right),{\nu }_{A}^{L}\left(z\right)\), and \({\nu }_{A}^{U}\left(z\right),\), respectively. Note that the intervals of \({\mu }_{A}\left(z\right)\) and \({\nu }_{A}\left(z\right)\) denote the degree of membership and non-membership, respectively, for each \(z \epsilon Z\) in the IVIFS \(A\), as shown below:

$$A = \left\{ {z,\left[ {\mu_{A}^{L} \left( z \right), \mu_{A}^{U} \left( z \right)} \right], \left[ {\nu_{A}^{L} \left( z \right), \nu_{A}^{U} \left( z \right)} \right]/z \in Z} \right\}$$
(3)

where \({0\le \mu }_{A}^{L}\left(z\right)\le {\mu }_{A}^{U}\left(z\right)\le 1\) and \({0\le \nu }_{A}^{L}\left(z\right)\le {\nu }_{A}^{U}\left(z\right)\le 1\). Moreover, for each IVIFS \(A\), where \(z\in Z\), there exists an IVIF index \({\pi }_{A}\left(z\right)\) that denotes the range of the hesitancy membership degree of the element \(z \epsilon Z\) in \(A\), as proposed by Park et al. (2008) and can be calculated as follows.

$$\pi_{A} \left( z \right) = \left[ {\pi_{A}^{L} \left( z \right), \pi_{A}^{U} \left( z \right)} \right] = \left[ {1 - \mu_{A}^{U} \left( z \right) - \nu_{A}^{U} \left( z \right),1 - \mu_{A}^{L} \left( z \right) - \nu_{A}^{L} \left( z \right)} \right]$$
(4)

For convenience, an IVIFN is denoted by α. Its formal definition was proposed by Xu and Yager (2006) as \(\alpha =\left(\left[{\mu }_{\alpha }^{L}, {\mu }_{\alpha }^{U}\right], \left[{\nu }_{\alpha }^{L}, {\nu }_{\alpha }^{U}\right], \left[{\pi }_{\alpha }^{L}, {\pi }_{\alpha }^{U}\right]\right)\). It is beyond the scope of this work to provide a detailed explanation of the rules, laws and operations of IVIFN. Interested readers are encouraged to consult Park et al. (2008), Büyüközkan and Göçer (2018), Garg and Kumar (2019)

3.2 Linguistic variables

One advantage of IVIFS with respect to other techniques is the possibility of elaborating reasoning assessments, which are given in terms of linguistic variables (e.g., “very good”, “good” and “medium good”) (Wang et al. 2015). Numerous studies have documented that linguistic variables are more suitable for describing real-world MCDM problems that involve vague, imprecise or uncertain information (Garg and Kumar 2019). The linguistic terms used in this work and their corresponding IVIFNs are provided in Table 2. Linguistic variables are efficient and suitable tools for denoting different activities in the real-world that can be assessed by the qualitative judgements of DMs rather than quantitative judgements (Zadeh 1975; Herrera et al. 1997). Linguistic variables are more appropriate for describing real-world MCDM problems involving imprecise and uncertain information (Garg and Kumar 2019). These linguistic variables can be expressed by IVIFNs. Many scholars have concluded that the process of collecting DMs’ judgements for MCGDM problems with a high degree of uncertainty via linguistic contexts leads to more reliable results (Ben-Arieh and Chen 2006; Rodríguez et al. 2013).

Table 2 Linguistic terms used to evaluate the importance of the DMs

Note that the IVIFS introduced in Table 2 are utilized here to assess the importance of each DM. The values were proposed by Büyüközkan and Göçer (2018). The linguistic terms’ IVIFNs and reciprocal IVIFNs are presented in Table 3. The relations shown in Table 3 are useful for calculating the relative importance and associated weight for each dimension and its indicators. The methodology proposed by Abdullah and Najib (2016) is replicated here.

Table 3 Linguistic terms used for the importance weighting

Based on the linguistic preferences in the IVIFN and the reciprocal information presented in Table 3, the proposed IVIF-AHP is shaped to prioritize the dimensions and their indicators based on their importance for COVID post-lockdown strategies. Next, the DSS methodology is introduced.

3.3 IVIF DSS methodology

Here, we introduce the DSS methodology, which includes eight steps and is based on the IVIF-AHP framework, as illustrated in Fig. 2. These steps are similar to those used for the traditional AHP. However, among the most important differences are the utilization of linguistic variables for collecting DM information and the application of an interval-value measurement scale. The concept of entropy is utilized for assigning weights to each DM according to their respective backgrounds and expertise. The scale is revised by adding interval-value hesitation degrees. These concepts are explained in detail below.

Fig. 2
figure 2

Flowchart of the methodology

3.3.1 Defining the goal, dimensions and indicators

The dimensions and related indicators are defined based on the operational definitions of the variables. Note that the set of indicators aligns with the main objective of the research. A general definition of the problem is provided below:

  • Let \({C}_{i}\) be a finite set of \(n\) dimensions with a weight vector\([{W}_{1}, {W}_{2}, . . . ,{W}_{n}]\), where \({W}_{i}>0, i=\mathrm{1,2},\dots n, {\sum }_{i=1}^{n}{W}_{i}=1\)

  • For each \(i\)-th dimension, let \({C}_{ij}\) be a finite set of \(m\) indicators associated with a local weight vector\([{w}_{1}, {w}_{2}, . . . ,{w}_{m}]\), where\({w}_{j}>0, j=\mathrm{1,2},\dots m, {\sum }_{j=1}^{m}{w}_{j}=1\). Here, \({C}_{ij}\) denotes the \(j\) th indicator of the main \(i\) th dimension.

  • Let \(E = \{{e}_{1}, {e}_{2}, . . . , {e}_{K}\}\) be the same set of \(k\) DMs and \(E\) be associated with a weight vector denoted as\([{\lambda }_{1}, {\lambda }_{2}, . . . ,{\lambda }_{k}]\), where\({\lambda }_{k}>0, k=\mathrm{1,2},\dots K, {\sum }_{k=1}^{K}{\lambda }_{k}=1\).

3.3.2 Determining the group of DMs

A GDM is defined as two or more individuals interacting to choose from a set of alternatives. Decisions made by groups are believed to be more robust since they are not attributable to a single individual but rather to the group (Merigó and Casanovas 2011). For real-life situations and problems, evidence shows that collaborative decision-making is one of the most important sources of creativity and accuracy (Ye 2011). Evidence also shows that groups with more than 6 individuals complicate the process (Abdullah and Najib 2016). Therefore, five experts participated in this study.

3.3.3 Computing the weights of the DMs

To develop a measure of the importance of each DM in the group in terms of linguistic expressions, three main criteria were considered. (1) Experience refers to the number of years that each DM has been in the professional field since the date of his/her first job. (2) Knowledge refers to the area in which the DM is specialized, and (3) responsibilities indicates whether the DM supervises other people or staff.

\({e}_{k} (k=\mathrm{1,2},\dots K)\) represents the participants. Then,\(d^{\left( k \right)} = \left( {\left[ {\left( {\mu_{d}^{L} } \right)^{\left( k \right)} , \left( {\mu_{d}^{U} } \right)^{\left( k \right)} } \right], \left[ {\left( {\nu_{d}^{L} } \right)^{\left( k \right)} , \left( {\nu_{d}^{U} } \right)^{\left( k \right)} } \right], \left[ {\left( {\pi_{d}^{L} } \right)^{\left( k \right)} , \left( {\pi_{d}^{U} } \right)^{\left( k \right)} } \right]} \right)\). This measure also includes the uncertainty preference rate for the importance option provided by each DM. Note that the equivalences provided in Table 2 are applied here to obtain numerical values for each\({d}^{\left(k\right)}\), where k = 1,…,K represents the number of DMs participating in the study. Next, a fuzzy entropy measure is applied for each IVIFN to estimate the DM weight vector (denoted as\([{\lambda }_{1}, {\lambda }_{2}, . . . ,{\lambda }_{K}]\)) (Ye 2010). Hence, the fuzzy entropy measure of\({d}^{\left(k\right)}\), denoted as\(EN\left({d}^{\left(k\right)}\right)\), is calculated as follows:

$$EN\left( {d^{\left( k \right)} } \right) = \left[ {\cos \left( {\frac{{\pi \left[ {1 + \left( {\mu_{d}^{L} } \right)^{\left( k \right)} + p\left( {\left( {\mu_{d}^{U} } \right)^{\left( k \right)} - \left( {\mu_{d}^{L} } \right)^{\left( k \right)} } \right) - \left( {\nu_{d}^{L} } \right)^{\left( k \right)} - q\left( {\left( {\nu_{d}^{U} } \right)^{\left( k \right)} - \left( {\nu_{d}^{L} } \right)^{\left( k \right)} } \right)} \right]}}{4}} \right) + \cos \left( {\frac{{\pi \left[ {1 - \left( {\mu_{d}^{L} } \right)^{\left( k \right)} - p\left( {\left( {\mu_{d}^{U} } \right)^{\left( k \right)} - \left( {\mu_{d}^{L} } \right)^{\left( k \right)} } \right) + \left( {\nu_{d}^{L} } \right)^{\left( k \right)} + q\left( {\left( {\nu_{d}^{U} } \right)^{\left( k \right)} - \left( {\nu_{d}^{L} } \right)^{\left( k \right)} } \right)} \right]}}{4}} \right) - 1} \right] \times \frac{1}{\sqrt 2 - 1}$$
(5)

where \(0\le EN\left({d}^{\left(k\right)}\right)\le 1\) and \(p, q\epsilon [\mathrm{0,1}]\) are two fixed numbers. Then, the entropy weights of the \(k\)-th DM are determined as follows:

$$\lambda_{k} = \frac{{1 - EN\left( {d^{\left( k \right)} } \right)}}{{K - \mathop \sum \nolimits_{k = 1}^{K} EN\left( {d^{\left( k \right)} } \right)}}$$
(6)

where \({\lambda }_{k}\epsilon [\mathrm{0,1}], k=\mathrm{1,2},\dots K, {\sum }_{k=1}^{K}{\lambda }_{k}=1\).

3.3.4 Obtaining the individual pairwise comparison matrices

Each DM is asked to provide his/her judgement regarding each dimension and the related indicators through individual pairwise comparison evaluation matrices. Linguistic terms, which are shown in Table 3, are used to evaluate the dimension and indicators. The collected data were transformed into IVIFNs, as shown in Table 3. Next, the corresponding individual pairwise comparison matrix for the \(k\)th DM used to evaluate the dimensions can be constructed as follows:

$$T_{iy}^{\left( k \right)} = \left[ {t_{iy}^{\left( k \right)} } \right]_{n \times n} = \left[ {\begin{array}{*{20}c} {\begin{array}{*{20}c} - \\ {t_{21} } \\ \end{array} } & {\begin{array}{*{20}c} {t_{12} } \\ - \\ \end{array} } & {\begin{array}{*{20}c} {\begin{array}{*{20}c} \ldots & {t_{1n} } \\ \end{array} } \\ {\begin{array}{*{20}c} \ldots & {t_{2n} } \\ \end{array} } \\ \end{array} } \\ \vdots & \vdots & {\begin{array}{*{20}c} { \vdots } & \vdots \\ \end{array} } \\ {t_{n1} } & {t_{n2} } & {\begin{array}{*{20}c} \ldots & - \\ \end{array} } \\ \end{array} } \right]^{\left( k \right)}$$
(7)

where \(t_{iy}^{\left( k \right)} = \left( {\left[ {\mu_{{t_{iy} }}^{L} , \mu_{{t_{iy} }}^{U} } \right], \left[ {\nu_{{t_{iy} }}^{L} , \nu_{{t_{iy} }}^{U} } \right], \left[ {\pi_{{t_{iy} }}^{L} , \pi_{{t_{iy} }}^{U} } \right]} \right)^{\left( k \right)}\)\(, i,y=\mathrm{1,2},\dots ,n;i\ne y\), and \({{t}_{iy}}^{\left(k\right)}\) denotes the evaluation importance degree of dimension \(i\) over dimension \(y\) according to DM\(k\). In the same way, the individual pairwise comparison matrices for the \(j\) th indicator with respect to the main \(i\) th dimension can be obtained.

3.3.5 Calculating the aggregated IVIFNs and obtaining the group evaluation IVIFNs

Interval-valued intuitionistic fuzzy weighted averaging (IIFWA) is applied to aggregate the evaluation matrix \({{T}_{iy}}^{\left(k\right)}= {\left[{{t}_{iy}}^{\left(k\right)}\right]}_{n\times n}\) with \({T}^{k}= {\left[{{t}_{i}}^{\left(k\right)}\right]}_{n\times 1}\) for k = 1,…,K. The IIFWA operator inputs the vector of the weights of the DMs, as illustrated in formula (8).

$${{t}_{i}}^{\left(k\right)}=\mathrm{IVIFWA}\left({{t}_{i1}}^{\left(k\right)}, {{t}_{i2}}^{\left(k\right)},., {{t}_{in}}^{\left(k\right)}\right)=$$
$$\left( {\left[ {1 - \mathop \prod \limits_{y = 1}^{n} \left( {1 - \left( {\mu_{{t_{iy} }}^{L} } \right)^{\left( k \right)} } \right)^{{\lambda_{k} }} ,1 - \mathop \prod \limits_{y = 1}^{n} \left( {1 - \left( {\mu_{{t_{iy} }}^{U} } \right)^{\left( k \right)} } \right)^{{\lambda_{k} }} ,} \right], \left[ {\mathop \prod \limits_{y = 1}^{n} \left( {1 - \left( {\nu_{{t_{iy} }}^{L} } \right)^{\left( k \right)} } \right)^{{\lambda_{k} }} ,\mathop \prod \limits_{y = 1}^{n} \left( {1 - \left( {\nu_{{t_{iy} }}^{U} } \right)^{\left( k \right)} } \right)^{{\lambda_{k} }} } \right]} \right)$$
(8)

where \({{t}_{i}}^{\left(k\right)}\) represents the importance degree for the i-th dimension, where \(i=1,\dots ,n\).

Additionally, for the kth DM, the IVIFN is given by \(t_{i}^{\left( k \right)} = \left( {\left[ {\mu_{{t_{i} }}^{L} , \mu_{{t_{i} }}^{U} } \right], \left[ {\nu_{{t_{i} }}^{L} , \nu_{{t_{j} }}^{U} } \right], \left[ {\pi_{{t_{i} }}^{L} , \pi_{{t_{i} }}^{U} } \right]} \right)^{\left( k \right)}\). The IVIF values for the \(j\)th indicators with respect to their respective i-th dimension are computed in similar form. Finally, the IVIF values, which are based on the DMs’ evaluations, are aggregated into a unique main dimension by utilizing the IVIFWA operator, as shown in the following expression (Xu and Cai 2012):

$$t_{i} = {\text{IVIFWA}}\left( {t_{i}^{\left( 1 \right)} , t_{i}^{\left( 2 \right)} ,., t_{i}^{\left( K \right)} } \right) = \left( {\left[ {1 - \mathop \prod \limits_{k = 1}^{K} \left( {1 - \left( {\mu_{{t_{i} }}^{L} } \right)^{\left( k \right)} } \right)^{{\lambda_{k} }} ,1 - \mathop \prod \limits_{k = 1}^{K} \left( {1 - \left( {\mu_{{t_{i} }}^{U} } \right)^{\left( k \right)} } \right)^{{\lambda_{k} }} ,} \right], \left[ {\mathop \prod \limits_{k = 1}^{K} \left( {1 - \left( {\nu_{{t_{i} }}^{L} } \right)^{\left( k \right)} } \right)^{{\lambda_{k} }} ,\mathop \prod \limits_{k = 1}^{K} \left( {1 - \left( {\nu_{{t_{i} }}^{U} } \right)^{\left( k \right)} } \right)^{{\lambda_{k} }} } \right]} \right)$$
(9)

where \({\lambda }_{k}\epsilon [\mathrm{0,1}], {\sum }_{k=1}^{K}{\lambda }_{k}=1\) and \({t}_{i}\) indicates a group decision IVIF value for the i-th dimension, where i = 1,…,n and is represented as \(t_{i} = \left( {\left[ {\mu_{{t_{i} }}^{L} , \mu_{{t_{i} }}^{U} } \right], \left[ {\nu_{{t_{i} }}^{L} , \nu_{{t_{i} }}^{U} } \right], \left[ {\pi_{{t_{i} }}^{L} , \pi_{{t_{i} }}^{U} } \right]} \right)\). The same procedure is applied to obtain group decision IVIFNs for the \(j\) th indicator with respect to the i-th dimension.

3.3.6 Obtaining consistency ratio (CR)

The CR is useful for evaluating the extent to which each DM assessment is reliable and comprises acceptable measures of variation. Using the dimensions and their indicators as inputs, the CR is calculated. The random index (RI), for which the respective values are provided in Table 4, is applied here to compute the CR. The formula proposed by Saaty (1990) is utilized to calculate the CR, as shown below.

$$CR = \frac{{RI - \frac{{\sum \pi_{{t_{i} }}^{U} }}{n}}}{n - 1}$$
(10)

where n represents the number of dimensions and \({\pi }_{{t}_{i}}^{U}\) indicates the upper hesitancy value. The CRs for all indicators included in this analysis are calculated in similar form. According to Saaty (1990), a CR is considered acceptable if it has a value of 0.10 or lower (CR ≤ 0.10). If the CR is higher than 0.10 (CR > 0.10), the DM’s assessment is considered inconsistent and therefore should be either revised or dropped.

Table 4 Random indexes

3.3.7 Calculating the weights of the dimensions and indicators

The weights of the dimensions and their indicators are objectively determined by using the IVIF values as inputs (see Table 3). The dimension weights, denoted by \({W}_{i}\), where \({W}_{i}\epsilon \left[\mathrm{0,1}\right]; i=\mathrm{1,2},\dots n; {\sum }_{i=1}^{n}{W}_{i}=1\), are obtained using the equations proposed by Büyüközkan et al. (2020), as illustrated in (11) and (12).

$$W_{i} = \frac{{1 - H\left( {c_{i} } \right)}}{{n - \mathop \sum \nolimits_{i = 1}^{n} H\left( {c_{i} } \right)}}$$
(11)

where

$$H\left( {c_{i} } \right) = 1 - \frac{{\mathop \sum \nolimits_{k = 1}^{K} \frac{{\lambda_{k} \left[ {\mu_{{t_{i} }}^{L} + \mu_{{t_{i} }}^{U} } \right]}}{2}}}{{\sqrt {\mathop \sum \nolimits_{k = 1}^{K} \frac{{\lambda_{k} \left[ {\left( {\mu_{{t_{i} }}^{L} } \right)^{2} + \left( {\mu_{{t_{i} }}^{U} } \right)^{2} + \left( {\nu_{{t_{i} }}^{L} } \right)^{2} + \left( {\nu_{{t_{i} }}^{U} } \right)^{2} } \right]}}{2}} }}$$
(12)

For each \(i\)th dimension (\(i=\mathrm{1,2},\dots n\)), the local weights for indicators \({C}_{ij}\), denoted by \({w}_{j}\), where (\({w}_{j}\epsilon \left[\mathrm{0,1}\right];\, j=\mathrm{1,2},\dots m; {\sum }_{j=1}^{n}{w}_{j}=1)\), are obtained by recursively utilizing formulas (11) and (12) for each dimension.

3.3.8 Ranking the dimensions and indicators

This last step comprises three types of results. First, the overall relative importance of each dimension (denoted by \({W}_{i}\)) is calculated. The local weight for each indicator allocated to the i-th dimension is calculated (identified as \({w}_{j}\)). The global weights that represent the relative importance of all indicators with respect to the goal are then calculated.

The abovementioned IVIF-AHP framework is used for an original application of linguistic scales and three types of IFS, namely, membership, non-membership and hesitation, to rank the relevant COVID-19 post-lockdown indicators. By using this methodology, the comparison matrix judgements typical of AHP problems yield more accurate results, as shown in the following section.

4 Results

It is relevant for governments and DMs to evaluate the importance of the set of indicators introduced in Table 1 when designing a post-COVID-19 lockdown reopening strategy. Based on data collected from a group of DMs, a practical case study is presented in this section. To ensure the explanation is “easy-to-follow”, this illustrative case study follows the same eight-step methodology presented in the last section.

4.1 Defining the goal, dimensions and indicators

The structure of the IVIF-AHP problem consists of four main dimensions: economic growth, environmental protection, societal well-being and individual well-being. Similarly, there are 16 indicators that are hierarchically allocated to these dimensions. The overall framework is presented in Fig. 3. The operational definitions of the indicators, which includes a detailed explanation of each concept based on a review of the literature, is provided in Appendix 1. The next step consists of determining the group of experts that will participate in the study.

Fig. 3
figure 3

The hierarchical model of the indicators that should be considered for a post-COVID-19 lockdown reopening strategy

4.2 Determining the group of DMs

A committee consisting of five experts,\({e}_{k} (k=\mathrm{1,2},\dots 5)\) who were knowledgeable in areas such as the decision sciences, business, engineering and computer sciences was invited to participate. The group included professionals from both academia and industry. All participants were highly experienced in assessing indicators and working with multi-criteria decision analysis problems. The participants were asked to provide their assessment by responding to an electronic questionnaire that was designed for these purposes and featured the typical hierarchical structure of AHP problems. The academic background and profile of each participant is presented in Table 5.

Table 5 Academic backgrounds and profiles of the DMs

4.3 Computing the weights of the DMs

Considering the linguistic terms and their respective IVIFNs (see Table 1), the weights of the DMs \({e}_{k} (k=\mathrm{1,2},\dots 5)\) were calculated. First, the DMs were asked to evaluate their knowledge of the indicators that might affect reopening strategies following COVID-19 lockdowns. Later, the scale presented in Table 2 was used for these purposes and adapted to our electronic questionnaire. The weight of each DM in terms of the linguistic terms and IVIFN are presented in Table 6. Note that formulas (5) and (6) were applied to calculate these weights.

Table 6 Importance ratings of the DMs in terms of the linguistics terms and related weights

4.4 Obtaining the individual pairwise comparison matrices

By using the linguistic terms shown in Table 3, DMs assessed the importance of each dimension with respect to the others and then the relative importance of indicators. The responses provided by the DMs regarding each comparison is provided in Appendix 2. The linguistic terms reported by the DMs were later transformed to IVIFNs by following the equivalences shown in Table 3. Note that formula (7) was utilized to generate the pairwise comparison matrices. A total of five matrices were computed (see Appendix 2). While the first matrix referred to the dimensions, the remaining four matrices related the same number of indicator groups.

4.5 Calculating the aggregated IVIFNs and obtaining the group evaluation IVIFNs

By employing the IVIFNs, the equivalences, dimensions and indicators are aggregated with DMs’ respective weights through the application of formula (8). The output of this process in terms of the dimensions is provided in Table 7. Table 7 provides only the aggregated IVIFNs of the main dimensions. Later, the same calculation is replicated to aggregate the indicator evaluations into their respective dimensions by considering the IVIFN equivalences.

Table 7 Aggregated IVIFNs for the main dimensions

The IVIFWA operator, previously introduced in formula (9), is used here. Table 8 illustrates the outputs regarding the dimensions and four groups of indicators.

Table 8 Group evaluation of the IVIFNs for the main dimensions and the corresponding indicators

4.6 Obtaining the CR

The CRs for the dimensions group evaluation IVIFNs were calculated. By using formula (10), the CR for the dimensions can be obtained as follows.

$$CR = \frac{{0.90 - \frac{{\left( {0.73 + 0.68 + 0.65 + 0.66} \right)}}{4}}}{3} = 0.074$$

Similarly, the CRs were obtained for all IVIFNs; the CRs were 0.067, 0.095, 0.092 and 0.074 for C1, C2, C3 and C4, respectively. Note that all the CR values were below 0.10; therefore, all aggregation matrices are consistent.

4.7 Calculating the weights for the dimensions and associated indicators

By applying formulas (11) and (12), the weights for the dimensions and indicators are obtained. Local weights, which refer to the relative importance to indicators within respect to its groups, are calculated in similar form. Based on the above, three types of weights are provided. The first weight is related to dimension importance, followed by local and global indicator weights. The weight values of all dimensions and the local weight values of the indicators included in the study are shown in Table 9. The global weights are provided in Table 10.

Table 9 Weights of the main dimensions and associated indicators
Table 10 Global weights of the indicators

4.8 Ranking the dimensions and indicators

The ranking of the dimensions, based on their weights, is carried out. “Societal well-being” (C3) is the most important dimension, followed by “Economic growth” (C1), “Individual well-being” (C4) and “Environmental protection” (C2).

In relation to the global weights, the five most important indicators are “Reduction in crime and violence” (C32), “Reduction in unemployment” (C32), “Retail sales (C13), “Freshwater source preservation” (C23) and “Nutrition” (C23). Note that “Societal well-being” and “Environmental protection” are the two dimensions with four indicators. The complete ranking of all indicators with respect to the global weights is provided in Table 10.

5 Discussion

A DSS is a tool for collecting, organizing and analysing information to solve decision-making problems. A DSS can be used for either single or GDM by ensuring the results are more robust when dealing with complex or unstructured scenarios. In this paper, a DSS is proposed to help stakeholders by prioritizing indicators for designing and implementing a post-COVID-19 strategy. Considering that COVID-19 pandemic and related lockdowns are complex scenarios, mainly featuring vagueness and uncertainty, the proposed tool is based on the IFS AHP (Abdullah and Najib 2016). By distinguishing among membership, non-membership and hesitancy, our DSS is capable of more accurately calculating the DM’s assessments. Four dimensions are considered to shape the main framework of the problem. International organizations, such as the United Nations, World Bank and International Labour Organization, among others, widely utilize these dimensions-indicators for ranking, comparing or investigating countries (FAO 2021b; ILO 2021; IMF 2021; OECD 2021a; UN 2021; WHO 2021). Each dimension comprises four indicators. Our DSS was utilized by DMs from industry and academia. The assessment made by the group of DMs yielded three types of results: dimension importance and the prioritization of the global and local indicators.

“Societal well-being” was the most important dimension, followed by “Economic growth”, “Individual well-being” and “Environmental protection”. Addressing poverty by guaranteeing access to material resources to help individuals achieve acceptable levels of quality of life and ensuring that worthy sources of employment are available in a clean and safe environment, should be key aspects of post-COVID-19 strategies. According to Bick and Blandin (2020), the unemployment rate in the U.S. increased to 20% by the beginning of May 2020. Martin et al. (2020) state that unemployment is among the most devastating consequences of policies that restrict movement and lockdowns, since unemployed individuals are mainly focused on meeting their primary needs (e.g., food, housing and clothing), and therefore, the risk of being infected is secondary. Successful COVID-19 prevention policies are inversely related to unemployment. Ensuring socioeconomic stability and employment confidence after lockdowns must be a priority for governments when post-COVID lockdown policies are designed and implemented (Tran et al. 2020).

Research provides evidence of the negative effects of lockdowns on individuals’ mental health. The negative effects are manifested as anxiety, depression or eating disorders that eventually lead to several types of violence (Piquero et al. 2020). On the other hand, crime and violence events are products of these circumstances. The likelihood that crime or violence occurs is given by the convergence of the perpetrator and victim in the absence of surveillance (Estévez-Soto 2020). Therefore, lockdowns lead to a decrease in robberies in public places but an increase in several types of domestic violence, since the perpetrators and victims spend more time at home. A study reported that from March to May 2020, there was a substantial decrease in robberies and residential burglaries but a significant increase in domestic violence in Los Angeles and Indianapolis (Mohler et al. 2020). Another study in Mexico City reported a significant increase in the number of calls received by the Violence Against Women (VAW) help line due to the lockdown. This increase is positively correlated with the number of domestic violence events indicated by official statistics (Estévez-Soto 2020).

As reported in the literature, crime, violence and unemployment are among the variables most negatively affected by lockdowns (Mohler et al. 2020). Our proposed DSS is consistent with this evidence, in the sense that governments and policy makers should consider these variables as priorities when the lockdowns are lifted. Note that the dimension “Environmental protection” appears last in terms of importance. This does not mean that protecting the environment is not relevant or unimportant, quite the opposite. Based on a general framework composed of conflicting and negatively correlated indicators, the DSS helps DMs focus their attention on the indicators that should be prioritized when the lockdowns are lifted. By differentiating between immediate priorities (Societal well-being) and long-term priorities (Environmental protection), the system helps policy makers by strategically focusing their attention on priorities, given the complex and unstructured situation.

6 Conclusions

By the end of January 2021, the number of positive COVID cases worldwide was around 108, 148,755 and 2,267,768 deaths (JHU 2021). Humankind is facing the one of the biggest economic, social and public health crises. Besides the irreparable deaths and the related suffering to relatives, the pandemic is bringing devastating effects on income and employment of millions of families around the world, and therefore deteriorating the well-being on individuals and societies. As commonly occurs on this type of phenomena, the poorest part of the population is the most affected since lockdowns and curfews lead to higher levels of inequality, poverty and setbacks on nutrition, well-being and economic development. Since the emergence of the virus until January 2021, several countries passed from believing the pandemic was under control to face new lockdowns, which forced governments and policy makers to rethink strategies and actions. Although, plenty of scientific evidence related to key aspects of the pandemic emerged on the past few months, at this moment it is not possible to carry out an accurate assessment of how effective governments' strategies and actions have been, since the pandemic is not over.

Under this perspective, this study makes valuable contribution by making available to governments, policy makers and stakeholders a flexible and easily adaptable methodology, which can be replicated to assess future stages of the pandemic. Since the proposed framework is based on indicators introduced by organizations as World Health Organization (WHO), International Monetary Found (IMF) and Food and Agriculture Organization (FAO), stakeholders can find here a scientific-based approach for conducting further research related with the pandemic evolution.

In this work, we do not attempt to set the foundations for designing post-COVID lockdown strategies. Our goal is much more modest. Considering that the COVID pandemic is a phenomenon with multiple variables, we contribute to the discussion of how this complex problem should be approached and highlight additional elements for analyses. Given that the DMs participating in this study are not directly involved in public policy making, our results are more illustrative than prescriptive. Considering this framework is based on an approach that takes generic indicators proposed by international organizations, country level implementations should be carefully studied, since each country has peculiarities and not all proposed indicators might be suitable for all countries. Thus, we attempt to put at the hand of governments, policymakers and stakeholders, consistent and evidence-based tools to evaluate difficult choices and avoid further deterioration in the economy, employment and wellbeing for societies and individuals. Additionally, the present work illustrates the suitability of decision science tools for tackling complex and unstructured problems, such as the COVID pandemic, and consequently make an original contribution in that direction. Additional research involving expert public policy makers in the DSS process is needed to produce helpful results for society. Further research related to the present pandemic or other pandemics that might arise in the future can find in this work a set of general guidelines through which decision-making processes can be conducted within a clear and structured framework.