key: cord-1029938-u79q19zw authors: Babaee Tirkolaee, Erfan; Goli, Alireza; Ghasemi, Peiman; Goodarzian, Fariba title: Designing a sustainable closed-loop supply chain network of face masks during the COVID-19 pandemic: Pareto-based algorithms date: 2021-12-15 journal: J Clean Prod DOI: 10.1016/j.jclepro.2021.130056 sha: 057523d165bce464494f167e8926ad7b0168c18f doc_id: 1029938 cord_uid: u79q19zw This study develops a novel mathematical model to design a sustainable mask Closed-Loop Supply Chain Network (CLSCN) during the COVID-19 outbreak for the first time. A multi-objective Mixed-Integer Linear Programming (MILP) model is proposed to address the locational, supply, production, distribution, collection, quarantine, recycling, reuse, and disposal decisions within a multi-period multi-echelon multi-product supply chain. Additionally, sustainable development is studied in terms of minimizing the total cost, total pollution and total human risk at the same time. Since the CLSCN design is an NP-hard problem, Multi-Objective Grey Wolf Optimization (MOGWO) algorithm and Non-Dominated Sorting Genetic Algorithm II (NSGA-II) are implemented to solve the proposed model and to find Pareto optimal solutions. Since Meta-heuristic algorithms are sensitive to their input parameters, the Taguchi design method is applied to tune and control the parameters. Then, a comparison is performed using four assessment metrics including Max-Spread, Spread of Non-Dominance Solution (SNS), Number of Pareto Solutions (NPS), and Mean Ideal Distance (MID). Additionally, a statistical test is employed to evaluate the quality of the obtained Pareto frontier by the presented algorithms. The obtained results reveal that the MOGWO algorithm is more reliable to tackle the problem such that it is about 25% superior to NSGA-II in terms of the dispersion of Pareto solutions and about 2% superior in terms of the solution quality. To validate the proposed mathematical model and testing its applicability, a real case study in Tehran/Iran is investigated as well as a set of sensitivity analyses on important parameters. Finally, the practical implications are discussed and useful managerial insights are given. Recent Coronavirus disease is an infectious disease that was first found in December 2019 in Wuhan, China, while mortality rates from the disease have been reported to be up to 5%, often transmitted by tiny respiratory droplets from one person to another (Mardani et al., 2020a; . Statistics show that as of 10.25.2020, the number of infected people with this deadly virus in Iran, the USA, China, Spain, and Germany is 526490, 8296249, 85659, 982723, and 357538, respectively. The number of deaths in these countries is 30123, 2223730, 4634, 33775, and 9838, respectively (World Health Organization, 2020) . The use of masks (N95, KN95, and surgical masks) is one of the preventive measures that can limit the prevalence of some respiratory viral diseases containing the COVID-19 (Nilashi et al., 2020; Kairon and Bhattacharyya, 2020; . Most people in the world use face masks at the same time as the coronavirus outbreak because masks can prevent the transmission of large droplets that are released into the air after sneezing and coughing, so face masks are useful (Cheng et al., 2020; Omori et al., 2020; Ullah et al., 2020; Wang et al., 2020) . In addition, the mask prevents the person from touching their nose and mouth with their hands, thus reducing their chances of contracting the virus distribution, and disposal centers are established and the amount of stored masks in each center directly affects the costs of the supply chain and, more importantly, the health of the individual (Mollenkopf et al., 2020) . For example, the production of the Shanghai-GM-Wuling (SGMW) car stopped production and resumed mask production in the COVID-19 situations because of decreasing vehicle demand and increasing mask demand (Betti and Ni, 2020) . Greenness, circularity and bioeconomy bring global sustainability concepts into the mainstream (D' Amato et al., 2017) . Although they are interrelated towards sustainable utilization of natural resources within economic cycles, their dependencies, relations, advantages/disadvantages and best approaches have not been discussed enough in the literature, particularly referring to the public sector such as healthcare. During the COVID-19 pandemic, sustainable development of different supply chain networks becomes doubly important (Barbier et al., 2020; Tirkolaee et al., 2021) . The growing importance of a sustainable supply chain design in today's turbulent world is not hidden from everyone (Mardani et al., 2020b) . The supply chain of face masks as one of the most critical sanitary equipment is highly dependent on these requirements. It is obvious that a sustainable supply chain can minimize The development of a multi-objective, multi-period, multi-echelon, and multi-commodity mathematical model considering location, supply, production, distribution, collection, quarantine, recycling, and burial decisions are some of the novelties of this research. In addition, considering the proposed sustainability in the proposed network for a real case study and simultaneous attention to the three pillars of sustainability including economic, environmental, and social aspects are among the other innovations of this research. The considered location in this research is as a discrete parameter. The location is such that the optimal points are selected from the pre-defined candidate points. The selection of optimal locations for the establishment of facilities will reduce transportation costs and as well as the resulting environmental pollution. The framework of the rest sections is explained as follows. Section 2 provides a survey on similar research studies in the literature. Section 3 investigates the problem and formulates the presented mathematical model. Our proposed solution methodology is provided in Section 4. Section 5 proposes the validation of the presented methods. A real case study is then investigated in Section 6. Eventually, the discussion, concluding remarks, and outlook of the study are given in Section 7 A variety of helpful research studies can be found in the literature, concentrating on CLSCN design (e.g., see Polo et al., 2019; Sazvar et al., 2021) . However, in this section, the most relevant studies related to the supply chain of medical items are reviewed in brief. A resilient-sustainable model was proposed by Zahiri et al. (2017) to develop a supply chain of pharmaceutical distribution and production under uncertainty. They utilized possibilistic-stochastic programming to tackle the uncertainty of parameters. The establishment of production centers using built-in primary technology levels and considering intermodal freight transportation system was their main contribution. Finally, the proposed model was solved using the Pareto-based lower bound approach for the case study and the outcomes represented the proper performance of the mathematical model. Savadkoohi et al. (2018) designed a medical equipment supply chain network under uncertainty. The possibilistic programming method was utilized to cope with uncertainty. The main purpose of the presented model was to minimize supply chain costs in Iran. The location of the distribution and production centers and determining the inventory of the centers according to the perishability of medicines and equipment were among their strategic decisions. Xu et al. (2018) addressed the decision and coordination in the dual-channel supply chain network arising out of low-carbon preference and channel substitution. They developed the decisionmaking models of the centralized and decentralized supply chain network, which contain one manufacturer and one retailer. Additionally, they designed an improved revenue-sharing contract to effectively coordinate the manufacturer and retailer. De Melo Sathler et al. (2019) presented a model for the location of the medical equipment DCs and the allocation of the medical equipment to medical centers. Their main aim was to maximize the level of patient satisfaction in the face of resource scarcity. They considered a case study in Brazil and the results represented that with 4% more specialist staff, 15% more medical equipment can be delivered to medical centers. Nasrollahi and Razmi (2019) developed a medicine supply chain involving manufacturers, distributors, hospitals, and patients. They presented a multi-period model for determining the flow of medicines and transferred medical equipment between centers. The main purpose of their research was to clarify the demand coverage along with minimizing costs. Eventually, the proposed model was solved using two approaches including NSGA-II and Multi-Objective Particle Swarm Optimization (MOPSO). Weraikat et al. (2019) presented a mathematical model for managing the distribution of medical facilities. In addition, they considered sustainability along with the perishable products. The Vendor-Managed Inventory (VMI) approach was utilized to minimize unsatisfied demand. Moreover, the Monte-Carlo simulation approach was employed to cope with the uncertainty of the presented model. The outcomes demonstrated the suitable efficiency of the proposed model. Roshan et al. (2019) formulated a supply J o u r n a l P r e -p r o o f chain network for the distribution of medical equipment under crisis conditions. They designed a multiobjective mathematical model for minimizing distribution and production costs and maximizing unsatisfied demand. In order to increase social responsibility, they considered the maximization of employment rates in established centers. Finally, they reported the proposed sustainable model for the Seattle/USA case study. The main challenges, limitations, opportunities and innovations of solid waste management were discussed by Sharma et al. (2020) during and post COVID-19 pandemic. They pointed out the high importance of efficient medical and hazardous waste (e.g., face masks, gloves and other biomedical waste) management. Sangkham (2020) estimated the number of face masks and the number of mask disposal centers during the COVID-19 outbreak. For this purpose, centers such as hospitals, mask manufacturers, and landfills in Asia were surveyed. The parameters to estimate the number of required masks included the number of population and the acceptance rate of the mask by the community. The acceptance rate of the mask by the community was estimated at 80%. Finally, the number of required masks and the number of disposal centers were estimated for 54 Asian countries including Iran. Klemeš et al. (2020) studied to minimize the effects of plastic waste and waste paper during the outbreak of the COVID-19. In addition, they investigated the effects of pandemics on the amount of plastic waste on the effects of pandemics on environmental pollutants, including carbon dioxide. They took into account a real case study including hospitals, landfills, and households. The obtained results demonstrated that in pandemic conditions, as the amount of released pollutants by vehicles decreases, the amount of plastic pollutants increase sharply, which results in increasing the number of landfills and job creation. Alizadeh et al. (2020) designed a direct and reverse medicine supply chain according to biological risks. In the direct direction, warehouses and DCs are located, and in the reverse direction, medical waste disposal centers are located. The goal of their suggested model was maximizing the difference between revenues and costs and minimizing the time distance between clinics and sterilization. The case study was examined in Tehran/Iran and the Bounded De Novo Programming method was utilized to solve it. Goodarzian et al. (2020) suggested a mathematical model for production and distribution in the pharmaceutical supply chain network. In their proposed model, the costs of ordering, purchasing, and production were considered as a fuzzy parameter and robust fuzzy technique was employed to deal with uncertain parameters. Then, the multi-modal transportation system was considered as one of the contributions. To minimize transportation costs and to minimize medicine delivery time were the most important objectives of their paper. The outcomes represented the suitable efficiency of their suggested model. Kargar et al. (2020a) examined a reverse supply chain for hospital waste disposal management. Therefore, they presented a multi-period and multi-item stochastic model with three objectives including decreasing transportation costs, inventory holding, and established disposal centers along with minimizing all stored hospital waste. Finally, the proposed model was solved by fuzzy goal programming and the outcomes indicate the proper performance of the model. In another research, Kargar et al. (2020b) proposed a mathematical model to optimally deal with hospital waste disposal, such as masks and face shields within a reverse logistics network during the COVID-19 pandemic. The allocation of hospitals to disposal centers along with the location of disposal centers was one of the decisions that did in their research. Their main goal was to minimize location and transportation costs while minimizing uncollected waste from the centers. A revised version of Multi-Choice Goal Programming (MCGP) approach was used to find the best solutions. Franco and Alfonso-Lizarazo (2020) introduced a simulation-optimization model for the distribution of medicines and medical equipment. They considered a bi-level supply chain network including hospitals and pharmacies. Moreover, two stochastic models were proposed for the considered problem. The first model considered the expiration date of medicines along with the level of service to patients and the second model considered inventory control and allocation of medical equipment to hospitals. Finally, the proposed J o u r n a l P r e -p r o o f models were solved using the ε-constraint approach and the results indicated a 16% reduction in supply chain costs. Recently, Vanapalli et al. (2021) provided strategies and challenges for effective plastic waste management in the COVID-19 pandemic. They addressed the inefficiencies and inadequacies of the current waste management system to cope with the raised dependence on plastic. Then, these problems aggravated its leakage and mismanagement into the environment, then triggered a novel environmental crisis. Moreover, a futuristic review was conducted by Purnomo et al. (2021) on the possible thermochemical treatments for the COVID-19 related medical waste. They also examined the possible conversion to fuels by analyzing the characteristics of waste including face masks, clothes, goggles, gloves and disinfectant containers. Mosallanezhad et al. (2021) designed a disaster relief supply chain network for personal protection equipment in the COVID-19 condition. They formulated a multiobjective, multi-product, and multi-period model to minimize total cost and shortage considering the COVID-19 outbreak. To solve their model, they used four well-known multi-objective metaheuristics. They suggested the Taguchi method to tune the algorithm's parameters and compare the proposed algorithms based on assessment metrics. Goodarzian et al. (2021a) developed an integrated sustainable medical supply chain network in the COVID-19 pandemic. They designed a production-distributioninventory-allocation-location problem for the multi-objective, multi-level, multi-product, and multiperiod model. To solve their model, hybrid methods based on meta-heuristic algorithms were developed. In addition, they suggested the response surface method to tune the algorithms' parameters. Finally, a real case study was considered to validate their proposed model and show the performance of their proposed methods. Xu et al. (2021) investigated the effect of the COVID-19 epidemic on port performance in China. They considered a panel regression model to analyze the effect on port efficiency. Finally, 14 major ports in China were chosen as a case study. Goodarzian et al. (2021b) proposed a new sustainable-resilience healthcare network under uncertainty to handle the COVID-19 epidemic. They formulated and designed an optimization-simulation mixed-integer linear programming model for a multi-period, multi-product, multi-objective, and multi-echelon problem. To cope with uncertain parameters, stochastic chance-constraint programming was proposed. They developed three hybrid meta-heuristic algorithms to solve their model. Finally, they suggested a real case study to show the efficiency of their developed model. Pedroza-Gutiérrez et al. (2021) explained adaptive governance and coping strategies in the Yucatan Peninsula coasts facing COVID-19. They considered the COVID-19 crisis as a social driver of change with new and different forms of governance. Also, for coastal areas, a multi-level governance strategy to recover from the COVID-19 epidemic and economic crises was needed. Additionally, coastal communities had difficulties demonstrating non-fisheries-based or non-recreate-based resilience. Milanes et al. (2021) presented forty-three positive and negative practices and actions on beaches during the COVID-19 lockdown. In addition, they designed twenty-seven new recommendations for coastal planning and beach management, of which several recommendations were validated in the Republic of Cuba. Finally, experts in beach and coastal management from two Ibero-American networks were consulted. Goodarzian et al. (2021c) proposed a new multi-objective optimization approach for the pharmaceutical supply chain network design to minimize the total cost and the delivery time of pharmaceutical products while maximizing the reliability of the transportation system. Also, a new mixed-integer non-linear programming model was developed for the productionallocation-distribution-inventory-ordering-routing problem. Three new heuristics (H-1), (H-2), and (H-3) and two new meta-heuristic algorithms, namely, an Improved Social Engineering Optimization (ISEO) and Hybrid Firefly and Simulated Annealing Algorithm (HFFA-SA) were developed to solve their model. Finally, their proposed model was evaluated through extensive simulation experiments by analyzing different criteria. According to the examined papers, the novelties of this study are discussed as follows: i. Developing a novel multi-objective, multi-period, multi-level, multi-product MILP model to design a sustainable mask CLSCN during the COVID-19 pandemic, ii. Addressing locational, supply, production, distribution, collection, quarantine, recycling, reuse and disposal decisions within the proposed network, iii. Studying the sustainability in mask CLSCN design problem in order to minimize the total human risk as much as possible and besides minimizing the total cost and total environmental pollution, iv. Developing two multi-objective meta-heuristic algorithms of Multi-Objective Grey Wolf Optimization (MOGWO) and NSGA-II to tackle the complexity of the problem, v. Enhancing the performance of the algorithms using the Taguchi design method, vi. Evaluating the performance of the algorithms using well-known metrics, problem instances and a statistical test, vii. Validating the proposed model and solution methods using a real case study in Tehran/Iran. In this section, a sustainable mask CLSCN is modeled to deal with the production, distribution, and recycling processes of the face masks in the COVID-19 pandemic condition. In the forward supply chain structure, different levels of suppliers, factories, DCs, and customer centers are defined. On the other hand, in the reverse supply chain structure, recycling centers, collection centers, quarantine centers, and disposal centers are considered. In addition to raw materials obtained from suppliers, factories produce different types of masks (N95, KN95, and surgical masks). With the help of the DCs, these products are transferred to customer centers including hospitals and infirmaries. A major part of demand at customer centers is then collected using sanitary packaging and transported to collection centers, then they are taken to the quarantine centers and placed there for four days. They are then crushed into small pieces and exposed to ultraviolet light to guarantee that they are thoroughly disinfected before the recycling process begins. In the next step, the disinfected masks are transported to recycling centers. Here, there are two output flows from recycling centers to disposal centers and suppliers. In forward logistics, since many products are dispatched and distributed among different levels, a fleet of vehicles with different sizes are used to handle the transportation process, but in reverse logistics, since fewer products are transported, a fleet of mini vehicles are used according to the economical aspect. To address sustainable development, three objective functions are studied including social, economic, and environmental impacts. The 1st objective is minimizing the whole cost of the supply chain containing the establishment costs, transportation costs, and employment costs of vehicles. The 2nd objective tries to minimize the whole pollution emission throughout the supply chain. Finally, the 3rd objective seeks for minimizing the infection risk imposed on the population around newlyestablished sites. To better understand the components of the proposed CLSCN, Figure 1 represents a schematic view. The assumptions of the presented model are explained as follows:  Multiple levels of suppliers, factories, DCs, customer centers, disposal centers, recycling centers, collection centers, and quarantine centers are considered to configure the proposed closed-loop supply chain,  Locational decisions are made at the levels of factories, DCs, disposal centers, collection centers, quarantine centers, and recycling centers at the beginning of the time horizon,  Operations and transportations between different levels and operational activities at different facilities are the origins of pollution emission,  Two types of transportation systems are defined including transportation trucks and recycling mini trucks,  All produced masks at factories should be transported to DCs in each time period,  Facilities located on reverse logistics flow yield infection risks,  Parameters are deterministic,  Capacity of facilities is limited,  A time horizon is regarded,  Four raw materials and three final products are needed to be supplied,  Shortage in the main market (customer centers) is not allowed. Now, the proposed mathematical model is given as follows: Set of suppliers ( ∈ ), Set of factories ( ∈ ), Set of DCs ( ∈ ), Set of customer centers ( ∈ ), J o u r n a l P r e -p r o o f Set of collection centers ( ∈ ), Set of quarantine centers ( ∈ ), Set of recycling centers ( ∈ ), Set of disposal centers ( ∈ ), Set of raw materials ( ∈ ; | | = 4), Set of masks ( ∈ ; | | = 3), Set of disposal centers ( ∈ ), Set of transportation trucks ( ∈ ), Set of recycling mini trucks ( ′ ∈ ′), ′ Set of planning periods ( ∈ ). Fixed establishment cost of factory , Fixed establishment cost of DC , Fixed establishment cost of collection center , Fixed establishment cost of quarantine center , Fixed establishment cost of recycling center , Fixed establishment cost of disposal center , Fixed cost of using transportation truck , ′ ′ Fixed cost of using mini recycling truck ′, Unit purchasing cost of material from supplier , Unit production cost of mask by factory , Unit processing cost of mask by DC , Unit processing cost of mask by collection center , Unit processing cost of mask by quarantine center , Unit processing cost of mask by recycling center , Unit processing cost of mask by disposal center , Unit cost of the transportation of material from supplier to factory , Unit cost of the transportation of mask from factory to DC , Unit cost of the transportation of mask from DC to customer center , Unit cost of the transportation of mask from customer center to collection center , Unit cost of the transportation of mask from collection center to quarantine center , Unit cost of the transportation of mask from quarantine center to recycling center , Unit cost of the transportation of mask from recycling center to disposal center , Unit cost of the transportation of mask from recycling center to supplier , Unit consumption coefficient of material to produce mask , Demand of mask at customer center in time period , Return rate of used masks to be transported from customer center to collection centers in time period , Conversion rate of material to be transported to suppliers from recycling center in time period , Disposal rate of masks transported to recycling center in time period , Capacity of supplier to provide material in time period , Capacity of factory to produce mask in time period , Capacity of DC to process mask in time period , Capacity of collection center to process masks in time period , J o u r n a l P r e -p r o o f Capacity of quarantine center to process masks in time period , Capacity of recycling center to process masks in time period , Capacity of disposal center to process masks in time period , Capacity of supplier to store recycled material transported by recycling centers in time period , Capacity of transportation truck , ′ ′ Capacity of mini recycling truck ′, Unit pollution emission for transporting material from supplier to factory , Unit pollution emission for transporting mask from factory to DC , Unit pollution emission for transporting mask from DC to customer center , Unit pollution emission for transporting mask from customer center to collection center , Unit pollution emission for transporting mask from collection center to quarantine center , Unit pollution emission for transporting mask from quarantine center to recycling center , Unit pollution emission for transporting mask from recycling center to disposal center , Unit pollution emission for transporting mask from recycling center to supplier , Unit pollution emission for providing material at supplier , Unit pollution emission for producing mask at factory , Unit pollution emission for processing mask at DC , Unit pollution emission for processing masks at collection center , Unit pollution emission for processing masks at quarantine center , Unit pollution emission for processing masks at recycling center , Unit pollution emission for processing masks at disposal center , Population size around collection center , Population size around quarantine center , Population size around recycling center , Population size around disposal center . J o u r n a l P r e -p r o o f (1) , , , ∀ ∈ , ∈ , ∈ , ∈ , ∈ , ∈ , ∈ , ∈ , ∈ , ′ ∈ ′ , , ′ , ′ , ′ , , ′ ≥ 0, , ∈ + , ∀ ∈ , ∈ , ∈ , ∈ , ∈ , ∈ , ∈ , ∈ , ∈ , ′ ∈ ′ , ∈ , ∈ . The first objective function (1) shows minimizing the total cost of the CLSCN design including the fixed establishment costs, fixed employment costs of trucks, transportation costs, and operational costs at different facilities. The second objective function (2) indicates minimizing the whole pollution emission incurred by transportation and operational processes. The third objective function (3) states minimizing the total human risk imposed by the establishment of facilities such as the number of people around collection centers, quarantine centers, recycling centers, and disposal centers is minimized. Eqs. (4)-(10) represents the capacity constraints of suppliers, factories, DCs, collection centers, quarantine centers, recycling centers, and disposal centers, respectively. Furthermore, decisions on the facility location at different levels can be made according to these equations. Eqs. (12)-(19) indicate the capacity constraints of trucks employed by suppliers, factories, DCs, collection centers, quarantine centers, recycling centers, and disposal centers, respectively. Eq. (20) states the flow balance from suppliers to factories considering consumption rates of raw materials to produce masks. Eq. (21) In this research, MOGWO and NSGA-II are applied to tackle the complexity of the problem. The main reasons for choosing MOGWO and NSGA-II algorithms can be stated in two parts. I. The purpose of this study is to examine the performance of novel algorithms in the field of supply chain network design. Among the emerging multi-objective algorithms, MOGWO algorithm is chosen because of its strong structure in terms of both intensification and diversification (Mirjalili et al., 2016) . II. To evaluate the performance of the proposed MOGWO algorithm, comparisons must be made with one of the strongest and well-known algorithms. In this regard, studies have demonstrated that the NSGA-II algorithm is one of the best multi-objective meta-heuristic algorithms . Therefore, in this research, MOGWO and NSGA-II algorithms are implemented and their results are evaluated both in terms of quality and dispersion of Pareto solutions using several performance metrics. Gray Wolf Optimization (GWO) algorithm is a meta-heuristic algorithm inspired by the hierarchical structure and social behavior of gray wolves while hunting. This algorithm is population-based, has a simple process, and can be easily generalized to large-scale problems (Makhadmeh et al., 2021) . The GWO was introduced by Mirjalili et al. (2014) . Gray wolves are regarded as apex predators, which are at the top of the food chain pyramid. Gray wolves prefer to live in a group, each group has an average of 5-12 members. All members of this group have a very precise hierarchy of social domination and have specific tasks. In each herd of wolves for hunting, there are 4 degrees, which is modeled as a pyramidal structure including alpha ( ), beta ( ), delta ( ), omega ( ). Leader wolves are called alpha groups, which can be male or female. These wolves dominate the herd. Beta wolves help alpha wolves in the decision-making process and are also prone to be chosen instead. Delta wolves are lower than beta wolves including older wolves, predators, and baby care wolves. Omega wolves show the lowest rank in the hierarchy that has the least rights over the rest of the group. After all, they eat and do not participate in the decision-making process. In GWO, the most appropriate solution is taken into account as alpha, and the second and third appropriate solutions are named beta and delta, respectively. The rest of the solutions are considered omega. In the GWO, hunting is driven by alpha, beta, and delta. Solution omega follows these three wolves. To model the circumvention behavior, Eqs. (30) and (31) Grey wolves have the ability to detect the position of prey and bypass them. Hunting is guided by Alpha. Delta and Beta may also occasionally participate in hunting. Moreover, in an absolute search space, we have no solution for the optimal position (hunting). To mathematically simulate the hunting behavior of grey wolves, we suppose that Alpha (the best candidate solution), Delta, and Beta are sufficiently aware of the potential hunting position. So, the first three obtained solutions are saved and other search agents are obliged to update their position based on the position of the best search agents. This operation is performed based on Eqs. (34)-(36): ( + 1) = 1 + 2 + 3 3 . In GWO, the search process starts by creating a random population of grey wolves (candidate solutions). In the iteration time, Beta, Delta, and Alpha wolves estimate the probable hunting position. Each candidate solution updates its distance with the prey. Parameter is decreased from 2 to 0 to boost the detection and attack process to prey. If | |>1, candidate solutions diverge; else if | |<1, candidate solutions converge. In the MOGWO, the general framework of the grey wolf algorithm is implemented. The variation is that each iteration of the non-dominated solutions is separated in the manner presented in Sub-section 4.2.1. The pseudo-code of the MOGWO is represented in Figure 2 . To rank a set of solutions and place them on various fronts in terms of the degree of non-dominated, first, the following two parameters are computed for each solution; , which is actually the number of solutions that have overcome the solution , and , which is the set of solutions that have been dominated by the solution . All solutions with =0 are ranked first and placed in the first front. Then, for every solution of with rank one, every member ( ) of the set is met and its is reduced by one unit. Thus, if = 0, then the solution takes second place and is in the second front. In the same way, the next fronts are formed. Figure 3 shows the pseudo-code for how to fast non-dominated sorting (FNS) in the proposed MOGWO. NSGA-II algorithm is one of the most widely-utilized and robust algorithms for treating multi-objective optimization problems and has an acceptable efficiency (Onan et al., 2015; Hu et al., 2020) . Srinibas and Deb (1995) suggested the NSGA optimization method for solving multi-objective optimization problems. In general, the structure of the NSGA-II is given in Figure 4 . In the first step of a multi-objective genetic algorithm, several initial solutions must be generated. Initial solutions must be randomly selected to cover the set of justified spaces. Therefore, random number generation functions are used in software environments. Fitness means that one solution is superior to another. In optimization problems, the superiority of one solution over another is specified by the objective function of the problem. Therefore, the fitness function will be the same as the objective function of the mathematical model. If the objective function is of the minimum sort, the fitness function is defined as − or 1/ ( represents the formula of the objective function), but if the objective function is of the maximum sort, the concept of fitness is exactly the same as the objective function value. It is first necessary to define the concept of "dominance": I. dominates all members of set . The members of set shows points that are at least in one sense worse than and in no way better than . II. is dominated by all members of set . Here, is not better than any of the members in set and is worse in at least one respect. However, the members in sets and are better in one sense, so we cannot choose one. In this step of NSGA-II, between the set of solutions, we select and separate the ones that are nondominated to the other solutions are selected and separated. These non-dominated solutions are essentially the Pareto front or the front of efficiency found by the genetic algorithm. This Pareto front is improved in each iteration to report the best Pareto front by the NSGA-II at the end of the genetic algorithm. The crossover operation in the NSGA-II shows exactly the same as the crossover operation in the singleobjective genetic algorithm. In this section, the solutions that were not non-dominated are separated, then composed in pairs according to a predetermined mechanism (tournament, random, or roulette wheel), and generate new solutions. This composition of solutions is done in methods such as uniform, single-point, two-point, etc. The aim of the crossover is to compose the features of various solutions with the aim of generating better solutions. The mutation operation in the NSGA-II is exactly the same as the mutation operation in the singleobjective genetic algorithm. At this stage, some solutions are chosen at random and some of their cells are modified at random. The aim of the mutation is to make sudden alters in the solutions in order to escape the trap of optimal local. Each solution of the studied problem consists of two different structures. The first structure is for location decision variables and the second structure is for transport decision variables in the presented network. In the first structure, a set of numbers between 0 and 1 is presented and numbers less than 0.5 mean no selection/establishment and numbers above 0.5 mean selection/establishment. For example, considering 5 potential points for the establishment of a DC, Figure 5 shows an example of the first structure for the location of the DCs. According to Figure 5 , DCs 1, 2, and 4 are established. This structure is repeated for all network facilities (transport centers and equipment) that must be decided about selected/established. In the second structure, which determines the volume of transmission between different bases and destinations, a matrix with values between 0 and 1 is used. The number of rows in this matrix is equal to the number of bases and the number of columns is equal to the number of destinations. It should be noted that these matrices are repeated for each time period. For example, if we have 3 distributors and 4 customers, an example of the solution is indicated in Figure 6 . In order to decode the above matrix, several steps are performed sequentially. First, according to the generated random numbers from large to small, we start to meet the needs of the customer until all the needs of the customer are met. The first step is represented in Figure 7 . According to the selected random number, the amount of customer demand number 4 should be answered as much as possible first. Therefore, the amount of 200 units from distributor number 2 is allocated to customer 4. In the next step, a larger random number is selected from the available random numbers, the example of which is indicated in Figure 8 . According to the selected random number, the amount of 120 units is transferred from distributor number 3 to customer 2 and the capacity and demand values are updated (according to Figure 9 ). This process continues until the end of full coverage of customer demand. Figure 10 shows how full coverage of customer demand by distributors. Therefore, customer demand was met in this way. Finally, the sent quantities from the distributor to the customer are presented in Figure 11 . This designed structure for different types of bases, destinations, and in different periods as well as for different types of transportation systems are repeated to finally provide a complete solution to the closed-looped supply chain network design problem. In this section, the suggested indicators are usually utilized to compare the performance of metaheuristic algorithms that are explained in the next sub-sections. This index is utilized to compute the amplitude of the optimal Pareto front solutions gained by the algorithm. The larger the value of this index, the higher the priority of the algorithm. The value of this index is computed using Eq. (37): where Min and Max represent the minimum and maximum values of the objective function between all obtained non-dominated solutions from the algorithm, respectively. Using this metric, a close distance is obtained among the obtained non-dominated solutions and the ideal point. The smaller the value of this index, the higher the priority of the algorithm. The MID is calculated utilizing Eq. (38). where 1 and 2 are considered as the largest and the smallest value between non-dominated solutions, respectively. This index also called the spread metric, is utilized to calculate the degree of variation of Pareto solutions. The larger the value of the SNS, the higher the priority of the algorithm. The value of the SNS is computed using Eq. (39). In Eq. (39), represents the number of non-dominated solutions and also the value of C i is calculated using Eq. (40): where Eq. (40), 1 and 2 are the values of the first and second objective functions, respectively, for the non-dominated solution . The NPS is employed to calculate the number of the gained non-dominated solutions by the proposed algorithm. The larger the NPS metric, the higher the priority of the algorithm. Here, the Taguchi design approach is used to tune the parameters of the algorithms efficiently. To design experiments in the proposed algorithms, firstly, three different levels are defined for its parameters. And then the predefined tests are performed in this algorithm utilizing random test problems in mediumsized. The proposed values for the parameters of the MOGWO and NSGA-II algorithms are shown in Table 1 . Then, with Taguchi L9 design, different experiments were created and the MOGWO and NSGA-II algorithms were implemented for each of them. The implementation results are reports in Table 2 . Therefore, by presenting these outputs by using MINITAB software, the S/N ratios of the proposed algorithms are displayed in Figures 12 and 13 . Hence, according to the presented results in Figures 12 and 13 , the best value of each parameter is determined and other examples are executed with these values of the algorithm parameters. Table 3 presents the optimal values of the parameters. In this part, to validate the presented problem, the efficiency of the two algorithms was investigated. According to the superiority of NSGA-II in the MID metric is very slight, so in general, the MOGWO algorithm can be introduced as the superior algorithm. Since this analysis may not be complete and sufficient, thirty numerical problems in different sizes are generated and comparison metrics were calculated for both algorithms. It should be noted that the size of the test problems is divided into three groups: small, medium and large. The basis of this division is based on modeling from (Pahlevan et al., 2021) . The input information for these thirty problems is presented in Table 4 . Also, the results related to different assessment metrics for the MOGWO algorithm are presented in Table 5 and for the NSGA-II algorithm in Table 6 . Table 4 . Information of the problem instances. #Suppliers #Factories #Collection centers #DCs #customers #Quarantine centers #Recycle centers #Disposal centers #Periods #Trucks #Mini trucks Small scale 1 1 2 2 1 4 1 1 1 1 10 5 2 1 3 2 2 5 1 1 2 1 10 5 3 2 3 3 3 6 1 2 3 1 10 5 4 2 4 3 4 7 2 2 3 2 15 8 5 3 4 4 5 8 2 2 4 2 15 8 6 3 5 4 7 9 2 3 4 2 15 8 7 4 5 5 9 10 3 3 5 3 20 10 8 4 6 5 10 11 3 3 5 3 20 10 9 5 6 6 11 12 3 4 6 3 20 10 10 5 7 6 12 13 4 4 6 4 25 13 Median Scale 11 6 8 7 13 15 6 4 7 6 30 15 12 7 10 9 15 20 8 5 7 6 35 18 13 8 12 11 17 25 10 5 7 8 40 20 14 9 14 13 19 30 12 5 8 8 45 23 15 10 16 15 21 35 14 6 8 10 50 25 16 11 18 17 23 40 16 6 8 10 55 28 17 12 20 19 25 45 18 6 9 12 60 Problem size #P #Suppliers #Factories #Collection centers #DCs #customers #Quarantine centers #Recycle centers #Disposal centers #Periods #Trucks #Mini trucks 18 13 22 21 27 50 20 7 9 12 65 33 19 14 24 23 29 55 22 8 9 14 70 35 20 15 26 25 30 60 24 9 10 14 75 38 Large Scale 21 20 30 30 40 70 26 10 10 16 80 40 22 25 40 35 50 80 28 11 10 18 90 It is clear that in Tables 5 and 6 , both algorithms have been able to solve the optimal solutions of thirty problems and the process of the presented values can be seen quite logically. In the following, the results related to each metric are analyzed graphically. The obtained results of comparing the two algorithms MOGWO and NSGA-II based on the DM metric in Tables 5 and 6 , and also Figure 14 shows that the average value of DM for the MOGWO shows 583.07 and the NSGA-II indicates 468.12. As shown in Figure 14 , in all the solved problems, the MOGWO algorithm has obtained more than the DM metric. Overall, the MOGWO algorithm performed about 25% better than the NSGA-II in terms of DM metric. The results of comparing the two algorithms MOGWO and NSGA-II based on the MID metric in Tables 5 and 6 and also Figure 15 show that the average value of MID for the MOGWO shows 0.434 and the NSGA-II is equal to 0.424. Also, as shown in Figure 15 , the two algorithms have obtained very close values in the MID assessment metric. In total, the MOGWO algorithm performed about 2% more than the NSGA-II in the MID. In terms of the MID, this indicates the very slight superiority of the NSGA-II. The results of comparing the two algorithms MOGWO and NSGA-II based on the SNS metric in Tables 5 and 6 as well as Figure 16 show that the average SNS value for the MOGWO indicates 732.93 and the NSGA-II shows equal to 613.77. Also, as illustrated in Figure 16 , in all the solved problems, the MOGWO algorithm has obtained more than the SNS metric. The MOGWO algorithm shows about 19% superiority against the NSGA-II in terms of the SNS metric. The results of comparing the two algorithms MOGWO and NSGA-II based on the NPS metric in Tables 5 and 6 as well as Figure 17 show that the average NPS value for the MOGWO displays 7.83 and the NSGA-II represents equal to 6.87. Also, as represented in Figure 17 , in all the solved problems, the MOGWO algorithm has obtained more than the NPS metric. All in all, the MOGWO algorithm performed about 13% better than the NSGA-II in terms of the NPS metric. To further assess the indicators, a paired t-test comparison between algorithms was performed. This test is performed separately for each indicator. The confidence level is 95%. Details of the results are provided in Table 7 . It is clear that in Table 7 , the Sig value for the DM, SNS, and NPS assessment metrics is approximately equal to zero. Accordingly, the hypothesis test related to them is meaningful. This meaningful interpretation shows that the difference between the MOGWO and NSGA-II is also statistically proven and therefore the superiority of the MOGWO algorithm in terms of these assessment metrics is confirmed. Regarding the MID metric, the sig value is 0.058, which is more than 0.05, and therefore the significance of the difference among the two algorithms in terms of the MID metric is not confirmed, and it can be said that these two methods have the same average performance. Overall, it can be argued that the MOGWO meta-heuristic algorithm, with its new and more powerful tools than the NSGA-II algorithm, can provide a better and stronger Pareto front on multiobjective optimization problems. Tehran is the most populous city in Iran and one of the densest cities in Asia. The city, with a population of about 9 million people has faced many problems suffered great human and financial losses in the COVID-19 outbreak. The lack of masks is one of the main problems in this city and the supply chain management of the production, distribution, recycling, and disposal of masks in this city seems necessary. Therefore, in this study, the city of Tehran has been selected as a case study. Figure 18 shows a case study map. As it is clear that, 6 suppliers, seven manufacturers, 6 DCs, 8 customer centers, 4 collection centers, 7 quarantine centers, 8 recycling centers, and 9 disposal centers are considered. Table 8 indicates the fixed cost of established quarantine centers and factories in dollars. For example, the fixed cost of an established Evin DC is $250,000. Then, the transportation costs from collection centers to quarantine centers for different types of masks in dollars are reported in Table 9 . For instance, the cost of transportation from Sabalan Collection Center to Shariati Quarantine Center for N95, KN95 masks, and surgical masks is $0.2, $0.4, and $0.3, respectively. Therefore, the raw material capacity for suppliers and the demand for different types of masks in the customer centers in kilograms are displayed in Tables 10 and 11, respectively. It is clear that in Table 12 , in Pareto optimal solutions, the values of the second objective function decreased as the first objective function value increased. In other words, the worse the solutions are in terms of the first objective function, the better the second function. The same is true of the second and third objective functions. The better the values of the second objective function, the worse the values of the third objective function. Therefore, none of the solutions is completely superior to the other and all the solutions are non-dominated. In order to examine the outputs of MOGWO in more detail, a solution from this algorithm is presented graphically in Figure 20 . As can be seen in Figure 20 , the obtained result from the MOGWO algorithms based on all the assumptions and conditions of the mathematical model are considered and are therefore sufficiently accurate. In this sub-section, the effect of one of the most important and significant parameters of the mathematical model on the Pareto front of the problem is investigated. In this regard, the raised problem in the validation section is examined. In the following, the value of the demand parameter fluctuates between -20% to + 20%. For example, if the demand for a customer is equal to 10, the demand value between 8 and 12 is checked for that. Then, the values of Pareto solutions are reported. The outcomes are reported in Table 13 . In order to better analyze the results, the average values of each of the objectives are reported in Table 14 Sensitivity analyses indicate that raising the amount of demand causes raising each of the objective functions. In terms of the total cost, this trend is broken linearly, as well as in terms of total environmental pollution. However, in terms of human risk, there is a non-linear trend that with increasing demand, the slope of raising the third objective function value has also increased. Therefore, customer demand is a very important and key parameter that can have drastic changes in any of the objective functions. Central governments and municipalities are mainly responsive for defining and enacting public strategies to move towards sustainable development during pandemics, especially in line with healthcare sectors as the frontline organizations to restrain destructive impacts as much as possible. The supply chain design of sanitary equipment such as face mask is one of the most significant issues in this regard. Evaluating the effects of policy actions within such a supply chain is in need of high-quality data, indicators and measures to surveil the chain performance considering the use of raw materials, production and distribution of final products, the management of waste flows, disposal, recycling and circularity practices. Accordingly, this study tried to configure a sustainable CLSCN of face masks during the COVID-19 outbreak using a novel multi-objective MILP for the first time. To this end, a multi-level, multi-period, multi-product network was considered in order to address the locational, supply, production, distribution, collection, quarantine, recycling, reuse, and disposal decisions. The objectives were to simultaneously minimize the total cost, environmental pollution and human risk caused by the establishment of facilities. To tackle the problem efficiently, two meta-heuristic algorithms of MOGWO and NSGA-II were utilized, which were calibrated optimally using the Taguchi design method. Then, four assessment metrics of MID, DM, SNS, and NPS were applied to assess the performance of the suggested algorithms. The outcomes showed that the average of the MID for MOGWO and NSGA-II is 0.4343 and 0.4244, respectively, and also the average of the DM for the MOGWO and NSGA-II was 583.0724 and 468.1279, respectively. Furthermore, the average of the SNS for MOGWO and NSGA-II was 732.9315 and 613.7712, respectively. The MOGWO outperformed the NSGA-II for 25%, 2%, 19%, and 13% in terms of the DM, MID, SNS, and NPS metrics, respectively. Accordingly, the outcomes of the t-test with a 95% confidence interval indicated the superiority of the MOGWO algorithm. In the next stage, a real case study was investigated in Tehran/Iran using the MOGWO algorithm. The optimal policy was determined and discussed accordingly. Furthermore, the outcomes of the sensitivity analyses demonstrated that as demand increases, the total cost increase with a linear behavior approximately. Similarly, with increasing demand, the amounts of environmental pollution and human risk also increase. As one of the most important management insights that can be considered in this research is providing several robust solutions for decision-makers in order to facilitate the process of decision-making. In this study, some intelligent methods were utilized to solve the problem that can ultimately provide several optimal solutions. This feature helps supply chain managers to choose and implement the most appropriate solutions based on their preferences which cannot be always considered in such optimization models. Here, to be more specific and according to the sensitivity of human risk and the possibility of human losses due to the establishment of facilities, it is suggested to managers and decision-makers to specify the potential locations of facilities as much as possible in areas with less population. This may increase transportation costs, but it seems necessary according to the priority of human risk. Moreover, air pollution can double the effects of COVID-19, especially in Tehran. Therefore, managers and decision-makers are advised to use hybrid vehicles that have less environmental pollution. Moreover, the optimal location of facilities can greatly reduce environmental pollution. In addition, supply chain costs can be reduced by using recycling operations. Finally, managers are recommended to focus on recycling operations due to high production costs. The main limitations of the current paper include the following:  Usually, there is no specific database of transportation costs, so in this study, drivers' assessments were employed to estimate transportation costs,  The demand amount was just estimated based on the report of the Ministry of Health and based on the report of the medical centers,  Uncertainty was not addressed in the problem. Finally, suggestions for future works are as follows:  Estimating the demand parameter using simulation methods to increase the accuracy,  Providing other aims such as minimizing the shortage of products and maximizing the reliability of the network to make the model more compatible with sustainable development,  Employing uncertainty approaches such as fuzzy programming (Das et al., 2020; Tirkolaee et al., 2021) and robust optimization (Özmen et al., 2017; Tirkolaee et al., 2020a) , grey systems (Roy et al., 2017; and stochastic optimal control (Temoçin and Weber, 2014; to deal with the uncertain nature of the problem,  Integrating routing and inventory decisions in the proposed mask CLSCN for optimizing the problem more practically (Alinaghian et al., 2021) . An augmented Tabu search algorithm for the green inventory-routing problem with time windows. Swarm and Evolutionary Computation Surgical mask partition reduces the risk of noncontact transmission in a golden Syrian hamster model for Coronavirus Disease 2019 (COVID-19) Clinical progression of patients with COVID-19 in Shanghai The role of community-wide wearing of face mask for control of coronavirus disease 2019 (COVID-19) epidemic due to SARS-CoV-2 A comparative analysis of sustainability avenues Integration of facility location and equipment allocation in health care management Application of Type-2 Fuzzy Logic to a Multiobjective Green Solid Transportation-Location Problem With Dwell Time Under Carbon Tax, Cap, and Offset Policy: Fuzzy Versus Nonfuzzy Techniques The prospects of waste management in the hospitality sector post COVID-19. Resources, Conservation and Recycling, 105272 A multi-objective pharmaceutical supply chain network based on a robust fuzzy model: A comparison of meta-heuristics Multi-objective based scheduling algorithm for sudden drinking water contamination incident. Swarm and Evolutionary Computation Comparative study of variational quantum circuit and quantum backpropagation for COVID-19 outbreak predictions A reverse supply chain for medical waste: A case study in Babol healthcare sector Reverse logistics network design for medical waste management in the epidemic outbreak of the novel coronavirus (COVID-19) Minimising the present and future plastic waste, energy and environmental footprints related to COVID-19 A novel extended approach under hesitant fuzzy sets to design a framework for assessing the key challenges of digital health interventions adoption during the COVID-19 outbreak Evaluation of green and sustainable supply chain management using structural equation modelling: A systematic review of the state of the art literature and recommendations for future research Multi-objective grey wolf optimizer: a novel algorithm for multi-criterion optimization A transformative supply chain response to COVID-19 Disaster relief supply chain design for personal protection equipment during the COVID-19 pandemic A mathematical model for designing an integrated pharmaceutical supply chain with maximum expected coverage under uncertainty Intelligent Recommender Systems in the COVID-19 Outbreak: The Case of Wearable Healthcare Devices An evolutionary multi-objective optimization approach to disaster waste management: A case study of Istanbul Robust optimization in spline regression models for multi-model regulatory networks under polyhedral uncertainty Sustainable supply chain network design using products' life cycle in the aluminum industry Adaptive governance and coping strategies in the Yucatan Peninsula coasts facing COVID-19. Ocean and Coastal Management, 212, 105814 Multi-objective two-stage grey transportation problem using utility function with goals Designing a sustainable closed-loop pharmaceutical supply chain in a competitive market considering demand uncertainty, manufacturer's brand and waste management Face mask and medical waste disposal during the novel COVID-19 pandemic in Asia. Case Studies in Chemical and Environmental Engineering, 2, 100052 Challenges, opportunities, and innovations for effective solid waste management during and post COVID-19 pandemic. Resources, conservation and recycling A robust green location-allocation-inventory problem to design an urban waste management system under uncertainty Multi-objective optimization for the reliable pollution-routing problem with cross-dock selection using Pareto-based algorithms Challenges and strategies for effective plastic waste management during and post COVID-19 pandemic Improving sustainability in a two-level pharmaceutical supply chain through Vendor-Managed Inventory system Coronavirus disease 2019 (COVID-19): situation report Decision and coordination in the dual-channel supply chain considering cap-and-trade regulation The effect of COVID-19 pandemic on port performance: Evidence from China. Ocean and Coastal Management Expert Recommendations for Tracheal Intubation in Critically III Patients with Noval Coronavirus Disease Table 11 . Demand distribution for the masks at customer centers (Kg). Valiasr N95 1100 1000 850 1200 1000 900 1400 1800 KN95 1500 1250 1000 1000 1600 1350 1500 2000 Surgical 500 800 300 750 500 650 450 700 After optimizing this problem with MOGWO algorithms, the Pareto front of each is obtained, which is shown in Figure 19 . Moreover, the set of Pareto solutions are presented in Table 12 . ☒ The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.☐The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: J o u r n a l P r e -p r o o f