Investigating a Dual-Channel Network in a Sustainable Closed-Loop Supply Chain Considering Energy Sources and Consumption Tax

This paper proposes a dual-channel network of a sustainable Closed-Loop Supply Chain (CLSC) for rice considering energy sources and consumption tax. A Mixed Integer Linear Programming (MILP) model is formulated for optimizing the total cost, the amount of pollutants, and the number of job opportunities created in the proposed supply chain network under the uncertainty of cost, supply, and demand. In addition, to deal with uncertainty, fuzzy logic is used. Moreover, four multi-objective metaheuristic algorithms are employed to solve the model, which include a novel multi-objective version of the recently proposed metaheuristic algorithm known as Multi-Objective Reptile Search Optimizer (MORSO), Multi-Objective Simulated Annealing (MOSA), Multi-Objective Particle Swarm Optimization (MOPSO), and Multi-Objective Grey Wolf (MOGWO). All the algorithms are evaluated using LP-metric in small sizes and their results and performance are compared based on criteria such as Max Spread (MS), Spread of Non-Dominance Solution (SNS), the number of Pareto solutions (NPS), Mean Ideal Distance (MID), and CPU time. In addition, to achieve better results, the parameters of all algorithms are tuned by the Taguchi method. The programmed model is implemented using a real case study in Iran to confirm its accuracy and efficiency. To further evaluate the current model, some key parameters are subject to sensitivity analysis. Empirical results indicate that MORSO performed very well and by constructing solar panel sites and producing energy out of rice waste up to 19% of electricity can be saved.


Introduction
In recent years, Agricultural Supply Chain (ASC) management has been the focus of many researchers due to its limited shelf life and the variety of demand and cost [1]. Moreover, the concept of sustainability in the logistics of these products with the goal to consider the environmental, economic, and social dimensions has been simultaneously set forth as another significant issue in the world [2,3]. Sustainable agriculture provides a potential solution to enable agricultural systems to feed a growing population within the changing environmental conditions [4]. The agricultural products industry involves the most practical and challenging applications in this field. Among such products, rice is a product consumed as the main part of meals in most parts of the world. Over half of the world population feeds on this basic grain, and this product plays a critical role in boosting the economies of the producing nations [8]. In 2018, global rice production was 782 million tons, 50% of which belonged to China and India, as depicted in Figures 2 and 3 (www.fao.org (accessed on 3 March 2022)). The agricultural products industry involves the most practical and challenging applications in this field. Among such products, rice is a product consumed as the main part of meals in most parts of the world. Over half of the world population feeds on this basic grain, and this product plays a critical role in boosting the economies of the producing nations [8]. In 2018, global rice production was 782 million tons, 50% of which belonged to China and India, as depicted in Figures 2 and 3 (www.fao.org (accessed on 3 March 2022)).
In addition to consuming rice as a main meal, the bran is widely utilized in the production of toiletry and pharmaceutical industries. Rice straw can also be used in the production of bio-energy and compost [9] (see Figures 4 and 5).
Most of the rice-producing countries suffer heavy and frequent losses in crop harvest because of poor maintenance of technologies, inefficient supply chains, and farmers' incapability to market the products. In research conducted by the World Bank and FAO, it has been claimed that, on average, around 8% to 26% of the rice produced in developing nations is destroyed due to poor facilities and post-harvest problems (www.fao.org (accessed on 3 March 2022)).  In addition to consuming rice as a main meal, the bran is widely utilized in the production of toiletry and pharmaceutical industries. Rice straw can also be used in the production of bio-energy and compost [9] (see Figures 4 and 5).       Most of the rice-producing countries suffer heavy and frequent losses in crop harvest because of poor maintenance of technologies, inefficient supply chains, and farmers' incapability to market the products. In research conducted by the World Bank and FAO, it has been claimed that, on average, around 8% to 26% of the rice produced in developing nations is destroyed due to poor facilities and post-harvest problems (www.fao.org (accessed on 3 March 2022)).
In addition, the rice-processing industry is heavily dependent on energy consumption. With environmental concerns such as greenhouse gas emissions, there is an urgent need to switch to non-carbon, renewable, and clean energy sources such as solar and biomass energy for rice processing [10]. Recent studies revealed that rice waste such as the husk and straw has great potential for electricity generation [8]. Instead of being disposed, rice waste that is properly used can bring about positive environmental effects (see Figures 6 and 7). Selecting an optimal combination out of traditional and renewable energy sources for rice production and processing can slightly lower some of the energy supply problems in this area [11]. In addition, the rice-processing industry is heavily dependent on energy consumption. With environmental concerns such as greenhouse gas emissions, there is an urgent need to switch to non-carbon, renewable, and clean energy sources such as solar and biomass energy for rice processing [10]. Recent studies revealed that rice waste such as the husk and straw has great potential for electricity generation [8]. Instead of being disposed, rice waste that is properly used can bring about positive environmental effects (see Figures 6 and 7). Selecting an optimal combination out of traditional and renewable energy sources for rice production and processing can slightly lower some of the energy supply problems in this area [11]. Most of the rice-producing countries suffer heavy and frequent losses in crop harvest because of poor maintenance of technologies, inefficient supply chains, and farmers' incapability to market the products. In research conducted by the World Bank and FAO, it has been claimed that, on average, around 8% to 26% of the rice produced in developing nations is destroyed due to poor facilities and post-harvest problems (www.fao.org (accessed on 3 March 2022)).
In addition, the rice-processing industry is heavily dependent on energy consumption. With environmental concerns such as greenhouse gas emissions, there is an urgent need to switch to non-carbon, renewable, and clean energy sources such as solar and biomass energy for rice processing [10]. Recent studies revealed that rice waste such as the husk and straw has great potential for electricity generation [8]. Instead of being disposed, rice waste that is properly used can bring about positive environmental effects (see Figures 6 and 7). Selecting an optimal combination out of traditional and renewable energy sources for rice production and processing can slightly lower some of the energy supply problems in this area [11].  Besides supplying energy sources, another important issue that decision-makers and managers may run into when designing a rice supply chain network is the uncertainty of decision making. Uncertainty influences the composition and coordination of the rice supply chain and may take place in several ways, such as the uncertainty in supply and demand. Under real-world conditions, rice supply and customer demand are some of the most common uncertain parameters due to the lack of sufficient information. The high uncertainty level in the supply chain brings about a challenging situation for predicting the future. Fuzzy logic is one of the approaches that is effectively modeling the parameter under uncertainty. The results indicate that fuzzy optimization Besides supplying energy sources, another important issue that decision-makers and managers may run into when designing a rice supply chain network is the uncertainty of decision making. Uncertainty influences the composition and coordination of the rice supply chain and may take place in several ways, such as the uncertainty in supply and demand. Under real-world conditions, rice supply and customer demand are some of the most common uncertain parameters due to the lack of sufficient information. The high uncertainty level in the supply chain brings about a challenging situation for predicting the future. Fuzzy logic is one of the approaches that is effectively modeling the parameter under uncertainty. The results indicate that fuzzy optimization is a proper option for decision making under uncertainty due to its potential to deal with uncertainty [7].
Pursuant to the mentioned subjects, in this study, optimizing the dual-channel CLSC network of rice is investigated while considering energy sources. Hence, a Mixed Integer Linear Programming (MILP) model is formulated for optimizing the total costs, the quantity of pollutants, and the number of job opportunities created in the proposed logistics network under uncertainty over the cost of purchase and supply and demand. The model aims to locate the proper sites for the construction of solar panels and bio-refineries and to determine the optimal product flow and its waste among the network facilities and the inventory level during each period while considering sustainability factors. Furthermore, fuzzy logic is applied for dealing with uncertainty. After that, four multi-objective metaheuristic algorithms, i.e., MORSO, MOSA, MOPSO, and MOGWO, are used to solve the model, and their results and performance are compared based on several criteria. Finally, the programmed model is implemented on a real case study in Iran in order to verify its accuracy and efficiency. The novel aspects reported in this paper are mainly as follows: Considering the energy sources are rarely observed in the ASC optimization models, in this paper, a sustainable double-channel CLSC network for rice product considering energy sources and tax on energy consumption is designed. Then, a new MILP model is formulated for optimizing the total costs, the quantity of pollutants, and the job opportunities created in the rice CLSC in a fuzzy environment. After that, to solve the model and find the Pareto solutions, a new multi-objective version of the recently released algorithm termed Multi-Objective Reptile Search Optimizer (MORSO) is employed to solve the proposed model in high dimensions. Next, its results and performance are compared to MOSA, MOGWO, and MOPSO based on criteria such as MS, SNS, MID, and CPU time.
The present paper is organized into six sections. In the continuation, the literature is reviewed in Section 2. The problem and the mathematical model are stated in Section 3. The solution approach is addressed in Section 4. The proposed model's accuracy is verified by implementing it onto a real case study in Iran in Section 5. Moreover, the metaheuristic algorithm parameters are tuned and the results of solving the model are analyzed and compared. For further analyzing the model, sensitivity analysis is run on the key parameters. The conclusions and suggestions for research development are presented in Section 6.

Literature Review
As stated before, producing rice is both economically and socially critical for most riceproducing countries. However, there are few studies related to the logistics optimization of this product. Thus, in this section, the studies performed on supply chain optimization of other agricultural products are reviewed in order to discover the existing research gap. Therefore, the published papers related to this area are reviewed in two sub-sections called ASC optimization and ASC sustainability.

ASC Optimization
As mentioned in the prior sections, agriculture is one of the most salient sectors influencing the economies of many countries. As a result, many researchers have made efforts to enhance the logistics management of agricultural products using novel methods such as operational research. One of the primary studies on programming models in ASC was conducted by [12]. They conducted research on various agricultural products, both perishable and non-perishable, and vegetables. The authors of [13] addressed the production and distribution of perishable products and presented a mathematical model for optimizing the product's quality and freshness. The author of [14] presented an integrated model for investment-related decisionmaking in the fruit and vegetable industry. The authors of [15] proposed a transportation programming model for a fruit supply chain in which a fruit logistics center was established by various storage centers to meet demand during the non-harvest seasons. The authors of [16] proposed a mathematical model for optimizing the cost of purchasing, transporting, and storing fresh agricultural products. They tested the performance of their proposed model on a real case study in an apple juice factory.
The authors of [17] presented a multi-period, non-linear mathematical programming model for optimizing the transportation and operational costs in the logistics system of wheat. The authors of [18] developed a new mathematical model with the goal of reducing the costs of citrus CLSC and maximizing responsiveness to customers' demand in forward and reverse flow. The work in [19] designed a logistics network for a wheat supply chain and proposed a mathematical programming model for minimizing the total costs under demand uncertainty. They used robust optimization to deal with uncertainty.
The authors of [20] presented a non-linear mathematical model for optimizing the costs and the total profit of an agricultural logistics network. They employed several metaheuristic algorithms such as NSGA-II, MOICA, and MOPSO. In other research, the authors of [21] presented a mathematical model for optimizing the costs in a rice supply chain network. They applied some well-known metaheuristic algorithms for solving the proposed model. Their derived results denoted the proposed model and their solution methods as being valid, practical, and effective.
The authors of [22] developed a novel CLSC network for walnut product. Then, they formulated an MILP model for minimizing the total costs of the proposed network. In addition, they used several metaheuristic algorithms for solving the presented model. Their results indicated their solution approach and the model being valid. Developing a CLSC network for the avocado industry was another study in [23]. An MILP mathematical model was formulated for optimizing the total costs and the number of the job opportunities created in the proposed network. They employed exact methods to solve the presented model. The authors of [24] developed a multi-echelon supply chain network for the sugarcane industry. Following that, they formulated an MILP mathematical model to minimize the total network costs, and due to its NP-hard nature, they applied some hybrid metaheuristic algorithms to solve it. Their results revealed the efficiency of the model and the solution approach.

ASC Sustainability
In recent years, the sustainability concept in supply chains, especially the logistics of agricultural products, has appealed to many researchers. The environmental and social effects of producing and consuming agricultural products have led the decision-makers toward sustainable management, and several related studies have been performed by researchers. For instance, the authors of [25] presented a two-objective mathematical model to minimize demand overruns and to maximize profits in a pear logistics network. They used a lexicographic method to solve their model. The authors of [26] presented a sustainable, multi-objective MILP mathematical model for optimizing the CLSC network of mushroom production by considering a variety of recycling technologies. Their results revealed that using technology increases the overall profit by up to 12% and reduces emissions by about 28%. The authors of [27] presented a mathematical model for a wheat supply chain network, pursuing the goal of optimizing the total costs (including the fixed costs for selecting the suppliers and locating the warehouses and the variable costs for maintenance, transportation, and production). They employed GAMS software for solving the model, the results of which indicated that their proposed model can be used in making decisions about wheat import and distribution.
The authors of [28] proposed a multi-objective MILP model for optimizing wheat production in Spain while considering the environmental impacts.
The authors of [29] analyzed the agro-food supply chain design, focusing on the sustainable dimensions in their multi-objective model. The work in [30] presented a multiobjective non-linear programming model for the sustainable supply chain of perishable agro-food products produced through organic and non-organic methods. Their goal was to strike a balance between the production and consumption of the organic and non-organic products and reduce the costs, lower environmental degradation, and increase the levels of consumer health, wherein they used the Epsilon constraint method to solve the model.
The work in [31] presented a mathematical programming model with the goal of optimizing the costs in an apple supply chain and investigated the environmental effects in their network.
The work in [2] designed an MILP model to optimize the costs, meet the demand, and reduce CO 2 emissions in a citrus CLSC network. Moreover, they employed three metaheuristic algorithms called MOTGA, NRGA, and NSGA-II. Their results and analysis suggested MOTGA outperformed in this regard.
The authors of [32] prepared a multi-objective mathematical programming model for optimizing the costs, water consumption, and the number of job opportunities created in a wheat logistics network under demand uncertainty. Next, they used a simulation approach for estimating the demand. The authors of [33] presented a bi-objective mathematical model for optimizing the profit and the quantity of the pollutants in a green pistachio supply chain. In their model, the demand and costs were assumed as uncertain. They used a robust possibilistic programming model to deal with uncertainty. In addition, they employed the Epsilon constraint method to solve the model. The work in [27] proposed a multi-objective MILP model to optimize the costs and the quantity of the pollutants emitted in a date supply chain. They used the LP-metric method and some metaheuristic algorithms to solve the model.

Research Gap and Motivation for Research
In the current study, in order to come across the study gap, we analyzed several papers about ASC optimization. A brief review of some related studies is given in Table 1. Briefly, by investigating the mentioned papers in the table above, it can be concluded that few papers have been published about rice logistics optimization. Furthermore, the sustainability dimensions in ASC optimization have rarely been considered concurrently. The agricultural sector, on the one hand, has enormous environmental effects and plays a critical role in climatic changes, water scarcity, land degradation, deforestation, and other processes, and on the other hand, the income of many rural people depends on it and the sector provides several job opportunities for them. Thus, considering the economic, environmental, and social factors simultaneously and creating a trade-off among them greatly help to improve the ASC network, particularly for rice, an issue that has rarely been analyzed. Besides the sustainability concept, uncertainty has seldom been considered in the optimization models in the agricultural sector. Undoubtedly, the supply chain of agricultural products is challenging, with various uncertainties in the real world, and therefore it is necessary to employ proper approaches for modeling under uncertain conditions for successful optimization and effective decision making in the supply chain.
Furthermore, considering the energy sources is rarely observed in ASC optimization models. There is no doubt that the processing of agricultural products, especially rice, in various processes, including its own processing, requires electrical energy. Processing stages are significantly energy-consuming, which exerts environmental effects. Thus, it urgently necessitates employing renewable and clean energy sources for rice processing. As a result, selecting an optimal combination of traditional and renewable energy used for processing can boost energy efficiency and mitigate the environmental consequences. To the best of our knowledge, the present study is one of the first papers that introduces a mathematical model with the goal of optimally using rice waste and properly locating solar panel sites and bio-refineries for minimizing the fixed and operational costs intending to generate electricity. The asterisk means the model is single period or multi period.
Moreover, when investigating the relevant studies, the traditional and online channels in the ASC network design have rarely been analyzed. Today, with the expansion of E-commerce, most consumers are increasingly becoming fond of purchasing agricultural products online. The online purchase of agricultural products has attracted many researchers due to advantages such as more competitive prices, timely purchase, easy product delivery, and fast product delivery [36].
In order to bridge the research gap in this study, we introduce a dual-channel CLSC network, including the traditional and online purchasing of rice, while taking into account the energy sources. Then, we present a three-objective mathematical model for optimizing the costs, the quantity of the emitted pollutants, and the number of job opportunities created in the rice supply chain network in a fuzzy environment. Next, the designed model is solved via some multi-objective metaheuristic algorithms and their results and performance are compared. Furthermore, using a case study in Iran, the accuracy of the programmed model is surveyed. In the following sections, the problem model and its solution method are presented and described.

Problem Statement
In this section, the logistics network and the proposed mathematical model for optimizing the objectives are described. The proposed multi-period network in Figure 8 involves the producers (farmers), the distribution centers, the solar panel sites, the bio-refineries, the rice factories, the warehouses, the markets (retailers), the customers, and the toiletry and pharmaceutical industries. In this network, the farmers harvest rice in two periods and transfer the unprocessed rice (paddy) from farms to the distribution centers. Moreover, rice harvest-produced waste, i.e., the straw, is transferred to the recycling centers and bio-refineries to produce compost and convert it into electricity, respectively. The paddy from the farmers is transferred to the distribution centers in two months because the harvest is completed in two months at a maximum level. Furthermore, the maximum shelf life of rice in the distribution centers is 6 months and, from there, it is transferred to the rice factories for processing. Figure 9 depicts the production and processing stages of rice. As rice processing requires a lot of electricity, some points have been considered as potential sites for constructing solar panels and bio-refineries with the goal of supplying energy.  Moreover, rice processing-induced waste, i.e., rice husk, is transferred to the biorefineries due to its potential to generate electricity. Ultimately, both the processed rice and the bran are dispatched from the factory to both the markets and directly to the consumers to satisfy their demand. As observed in Figure 8, a dual channel was introduced in the proposed supply chain network in which the customers can purchase rice and bran both online and directly from the warehouse and from the market in a traditional manner. Some rice bran can also be sent to the toiletry and pharmaceutical industries. It is worth noting that the location of all facilities except for the new distribution centers, the new recycling centers, the solar panel sites, and the bio-refineries, is fixed. The goal behind this model is to appropriately locate the construction of solar panels, the biorefineries, and the recycling centers to create a balance among the total cost, the quantity of the pollutants, and the number of the job opportunities created by considering the energy sources in the network so that the consumers' demand is met. The indicators, the parameters, and the variables of decision making are presented in the Appendix A, and here, the mathematical model is explained. The proposed model's assumptions are the following: • Transportation costs between the network facilities are consistent with the distance.

•
Processing centers have limited storage capacity.  Moreover, rice processing-induced waste, i.e., rice husk, is transferred to the biorefineries due to its potential to generate electricity. Ultimately, both the processed rice and the bran are dispatched from the factory to both the markets and directly to the consumers to satisfy their demand. As observed in Figure 8, a dual channel was introduced in the proposed supply chain network in which the customers can purchase rice and bran both online and directly from the warehouse and from the market in a traditional manner. Some rice bran can also be sent to the toiletry and pharmaceutical industries. It is worth noting that the location of all facilities except for the new distribution centers, the new recycling centers, the solar panel sites, and the bio-refineries, is fixed. The goal behind this model is to appropriately locate the construction of solar panels, the biorefineries, and the recycling centers to create a balance among the total cost, the quantity of the pollutants, and the number of the job opportunities created by considering the energy sources in the network so that the consumers' demand is met. The indicators, the parameters, and the variables of decision making are presented in the Appendix A, and here, the mathematical model is explained. The proposed model's assumptions are the following: • Transportation costs between the network facilities are consistent with the distance.

•
Processing centers have limited storage capacity. Moreover, rice processing-induced waste, i.e., rice husk, is transferred to the biorefineries due to its potential to generate electricity. Ultimately, both the processed rice and the bran are dispatched from the factory to both the markets and directly to the consumers to satisfy their demand. As observed in Figure 8, a dual channel was introduced in the proposed supply chain network in which the customers can purchase rice and bran both online and directly from the warehouse and from the market in a traditional manner. Some rice bran can also be sent to the toiletry and pharmaceutical industries. It is worth noting that the location of all facilities except for the new distribution centers, the new recycling centers, the solar panel sites, and the bio-refineries, is fixed. The goal behind this model is to appropriately locate the construction of solar panels, the bio-refineries, and the recycling centers to create a balance among the total cost, the quantity of the pollutants, and the number of the job opportunities created by considering the energy sources in the network so that the consumers' demand is met. The indicators, the parameters, and the variables of decision making are presented in the Appendix A, and here, the mathematical model is explained. The proposed model's assumptions are the following:

•
Transportation costs between the network facilities are consistent with the distance. • Processing centers have limited storage capacity. • Deficiency cost has not been considered. • Rice production capacity in factories is also limited. • Excess electrical energy should be connected to the mains (grid) electricity.

Problem's Model
After describing the offered problem and assumptions, the expanded multi-objective mathematical model is formulated as the following: The model's objective function z 1 related to the economic aspect of sustainabilty maximizes the total cost. This amount of the sum of the cost refers to the consumers' rice purchase (1a), the cost of establishing new facilities (1b), transportation cost (1c), the cost of production for farmers, the cost of maintaining (storing) the product, and the cost of production for the rice factories (1d), and the cost of energy tax (1e). Objective function z 2 related to the environmental aspect of sustainability minimizes the total amount of CO 2 emitted from constructing the new facilities (2a) and the vehicles traveling among the network facilities (2b). The vehicle has a certain capacity (capv) that is based on the quantity of the product transported in each stage (for instance, n for the first time (n = Xa ijt / capv )), the number of transportation times is determined; this vehicle needs a certain amount of fuel per k(ϕ) for transportion, and when this fuel is used by the vehicle, a certain amount of CO 2 per km (RF) is emitted.
Objective function z 3 is associated with the third aspect of sustainability and maximizes the fixed and variable job opportunities. The fixed job opportunities are independent of production quantity, but the variable job opportunities are dependent on the production quantity. This objective function consists of three segments: the 1st one optimizes the fixed job opportunities created by opening the new facilities; the 2nd one optimizes the variable job opportunities of rice production in the farms; the 3rd one optimizes the variable job opportunities of rice production in the factories. Our target behind presenting this model was to strike a balance among the costs, the quantity of the pollutants, and the number of the fixed and variable job opportunities. This model's constraints are the following: Constraints: Constraint (4) states that the amount of rice plant harvested by the producers is the maximum of their production capacity. Constraint (5) displays the equilibrium in the distribution centers and states that the sum of the input product to each distribution center is larger than or equal to the output product. Constraint (6) balances the inventory of the unprocessed rice in the distribution centers. Constraint (7) states that the condition for dispatching the product to the distribution centers is the establishment of that center. Constraint (8) denotes that the quantity of the unprocessed rice transferred to each factory has to be less than its storage capacity. Constraints (9) and (10) show that the demand for the rice required by the customers has to be met through traditional and online shopping. Constraints (11) and (12) state that the rice bran demanded by the pharmaceutical and toiletry industries has to be satisfied. Constraint (13) indicates that the quantity of the produced product has to be less than the production capacity of each factory. Constraint (14) demonstrates that the amount of rice or bran produced in factories is gained as a fraction of paddy.
Constraints (15) and (16) state that the total bran and rice produced at each factory has to be delivered to the consumers. Constraints (17) and (18) display the storage capacity of the collection centers. Constraint (19) shows that the quantity of the unprocessed rice transported to the distribution centers has to be equal to the amount of the produced rice plant multiplied by its conversion rate to paddy. Constraints (20) and (21) demonstrate that the energy required for processing rice is supplied from a renewable energy source if it is opened. Constraint (22) displays that rice straw is transported to the recycling centers if that center is opened. Constraint (23) states that rice straw and husk are dispatched to the bio-refineries if the center is opened. Constraint (24) shows that the quantity of rice transported to each center has to be less than its storage capacity. Constraint (25) displays the compost demanded by the farmers that has to be met during each period. Constraint (26) demonstrates that the amount of rice straw transferred to the bio-refineries and the recycling centers is a fraction of the quantity of the produced rice plant.
Constraint (27) states that the amount of compost produced in recycling centers is equal to the amount of the received rice straw multiplied by its conversion rate to compost. Constraint (28) demonstrates that the energy generated from rice waste in each bio-refinery is equal to the amount of the received rice straw and husk multiplied by their conversion rate into energy. Constraint (29) indicates that the amount of rice husk transported to the bio-refinery is limited and equal to the amount of the received paddy multiplied by its conversion rate to husk.
Constraint (30) indicates that the quantity of the rice and the produced bran in each factory is bound to the available energy. Constraint (31) shows that the energy generated from each source is bound to the energy generation capacity, and Constraint (32) balances the flow in each market. Finally, Constraint (33) depicts the model's type of decision variables.

Uncertainty Model
In the proposed model, the parameters of purchase cost, supply, and demand are uncertain. In this section, using the fuzy logic proposed by Jimenez's method [37], we introduce the uncertain counterpart of the proposed model. As Jimenez's method does not impose additional inequality constraints, it is efficient in solving fuzzy decision problems [7]. Pursuant to this method, the fuzzy parameters presented in Equations (34) and (35) can be represented by a set of triangular fuzzy numbers for the optimistic, realistic, and pessimistic cases. Regarding Jimenez's method, the Expected Interval (EI) and the Expected Value (EV) of a fuzzy parameter can be estimated as given in the following: The fuzzification of a parameter based on Jimenez's method can be presented as follows. Suppose a fuzzy parameter of demand in the optimistic, realistic, and pessimistic state including the values c w =3, 4, and c v =6 units. In Jimenez's method, the expected value of the fuzzy parameter can be estimated as the following: To convert the set of constraints like a i X ≥ b i , i = 1, 2, . . . I, the following convertion is performed: If the uncertainty results in the model's infeasibility, the following set of equations can be accepted: Furthermore, when the constraints become non-fuzzy, the membership function proposed by Torabi-Hassani (TH) in [37] is used according to Equation (38): Moreover, when we intend to minimize the objective function, the following membership function is employed: where the positive ideal solutions (a-PIS) and the negative ones (a-NIS) of the objective function (Z) at the feasibility level (a) are considered for the total cost's objective function.
Thus, the fuzzy model equivalent to the original problem's model can be formulated as follows: Problem's model: Equations (2) and (3) Equations (11) and (33).

Solution Approach
In the current study, a multi-objective MILP model is proposed to create a balance among the cost, the quantity of the pollutants, and the job opportunities in a rice logistics network. In order to solve this model, four metaheuristic algorithms, including a multiobjective version of the reptile search optimizing algorithm called MORSO and three popular algorithms, namely MOSA, MOPSO, and MOGWO, with priority-based encoding are used as the proposed approach. Moreover, the LP-metric method is applied to evaluate the performance of the aforementioned algorithms.
The algorithms of MOPSO and MOSA have been further explained in [21] and MOGWO has been given in [38]. Therefore, the explanation of these algorithms is taken for granted and the MORSO algorithm is described separately in a sub-section. In this section, we explain the proposed solution algorithm in the sub-section Encoding and Decoding and indicate how to satisfy the model's constraints using this approach. Then, in the continuation, the criteria for comparing the algorithms are described.

Encoding and Decoding
In this research, the priority-based encoding method developed by [39] is used to display the initial solution. Considering the encoding and decoding design and method, a small size example is needed to meet the constraints. Here, in Figure 10, the initial solution's structure is illustrated.

Solution Approach
In the current study, a multi-objective MILP model is proposed to create a balance among the cost, the quantity of the pollutants, and the job opportunities in a rice logistics network. In order to solve this model, four metaheuristic algorithms, including a multi-objective version of the reptile search optimizing algorithm called MORSO and three popular algorithms, namely MOSA, MOPSO, and MOGWO, with priority-based encoding are used as the proposed approach. Moreover, the LP-metric method is applied to evaluate the performance of the aforementioned algorithms.
The algorithms of MOPSO and MOSA have been further explained in [21] and MOGWO has been given in [38]. Therefore, the explanation of these algorithms is taken for granted and the MORSO algorithm is described separately in a sub-section. In this section, we explain the proposed solution algorithm in the sub-section Encoding and Decoding and indicate how to satisfy the model's constraints using this approach. Then, in the continuation, the criteria for comparing the algorithms are described.

Encoding and Decoding
In this research, the priority-based encoding method developed by [39] is used to display the initial solution. Considering the encoding and decoding design and method, a small size example is needed to meet the constraints. Here, in Figure 10, the initial solution's structure is illustrated.  Figure 10. Schematic design of proposed arrays.
Suppose that the number of the producers, the distribution centers, the rice factories, the pharmaceutical industries, the toiletry industries, the markets, the customers, the solar panel sites, the bio-refineries, and the recycling centers are 2, 1, 1, 3, 1, 2, 1, 2, 2, and 1, respectively. The proposed array is a matrix with six rows and with the number of columns as 3 * + 2 * + 3 * + + + 2 * + 2 * + 2 * + . The cells of the 1st sub-segment are filled with random numbers between 0 and 1. In the next stage, as seen in Figure 11, the cells of the 1st sub-segment are sorted by their priority. Sorting the numbers is performed separately for each segment. Regarding the encoding of the 1st segment, whose steps are displayed in Table 2, Constraints (4)-(6) can be satisfied. Moreover, using the 2nd and 3rd segments' encoding displayed in Tables A1 and A2 in Appendix A, Constraints (17) and (18) are satisfied. In addition, Suppose that the number of the producers, the distribution centers, the rice factories, the pharmaceutical industries, the toiletry industries, the markets, the customers, the solar panel sites, the bio-refineries, and the recycling centers are 2, 1, 1, 3, 1, 2, 1, 2, 2, and 1, respectively. The proposed array is a matrix with six rows and with the number of columns as 3 * i + 2 * j + 3 * f + u + k + 2 * m + 2 * b + 2 * o + u. The cells of the 1st sub-segment are filled with random numbers between 0 and 1. In the next stage, as seen in Figure 11, the cells of the 1st sub-segment are sorted by their priority. Sorting the numbers is performed separately for each segment. Regarding the encoding of the 1st segment, whose steps are displayed in Table 2, Constraints (4)-(6) can be satisfied. Moreover, using the 2nd and 3rd segments' encoding displayed in Tables A1 and A2 in Appendix A, Constraints (17) and (18)       For t = 1:T Inputs: I = set of producers J = set of distribution centers Ca(i, t) = production capacity of producer i in period t D(j, t) =capacity V(L + N) =encode Dis(i, j) = Distance between nodes Outputs: Xaloc(i, j, t) =amount in period t Em(i, j, t) =Amount of CO 2 emission caused by transferring quantity between node i and in period t W(j) = binary variable shows the distribution centers j is opened

Reptile Search Algorithm
The Reptile Search Optimizer (RSO) proposed by [40] is a population-based algorithm that simulates the hunting mechanisms and social behavior of crocodiles in the wild. The algorithm is inspired by some of the key features of crocodile behavior, including encircling the prey and the coordination of crocodiles during an attack. The mechanisms of encircling and hunting prey are mathematically modeled as follows.

Encircling Phase
In this stage, the search space for finding a better solution is analyzed based on two main strategies including high walking and belly walking. Selecting the strategy depends on the number of iterations so that as long as t ≤ T 4 , the high walking strategy is selected, and as long as > T 4 , t ≤ T 2 , the belly walking strategy is chosen. To update the crocodiles' position in the exploration stage in each iteration, the following equations are applied: where Best j (t) stands for the jth position in the best-obtained solution in iteration t and rand is the random number between 0 and 1 and T is the maximum number of iterations., which is calculated according to Equation (46), is a hunting operator for the position in the ith solution. The parameter β is responsible for controlling the exploration accuracy in each iteration, which is equal to 0.1. The reduction function, calculated according to Equation (47), reduces the search space. r1 is a random value between the range [1, N], and x r1,j indicates the position of the ith solution. N is the number of the solutions. Evolutionary Sense (ES (t)) is also a probability ratio of randomly decreasing values between 2 and −2 throughout the number of iterations, which is calculated using Equation (48).
In Equation (47), ε is a small value and r2 is a random number between [1, N]. In addition, r3 denotes a random integer between and 1. P i,j is the percentage difference between the jth position of the best-obtained solution and the jth position of the current solution, which is calculated using Equation (49): where, in Equation (49), M(x i ) is the average positions of the ith solution, calculated using Equation (50). UB j and LB j are the upper and lower boundaries of the position, respectively. Finally, α is a parameter, which controls the exploration accuracy during iterations, which is fixed equal to.

Hunting Simulation
The goal behind this mechanism is to escape being trapped in the optimal local points, which is based on two strategies during hunting, coordination and cooperation. Like the encircling mechanism, selecting the strategy depends on the number of iterations so that when, crocodiles select a hunting coordination strategy; otherwise, the cooperation hunting strategy is selected, which is conditioned by t > 3T 4 , t ≤ T. To update the position of the crocodiles in this phase in each iteration, the following equations are applied: where the values of each parameter are similar to what has been pointed out in the previous part. Consequently, the computational complexity of the proposed RSA is as follows: where T is the number of iterations, N presents the number of used solutions, and D presents the solution size. The pseudo-code of the RSO algorithm is depicted in Table A3 in the Appendix A.

Multi-Objective Reptile Search Optimization Algorithm
In this sub-section, a multi-objective version of the RSO known as MORSO is proposed in order to solve the proposed multi-objective model. Like other multi-objective metaheuristic algorithms, we deal with concepts such as archive, grid approach, and leader selection in this paper.

Archive and Grid Approach
Archive is in charge of saving, controlling, and retrieving the optimal achieved Pareto solutions. During each iteration, the position of the search agents is updated based on the mechanism of the RSO algorithm and the new obtained solutions are compared with the archive members. If the new solution dominates one of the archive members, it will substitute with it. The archive has a limited capacity, and to delete a solution, the grid approach is employed to select one of the members in the most crowded part of the archive and set it aside.

Selecting a Leader
For selecting a leader, one of the best-obtained optimal solutions in the archive is selected by the roulette wheel method, and other search agents update their position in order to attack the prey. The flowchart and the pseudo-code of the MORSO algorithm are illustrated in Figure 12 and Table 3, respectively.

Archive and Grid Approach
Archive is in charge of saving, controlling, and retrieving the optimal achieved Pareto solutions. During each iteration, the position of the search agents is updated based on the mechanism of the RSO algorithm and the new obtained solutions are compared with the archive members. If the new solution dominates one of the archive members, it will substitute with it. The archive has a limited capacity, and to delete a solution, the grid approach is employed to select one of the members in the most crowded part of the archive and set it aside.

Selecting a Leader
For selecting a leader, one of the best-obtained optimal solutions in the archive is selected by the roulette wheel method, and other search agents update their position in order to attack the prey. The flowchart and the pseudo-code of the MORSO algorithm are illustrated in Figure 12 and Table 3, respectively.

Evaluation Indices of Algorithms' Performance
In this research, five indices were used to evaluate the algorithms' performance. These metrics measure different criteria, which can be listed as follows.

Number of Pareto Solutions (NPS)
In this criterion, the number of Pareto solutions is computed. Every method with a higher NPS criterion is better. The ideal state in this method is that a higher number of Pareto solutions get distributed more evenly in the possible space.

Mean Ideal Distance (MID)
This index is used to calculate the distance between Pareto solutions. Considering Equation (52), the lower this index, the higher the algorithm's performance. In the current study, the ideal point equals the minimum of each of the objective functions. This index is calculated as shown by Equation (52).

Maximum Spread (MS)
This criterion measures the spread of non-dominated solutions. The more spread out the non-dominated solutions are, the larger this index will be, and the higher this value, the more appropriate. This criterion is estimated by Equation (53)

Spread of Non-Dominance Solutions (SNS)
This index is used to check the variety of non-dominance solutions and is calculated by Equation (54). The higher this index, the better the algorithm performs. In this equation, F1 i , F2 i , and F3 i are the values of the 1st, 2nd, 3rd objective functions for the non-dominance ith solution.

Computational Time Index (CPU Time)
In large-scale problems, one of the most critical indices is the CPU time, and the lower this value, the lower the efficiency of the algorithm.

Validation and Analysis of Results
In this section, the accuracy of the programmed model is analyzed by implementing it onto a real case study. Then, in another sub-section, the parameters of the proposed algorithms are tuned to come up with better results. In the third sub-section, to further evaluate the model, sensitivity analysis is performed on some key parameters.

Case Study
In this section, the accuracy of the programmed model is analyzed by implementing it onto a real case study in Iran, i.e., in Mazandaran Province. At the moment, over 300,000 hectares are under rice cultivation in Iran and this province is the largest producer of rice in Iran [21]. In this study, in order to collect the data, some of the province-based cities were considered as the position of the network facilities; Figure 13 depicts the position of the cities. Moreover, ten sample problems were generated based on the number of the network facilities for evaluating the proposed model's efficiency, as seen in Table 4. As the information collected from the farmers indicated, the conversion rate of rice plant to paddy is about 0.8. Therefore, the conversion rate of rice to straw is about 0.2. In addition, the conversion rate of paddy to processed rice is about 0.64, and the conversion rates of paddy to bran and to husk are about 0.26 and 0.1, respectively. In bio-refineries, around 26.7 kWh of electricity can be supplied per kg of rice straw and husk. The power generation capacity of solar panel sites ranges from 200 to 400 MWh. number of the network facilities for evaluating the proposed model's efficiency, as seen in Table 4. As the information collected from the farmers indicated, the conversion rate of rice plant to paddy is about 0.8. Therefore, the conversion rate of rice to straw is about 0.2. In addition, the conversion rate of paddy to processed rice is about 0.64, and the conversion rates of paddy to bran and to husk are about 0.26 and 0.1, respectively. In bio-refineries, around 26.7 kWh of electricity can be supplied per kg of rice straw and husk. The power generation capacity of solar panel sites ranges from 200 to 400 MWh.   7  65  24  11  24  16  12  28  15  25  24  22  8  72  30  14  28  20  15  35  16  31  30  26  9  80  34  17  32  24  19  42  18  37  39  34  10  82  38  20  36  28  21  49  20  42  40  39 The transportation costs among the facilities correspond with the distance, which is provided in Table A4 in Appendix A. In addition, a nine-ton truck is used to transport the products. This vehicle uses R = 0.0832 fuel per k of transportation and its fuel-emitted CO 2 rate (ϕ) is 3.15. The values of these two parameters have been extracted from [38]. Some data are also generated randomly, as demonstrated in Table A5 in Appendix A.

Results Analysis
Firstly the proposed algorithms are tuned to achieve better results. After that, the results from solving the model are analyzed, and for further model evaluation, sensitivity analysis is performed on parameters such as the energy supply capacity and customers' demand.

The Tune of the Algorithm's Parameters
There are several methods for tuning the metaheuristic algorithms' parameters with the goal of improving their performance, out of which the Taghuchi method is utilized in the present paper. Taguchi developed new statistical concepts and combined and established particular groups of orthogonal arrays to present the tests. This method classifies a group of factors based on orthogonal arrays into two main parts, including control and noise factors, and while maximizing the effect of the control factors, it minimizes that of the noise factors according to the following equation [41].
where Y stands for the solution value and n is the number of the orthogonal arrays. Moreover, Equation (56) presents the selected responses in this study. Two main concepts related to the solution are convergence and diversity. The MID metric measures the convergence of the solutions and the variety of Pareto solutions is gained by the MS metric [2].
It is necessary to determine the parameters we need to tune in every algorithm to perform the Taguchi test. In the present study, three levels were considered for all parameters given in Table 5. The MOGWO algorithm has three factors, and MOPSO, MOSA, and MORSO algorithms have four factors. By performing the experiment in Minitab software, L9 orthogonal arrays are proposed. The orthogonal arrays of each algorithm and their derived results are illustrated in Tables A6-A9 in Appendix A. Moreover, the S/N plots for the mentioned algorithms are displayed in Figures 14-17, by which the highest level is the best level for each algorithm. As seen in the plots, for example, in the simulated annealing algorithm, the best level for the maximum 3rd level iteration factor is the 3rd one, the best level for the population number factor is the 3rd level, for Phi 1 factor this is the 2nd level, and for Phi 2 factor, the best level is the 1st level.  The orthogonal arrays of each algorithm and their derived results are illustrated in Tables A6-A9 in Appendix A. Moreover, the S/N plots for the mentioned algorithms are displayed in Figures 14-17, by which the highest level is the best level for each algorithm. As seen in the plots, for example, in the simulated annealing algorithm, the best level for the maximum 3rd level iteration factor is the 3rd one, the best level for the population number factor is the 3rd level, for Phi1 factor this is the 2nd level, and for Phi2 factor, the best level is the 1st level.

Main Effects Plot for SN ratios Data Means
Signal-to-noise: Smaller is better

Analysis of Results
We designed 10 sample problems for evaluating the proposed model's efficiency. For this purpose, a computer with 4 GB of RAM and 2.2 GHz CPU was employed and encoding the proposed models and algorithms was run in MATLAB software. The first to fourth sample problems were solved by LP-metric and the metaheuristic algorithms and the rest were solved by metaheuristic algorithms regarding the model being an NPhard.
To evaluate the aforementioned metaheuristic algorithms, a one-to-one comparison was performed as the criteria indicated in Section 4. The top performance is based on

Main Effects Plot for SN ratios Data Means
Signal-to-noise: Smaller is better

Analysis of Results
We designed 10 sample problems for evaluating the proposed model's efficiency. For this purpose, a computer with 4 GB of RAM and 2.2 GHz CPU was employed and encoding the proposed models and algorithms was run in MATLAB software. The first to fourth sample problems were solved by LP-metric and the metaheuristic algorithms and the rest were solved by metaheuristic algorithms regarding the model being an NPhard.
To evaluate the aforementioned metaheuristic algorithms, a one-to-one comparison was performed as the criteria indicated in Section 4. The top performance is based on

Main Effects Plot for SN ratios Data Means
Signal-to-noise: Smaller is better

Analysis of Results
We designed 10 sample problems for evaluating the proposed model's efficiency. For this purpose, a computer with 4 GB of RAM and 2.2 GHz CPU was employed and encoding the proposed models and algorithms was run in MATLAB software. The first to fourth sample problems were solved by LP-metric and the metaheuristic algorithms and the rest were solved by metaheuristic algorithms regarding the model being an NPhard.
To evaluate the aforementioned metaheuristic algorithms, a one-to-one comparison was performed as the criteria indicated in Section 4. The top performance is based on

Main Effects Plot for SN ratios Data Means
Signal-to-noise: Smaller is better

Main Effects Plot for SN ratios Data Means
Signal-to-noise: Smaller is better Figure 17. Signal to noise plot of MORSO.

Analysis of Results
We designed 10 sample problems for evaluating the proposed model's efficiency. For this purpose, a computer with 4 GB of RAM and 2.2 GHz CPU was employed and encoding the proposed models and algorithms was run in MATLAB software. The first to fourth sample problems were solved by LP-metric and the metaheuristic algorithms and the rest were solved by metaheuristic algorithms regarding the model being an NP-hard.
To evaluate the aforementioned metaheuristic algorithms, a one-to-one comparison was performed as the criteria indicated in Section 4. The top performance is based on lower MID and CPU time criteria and higher NPS, MS, and SNS. Table 6 displays the results of this comparison. As seen in the table, metaheuristic algorithms have less execution time compared to the LP-metric method. As the number of sample problems increases, the execution time of the model with the LP-metric method should increase, which indicates that the model is NP-hard. Therefore, we need meta-algorithms in order to solve the model in a reasonable time and reduce its complexity [42,43]. In order to compare the algorithms in terms of performance, ANOVA was used. The intervals plot (at confidence level 95%) are sketched pursuant to the data in Table 6 separately for each of the algorithms and according to each criterion in Figures 18-22. As perceived from these plots, the MORSO algorithm performed better than other algorithms in terms of NPS, MS, and MID criteria. These plots indicate that MOSA is faster than other algorithms. In terms of SNS criteria, all three algorithms performed similarly close to each other. Thus, it can be concluded that the MORSO algorithm performed better than other algorithms. An example of non-dominance solutions for the first sample problem for the mentioned metaheuristic algorithm is shown in Figures 23-26. close to each other. Thus, it can be concluded that the MORSO algorithm performed bet ter than other algorithms. An example of non-dominance solutions for the first sampl problem for the mentioned metaheuristic algorithm is shown in Figures 23-26.    close to each other. Thus, it can be concluded that the MORSO algorithm performed bet ter than other algorithms. An example of non-dominance solutions for the first sampl problem for the mentioned metaheuristic algorithm is shown in Figures 23-26.    close to each other. Thus, it can be concluded that the MORSO algorithm performed bet ter than other algorithms. An example of non-dominance solutions for the first sampl problem for the mentioned metaheuristic algorithm is shown in Figures 23-26.

Sensitivity Analysis
In this sub-section, for further evaluation of the proposed model and the top algo rithm, the sensitivity analysis is conducted in two states. It is worth mentioning that sen sitivity analysis was performed for the first sample problem. In the first state, the sensi tivity analysis is performed on the parameter of capacity of each energy source, and in the second state, on the parameter of customers' demand for the processed rice.

Sensitivity Analysis
In this sub-section, for further evaluation of the proposed model and the top algo rithm, the sensitivity analysis is conducted in two states. It is worth mentioning that sen sitivity analysis was performed for the first sample problem. In the first state, the sensi tivity analysis is performed on the parameter of capacity of each energy source, and in the second state, on the parameter of customers' demand for the processed rice.

Sensitivity Analysis
In this sub-section, for further evaluation of the proposed model and the top algo rithm, the sensitivity analysis is conducted in two states. It is worth mentioning that sen sitivity analysis was performed for the first sample problem. In the first state, the sensi tivity analysis is performed on the parameter of capacity of each energy source, and in the second state, on the parameter of customers' demand for the processed rice.

Sensitivity Analysis
In this sub-section, for further evaluation of the proposed model and the top algorithm, the sensitivity analysis is conducted in two states. It is worth mentioning that sensitivity analysis was performed for the first sample problem. In the first state, the sensitivity analysis is performed on the parameter of capacity of each energy source, and in the second state, on the parameter of customers' demand for the processed rice.

Sensitivity Analysis on Energy Source Capacity Parameter
In this state, it is assumed that the capacity of each energy source would decrease or increase consistent with the conditions listed in Table 7.
First condition: It is assumed that all electricity required for product processing should be supplied by the mains electricity and not by the renewable energy sources.
Second condition: It is assumed that the capacity of mains electricity would decrease by 20% and the capacity of renewable energy sources would increase by 30%.
Third condition: It is assumed that the capacity of the mains electricity would decrease by 30% and the capacity of renewable energy sources would increase by 40%. Having solved the model under these conditions by the LP-metric method, the results are provided in Figures 27-29. The first condition states that the utilization of the mains electricity would increase by about 19%. The second and third conditions state that with the renewable sources' capacity increase, the utilization of mains electricity would decrease by about 7% and 9%, respectively.
In this state, it is assumed that the capacity of each energy source would decrease or increase consistent with the conditions listed in Table 7.

Bio-Refinery Solar Panels Mains Electricity
Bio-Refinery Solar Panels Mains Electricity Condition First condition: It is assumed that all electricity required for product processing should be supplied by the mains electricity and not by the renewable energy sources.
Second condition: It is assumed that the capacity of mains electricity would decrease by 20% and the capacity of renewable energy sources would increase by 30%.
Third condition: It is assumed that the capacity of the mains electricity would decrease by 30% and the capacity of renewable energy sources would increase by 40%. Having solved the model under these conditions by the LP-metric method, the results are provided in Figures 27-29. The first condition states that the utilization of the mains electricity would increase by about 19%. The second and third conditions state that with the renewable sources' capacity increase, the utilization of mains electricity would decrease by about 7% and 9%, respectively.

Sensitivity Analysis on Demand
In the second state, sensitivity analysis was carried out on the traditional shopping consumers' demand parameter. Sensitivity analysis is performed under five states in which the demand parameter decreases and increases. Suppose that consumers' demand increases by 10% and 20% in two states and decreases by the same level in two states. The third state is consistent with the base state. After solving the model under all conditions, the results are displayed in Figures 27-29.

Sensitivity Analysis on Demand
In the second state, sensitivity analysis was carried out on the traditional shopping consumers' demand parameter. Sensitivity analysis is performed under five states in which the demand parameter decreases and increases. Suppose that consumers' demand increases by 10% and 20% in two states and decreases by the same level in two states. The third state is consistent with the base state. After solving the model under all conditions, the results are displayed in Figures 27-29.  According to these figures, as the demand increases, all three objective functions of the model also increase, and vice versa. In other words, if customer demand for rice increases, the costs, pollutants, and job opportunities created in the supply chain network will also increase. Due to the increasing demand for rice, farmers are forced to produce more crops. As a result, production costs and energy consumption taxes will also increase, and, therefore, the first objective function will increase. In addition, more product will be moved between facilities and more pollutants will be produced. Therefore, the second objective function will also increase. Finally, due to the production of more crops, more job opportunities will be created in rice farms. Thus, the third objective function will also increase. Therefore, we can conclude that demand and all three objective functions of the model have a direct relation.

Conclusions, Managerial Insight, and Future Works
At the present time, rice feeds over half of the world population, which is globally crucial in our food systems. Establishing more efficient and sustainable rice value chains could enable the UN to achieve the Sustainable Development Goal of Zero Hunger by 2050. In order to achieve this goal, it is essential to properly manage supply chains and create a balance between the supply and demand of rice, which can solve many supply problems of this product.
In the present study, a dual-channel CLSC network of a sustainable supply chain was designed for a rice plant while considering energy sources in an uncertain environ-  According to these figures, as the demand increases, all three objective functions of the model also increase, and vice versa. In other words, if customer demand for rice increases, the costs, pollutants, and job opportunities created in the supply chain network will also increase. Due to the increasing demand for rice, farmers are forced to produce more crops. As a result, production costs and energy consumption taxes will also increase, and, therefore, the first objective function will increase. In addition, more product will be moved between facilities and more pollutants will be produced. Therefore, the second objective function will also increase. Finally, due to the production of more crops, more job opportunities will be created in rice farms. Thus, the third objective function will also increase. Therefore, we can conclude that demand and all three objective functions of the model have a direct relation.

Conclusions, Managerial Insight, and Future Works
At the present time, rice feeds over half of the world population, which is globally crucial in our food systems. Establishing more efficient and sustainable rice value chains could enable the UN to achieve the Sustainable Development Goal of Zero Hunger by 2050. In order to achieve this goal, it is essential to properly manage supply chains and create a balance between the supply and demand of rice, which can solve many supply problems of this product.
In the present study, a dual-channel CLSC network of a sustainable supply chain was designed for a rice plant while considering energy sources in an uncertain environ- According to these figures, as the demand increases, all three objective functions of the model also increase, and vice versa. In other words, if customer demand for rice increases, the costs, pollutants, and job opportunities created in the supply chain network will also increase. Due to the increasing demand for rice, farmers are forced to produce more crops. As a result, production costs and energy consumption taxes will also increase, and, therefore, the first objective function will increase. In addition, more product will be moved between facilities and more pollutants will be produced. Therefore, the second objective function will also increase. Finally, due to the production of more crops, more job opportunities will be created in rice farms. Thus, the third objective function will also increase. Therefore, we can conclude that demand and all three objective functions of the model have a direct relation.

Conclusions, Managerial Insight, and Future Works
At the present time, rice feeds over half of the world population, which is globally crucial in our food systems. Establishing more efficient and sustainable rice value chains could enable the UN to achieve the Sustainable Development Goal of Zero Hunger by 2050. In order to achieve this goal, it is essential to properly manage supply chains and create a balance between the supply and demand of rice, which can solve many supply problems of this product.
In the present study, a dual-channel CLSC network of a sustainable supply chain was designed for a rice plant while considering energy sources in an uncertain environment. The network included the producers, the distribution centers, the rice factories, the pharmaceutical industries, the toiletry industries, the markets, the customers, the solar panel sites, the bio-refineries, and the recycling centers. The grid also supplied electricity needed to process rice through bio-refineries and solar panel sites. In this network, the required electrical energy for rice processing was supplied through bio-refineries and the solar panel site as well. Regarding how rice waste as biomass has the potential to generate energy, this study addressed this issue and proposed a mathematical model for optimally using rice waste for energy generation.
Then, the total costs, the quantity of the emitted pollutants, and the fixed and variable job opportunities created in the proposed network were optimized via an MILP model under the uncertainty of cost, supply, and demand. The goal of the proposed model was to determine the optimal quantity of rice production and distribution and its waste among the network facilities.
The costs include the total transportation costs among the network facilities, the production costs for the farmers and rice factories, the product maintenance costs for the distribution centers, and the fixed costs for reopening new distribution centers, the bio-refineries, and the solar panel sites. In the present research, the parameters cost, supply, and demand were assumed uncertain, and fuzzy logic was applied to deal with uncertainty. To solve the proposed model, MOPSO, MOSA, and MOGWO algorithms and a new multi-objective version of the Reptile Search Optimizer called MORSO were used. Then, the mentioned algorithms were validated using LP-metric in small-sized samples. Furthermore, their results and performance were compared based on criteria such as MS, SNS, NPS, MID, and CPU time. Moreover, for confirming the model's validity, ten sample problems in different sizes were designed. In addition, for coming up with the best performance, the parameters of all three algorithms were tuned by the Taguchi method. After solving the model, the derived results were evaluated and these algorithms' performances were analyzed by ANOVA through the interval plots at the confidence level of 95%. A significant statistical difference was observed in terms of the performances of these algorithms, and according to the statistical tests, the MORSO algorithm performed better than other algorithms in terms of NPS, MS, and MID criteria. Furthermore, for more model evaluation, the sensitivity analysis was performed on the key parameters. The numerical results indicated that it is possible to save up to 19% of electricity consumption by constructing the solar panel sites and producing energy from the rice waste.
Considering management, the present study's results can benefit the relevant managers and the countries producing rice. The present study's findings will be highly useful and involve remarkable applications for decision making on opening solar panel sites and bio-refineries for supplying the required electrical energy for rice processing through focusing on the sustainable dimensions. The related results and findings have been provided. The findings of this research in the areas of clean energy consumption and job opportunities can be stated as follows: • Job opportunity: Agricultural growth has led to increasing the productivity and income of small and marginal farmers and raising the employment and wages of workers. Consequently, it is critical to consider this dimension of sustainability in supply chain optimization. The proposed model pursues the goal of raising job opportunities and its results revealed that increasing rice production could boost the intended goal. • Clean energy consumption: Using solar energy as a source of clean energy has gained importance considering the non-renewable nature of fossil sources such as oil and gas. Solar energy and biomass are viewed as a source of clean energy and the cost of generating electricity by them is less than that generated by fossil fuels. Moreover, they emit less pollution and fewer greenhouse gases. Regarding how rice straw and husk have the potential to generate energy, the current research has dealt with this matter and a mathematical model has been formulated to make strategic decisions about the construction of solar panel sites and bio-refineries and optimally using rice waste for energy generation.

Limitations of the Current Study
The current research has many limitations and it seems that different directions can be considered for its development. To develop this study in the future, different methods of rice production and the capacity of energy sources could be considered. Moreover, water resources are of crucial importance in rice production, thus sustainability should be promoted by considering water resources issues in future works. Moreover, considering other uncertainty approaches such as stochastic, probability, and possibilistic is one of the issues that could be emphasized by the researchers in this field and could be incorporated to improve the model. Furthermore, solving the proposed model by the heuristic methods and other metaheuristic algorithms and comparing their results might be interesting. Finally, integrating the proposed model with topics such as the Internet of Things and Industry 4.0 could also be considered by researchers.

Conflicts of Interest:
The authors declare that they have no conflict of interest.

Abbreviations Sets and Indicators: I
Set of producers i ∈ I J 1 Set of existing distribution centers j 1 ∈ J 1 J 2 Set of new distribution centers j 2 ∈ J 2 J Set of all distribution centers j ∈ Fixed cost for launching recycling distribution center o 2 f cb b Fixed cost for launching bio-refinery b prm pmt Cost of purchasing each kg of p-type product from market m (traditional purchase) during period t pr f p f t Cost of purchasing each kg of p-type product from factory f (online shopping) during period t pcap f p f t Production capacity of p-type product in factory f during period t cpa i Cost of production per kg of rice plant for farmer i ctax r Tax on consumption per unit of electricity supplied from source r cpt p f r Production cost per kg of p-type product in factory f, whose electricity is supplied from source r chj jt Cost of per kg unprocessed rice maintenance at distribution center j during period t CT Cost of transporting by vehicle per kilometer disa ij Distance between producer i and distribution center j disb j f Distance between distribution center j and rice factory f disc f v Distance between rice factory f and pharmaceutical industry v disd f k Distance between rice factory f and toiletries industry k dise f m Distance between rice factory f and market m dis f f c Distance between rice factory f and customer c disg mc Distance between market m and customer c dish io Distance between producer i and recycling center o disi oi Distance between recycling center o and producer i disj ib Distance between producer i and bio-refinery b disk f b Distance between rice factory f and bio-refinery b cap it Product's production capacity by producer i during period t capr rt Max energy level supplied from source r during period t dct pct Demand for p-type product by consumer c in traditional purchase during period t dco pct Demand for p-type product by consumer c in online shopping during period t dci it Demand for compost by farm i during period t dv vt Demand for rice bran by pharmaceutical industry v during period t dk kt Demand for rice bran by toiletries industry k during period t ER Electrical energy required to process per kg of rice M A very large positive number η Paddy to husk conversion rate δ Rice husk to electrical energy conversion rate in bio-refineries λ Straw to electrical energy conversion rate in bio-refineries ψ Straw to compost conversion rate cpo o Cost of production per kg of compost at recycling center o πo o CO 2 quantity released from constructing recycling center o πs s CO 2 quantity released from constructing solar panel s πb b CO 2 quantity released from constructing bio-refinery b RF Fuel required for vehicle per travel per k capv Capacity of vehicle β Rice plant to paddy conversion rate 1 − β Rice plant to straw conversion rate θ p Paddy to product of type-p conversion rate

Decision Variables
XSe rt Amount of electrical energy that should be supplied from source r during period t Xhj jt Quantity of paddy rice stored in the distribution center j during period t Xa ijt Quantity of rice plant transferred from producer i to distribution center j during period t Xb j f t Quantity of paddy rice transferred from distribution center j to rice factory f during period t Xc f vt Quantity of rice bran transferred from rice factory f to pharmaceutical industry v during period t Xd f kt Quantity of rice bran transferred from rice factory f to toiletries industry k during period t Xe p f mt Quantity of type-p product transferred from rice factory f to market m during period t X f p f ct Quantity of type-p product transferred from rice factory f to customer c during period t Xg pmct Quantity of type-p product transferred from market r to customer c during period t

Xh iot
Quantity of compost transferred from producer i to recycling center o during period t Xi oit Quantity of compost transferred from recycling center o to producer i during period t Xj ibt Quantity of rice straw transferred from farm i to bio-refinery b during period t Xk f bt Quantity of rice bran transferred from rice factory f to bio-refinery b during period t Xq it Amount of rice plant produced by producer i during period t Xrh f prt Quantity of type-p product produced from rice factory f using source r during period t XSe rt Quantity of electrical energy that should be supplied from source r during period t Variables Zero and One: if recycling center o is established, 1 Appendix A   Table A1. The proposed priority-based decoding procedure of segment 2.
For t = 1 to T Inputs J = set of distribution centers F = set of rice factories D( f , t) = capacity of rice factories f in period t Ca(j, t) = capacity of distribution center j in period t V(J + F) = encode solution of period t Dis(j, f ) : distance between nodes -----------------------------Outputs: yaloc(j, f , t) = amount of shipments between node j and f in period t Em(j, f , t) = amount of CO 2 emission caused by transferring quantity by vehicle between nodes I NV(j, t) = amount of remained goods in distribution center j at period t ---------------Ca(j, t) = Ca(j, t) + I NV(j, t − 1) step1 = yaloc(j, f , t) = 0 j ∈ J, f ∈ F while ∑iCa(j,t) > 0 or ∑iD(f,t) > 0 step2: select the value of first column of first sub-segment J for j index) select the value of first column of first sub-segment F for f index step3 : yaloc(j, f , t) = min(Ca(j, t), D( f , t)) Update demands and capacities Ca(j, t) = Ca(j, t) − yaloc    Table A3. Pseudo-code of RSO algorithm.