Designing a Recycling Network for the Circular Economy of Plastics with Different Multi-Criteria Optimization Approaches

: A growing plastic production increases the pressure on waste management systems, which have to cope with greater volumes of plastic waste. Increased plastics recycling can reduce environmental impacts by lowering the need for primary plastics production and thus fossil resources demand. Current research is mainly focused on identifying environmentally friendly recycling technologies for different waste streams. However, recycling capacities must also be expanded to handle the waste generated. Therefore, this paper develops multiple exemplary multi-criteria optimization models to design an optimal recycling network. The models are deployed in a case study for plastic packaging waste in Europe for an advanced mechanical recycling process. We compare the different multi-criteria optimization approaches, how they balance environmental and economic aspects differently, and how this affects the recycling network design. Finally, we compare the optimization approaches and ﬁnd goal programming the most promising approach for recycling network design that ensures a balance between economic and environmental objectives.


Introduction
In 2019, global plastic production reached almost 370 million tons, and European plastic production reached nearly 58 million tons [1]. Current forecasts indicate that the global annual production of plastic will grow further [2], resulting in an increasing amount of plastic waste that waste management systems must handle. In Europe, almost 40% of plastic is used in the packaging sector, from which only 32.5% was collected for recycling [1]. Primary plastic is produced from crude oil, and the production process is energy demanding, resulting in an intensive environmental footprint [3]. Recycling plastic waste and establishing a circular economy for plastics is a promising way to become more independent of fossil-based resources and reduce the plastic industry's environmental footprint [4]. The European Union (EU) acknowledges this in its legislation and demands a recycling rate for plastic packaging waste of 50% in 2025 and 55% in 2030 [5]. Moreover, improving the circular economy of plastics is part of the current research focusing on the design for recycling [6], advanced mechanical recycling [7], or alternative recycling technologies such as chemical recycling [8].
Multiple studies have assessed current and evolving recycling technologies for plastic packaging waste and compared them regarding economic and environmental indicators. Neo et al. [9] provide life cycle assessments for end-of-life (EoL) plastic waste options in India and Indonesia, including mechanical recycling, co-processing in cement kilns, incineration, landfilling, open dumping, and open burning. Faraca et al. [10] provide an overview of the economic and environmental performance of simple and advanced mechanical recycling processes and a chemical recycling process for plastic waste in Denmark. Civancik-Uslu et al. [11] compare the environmental performance

Materials and Methods
This section provides a short overview of mathematical optimization approaches for designing a recycling network, introduces multi-objective decision-making, and describes the approaches implemented in the case study.

Designing Reverse Logistic Networks
Reverse logistic networks are designed to optimize the recovery or reuse of end-of-life resources in supply chains [17]. Recycling activities are usually a complex recovery process consisting of various collection, sorting, and processing stages that must be coordinated. Resources that can not be recycled must be disposed of by either incinerating or landfilling them. Part of the problem of designing reverse logistic networks can be solved with a respective facility location problem. Here, the objective is to find the optimal placement of facilities in a network to minimize costs, environmental burdens, or both. Mathematical models of facility location planning have been extensively studied [25]. Further extensions of the facility location problem such as capacities, stochastic modeling, or dynamic modeling are provided, e.g., by Nickel et al. [19].
Nickel et al. [19] also give an introduction and overview of multi-criteria optimization. Multi-criteria optimization allows the assessment of multiple objectives simultaneously. Multi-criteria decisions can be separated into Multi-Attribute Decision-Making (MADM) and Multi-Objective Decision-Making (MODM) [18]. With MODM, the optimal solution is determined by mathematical modeling, whereas with MADM, a solution is selected from a predefined set of alternatives [18]. This study focuses on MODM solution approaches, as the design of the recycling network has a high number of possible solutions [26].
MODM approaches describe the solution space with mathematical models, whereby all solutions are permissible that fulfill specified constraints [18]. Several quantifiable objective functions describe the decision-maker's objective system [27]. Constraints restrict the solution space of the problem. Solving MODM problems includes investigating conflicting objectives [28]. The ideal solution consists of the optimum value for each objective, while an efficient solution has an objective that cannot be improved without a simultaneous detriment of another objective [29]. These efficient solutions are also defined as Pareto-optimal or Pareto-efficient [30]. All Pareto-efficient solutions form the Pareto front of the problem. The computation of this Pareto front is resource-consuming but helps decision-makers find the preferred solution. Based on the preference order of the decision-maker, approaches exist that reduce the complexity of calculating one optimal solution [18]. According to Colapinto et al. [31], goal programming algorithms are popular, while Farahani et al. [28] rather propose lexicographic methods. In the following, these two approaches are introduced as well as implemented and compared in this study.

Lexicographic Optimization
The lexicographic optimization is directly based on the preference order of the decisionmaker. The objectives are considered sequentially in a predefined order determined by the decision-maker [19]. The solution procedure starts with the optimization of the objective with the highest priority, while subordinate objectives are not considered in this first step [19]. If there are multiple optimal solutions for the superordinate objective, the next subordinate objective is optimized based on this set of solutions [19]. This procedure is repeated until there is one solution left or all objectives have been optimized [19]. Objective suppression can occur since some objective functions could be disregarded [19].

Goal Programming
The method of goal programming assumes that the decision-maker has a target value for each objective value [19]. A solution is calculated by minimizing the deviation of the objective value from the target value [19]. Three main goal-programming approaches are differentiated by metric and are suitable for different types of optimization problems: lexicographic, weighted, and Chebyshev goal programming [32]. In the lexicographic goal-programming approach, a series of sequential optimizations are performed to solve the problem [32]. Due to the formulation with priority levels, and thereby its relation to lexicographic optimization, the approach is not pursued further here. In weighted goal programming, direct trade-offs of the undesirable variances are performed by a weighted, normalized single-target attainment function. The weights are assigned to the deviation variables that generally require a reasoned idea about the weighting. Another disadvantage of this approach is that the sum of the variances is minimized, which may result in an objective being weighted less if its variance has a minor impact on the objective function [32]. The min-max approach of Flavell [33], also referred to as Chebyshev Goal Programming (CGP), uses a maximum norm. The norm minimizes the maximum deviation from each goal rather than the sum of all deviations. Therefore, this approach is suitable when a good balance of all objectives is to be achieved. Furthermore, the CPG approach is the only goal programming approach of the three presented that finds solutions that are not at extreme points. Thus, the CPG approach is selected out of the three goal programming approaches to be compared to the lexicographic approach. Moreover, the CGP approach generally leads to a trade-off between objectives in contrast to the lexicographic approach, where prioritized objectives usually dominate the solution.

Case Study: Plastic Recycling Network for Europe
The approaches to designing a reverse logistics network are investigated based on a European case study to recycle post-consumer plastic packaging waste. In Europe, there are considerable differences in the recycling, landfilling, and energy recovery of plastic packaging waste between countries [34]. Despite these differences in handling plastic waste, this case study assumes that there are no regional differences in the implemented sorting and recycling steps and that the organization of the waste handling systems is the same throughout Europe. This assumption describes an idealized and standardized waste management system within Europe and a European circular economy. It is also justifiable since this study focuses on implementing a reverse logistic network for waste treatment with MODM approaches and comparing solution approaches regarding trade-offs between economic and environmental objectives. The assessed network includes waste generation within European regions, placement of recycling sites to treat the plastic packaging waste, and transportation of waste and the recycled material to processing and demand sites, respectively. Transportation is limited to direct and single supply; tour planning is not included. All materials are transported by lorry, other forms of transportation are neglected. Exports to other world regions are also not considered, as the vision of a European circular economy in which all European plastic waste is also treated in Europe is modeled.
The plastic packaging waste within the European regions is forwarded to a material recovery facility (MRF). Here, the different plastic types are separated, ground, washed, and pressed to pellets [10]. Separated metal scrap is forwarded to respective metal production processes (e.g., adhering aluminum), and remaining sorting residues are incinerated to recover electricity and heat [10]. The plastic pellets are input for plastic converters in Europe that produce new plastic products to close the material cycle. Costs and environmental burdens are associated with opening and operating MRFs, recycling plastic waste, and transporting waste and recycled plastic. The assessed indicators in this study are recycling costs and a subset of standard life cycle assessment indicators found in the literature: global warming potential (GWP), terrestrial acidification (TA), terrestrial eutrophication (ET), and human toxicity (HTc). The indicators chosen as data for the selected process are provided by Faraca et al. [10].

Mathematical Formulation
The case study models a static and deterministic facility location problem and aims to place two facility types of the same technology but with different waste handling capacities: a large central facility and a smaller decentral facility. Thus, it is a capacitated two-stage reverse logistics model with multi-sourcing and direct delivery without interaction between facilities. The model is formulated as a mixed-integer problem (MIP). Table 1 provides an overview of the problem's sets, parameters, and variables.    Four sets are distinguished: First, set S of regions with recyclable waste contains all regions (NUTS 2) where the plastic waste in Europe is generated. Second, R is a set of regions where MRFs can potentially be placed. There are no limitations where MRFs can be opened, so this set also contains every considered region. Third, set D contains all regions where demand for recycled material exists. Further information on recyclable wastes and demand per region are described in Section 3.2. Forth, a set of input capacity classes for the MRFs is introduced. We distinguish between two capacity classes, and the models are formulated so that only one kind of MRF can be opened in one region at a time. The waste generated and its demand for recycled material are captured within the respective parameters for each region (see Sections 3.2.1 and 3.2.1). The costs or environmental burdens of the transportation, recycling waste, and operating an MRF, as well as building an MRF, are parameters within the model (see Sections 3.

and 3.2.4).
Costs and environmental burdens of opening and recycling plastic waste vary depending on the capacity class. The recycling costs of the material are also region-specific, especially as labor costs and labor cost levels differ throughout Europe. An environmental impact differentiation based on the different regions (such as applied by Schiessl et al. [35], Müller et al. [36]) is not conducted as the same MRF with the same recycling process is opened throughout Europe. Decision variables of the model are the amount of transported material from supply locations to MRFs (x sr,i ) and from MRFs to the demand locations (x rd,i ), and the binary variable y r,i indicates if an MRF of capacity i is placed in a region. A cost-minimizing MIP can be formulated using the sets, parameters, and variables within Equations (1) to (8): This basic formulation of the model includes transport costs between s, r, and d, region-specific recycling costs at the MRF (includes recycling costs of the plastic waste and operating costs of the MRF), and costs associated with opening an MRF in a region in its objective function (1). The total costs of the recycling network are minimized.
The objective function (1) is subject to constraint (2) that all waste within a region must be handled within the network. Constraint (3) describes the material flow conservation of the network, defining that all the waste entering an MRF must be handled and associated products must be forwarded to the demand points. Therein, a plastics recovery rate of 67% of the input material is assumed [10]. Constraint (4) ensures that no demand point is oversupplied. Each MRF has a capacity (decentral or central) that must not be exceeded (Constraints (5)). Constraint (5) also ensurses that plastic waste is only forwarded to regions where an MRF has been placed. Constraint (6) ensures that only one facility is placed per region. Suppose this facility does not provide the capacity to handle all plastic waste in that region, then the remaining plastic waste must be exported to other European regions. Constraints (7) and (8) define the range of values of the decision variables.
In the following subsections, the one-dimensional formalization of the problem with a single objective is extended to a multi-objective linear program, including the assessment of environmental indicators.

Lexicographic Optimization
Additional objective functions for the assessed environmental indicators GWP, TA, ET, and HTc are formulated to include the minimization of environmental burdens in the model ( (9) to (12)). Each objective function follows the same structure as the cost objective function (1). Environmental burdens for transporting material between different stages of the network, an environmental burden associated with material recycling, and an environmental burden associated with the opening of MRFs of different capacities are summed up to a total environmental impact per indicator.
These objective functions are considered sequentially in a predefined order, which is assumed here to start with cost minimization. However, any other order is likewise possible. The objective functions are then considered in the order listed. The network's GWP objective (9) is minimized within all possible solutions with minimal recycling costs. Following that, the TA objective (10), ET objective (11), and HTC objective (12) are considered. In a slight variation of the lexicographic approach, we allow a deviation ( = 0.01) from the optimal objective function values to relax the lexicographic constraint and enable more flexibility. Thus, the feasible solution space can deviate from an optimal solution by 1%.

Goal Programming
As described above, CGP follows a satisfying approach, minimizing the maximum deviation from each goal and aiming at a balance of the achievement of each objective instead of suppressed objectives. Therefore, an adapted MIP is formulated with the Equations (13)- (25): The continuous, non-negative auxiliary variable dist is introduced in the CGP formulation. It measures the distance of each solution to the optimal values of the single-criteria optimization and should be minimized (13). Constraints (14) to (18) do not differ from the former model formulation. Constraints (19) to (23) represent the distance measurement, where the distances between objective values are normalized with a percentage normalization. Consequently, the variable dist is interpreted as the relative deviation (in percent) of the objective value from the optimal value. Constraints (24) and (25) define the range of values of the decision variables, including dist.

Input Data
The input data for the case study include the waste volume in different European regions, the demand for recycled plastic pellets for the production of plastic products, costs, as well as environmental burdens associated with the transport and recycling of the plastic waste and the opening of MRFs. Furthermore, the distances between the different regions of the network are given.

Plastic Waste Volume
In Europe, the plastic packaging waste generation per capita, as well as the residents for each NUTS 2 region are provided by Eurostat [34,37]. With this information, the European waste volume of post-consumer plastic packaging waste is regionalized, multiplying the number of residents of each region with the plastic packaging waste generation per capita. The calculation of the total plastic packaging waste is based on the data for the year 2019.
For six countries, no current data was available. In these cases, data from culturally and regionally consistent countries was adopted. Table 2 shows this relationship for countries with missing data.  Figure 1 provides an overview of the plastic waste volume throughout Europe. Some metropolitan regions such as Paris and Milan are clearly visible, while in south-western Europe, more plastic waste arises than in central and northern Europe.

Recycled Plastic Demand
The plastic demand per region is not available in official statistics. However, the total plastic demand of each European country for 2019 is available [1]. The total plastic demand of each country is assigned to its NUTS 2 regions based on their number of residents.
The plastic demand data are not available for every single European country. For six countries, the data are not available. Therefore, the per capita demand from culturally and regionally comparable countries was adopted (Table 3).  An overview of the regionalized plastic demand is provided in Figure 2. Comparing  Figures 1 and 2, the plastic demand is similar to the waste generation in the respective regions since both datasets are based on the number of residents within each NUTS 2 region. The recycling and operating costs are based on the advanced mechanical recycling process introduced by Faraca et al. [10]. They describe the mechanical recycling of the plastic waste in an MRF similar to Cimpan et al. [38]. The gross costs for recycling 1 Mg of post-consumer plastic waste include recycling the material and operating the MRF (Faraca et al. [10]). Based on a size degression with a coefficient of 0.68 (derived from Cimpan et al. [38] and displayed in OPEX and CAPEX) and adjusted to 2019, the recycling costs in a decentral MRF are 415 EUR/FU, and in a central MRF, they are 267 EUR/FU. According to [38], 65% of the recycling costs are associated with labor costs. Therefore, we split the recycling costs provided by Faraca et al. [10] into a part associated with labor and a fixed part. Labor costs are regionalized by data for the labor cost levels in Europe [39].
Transportation costs are based on an Ecoinvent 3.8 dataset [40] and are adjusted for inflation (base year 2019). Therefore, the transportation of 1 Mg for 1 km results in average transportation costs of 0.174 EUR. All costs are summarized in Table 4.

Environmental Burdens
Environmental burdens are associated with the opening of MRFs, the recycling process itself, and waste and recycled material transportation. The environmental burdens associated with the opening of an MRF are based on an Ecoinvent 3.8 dataset [41] and the Environmental Footprint (EF) method of the European Commission [42]. Impacts of opening MRFs are scaled with a size degression coefficient of 0.8 for both capacity classes. Furthermore, Faraca et al. [10] provide environmental burdens for the recycling of 1 Mg plastic waste that is used here. Consistent with the recycling and operating costs, the environmental impacts of recycling and operating are scaled with the size degression factor of 0.68. All values are gross values, and credits for substituting primary plastic are not considered. The Ecoinvent 3.8 transportation data are considered for environmental burdens of the transport [40]. All burdens for the considered categories GWP, TA, ET, and HTc are summarized in Table 4.

Distances
Each European region is represented by its point with the highest population density. For modeling simplicity, it is assumed that all waste and demand for plastic pellets accumulate at this point. Additionally, if an MRF is opened in a region, it is also located at this point. We choose the point of the highest population density since it is likely that waste accumulation, plastic demand, and recycling infrastructure are associated with this region. Transport distances between these points are calculated using the Euclidean norm. However, the assumption of direct transport is not applicable when considering overland routes. Therefore, a tortuosity factor of √ 2 is introduced to embrace real road conditions [43,44].

Implementation
The model was implemented in Python. The Python API of CPLEX was used for the solution of the described model. Python 3.8.5 in the Anaconda distribution is installed on the machine used (Intel ® Core™ i5-8265U at 1.60 GHz; 8.0 GB RAM). The integration of the CPLEX solver into Python is accomplished via the docplex library. When solving the formulated programs, a termination criterion with a selected time limit of 600 s is set for each optimization run to achieve comparability.

Results
The results of the case study are presented in this section. First, the optimal objective values for single objective optimization are presented. They are the basis for determining relative deviations between multi-criteria optimization solutions and individual optimality for each objective. Furthermore, they are needed to normalize the dist variable in the CGP approach. Second, the results of the lexicographic optimization are provided. We calculated two approaches: one without deviation from the optimal solution and one where a small deviation of = 0.01 is allowed. Third, the results of the introduced CGP approach are presented. A comparison and discussion of the results follow in Section 5.

Single Objective Optimization
Optimal solution values of the presented cost-minimizing optimization problem for all indicators, including cost, GWP, TA, ET, and HTc, as well as the relative shares of the recycling effort, the opening of the MRFs, and transport are presented in Table 5. GWP and HTC optimization are dominated by the recycling effort, with shares of 77% and over 99% in the optimal single objective solutions. Regarding costs, TA, and ET, there is a balance between the impact of the recycling process and opening MRFs. For all indicators, transportation of the waste and plastic pallets has a relatively small impact, with a maximum share of 8% within the TA indicator.

Lexicographic Optimization
The lexicographic optimization approach in this study considers the objective function in the following sequence: costs, GWP, TA, ET, and HTc. The lexicographic approach finds a solution that reaches the optimal value for total costs as the costs are the prioritized objective (see Table 6). Minimizing the GWP within the remaining solution space results in a deviation of 33% from the optimal solution obtained in the single objective optimization. As expected, the deviation further increases for the following two objectives to 64% (TA) and 257% (ET). For the HTc impact, there is only a deviation of 11% from its optimal value from the single objective optimization, as there is no strong conflict with the cost objective.   A slight deviation ( = 0.01 = 1%) from the lexicographic objective function values is introduced to soften the lexicographic constraint and enable more flexibility in the solution approach. The results for this approach are presented in Table 7 and show a reduction in the deviation from the single objective optimal value of the GWP, TA, and ET compared to the strict lexicographic approach above. The cost optimization uses = 0.01, leading to a reduction in the deviation for GWP, TA, ET, and HTc by 7, 13, 124, and 3 percentage points, respectively. Relaxing the optimality for the single objectives can lead to a more balanced result.  Figure 4 displays the recycling network for the lexicographic solution with = 0.01. Here, and compared to the strict lexicographic solution, fewer decentral MRFs with smaller capacities are opened while the number of central MRFs only increases by four. They are placed in Western Europe, replacing decentral MRFs in The Netherlands. The layout also decreases the transport distance.
A sensitivity analysis ( Figure 5) is performed to identify the effects of parameter variations on the objective values of the different indicators. In this case, all cost and environmental impact parameters of the recycling process, opening facilities, and transportation are varied separately within their corresponding parameter groups (recycling, opening, and transport).
The analysis reveals the following dependencies between total cost and the individual parameter groups: a 10% reduction of all recycling costs results in a 5% decrease in total costs, a 10% reduction of all opening costs causes a 4% decrease in total costs, while a 10% reduction of transportation costs has no influence on the total costs. These results and their linear relationship show that recycling costs significantly impact the cost target.
In contrast, the results of the environmental impact categories do not demonstrate a linear correlation as the objectives are developed based on the previous cost optimization. The sequential procedure leads to contingent shifts in the solution space for the environmental indicators that can allow for more flexibility and aid in aligning the results with each optimum. Despite these non-linear dependencies, certain trends regarding the impacts of parameter variations on GWP, TA, ET, and HTc can be identified.
Concerning GWP, the recycling process is found to have the most significant impact, causing a 7% decrease in total GWP when burdens are reduced by 10%. In contrast, varying the GWP impact for facility openings does not influence the total GWP at all and reducing the value associated with transportation by 10% causes only a 3% decrease in total GWP. Transportation does, however, considerably influence the TA and ET impact categories. A 10% decrease in the parameters is found to cause a 4% (TA) and 5% (ET) reduction of total impact, while a 10% increase results in a 6% (TA) and 8% (ET) increase. In contrast, the final environmental indicator, HTc, is shown to be mostly influenced by the recycling process, with transportation having no impact.
The analysis results indicate competing objectives where cost and environmental impacts of GWP, TA, and ET are concerned. The softened lexicographic approach is, however, inherently incapable of calculating a more balanced solution, which considers both environmental and economic indicators. For this reason, an alternative method, namely the CGP approach, is implemented and presented in the next chapter.

Goal Programming
Solving the introduced optimization problem with the CGP approach leads to objective values presented in Table 8 and a recycling network displayed in Figure 5. Here, many decentral MRFs are opened to avoid heavy transport activity. Transportation is limited to short distances between neighboring regions. Central MRFs with higher capacity are placed in regions with very high waste volumes when a decentral MRF's capacity is insufficient for handling all the plastic packaging waste in that region. Table 8 and Figure 6 show different results for both lexicographic approaches. Less central MRFs are placed throughout Europe. Instead, multiple decentral MRFs reduce the transport activities. The difference in results can also be seen in deviations of the objectives from their optimal single objective value. Within the CGP approach, the costs deviate from their optimal value by 23%, and the HTc impact deviates by 23% from its optimal value. However, compared to the lexicographic approaches, the impacts of GWP (13% vs. 26%), TA (5% vs. 51%), and ET (23% vs. 133%) deviate less from their optimal values. Overall, the CGP approach reduces the maximum deviation from 133% in the adapted lexicographic approach to 23%. The CGP approach reaches a balance between the objectives, satisfying them with the objective value dist = 0.23. A sensitivity analysis is performed to vary the recycling, opening, and transport group parameters (Figure 7), illustrating their effects on the total cost and environmental impacts of the recycling network. In contrast to the previously applied lexicographic approach, utilizing the CGP method allows for linear relationships between all parameters and indicators. The results indicate that the total costs of the recycling network are mainly influenced by the recycling costs themselves. A 10% reduction in recycling costs leads to a 6% overall decrease, while the same reduction in opening or transportation costs causes only a 3% and 1% in total cost decrease, respectively.
These numbers remain consistent for all indicators as the parameters are sorted into their respective groups, and both cost and environmental impacts are varied simultaneously within said groups. It can, therefore, be concluded that the recycling process has the most significant impact overall. In addition, the opening of a facility plays a more consequential role in the examined indicators than the transport of materials.
Overall, the results from both the CGP approach and lexicographic optimization show that a trade-off occurs between the cost and environmental indicators and all individual environmental indicators in the design of plastic recycling networks within a European circular economy.

Discussion
The discussion is separated into the case study discussion (Section 5.1), including the results, assumptions made, and the set time limit of 600 s for each optimization run, and into the discussion of the selected MODM approaches (Section 5.2).

Case Study
The case study results are calculated with a time limit of 600 s. The limitation of the computation time is introduced to compare both MODM approaches. The results of the one-dimensional optimization are only slightly affected since optimal solutions with MIP gaps of <0.28% were found in all optimization runs. The time limit of 600 s was used in the MODM approaches. MIP gaps of up to 5% are reached in the single optimization runs for the different subordinate objectives in the strict lexicographic formulation. By allowing a slight deviation from the optimal values ( = 0.01), the MIP gaps are reduced to a maximum of 2%. In the CGP approach, a MIP gap of 1% is achieved.
The lexicographic approach and the implemented goal programming lead to different network layouts. The lexicographic approach generates a network with more central MRFs opened. Here, the cost objective dominates the solution process. There are long transport distances, and fewer MRFs are opened, especially in less populated regions (Scandinavia, Bulgaria, Romania). The goal programming approach generates a network with fewer central MRFs and multiple decentral MRFs, reducing transport distances. An MRF is opened in nearly every region, and transport activities are limited to neighboring regions. For the environmental objectives GWP, TA, and ET, the focus is on minimizing transport, which conflicts with the cost objective. Due to this conflict of objectives, the different MODM approaches result in significantly different recycling network designs.
The assumptions of the case study and their impact on the results need to be critically discussed. The current model is limited to one stage, and only one type of facility is considered where the plastic waste is processed. The MRFs described by Faraca et al. [10] obtain the input material from recycling centers where the post-consumer waste is already processed, and the plastic fraction is forwarded. The sorting steps and efforts upstream of the MRF are not considered in the current model. Other waste handling options such as landfilling or energy recovery are not considered to reduce the complexity. Further, existing waste handling facilities are not considered and included in the design of the recycling network. It is assumed that the same type of MRF is placed throughout Europe. The incineration of sorting residues and recycling of metals is included in the recycling burdens [10], but it is assumed that all processes are part of the MRF, and no transport to other facilities is needed. The demand for the by-products (metals and energy) is not included in the location decision. The different plastic material fractions separated in the MRFs are considered as one material stream forwarded to the plastic converters since the plastic demand per European country is not available by polymer. The MRFs are limited to two capacity classes. Other capacities are neglected, and the optimization is limited to costs and the introduced environmental burdens. Transportation is limited to land, specifically road use by means of heavy trucks or lorries. Other modes such as water transport are neglected. Although they would be eligible within the model, both sensitivity analyses show the low influence of transport on the overall objectives, so modeling was refrained from.
Despite the assumptions made and the model's limitations, it can be concluded from the case study that the economic objective and the various environmental objectives do compete with each other. However, The MODM approaches show that a good balance between economic interests and environmental objectives can be achieved.
Future research could adopt the model for case studies in other regions, including additional constraints, products, or waste streams (e.g., post-production waste). The optimization could be expanded to other environmental burdens, such as other impact indicators from life cycle assessment or material or carbon recycling quotas. Additionally, a three-or general multi-stage model with more facility types would allow further differentiation of recycling/production steps. Moreover, the model definition could include the existing recycling or waste handling infrastructure. Furthermore, transportation efforts and distance calculation could be calculated more precisely. Further, regions are represented by their point of highest population density. Using data on the European road network and including locations on the network would allow a more precise consideration of transportation distances.
The model also does not consider importing and exporting plastic waste and plastic pellets to regions outside Europe. This does not reflect the current reality. However, the modeled recycling network can be considered a model for a European circular economy for plastics independent of plastic waste exports. The model's constraints reflect this.

Multi-Objective Decision-Making Approaches
The lexicographic approach considers the objectives in a predefined order determined by the decision-maker. In this case study, the objectives are considered in the following order: costs, GWP, TA, ET, and HTc. Other sequences of the objectives are not considered. However, both lexicographic formulations lead to the dominance of the cost objective. Therefore, it can be concluded that a different objective order results in a solution close to the solution of the single optimization of the first objective.
The modification of the lexicographic approach allowing a slight deviation from the optimal objective function values, changes the standard approach found in the literature. The slight deviation of ( = 0.01) is introduced to soften the strict lexicographic constraint and enable more flexibility in the solution approach. After adaption, the lexicographic approach achieves a more balanced solution than its strict formulation. However, the cost objective still dominates the solution, and the adapted version is not as balanced as the solution generated in the CGP approach. There is no variation of since the comparison between the lexicographic approach, and CGP is the focus of this study.
The lexicographic approach leads to higher deviations from the single objective solutions, while the CGP approach leads to a balance between the objectives. In general, the solutions show the expected characteristics mentioned in the literature. Therefore, CGP is the more suitable MODM approach for modeling a recycling network if economic and environmental trade-offs are to be considered in the design of the network and if economic and environmental indicators are to be considered equally. However, CGP is a satisfying approach, and no optimality of the individual goals is guaranteed.
The comparison of MODM approaches is limited to the lexicographic approach and goal programming. Other approaches, such as the scalarization approach, are not considered. In contrast to the CGP, the scalarization approach allows the calculation of optimality [19]. However, the scalarization approach needs weighting factors for every objective in the objective function. These weighting factors usually reflect the subjective opinion of decision-makers. Therefore, this study did not include the scalarization approach.

Conclusions
This study focuses on the recycling network planning for plastic packaging waste in Europe, including opening MRFs, transporting waste to the MRFs, and transporting plastic pellets from the MRFs to plastic converters. An MODM problem arises as costs and environmental burdens of GHG, TA, ET, and HTc are included in the network planning. We use a lexicographic approach, an adapted lexicographic approach with an = 0.01 allowing slight deviations from the optimal value, and a CGP approach to solve the problem. Comparing the results, the strict formulation of the lexicographical approach optimizes the main objective while subordinate goals are optimized in the remaining solution space. This leads to high deviations from the single objective optimizations' optimal values. Allowing a slight deviation from optimality leads to a more balanced solution; however, there are still high deviations from the single objective optimizations. The CGP approach leads to the most balanced solution since economic and environmental indicators are considered equally in the design process. Therefore, the CGP approach is most suitable for designing a recycling network if economic and environmental trade-offs are to be considered in the design of the network.
The high objective function value of the CGP in the formulated case study indicates that the economic and environmental objectives compete within the chosen model design and the made assumptions. For the cost objective, a relatively small number of central MRFs with high capacity are preferable to exploit economies of scale. Combining central and decentral MRFs to limit transport activities is favorable for the environmental indicators GWP, TA, and ET. This results in two different recycling networks following either the lexicographic or the CGP approach.
The models presented can be extended to gain a more detailed overview of the economic and environmental impacts of a European recycling network designed from scratch. The CGP approach should be pursued to ensure a balance between the different impact indicators, especially since they appear contradictory.