key: cord-0759780-tnv37qv4
authors: Asadi, Zeinab; Khatir, Mohammad Valipour; Rahimi, Mojtaba
title: Robust design of a green-responsive closed-loop supply chain network for the ventilator device
date: 2022-03-14
journal: Environ Sci Pollut Res Int
DOI: 10.1007/s11356-022-19105-1
sha: 2566afe352f26bd346d291454dc4b1843c5fdde1
doc_id: 759780
cord_uid: tnv37qv4

This study aims to investigate the closed-loop supply chain network design problem considering the environmental and responsiveness features. For this purpose, a multi-objective mathematical model is suggested that minimizes the carbon emissions and the total costs and maximizes the responsiveness of the system. Due to the dynamic space of the business environment, uncertainty is an integral part of the supply chain problem. Therefore, this research applies the robust possibilistic programming method to cope with uncertainty. Afterwards, since the research problem has a high level of the complexity, a hybrid solution approach based on a heuristic method and the meta-goal programming method is developed to solve the research problem in a reasonable time. Then, due to the importance of the ventilator device during the recent pandemic (COVID-19), this study considers this product as a case study. The main contribution of the current study is to design a green-responsive closed-loop supply chain network under uncertainty using a multi-objective robust possibilistic programming model, for the first time in the literature, especially in the medical devices industry. On the other side, the other contribution of this study is to develop an efficient hybrid solution method. The achieved results demonstrate the efficiency of the offered model and the developed hybrid method. Eventually, by carrying out sensitivity analysis, the impact of some of the critical parameters on the model is investigated. Based on the obtained results, an increase in the demand sizes leads to increasing the environmental damages and the total costs while reducing the responsiveness level. On the other side, an increase in the rate of return leads to an increase in all of the objective functions. Also, the achieved results show that when the capacity parameter is increased, the total costs are decreased, but the responsiveness and environmental impacts are increased. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11356-022-19105-1.

Nowadays, increasing the level of competition in the market space results in increasing the importance of the supply chain (SC) management problem. In general, a supply chain consists of all the elements that act directly or indirectly in meeting the demands of the customer (Chopra and Meindl 2007; Fallahpour et al. 2018) . SC management is the handling of the flow of materials and services and includes all the processes that convert raw materials into final products (Chopra and Meindl 2007) . Today, managers are well aware that in order to increase the productivity of their business and gain a competitive advantage, they need to optimally design their supply chain network (SCN) (Pal et al. 2016 (Pal et al. , 2021b Fallahpour et al. 2018; Nayeri et al. 2020; Diabat and Jebali 2021; Salehi-Amiri et al. 2021; Sazvar et al. 2021a) . Hence, the supply chain network design (SCND) problem is known as one of the important problems in today's competitive environment (Polat and Gungor 2021) .

In the traditional approach, researchers only considered the financial aspect of the SCND problem, but in recent years the concerns about the environmental impacts of the SC are increased which has led to attracting the attention of researchers to design the environmental-friendly SCN, which so-called green supply chain (GSC) Chen and Ye 2021; Fallahpour et al. 2021a, b; Hong et al. 2021; Usmani et al. 2021; Waqas et al. 2021a ). On the other side, increasing environmental pollutions has led to passing some regulations and rules that oblige companies to comply with environmental considerations. In general, minimizing the greenhouse gas emissions is one of the main concerns of the green view (Dogan and Seker 2016; Dogan and Turkekul 2016; Dogan and Inglesi-Lotz 2017; Alamroshan et al. 2021) . Although at the first glance considering the environmental issues may lead to increasing the total costs of the SC, it can increase the profits of the SC by providing a green image from the company and retention of customers for the long term. There are some aspects to consider greenness in the SCN problem such as reducing carbon emissions, reducing energy consumption, and using environmental-friendly raw materials (Pishvaee et al. 2014; Waltho et al. 2018; Gholizadeh and Fazlollahtabar 2020; Fallahpour et al. 2021b; Nayeri et al. 2021b) .

One of the concepts that has received more attention in recent years is reverse logistics (RL) (Govindan and Gholizadeh 2021) . The ability of RL to utilize the end-of-use (EOU) or end-of-life (EOL) products can reduce environmental pollutions, drastically (Afra and Behnamian 2021) . In general, when the forward SC and the reverse SC are integrated, the traditional SC converts to the closed-loop one (Rahimi and Fazlollahtabar 2018; Gholizadeh and Fazlollahtabar 2020) . In a closed-loop supply chain network (CLSCN), the demands of customers are satisfied by the forward SC, whereas the EOU and EOL products collected from consumers are handled by RL (Gholizadeh et al. 2020b; Sazvar et al. 2021b ). The efficient configuring of the CLSCN can lead to decreasing the environmental impacts of the SC by (I) utilizing the recycled components and materials to manufacturing products and (II) applying environmentally friendly transportation systems (Amin and Zhang 2013) .

Since customers are the most important element of each SC, satisfying their demands is considered as the main target of the SCN. In this regard, an important criterion has been introduced in the literature named responsiveness. In general, researchers have described responsiveness as the ability of an SC to meet the demands of customers over the planning horizon (Rabbani et al. 2018; Nayeri et al. 2021b) . For better understanding the importance of the responsiveness, see the following example. The market share of the Motorola company was decreased from 60% in 1994 to 16% in 2002 due to company's inability to respond to fluctuations and change customer demand (Roh et al. 2014) . This example demonstrates that ignoring the responsiveness measure can lead to irreparable financial loss. Hence, considering this criterion in the SC problem seems to be necessary. Several features such as increasing the service level, reducing the lead times, and providing various products are defined to incorporate responsiveness aspect into the SCN problem (Pishvaee and Rabbani 2011; Rabbani et al. 2018; Sabouhi et al. 2020; Nayeri et al. 2021a) .

The importance of medical devices has been drastically highlighted during the recent pandemic . Based on World Health Organization (WHO) , there are about 2,000,000 types in the market. Among them, one of the most widely used devices in the pandemic period is the ventilator. A ventilator is a device that allows air to enter the lungs and allows it to be expelled from the lungs. The ventilator performs breathing work for patients who have temporary or permanent respiratory problems. Motivated by the importance of the mentioned goods during the COVID-19 period, this research selected this product as a case study.

Based on the above-mentioned points, this study attempts to investigate the green-responsive CLSCN (GRCLSCN) design problem. For this purpose, a multi-objective mixed integer programming (MOMIP) model is offered to minimize the environmental damages and total costs while maximizing the responsiveness of the SCN. Since uncertainty is known as one of the major issues of the SCND problem, this study applies the robust possibilistic programming (RPP) method to tackle the uncertainty in the model's parameters. Afterwards, considering a case study, the offered MOMIP is solved by combining the meta-goal programming (MGP) method with a heuristic method. The framework of the current research has been depicted in Fig. 1 . The main advantages of the current paper comparing to the previous studies are (I) simultaneously consideration of the green and responsive features in the CLSCN design problem, (II) developing an efficient hybrid method to solve the proposed MOMIP, and (III) investigating the supply chain of the ventilator product as one of the most important medical devices in the recent pandemic . Therefore, improvement and novelty of the paper from the existing literature can be presented in the two aspects: (i) this is the first study that proposed a multi-objective robust possibilistic programming model to design a green-responsive supply chain network for the medical devices industry, and (ii) the current study developed an efficient hybrid solution approach based on the meta-goal programming method and a heuristic algorithm to solve the research problem. Based on the above discussions, the main innovation of the current research are as follows:

-This study simultaneously considers the greenness and responsiveness aspects in the closed-loop supply chain that never occurred in the previous research. In the green aspect, governments have paid special attention to the environmental factors that oblige the firms to comply with environmental issues (Breque et al. 2021 ). On the other side, considering the responsiveness aspect leads to gaining the competitive advantage for the firm that can increase the total profits (Nayeri et al. 2021b ).

-Since the closed-loop supply chain network design problem is an NP-hard one, it has a high level of complexity and needs a high CPU time for solving (Tiwari et al. 2016; Çağrı et al. 2018) . In this regard, to overcome this issue, this research aims at developing an efficient hybrid solution method based on the meta-goal programming, which is applied to convert the multi-objective model to a single one, and a heuristic algorithm, which is employed to tackle the complexity of the proposed model. -Since in the recent pandemic the importance of the ventilator device has become more highlighted and the demand for this product has drastically increased, this study selects this important product as a case study.

The current study aims at answering to the following research questions: 1) How to design a green and responsive closed-loop supply chain network using a multi-objective robust possibilistic programming model? 2) How to develop an efficient hybrid solution to solve the problem? 3) What is the behavior of the research problem according to a change in the crucial parameters such as demand? 4) How is the behavior of the proposed model in terms of robustness?

In the rest of this paper, the literature review is reported in "Literature review" section. "Problem definition" section presents the mathematical model. The FRO counterpart are developed in "Robust possibilistic programming counterpart" section. The hybrid solution method is explained in "Solving method" section. Case study and obtained results Ratio of the met demand over the potential demand Fig. 1 The framework of this study are provided in "Computational results" section. Finally, conclusions are presented in "Conclusions" section.

Recently, several study papers studied the SC management problem due to the importance of it (for example, see Alizadeh et al. 2020; De and Mahata 2020, 2021; Pourmehdi et al. 2020; De 2021) . In this research, according to the nature of this study, this section is dedicated to reporting some of the most related previous papers.

In the two last decades, regarding the importance of environmental aspects of the SC, several studies have been conducted to investigate the GSC design problem. In this field, Zohal and Soleimani (2016) developed an ant colony to configure a green CLSC (GCLSC) for the gold industry. In this regard, they proposed a MOMIP model to minimize the environmental damages and the total costs. The achieved outcomes indicated the efficiency of the offered model and the proposed method. Soleimani et al. (2017) suggested a MOMIP to design a GCLSC. Then, the authors applied fuzzy programming to tackle uncertainty in the model's parameters. Afterwards, they employed the genetic algorithm (GA) to handle the research problem. The results indicated the validity and performance of the suggested model and the offered algorithm. Rad and Nahavandi (2018) configured a GCLSC network considering the quantity discounts by proposing a MOMIP model. Their model aimed to minimize the environmental impacts and the total costs while maximizing customer satisfaction. The authors applied the LP-metric method to solve the research problem. Su et al. (2019) examined the GCLSC problem considering thirdparty recycling. To do this, the authors presented a two-stage model, and then compared the impact of decentralized decision-making and centralized decision-making on the pricing and returns strategies of each participant. Zhen et al. (2019) suggested a bi-objective programming model (BOPM) for configuring a GCLSC in order to create a trade-off between CO2 emissions and operating cost. They utilized the stochastic approach to tackle the uncertainty that existed in the model's parameters. Then, to solve the proposed model, the authors employed the Lagrangian relaxation method. Mardan et al. (2019) offered a BOPM to configure a multiperiod GCLSC for the wire-and-cable industry. The authors combined the LP-metric and the Benders decomposition methods to solve their model. Guo et al. (2020) offered a non-linear programming model to configure a GCLSC considering offline and online sales modes. Due to the complexity of the proposed model, the author employed the particle swarm optimization (PSO) and the GA algorithms to solve the research problem. Huang et al. (2020) designed a GCLSC network considering carbon emissions and cost control. They proposed a BOPM to minimize environmental damages and the total costs. Also, to deal with the uncertainty, the authors employed the stochastic programming method. Finally, the research problem was solved using the -constraint method. Mohtashami et al. (2020) combined the mathematical model and queuing system to design a GCLSC considering energy consumption. They considered different types of returned products that shipped from the customers to collection centers. Also, the authors applied both exact and meta-heuristics methods to solve the research problem. Sadeghi et al. (2020) suggested a MOMIP to design a GCLSC with the location-routing problem. It should be noted that the authors selected an automotive spare parts manufacturing company as a case study. Li et al. (2021) tried to design a CLSCN with consumer preference for green design. The authors emphasize the impacts of the consumer preference for green product design and the mandated targets on the strategies of the firms. Jian et al. (2021) studied the GCLSC problem under consideration of the manufacturer's fairness concern. Their model configures a GCLSC network with profit-sharing contract coordination fairness. The obtained results showed that the profit-sharing contract could improve the performance of the studied SC. Liu et al. (2021) offered a novel robust fuzzy model to configure a GCLSC with distribution ambiguity. The authors selected Coco-Cola Company as a case study to show the validity of their study. Mandal and Pal (2021) studied a game-theoretic method in a competitive SC considering a two-tier credit facility and imperfect production. It should be noted that the authors investigated the research problem under the environment of carbon emissions in which a manufacturer and two rival retailers compete with each other. Pal and Sarkar (2021) studied the strategies of a dual-channel GSC considering recycling under retailer promotional effort. In this regard, the authors employed the game theory approach to determine the optimal strategies. Ghomi-Avili et al. (2021) offered a bi-level mathematical model to design a multi-period and singleproduct SC network under the competitive environment. The authors employed the robust optimization approach to tackle uncertainty and used the Karush-Kuhn-Tucker (K-K-T) optimality conditions to solve the research problem. Sherif et al. (2021) designed a two echelon SCN by proposing a mathematical programming model for the battery manufacturing industry. The authors employed the simulated annealing algorithm to solve the research problem in the large-sized test problems.

As mentioned, responsiveness implies the capability of the SC to satisfy the demand of the customers. In recent years, several works tried to investigate the responsive SCND problem, which we report some of the most important of them in this section. Pishvaee and Rabbani (2011) proposed two MIP models to configure a responsive SCN considering direct shipment and prohibited one. They developed a graph theoretic heuristic method to handle the models. Baghalian et al. (2013) suggested a robust model to design a multi-product SC against disruptions under uncertainty in both sides of supply and demand. They used a path-based formulation to consider uncertainties in the supply side. A mathematical model is offered by Dubey et al. (2015) for designing an SCN with sustainability and responsiveness concepts. Their model aimed to minimize the total lead-time of collection and delivery activities and the total costs. They applied several types of robust optimization to tackle the uncertainty of the model's parameters. Martí et al. (2015) offered a BOPM to design a responsive SCN considering carbon policy and transportation modes. The main goal of their study was to create a balance between the responsiveness of SC and the carbon footprint. It should be noted that the authors considered the shortage of products and the transportation time as the responsiveness level. A robust non-linear programming model was offered by Mohammaddust et al. (2017) to design a SCN with responsiveness feature under the uncertain environment. They tried to minimize the lead-time of SC as the responsiveness measures. Rabbani et al. (2018) developed a multi-layer responsive SCN. The authors considered the lateral transshipments strategy in the research problem as the responsiveness criterion. On the other hand, they tried to minimize the total delivery time in the objective function. A scenario-based programming model is offered by Sabouhi et al. (2020) to design a responsive and resilient SCN in which the service level was defined as the responsiveness measure. The authors employed a multi-cut L-shaped method to solve the research problem. Nayeri et al. (2021a, b) examined the inventory-location problem in the responsive SCN under uncertainty. They suggested a mathematical model and solved it using several metaheuristic algorithms. Also, to demonstrate the performance of the developed algorithm, the authors solved several numerical examples. Aboolian et al. (2021) investigated a responsive make-toorder SCN. The authors considered the demand and service time as uncertain parameters and applied the stochastic programming method to tackle uncertainty. Finally, they solve the proposed model employing a metaheuristic algorithm. Azaron et al. (2021) configured a responsive SCN by proposing a stochastic multi-stage mathematical model. Their model aimed to maximize the total profit and minimize the traveling time as the responsiveness measure.

In recent years, due to the special attention of governments to environmental issues, managers of the firms and companies have been obliged to adopting environmental-friendly strategies for their businesses (Nayeri et al. 2020; Breque et al. 2021) . In this regard, several studies have been conducted to investigate the business environmental strategies. For example, Banerjee (2001) discussed about corporate environmental strategies and actions. They studied the environmental strategies of over 250 US firms. Handfield et al. (2005) interviewed with firms from several countries to develop a framework for environmental SC strategy decision-making. They found that some goals such as recycling and using environmental-friendly raw materials were on the list of agenda of these firms. Zhulavskyi et al. (2017) studied the environmental responsibility strategies of business. They provided the main criteria to identify the strategies for environmental responsibility of companies. Pedram et al. (2017) suggested the strategy of using new and repaired products with different prices that led to using the returned product and reducing environmental damages. Rahimi and Fazlollahtabar (2018) adopted the strategy of using EOL and EOU products to control reverse flow, reduce wastes, and decrease environmental damages. Nayeri et al. (2020) considered the manufacturing and recycling technologies with lower environmental damages but higher setup costs in order to reduce the carbon emissions generated in their supply chain. Mohammed et al. (2017) studied their problem under different carbon emissions policies (e.g., tax and cap-and-trade policies) to control the carbon emission of their SCN. Homayouni et al. (2021) tried to utilize vehicles with lower carbon emission that led to reducing environmental damages of their SCN. Nayeri et al. (2021b) attempted to employ welding technology with lower emissions and also transportation modes with lower pollution to decrease the environmental impacts of their business. Waqas et al. (2018a, b) investigated the crucial barriers for implementing the reverse supply chain in the manufacturing industries. In this way, the authors developed a two-stage method based on the structural equation modeling and the Delphi method and employed SPSS and AMOS softwares. Waqas et al. (2021a, b ) mapped a roadmap towards environmental performance, green innovation, and competitive advantage utilizing big data analytics. The main contribution of the mentioned work was introducing possible moderators and mediators among the relationship of big data analytics and environmental performance. Waqas et al. (2020) tried to model the reverse supply chain barriers in the manufacturing industry using interpretive structural modeling and MICMAC methods. Based on the obtained results, financial, technological, and organizational barriers were identified as dependent barriers. Waqas et al. (2018b) tried to explore the consumer willingness to select environmentally friendly options and also their attitude towards sustainable transportation in China using a norm activation model. The results demonstrated that self-enhancement and self-transcendence have negative and positive moderating impacts on the acceptability of sustainable transportation.

In today's competitive market space, the majority of the firms are struggling with high operating costs. In this regard, employing cost reduction strategies can help these firms to increase their profits, dramatically. Cost reduction strategies are effective methods or principles to enhancement of efficiency of operations. Cost reduction strategies can decrease the total costs of operations while enhancing productivity and allowing reallocation of resources (Bragg 2010) . One of the widely used cost reduction strategies is the supply chain management (Bragg 2010) . In this regard, one of the most important goals of the supply chain management is to minimize the total costs of SC. For example, by reducing the production costs, the company can refocus budgeted resources on expanding operations. This helps companies to align the goals and also increasing market share. In this way, for example, Chen and Wang (2012) studied the SC problem considering budget constraint in which the total costs of SC could not excess a certain value (budget). Nagurney et al. (2017) adopted the limited budget strategy in order to reduce the total costs of their supply chain. (Ha et al. 2017 ) utilized the information sharing strategy to reduce to costs of their business. Nayeri et al. (2021a) employed Markovian-based decision-making approach to reduce the inventory and shortage costs in their business.

As can be seen in the literature, although many papers investigated the green or the responsive SCND problems, there is a lack of studies that simultaneously considered green and responsive features especially considering reverse logistics. Indeed, no study developed a green and responsive closedloop supply chain network. However, as aforementioned, all of the mentioned concepts (i.e., greenness, responsiveness, and CLSC) are very important and can help the firms to gain a competitive advantage. In this regard, to cover these gaps, this study offers a MOMIP model for designing a GRCLSCN, which minimizes the environmental damages and the costs and maximizes the responsiveness of the SCN. Then, a robust possibilistic method is utilized to deal with the uncertainty. Finally, as the CLSC network design problem is considered as an NP-hard one (Tiwari et al. 2016; Çağrı et al. 2018) , a heuristic method is combined with the meta-goal programming approach to solve the MOMIP. Based on the mentioned points, the main contributions of the current research are listed as below:

-To the best of our knowledge, this is the first research that proposes a multi-objective robust possibilistic programming model to design a green-responsive closed-loop supply chain network design problem. -This research develops an efficient solution approach based on a heuristic method and the meta-goal programming method (HMGP) to solve the proposed model in a reasonable time. -This paper investigates the ventilator device as a case study due to its importance during the recent pandemic .

In this section, we describe the research problem and formulate the MOMIP. The CLSCN is that designed in this research includes four facilities in the forward chain involving suppliers, manufacturing sites (MSs), distribution centers (DCs), and customers. Also, it has three echelons in the reverse chain involving collection centers (CCs), remanufacturing centers (RCs), and disposal centers. The designed CLSCN is illustrated in Fig. 2 . The material flows in this network are as follows. At the outset, the raw materials are purchased from the suppliers and then shipped to the MSs. Afterwards, the original goods are manufactured in the MSs and shipped to the DCs for shipping to the demand points. Then, in the reverse logistics, the EOU/EOL products are collected from customers and sent to the CCs. In the CCs, the collected products are inspected and the percentage of them that capable to repair are shipped to the RCs. The remaining of them (that cannot be repaired) are sent to the disposal center. This study considers two types of products, the first one is the original product that is produced in the MSs, and the second one is the repaired product that is sent from the RCs to the DCs and then shipped to demand points. Also, this study considers different types of manufacturing technologies and transportation modes that each of which has different costs, carbon emissions, and capacity. In this research, to incorporate the green concept in the problem, two strategies are considered: (i) considering the revers logistic that leads to re-use the EOU/ EOL products and (ii) minimizing the carbon emissions due to production, remanufacturing, and transportation activities.

This section is devoted to present the MOMIP. In the following, the necessary notations of the MOMIP are presented. There are three main reasons for proposing the MOMIP: firstly, the reason for proposing a multi-objective model is that there are some conflicting targets (i.e., minimizing total costs, minimizing environmental impacts, and maximizing responsiveness) in the research problem. In general, when we have several conflicting objective functions, we utilize multi-objective models. On the other side, since in our problem there are different types of variables (continuous, binary, and integer), we would rather propose a mixed integer programming model. Also, there are several papers that proposed MOMIP to design SCNs and obtained good results (for example, see Paydar et al. 2017; Mamashli and Javadian 2020; Nayeri et al. 2020 Nayeri et al. , 2021b Pourmehdi et al. 2020 ) that shows the appropriate performance of the MOMIP models in the SCN design problem. Shortage of the original product p for customer k SH2 dk Shortage of the repaired product d for customer k Based on the above definitions, the proposed MOMIP can be formulated as follows. The first objective function (OF) (1) minimizes the total costs. In OF1, the first part is the fixed-cost (FC) of contracting with suppliers, the second part shows FC of opening MSs, the third part denotes the FC for DCs, part 4 represents the FC of CCs, and part 5 is the FC of RCs. On the other side, part 6 shows the purchasing costs of raw materials, part 7 denotes the manufacturing costs, part 8 is distribution costs, part 9 shows the collection costs, part 10 represents the remanufacturing costs, and parts 11 and 12 are the shortage costs. Eventually, the last part shows the transportation costs.

The second OF minimizes the total emissions, in which the first part shows production activities emissions, the second part is remanufacturing activities emissions, and the last part is related to transportation process emissions.

The third OF (relation (3)) maximizes the responsiveness measure. This measure has been defined based on the ratio of satisfied demand over the potential demand. In this regard, the sum of the quantity of the original and repaired products that are shipped to customers is divided by the total potential demand. This equation is inspired by Fattahi et al. (2017) and Sabouhi et al. (2020) , who calculated the responsiveness level of their supply chain based on the ratio of the met demand over the total demand. Constraint (4) calculates the quantity of raw materials that should be purchased from the suppliers, and relation (5) is the suppliers' capacity constraint.

Constraint (6) shows that the quantity of products shipped from MSs to DCs is equal to the quantity of products sent from DCs to demand points. Relation (7) Constraint (9) determines the quantity of the original product sent to customers and the quantity of the shortage for the original product. Also, constraint (10) measures the amount of repaired products shipped to the customers and the amount of the shortage for the repaired product. Relation (11) is the capacity constraint of DCs.

Constraint (12) shows the quantity of returned products shipped to CCs. On the other side, relation (13) demonstrates the capacity constraint for CCs. Equation (14) shows the quantity of products shipped from CCs to RCs. Relation (15) calculates the quantity of products sent from CCs to disposal centers.

Equations (16) and (17) 

Since the supply chain network design problem involves different decision levels especially long-term decisions, uncertainty has been considered as a big issue in this problem. Ignoring the uncertainty in the design phase may lead to achieving unrealistic results that can result in irreparable financial losses (Gholizadeh et al. 2020a) . Usually, in the SCND problem, some historical data are available but there is a lack of knowledge on input data that can be gathered from the experts (Nayeri et al. 2020) . The mentioned condition is called epistemic uncertainty that known as a common type of uncertainty in the SCND problem (Bairamzadeh et al. 2018) .

In this regard, one of the most powerful methods to tackle epistemic uncertainty is the robust possibilistic programming (RPP) method that widely used in previous related studies and showed good performance (see Pishvaee et al. 2012b, a; Hamidieh et al. 2017; Ghaderi et al. 2018; Nayeri et al. 2020; Rabbani et al. 2020; Habib et al. 2021) . Based on the mentioned reasons, we have employed the RPP approach to tackle uncertainty in this study. To provide more explanations about the RPP method, consider the following compact model.

In model (39), c , d , and d ′ are the trapezoidal fuzzy numbers with four critical point (i.e., c = (c 1 , c 2 , c 3 , c 4 ) ); f , A,A ′ , B, s , and N are the deterministic parameters. y is the binary variable and x is the continuous variable. The RPP counterpart for model (39) can be formulated as follows:

In model (40), the first and second sections of the objective function showed the expected value (average) for of the costs ( E[Z] ). The third section shows the optimality

robustness that has been incorporated by penalizing the difference between the worst case from the average. In this section, shows the optimality robustness penalty cost and Z max is the worst case that is calculated according to Eq. (41). The last section of the objective function demonstrates the feasibility robustness. In this regard, i represents the feasibility robustness penalty costs and m shows the satisfaction level.

Based on the above-mentioned points, the RPP counterpart model for the research problem can be formulated as follows. 

In this section, the solution procedure of this study is presented. To solve the proposed MOMIP, we have two main challenges: (a) the proposed model is a multi-objective one and (b) the research problem is an NP-hard (Tiwari et al. 2016; Asadi-Gangraj and Nayeri 2018; Çağrı et al. 2018) .

In this regard, to cope with these challenges, we develop an efficient hybrid method. In this way, we combine the metagoal programming method (to transform the MOMIP model to a single-objective model) with a heuristic method (to cope with the complexity of the problem). In the following, we briefly describe the mentioned approaches.

One of the popular approaches to solve multi-objective problems is goal programming (GP) Mamashli et al. 2021b; Razavi et al. 2021 ).

In the field of the SCND problem, different variants of the goal programming method is applied and they achieved good results (see Hocine et al. 2020; Mamashli and Javadian 2020; Nayeri et al. 2020; Pourmehdi et al. 2020; Homayouni et al. 2021 ; Jamali et al. Mamashli et al. 2021a ). In the current study, the meta-goal programming (MGP) method offered by Urıá et al. (2002) is applied to solve the proposed MOMIP. The main advantages of this method are its ability to obtain more balanced solutions and its flexibility to model the preferences of the decision-makers (DMs) comparing to the other variants of GP (Benítez-Fernández and Ruiz 2019). Suppose t i is the goal value for i-th objective function that determined by the decision-maker (DM), n i shows negative deviation and p i is positive deviation. Now consider expression (47) that denoted algebraic structure of the generic i-th goal. Also, there are three types of deviations which interested readers can be refer to (Urıá et al. 2002 ; Benítez-Fernández and Ruiz 2019) for more study.

Because the goals of this research are "less is better," the unwanted deviation is positive ones (Urıá et al. 2002) . Urıá et al. (2002) defined three types of meta-goals as follows:

• Type 1: As shown in relation (48), there is a certain bound ( Q 1 ) that the sum of unwanted percentage deviation variables cannot surpass it.

It should be noted that i is a preferential weight and to normalizing the deviation variables, them dividing among their corresponding target value. • Type 2: As shown in relation (49), there is a certain bound ( Q 2 ) that the maximum percentage deviation should not be greater than it.

where D shows the maximum weighted percentage deviation.

• Type 3: As shown in relation (50), there is a certain bound ( Q 3 ) that the percentage of unachieved goals should not be greater than it.where G i is a binary variable and M i is a large number. Due to above definitions, the MGP is modeled as follows:

(47)

The objective function of (51) can in turn take the weighted or the MinMax form which in this paper we consider weighted form.

Since the research problem is known as an NP-hard one, this study employs an efficient heuristic method that widely used in the SCND problem (see Kaur and Singh 2018; Homayouni et al. 2021; Shafiee et al. 2021) to handle the complexity of the research problem. The steps of this algorithm are listed below (Kaur and Singh 2018):

First step: Relax the proposed model from the binary constraints (convert the binary variables to continuous ones). Second step: Solve the relaxed model. Third step: Record all non-zero quantities of relaxedvariables which are obtained by solving the proposed model. Fourth step: Equalize each non-zero quantity of relaxedvariables to 1, and set them as constraints to the main model. Fifth step: Solve the created model optimally.

The pseudo-code of the developed hybrid method is depicted in Fig. 3 .

The COVID-19 pandemic has led to a significant disruption in the world. This virus, which has negative effects on the human respiratory system, has led to many deaths in the past 2 years. After the recent pandemic (COVID-19), the demand for medical devices has dramatically increased. One of the most widely used medical devices during the recent pandemic is

the ventilator device. On the other side, the logistics process of this product is adapted with the proposed CLSC in the current study. In this regard, in many cases when this product becomes EOU or EOL, it is collected and remanufactured for re-using. Hence, this research has selected the ventilator device as a case study. The main purpose of using a ventilator is to establish mechanical ventilation and support the lungs to supply oxygen to the patient. In fact, the ventilator does not cure the disease, but can only keep the patient alive to cure the underlying disease. Figure 4 depicts the ventilator and its application. The input data for solving the proposed MOMIP can be found in the Supplementary materials. It should be noted that the required data are gathered from the related literature (such as Hasani et al. 2015; Pedram et al. 2017; Nayeri et al. 2020and the opinion of the experts of the case study company).

In this section, the offered MOMIP are solved and results are presented. At the outset, to show the performance of the developed solution method, the problem is solved in different sizes and the achieved results by the developed hybrid method (HMGP) are compared with the traditional MGP (TMGP). To do this, two measures are utilized: (i) GAP and (ii) CPU time, which the gap is calculated based on the following relation.

where HMGP sol is the solution obtained by the proposed hybrid method and TMGP sol represents the solution obtained by the TMGP method. Table 1 shows the results of solving the proposed MOMIP by the two mentioned methods. As shown in Table 1 , the HMGP can achieve optimal/nearoptimal solutions on a better time than the TMGP. In this 

One of the most important parameters of the offered model is the demand (demand of original product and demand of repaired product). In this section, we have analyzed the behavior of the research problem according to changing the size of the demands. For this purpose, the problem has been solved with different values for the mentioned parameters and reported the achieved results in Fig. 7 . According to Fig. 7 when demand is increased, the environmental impact and the total costs are increased but the responsiveness of the SC (the third objective function) is decreased. In this regard, a 30% increase in demand from its base case value leads to a 24% increase in the total costs, a 10% increase in carbon emissions, and an 8% decrease in responsiveness. In addition, by 30% decreasing the demand size from its base case, Z1 and Z2 are respectively decreased by 22% and 8%, and Z3 is increased by 7%.

In this section, we have conducted a sensitivity analysis to examine the role of the capacity parameter. Figure 8 illustrates the behavior of the research problem according to changing the capacity parameter. As can be seen in this figure, when the capacity parameter is increased, the first objective function is decreased, but the second and third objective functions are increased. In this regard, by shifting the value of capacity from − 30 to + 30%, the total cost and environmental impacts are decreased by 52% and 12%, respectively. Also, this change leads to increasing the responsiveness of the supply chain by about 13%.

This section is devoted to analyzing the behavior of the studied system according to changing the rate of return of the products parameter. In this regard, the problem is solved under different values for the mentioned parameter and the results are depicted in Fig. 9 . Based on this figure, an increase in the rate of return leads to an increase in all of the objective functions.

In this way, when the rate of return is increases from 0.3 to 0.7, Z1 is increased by 12%, Z2 is increased by 10%, and Z3 is increased by 8%.

Here, the effect of changing parameter on the offered MOMIP is analyzed. The achieved results from sensitivity analysis have been presented in Table 2 . According to the results, changes in negligibly affect the objective functions. The minimum satisfaction level remains constant and equal to 1, as well.

This section is dedicated to investigating the impact of optimality robustness penalty cost on the results. Table 3 shows the result of sensitivity analysis. As can be seen in Table 3 , when increases from 0.3 to 0.9, the model tries to get values closer than to the Z max . In general, by increasing the Fig. 5 Comparison of the TMGP and HMFP according to solutions' quality. a The first objective function. b The second objective function. c The third objective function value of , the total cost is increased, too. This behavior is occurred due to a phenomenon, namely, "robustness price." It should be noted that the robustness price shows the costs that entered in the problem to face the uncertainties.

In this section, we try to provide some useful managerial insights obtained from this study. These implications have been listed below.

-This study has focused on the CLSCN configure problem with environmental and responsiveness features. Hence, this research can give a good view to the leaders of the SC to manage their business according to the two mentioned features. Also, this study helps the managers to see the impact of the uncertainty on the SC management problem. -Based on Fig. 7 , by increasing the demand parameter, the first objective function (costs) is increased, too. The main reasons for this behavior are that increasing the demand results in increasing the establishment, production, distribution, shortage, and transportation costs. One of the efficient strategies to cope with this condition is sub-contracting that can prevent the enhancement of the cost. - Fig. 7 showed that an increase in demand has resulted in increasing environmental damages. In order to increase the greenness of the SC, utilizing the TMs and technologies with fewer carbon emissions can be a good action. Although these technologies and transportation modes have higher setup costs and lead to enhancement of the costs, the leaders must know that using greener facilities leads to providing a green image from the company in the eyes of customers, which results in increasing the loyalty of customers. - Fig. 7 indicates that when the demand size is increased, the responsiveness of the SC is decreased. This point demonstrates that the shortage of products is increased by increasing the demand size. To tackle this condition, sub-contracting may be a good strategy that can manage the capacity temporarily and increase the met demand. -According to Fig. 9 , a decrease in the rate of return parameter results in decreasing the responsiveness of the SC. In this regard, considering incentive plans for customers to return the goods can be a good strategy for Fig. 9 The behavior of the research problem according to changing the rate of return parameter increasing the rate of return that leads to increasing the responsiveness of the supply chain.

The main theoretical contribution of the current research is combining the green and responsive aspects in the CLSCN that was ignored by previous studies. In this regard, we have designed a CLSC network and considered the carbon emissions as an environmental measure and the ratio of met demand over total demand as the responsiveness measure.

On the other side, this paper has investigated the problem under uncertainty and employed the robust possibilistic programming method to tackle uncertainty. To the best of our knowledge, this is the first paper that has proposed a multiobjective robust possibilistic model to design a GRCLSCN.

Comparing to the similar approach to tackle the epistemic uncertainty (e.g., fuzzy programming), the employed approach (robust possibilistic programming) has two main advantages as follows: (i) this method considers the riskaversion level of output decisions that can help decisionmakers , and (ii) this method has capability to finding optimal value for the satisfaction levels (Talaei et al. 2016 ). On the other hand, this study has developed an efficient solution method by integrating a heuristic method and the meta-goal programming method. Comparing to the other versions of the goal programming, metagoal programming, which has been applied in this study, has several advantages such as its ability to obtain more balanced solutions and its flexibility to model the preferences of the DMs (Nayeri et al. 2021b) . Also, the heuristic method that utilized in this research is an efficient one to solve the complex problems with binary variables (Kaur and Singh 2018; Homayouni et al. 2021; Mamashli et al. 2021a; Shafiee et al. 2021) . Hence, this hybrid approach can solve the multi-objective complex models, efficiently. Eventually, another theoretical contribution of this research is to investigate the supply chain of the ventilator as one of the most important medical devices during the COVID-19 pandemic.

In this section, we attempt to compare this work with the previous studies. On the other hand, regarding the obtained results, there are some similar behaviors between this research and the previous studies, which demonstrate the validation of the offered model. In this way, for instance, the results of this research showed that when the demand increases, the environmental damages and the total cost increase, too. This behavior was also observed in the previous related studies ( Sahebjamnia et al. 2018; Nayeri et al. 2020 Nayeri et al. , 2021b Pourmehdi et al. 2020; Sazvar et al. 2021a ). On the other side, the achieved results demonstrated that the responsiveness of the supply chain has decreased by enhancing the demand. This behavior was also occurred in Fattahi et al. (2017) , Sabouhi et al. (2020) , and Nayeri et al. (2021b, a) that verifies the achieved results of this work. Also, based on results, by increasing the rate of return, the environmental impacts and the total costs have increased. Similar behavior was observed in Pedram et al. (2017) and Nayeri et al. (2020) , which shows the validation of the proposed model.

This study addressed the GRCLSCN design problem under uncertainty. In this way, a MOMIP model is offered that minimized the environmental impacts and the total costs and also maximized the responsiveness of the supply chain. To tackle the uncertainty that existed in the nature of the supply chain problem, this research applied the robust possibilistic programming method. Then, a hybrid solution approach based on a heuristic method and the MGP method was developed to solve the proposed model. Since the demand for ventilator device is extremely increased during the COVID-19 pandemic, this study selected this product as a case study. Afterwards, the offered MOMIP was solved using the developed hybrid method and results are reported. The results demonstrated the performance and efficiency of the offered model and developed algorithm. Based on the results of the sensitivity analyses, by increasing demand the total costs and the environmental impacts were increased, but responsiveness was decreased. On the other side, an increase in the capacity of the facilities resulted in decreasing the total costs and environmental impacts while leading to increasing the responsiveness. One of the research deficiencies of this study is ignoring the global factors. However, since the top firms usually work at global scale, the future studies can add the global factors such as taxes and tariffs to the current model. Another research deficiency of the current paper is ignoring disruptions. In this regard, the future studies can add the resilience strategies such as backup suppliers to the current problem. Also, future studies can develop efficient metaheuristic algorithms to solve the suggested model.

The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s11356-022-19105-1.

Author contribution Zeinab Asadi: conceptualization, methodology, software, validation, original draft preparation, visualization. Mohammad Valipour Khatir: supervision, investigation, investigation, validation, writing-reviewing and editing.

Mojtaba Rahimi: conceptualization, methodology, writing-original draft preparation, visualization.

Responsive make-to-order supply chain network design

Lagrangian heuristic algorithm for green multi-product production routing problem with reverse logistics and remanufacturing

The green-agile supplier selection problem for the medical devices: a hybrid fuzzy decision-making approach

Forward and reverse supply chain network design for consumer medical supplies considering biological risk

A multi-objective facility location model for closed-loop supply chain network under uncertain demand and return

A hybrid approach based on LP metric method and genetic algorithm for the vehicle-routing problem with time windows, driver-specific times, and vehicles-specific capacities

Farhang Doost A (2021) Designing profitable and responsive supply chains under uncertainty

Robust supply chain network design with service level against disruptions and demand uncertainties: a real-life case

Modelling different types of uncertainty in biofuel supply network design and planning: a robust optimization approach

Corporate environmental strategies and actions. Manag Decis Benítez-Fernández A, Ruiz F (2019) A meta-goal programming approach to cardinal preferences aggregation in multicriteria problems

A robust two layer green supply chain modelling under performance based fuzzy game theoretic approach

Pollution sensitive global crude steel production-transportation model under the effect of corruption perception index

De Nul L, Petrides A (2021) Industry 5.0 -towards a sustainable, human-centric and resilient European industry

Forward supply chain network design problem: heuristic approaches

Influence of creating an oligopoly through government intervention to improve partner collaboration intentions in the context of green supply chains

Trade credit contract with limited liability in the supply chain with budget constraints

Supply chain management. Strategy, planning & operation

Solving an EOQ model under fuzzy reasoning

A production inventory supply chain model with partial backordering and disruption under triangular linguistic dense fuzzy lock set approach

Solution of an imperfect-quality EOQ model with backorder under fuzzy lock leadership game approach

Multi-product and multi-period closed loop supply chain network design under take-back legislation

Analyzing the effects of real income and biomass energy consumption on carbon dioxide (CO2) emissions: empirical evidence from the panel of biomass-consuming countries

Determinants of CO2 emissions in the European Union: the role of renewable and non-renewable energy

CO2 emissions, real output, energy consumption, trade, urbanization and financial development: testing the EKC hypothesis for the USA

The design of a responsive sustainable supply chain network under uncertainty

An intelligence-based model for supplier selection integrating data envelopment analysis and support vector machine

A hyper-hybrid fuzzy decision-making framework for the sustainable-resilient supplier selection problem: a case study of Malaysian palm oil industry

An integrated approach for a sustainable supplier selection based on Industry 4.0 concept

Responsive and resilient supply chain network design under operational and disruption risks with delivery lead-time sensitive customers

A robust fuzzy stochastic programming model for the design of a reliable green closed-loop supply chain network

A multi-objective robust possibilistic programming approach to sustainable switchgrassbased bioethanol supply chain network design

Robust optimization and modified genetic algorithm for a closed loop green supply chain under uncertainty: case study in melting industry

A robust fuzzy stochastic programming for sustainable procurement and logistics under hybrid uncertainty using big data

A closed-loop supply chain robust optimization for disposable appliances

Competitive green supply chain network design model considering inventory decisions under uncertainty: a real case of a filter company

Robust network design for sustainable-resilient reverse logistics network using big data: a case study of end-of-life vehicles

Research on green closed-loop supply chain with the consideration of double subsidy in e-commerce environment

Information sharing in competing supply chains with production cost reduction

A robust possibilistic programming approach toward animal fat-based biodiesel supply chain network design under uncertain environment

A robust possibilistic programming model for a responsive closed loop supply chain network design

Integrating environmental management and supply chain strategies

Robust closed-loop global supply chain network design under uncertainty: the case of the medical device industry

A fuzzy goal programming with interval target model and its application to the decision problem of renewable energy planning

A robustheuristic optimization approach to a green supply chain design with consideration of assorted vehicle types and carbon policies under uncertainty

Quality information acquisition and disclosure with green manufacturing in a closed-loop supply chain

A multi-objective stochastic programming model to configure a sustainable humanitarian logistics considering deprivation cost and patient severity

Decision-making and coordination of green closed-loop supply chain with fairness concern

Heuristic modeling for sustainable procurement and logistics in a supply chain using big data

Managing a closed-loop supply chain with take-back legislation and consumer preference for green design

A novel robust fuzzy mean-UPM model for green closed-loop supply chain network design under distribution ambiguity

A heuristicbased multi-choice goal programming for the stochastic sustainable-resilient routing-allocation problem in relief logistics

Sustainable design modifications municipal solid waste management network and better optimization for risk reduction analyses

Designing a sustainable-resilient disaster waste management system under hybrid uncertainty: a case study

Effects of green innovation and advertisement in an imperfect production-based competitive supply chain under two-tier credit facility

An accelerated benders decomposition algorithm for a bi-objective green closed loop supply chain network design problem

Carbon footprint and responsiveness trade-offs in supply chain network design

Developing lean and responsive supply chains: a robust model for alternative risk mitigation strategies in supply chain designs

Multi-period planning of closed-loop supply chain with carbon policies under uncertainty

A green closed loop supply chain design using queuing system for reducing environmental impact and energy consumption

A supply chain network game theory model of cybersecurity investments with nonlinear budget constraints

Goal programming-based post-disaster decision making for allocation and scheduling the rescue units in natural disaster with time-window

Multi-objective fuzzy robust optimization approach to sustainable closedloop supply chain network design

A robust fuzzy stochastic model for the responsive-resilient inventory-location problem: comparison of metaheuristic algorithms

A multi-objective fuzzy robust stochastic model for designing a sustainable-resilient-responsive supply chain network

Price, delivery time, and retail service sensitive dual-channel supply chain

Two-phase deteriorated supply chain model with variable demand and imperfect production process under two-stage credit financing

Two-echelon competitive integrated supply chain model with price and credit period dependent demand

Optimal strategies of a dual-channel green supply chain with recycling under retailer promotional effort

An engine oil closed-loop supply chain design considering collection risk

Integrated forward and reverse supply chain: a tire case study

A graph theoretic-based heuristic algorithm for responsive supply chain network design with direct and indirect shipment

Robust possibilistic programming for socially responsible supply chain network design: a new approach

An accelerated Benders decomposition algorithm for sustainable supply chain network design under uncertainty: a case study of medical needle and syringe supply chain

Correction to: WEEE closed-loop supply chain network management considering the damage levels of returned products

Scenario-based design of a steel sustainable closed-loop supply chain network considering production technology

Farrokhi-Asl H (2020) A hybrid robust possibilistic approach for a sustainable supply chain location-allocation network design

A graph theory-based algorithm for a multi-echelon multi-period responsive supply chain network design with lateral-transshipments

A novel multi-objective optimization model for integrated problem of green closed loop supply chain network design and quantity discount

Optimization of a closed loop green supply chain using particle swarm and genetic algorithms

A robust optimization model of the field hospitals in the sustainable blood supply chain in crisis logistics

Implementation of a responsive supply chain strategy in global complexity: the case of manufacturing firms

A multicut L-shaped method for resilient and responsive supply chain network design

A mixed integer linear programming model for designing a green closed-loop supply chain network considering location-routing problem

Sustainable tire closed-loop supply chain network design: hybrid metaheuristic algorithms for large-scale networks

Designing a sustainable closed-loop supply chain network for walnut industry

A capacity planning approach for sustainable-resilient supply chain network design under uncertainty: a case study of vaccine supply chain

Designing a sustainable closed-loop pharmaceutical supply chain in a competitive market considering demand uncertainty, manufacturer's brand and waste management

A robust multi-objective optimization model for inventory and production management with environmental and social consideration: a real case of dairy industry

Integrated optimization of transportation, inventory and vehicle routing with simultaneous pickup and delivery in two-echelon green supply chain network

Fuzzy multiobjective sustainable and green closed-loop supply chain network design

A green closed-loop supply chain coordination mechanism based on third-party recycling

A robust fuzzy optimization model for carbon-efficient closed-loop supply chain network design problem: a numerical illustration in electronics industry

A hybrid territory defined evolutionary algorithm approach for closed loop green supply chain network design

Meta-goal programming

Mapping green technologies literature published between 1995 and 2019: a scientometric review from the perspective of the manufacturing industry

Green supply chain network design: a review focused on policy adoption and emission quantification

Critical barriers to implementation of reverse logistics in the manufacturing industry: a case study of a developing country

Triggering sustainable firm performance, supply chain competitive advantage, and green innovation through lean, green, and agile supply chain practices

Big data analytics as a roadmap towards green innovation, competitive advantage and environmental performance

Modeling reverse logistics barriers in manufacturing industry of Pakistan: an ISM and MIC-MAC approach

Green and sustainable closed-loop supply chain network design under uncertainty

Social and environmental responsibility strategies of business

Developing an ant colony approach for green closed-loop supply chain network design: a case study in gold industry