Vehicles Allocation for Fruit Distribution Considering CO 2 Emissions and Decisions on Subcontracting

An important problem in rural-area supply chains is how to transport the harvested fruit to urban areas. Lowand medium-capacity vehicles are used in Colombia to carry out this activity. Operating them comes with an inherent cost and generates carbon emissions. Normally, minimizing operating costs and minimizing carbon emissions are conflicting objectives to allocate such vehicles efficiently in any of the supply chain echelons. We designed a multi-objective mixed-integer programming model to address this problem and solved it via the ε-constraint method. It includes decisions mainly about quantities of fruit to transport and store, types of vehicles to allocate according to their capacities, CO2 emission levels of these vehicles, and subcontracting on the collection process. The main results show two schedules for allocating the vehicles, showing minimum and maximum CO2 emissions. Minimum CO2 emissions scheme require subcontracting and the maximum CO2 scheme does not. Then, a Pareto frontier shows that CO2 emissions level are inversely proportional to total management cost for different scenarios in which fruit supply was modified.


Introduction
Fresh foods such as fruits are highly perishable, which means their shelf life is limited over time.In addition, unlike other products, their quality decreases continuously during the activities of the supply chain.Therefore, the process of dealing with food supply chain networks is complicated, with respect to other types of products that flow in this type of networks [1].These products are easily affected by environmental factors and rot in the process of circulation and flow through the supply chain [2].Given that perishable products reduce its value over time, one of the main objectives in the distribution of these products, especially in the food supply chains, is to attend to their freshness while delivering them to the demand areas.This objective creates a direct effect on the response capacity of the network.Considering these facts, choosing an advantageous distribution network is a key factor for the logistics system administrator.Therefore, food supply chains focus on the quality of the products and the minimization of shipping times or the maximization of the quality of the products in the delivery time, especially due to the challenges associated with seasonality, supply peaks, long delivery times and perishability [3].
In addition to transport costs and response capacity in the optimization of a logistics system, the sustainable characteristics of the network must also be considered [1].The importance of the environment is an issue that has been gaining strength for years in global production schemes.Due to the observation of the consequences of the use of non-environmentally friendly technologies, companies must go beyond social and volunteer environmental initiatives [4].One of the biggest challenges for sustainable business models is to achieve a scale of operations that is adequate to satisfy the quantity and depth of needs in their markets [5].The daily growth and development of distribution companies has increased the amount of greenhouse gas (GHG) emissions, especially carbon dioxide (CO 2 ).According to Cai et al. [6], there is a link between the consumption of clean energy, economic growth and CO 2 emissions.For Wang [7], it is necessary to create energies that do not emit or rarely emit greenhouse gases to favor a low carbon economy.To achieve this type of economy, the supply chain represents one of the main opportunities to improve the development of energy efficiency, even for small companies [8].In the supply chains, some alternative fuels such as biodiesel can be implemented, thus becoming drivers of the use of sustainable energy [9].Along these lines, impact assessments of carbon emission mitigation policies require appropriate emission models and indicators of transport logistics activity in supply chains [10].Therefore, there is a need to integrate environmentally sound options in the consumption of clean energy in the research and practice of supply chain management [11] to finally make them green and competitive, based on the reduction of ecological impact of core activities, cost reduction, quality, performance and energy efficiency [12].
Nearly 86% of CO 2 emissions are from fossil fuels [13] and 40% of GHG emissions up to 2010 were attributed to transport.Of these, carbon emissions from road transport represented 70% [14,15].Wang et al. [16] mentioned that supply chains are responsible for more than 20% of global GHG emissions.Therefore, the effects on the environment of vehicles distributing products must be considered to design a green supply distribution network.This can benefit the company, since it has been proved that companies that implement strategies to improve the environment can increase their product sales and market share [17].
Compared to a traditional supply chain, a chain designed to have lower carbon emissions faces additional challenges in its development.Examples of these problems can be high investment in reducing emissions when changing technologies or coordinating deliveries to maximize the supply chain's profits, especially if they are perishable products.The implementation of sustainability practices has spread, and studies have analyzed their impact on performance, that is, efficiency, flexibility and responsiveness [18,19].According to Bourlakis et al. [20], these practices are aimed at protecting the environment to respond to external pressures, which is why these chains need to be managed efficiently and ecologically to improve environmental and social aspects [21].
In the light of the United Nations Framework Convention on Climate Change, Colombia has voluntarily committed to reducing GHG emissions by 20% by 2030.Hence, the country is positioned as an ambitious and proactive country in these issues, which constitutes an advance towards the fulfillment of the Sustainable Development Goals (SDGs) [22].The country has estimated that land transport (light and freight vehicles) constitutes the main source of GHG emission due to the consumption of diesel and gasoline.The nationwide distribution of emissions is determined by the distribution of population and by the intensity of commercial and industrial activity.For the capital of the Republic, Bogotá, it has been estimated that the consumption of fossil fuels in transport generates 10.59 million tons of CO 2 .The mobility dynamics of the country's capital means that 45% of emissions are generated by land transport, specifically by the transport of cargo and passengers in public transport.In this case, the production of CO 2 , generated because of land transport in Bogotá, is a problem that affects the environment by increasing the greenhouse effect as well as affecting the cost of travel in land transport including the design of the transport network.This increase in emissions is due to the precarious transport infrastructure in remote areas branching from main roads.
The main advantage of considering the calculation of CO 2 emissions in the design of a fruit supply chain to Bogotá is the implementation of measures to reduce and compensate emissions in the design decisions on the collection routes of this product [23].These strategies to transport perishable products generate a higher complexity for the administrator of the supply chain because these products must be transported quickly towards the wholesale markets to avoid their deterioration.That perishability condition prevents producers in rural areas to store these products for a long time [24].To guarantee that products arrive in adequate conditions, it is necessary to identify production zones as well as access routes to schedule times for collection, design routes for transportation, identify collection points, systematize the information to deliver products on time, compare different types of transportation [25] and determine costs associated with reducing emissions released in the transport of perishable products.This article deals with a supply chain of collection and distribution of perishable products in rural areas near the city of Bogotá.The proposed solution method seeks to establish an optimal allocation of vehicles that minimizes the total cost of transport and the CO 2 emissions generated by the supply chain.The methodological proposal has been implemented with real data from an Association of fruit producers in the department of Cundinamarca, Colombia.We addressed a tactical-operational problem and solved it through a mixed-integer programming model.
In the framework of the described context, this article seeks to make contributions aimed at improving supply chains in the environmental and performance aspects for perishable products.The main intended contributions are listed below.Firstly, we analyze a supply chain for perishable products and its reduction of carbon emissions.This paper addresses the theoretical and applied study of the impacts generated by the reduction of carbon emissions in tactical-operational decisions.We accomplish this by studying five scenarios, which include different types of vehicles as well as different forms of transport service acquisition (personal vs. subcontracted vehicles).Different ways of acquiring the service are included in the modeling because it is a practice currently carried out by different fruit producing associations in the study area.Secondly, we analyze whether the supply chain considered in this article achieves lower costs and lower CO 2 emissions for different variations of fruit supply considering a tactical problem of operations that generates a weekly plan of allocation of vehicles.
The rest of the paper is organized as follows.Section 2 reviews the relevant literature.Section 3 presents the methodology used to solve the problem and the formulation of the mathematical model that represents it.Section 4 uses an applied example of a real case of fruit distribution to illustrate the theoretical results.Section 5 analyzes the results, regarding to impacts of variations in fruit supply under normal conditions and related to CO 2 emissions in distribution route planning.The conclusions are presented in the final section.

Literature Review
The research literature related to the work presented in this article has two aspects: studies of CO 2 reduction in the distribution of perishable products in Colombia and supply chain models with restrictions on the reduction of carbon emissions on a general level.

CO 2 Reduction in Supply Chains
The relationships among energy consumption, CO 2 emissions and macroeconomic variables, in particular the correlations among energy consumption, CO 2 emissions and GDP growth, have attracted many researchers [26].To investigate the link between clean energy consumption, economic growth and CO 2 emissions, Cai et al. [6] used an ARDL Bootstrap bound test in G7 countries.Their results show that, to achieve a low carbon economy, it is necessary to invest in the consumption of renewable energy.Helgesen et al. [27] linked a bottom-up energy system model (TIMES) and a top-down computable general equilibrium model (REMES) and indicated a reduction of CO 2 emissions from transport in Norway in 2030 to 50% compared to the CO 2 emissions produced in 1990.Dente and Tavasszy [10] suggested the use of the ton-kilometer indicator to measure logistics activity and emissions.In general, the ton-kilometer indicator is used to measure logistics activity, and emissions are calculated as a linear function of this indicator.
Supply chain management of perishable products is different from the supply chain management for articles that have a relatively longer shelf-life [28].
As perishable food supply chains become more complex, incidents of contamination increase [29].Some of these incidents are related to non-optimized management during supply chain processes [30].At present, mathematical modeling proposes to optimize certain processes for these supply chains and reduce the decay of perishable products, but mainly deal with local production, inventory, distribution and retailing of these products [31].Supply chains for perishable products such as fruits are often more complex and more difficult to administer since food products have a short shelf-life [32].Unfortunately, the chains of perishable products consume a lot of energy [33] to prevent these products from decaying, and increase the emission rates, especially CO 2 [29].
One of the strategies used to reduce the emission rates is that used by Yang et al. [34], who stated that it is more effective to establish larger carbon quotas penalizing those companies that pollute the most.The problem is that a large investment has to be made to improve those polluting technologies.Another alternative to reduce emissions is the joint consideration of delays in payments for those polluting companies [35].These authors showed that the optimization of the production rate, the duration of the delay in payments or the size of the batch provides the optimal solution for carbon emissions and the cost of a supply chain system.Bouchery et al. [36] compared management costs and carbon emissions resulting from a non-coordinated two-echelon serial economic order quantity model to that of the centralized solution; their results show that coordination allows to reduce both costs and emissions in a supply chain.

CO 2 Reduction in the Distribution of Perishable Products in Colombia
Small farmers who supply food to the city of Bogotá face different challenges that endanger their means of subsistence and, by extension, the food security of the capital of Colombia [37].There are climatic, social, economic and environmental aspects, among others threatening the wellbeing of producers and population in general.The environmental aspects are currently in the initial stage of study.While many different strategies are needed at the national and local levels to address the increase in GHG emissions, solutions at the local level can play an important role in mitigating environmental quality problems in low-and middle-income countries such as Colombia [38].Fruit transport operations in rural areas of Colombia are generally characterized by the absence of technology to improve the conditions in which they move, reducing the competitiveness of the country's agricultural sector [5,39].
According to Reina and Adarme [40], in Colombia, there are few studies in the field of distribution logistics of perishable products, and there are no investigations related to the initial phases of the distribution process or the reduction of CO 2 emissions.The German Federal Ministry for the Environment, Nature Conservation, Building and Nuclear Safety [41] established that the total number of motor transport vehicles in Colombia is composed of about 250,000 cars of different capacities and ages.From that fleet, about a third is over 30 years old and generates the largest proportion of CO 2 emissions and air pollution.It is estimated that this type of vehicles is used by small transporters in rural areas to take their products to large cities for sale.
In relation to studies about GHG emissions, Cuellar et al. [42] compared different sources of CO 2 in the means of transport in Bogotá.His main result showed that public transport is not the main source of emissions.Consequently, he recommended undertaking studies on other sources of pollution such as cargo trucks.Regarding the territorial analysis of the emissions, the results of Rios [43] revealed that the most polluting places in the city of Medellín have the lowest population densities and the lowest variety of land uses.Thus, CO 2 emissions seem to be emitted by old vehicles.Román et al. [44] determined the drivers of change in CO 2 emissions between 1990 and 2012 in Colombia.This study analyzed the measures implemented by the Colombian authorities to mitigate polluting emissions, including the substitution of fossil fuels.The results allowed for the conclusion that efficiency and energy-saving policies are effective and should focus on three sectors: industrial, transport and residential.
For the transport sector, Valenzuela et al. [45] developed a model to estimate the GHG emissions generated by that sector in 2010.In this model, national emissions were projected for the 2010-2040 period, using Kaya factors.With these factors, various options were evaluated to mitigate emissions through a cost-effectiveness analysis.The results show that the country must enter a stage of increase in gasoline fuel savings to contribute to a reduction in CO 2 generation [46], indicated from their analysis of information that the Clean Development Mechanism projects implemented in Colombia since 2000 showed a reduction of 0.78%.This result is relatively low compared to countries such as China and India.Cardenas et al. [47] analyzed the effect of tax on carbon emissions and indicated that it has not reduced the use of fossil fuels in the country.However, simulations on this type of initiative indicate that evidently the cargo transport sector in Colombia will reduce the intensity of carbon emissions as a result of the adoption of technological strategies such as replacement of engines.There can also be a significant reduction through reduced use of means of transport, clearly with an improvement in the dimension of integrated spatial development (PTCI_Space) based mainly on its connectivity through a wide network of roads and transport systems [48].
To date, no studies have analyzed the relationship between the transportation of perishable products from rural areas near the city of Bogotá and the GHG emissions released by the vehicles that carry out these operations.According to Valenzuela et al. [45], efforts should be made for perishable products to travel shorter distances and to help reduce carbon emissions, which are the main driver of climate change.

Distribution Models of Perishable Products with Carbon Emission Reduction Constrains
Although harder emission controls on vehicles have reduced harmful gas emissions, attention has shifted towards the growth of CO 2 emissions from the cargo transport sector [8,49].In recent years, the development and application of operations research models in supply chains of perishable products have attracted many researchers in different countries [1].For example, Soysal et al. [50] performed an analysis of a fresh tomato supply chain including distribution costs dependent on truck load for a comprehensive assessment of CO 2 emissions and fuel consumption.Volpe et al. [51] calculated CO 2 emissions for all stages of the production chain from nut-based products to the final point of sale; Yang et al. [52] analyzed the decay process of perishable foods and determined the optimal temperature of the chain of cold with minimal emissions; and Gallo et al. [53] conducted a study, using mixed integer linear programming, for the design of sustainable supply chains for the distribution of fresh apples and ice cream.The number of researchers that link product perishability and transport emissions has risen [54].Validi et al. [55] proposed a bi-objective model for the distribution of milk in Ireland.Their model minimizes the total cost and CO 2 emissions in distribution channels.The presented sustainable distribution process allows to the supply chain decision maker to reduce the total carbon emission of the transport involved in the entire distribution process, while optimizing the total costs.Bortolini et al. [56] studied a distribution network of products grown in Italy and distributed in several European markets.Their study considered three objectives: total operating costs, total CO 2 emissions and delivery time of the products.Their results showed that the multi-objective perspective is of great help since it allows to reduce CO 2 emissions with a limited cost increase and an admissible delivery time.The carbon reduction is equal to 9.6% and the increase of the operating cost is of 2.7% compared to the traditional distribution network configuration.Soysal et al. [57] analyzed the benefits of horizontal collaboration related to product perishability, CO 2 emissions from transport operations and logistics costs in the Inventory Routing Problem (IRP) with multiple suppliers and customers.Through the development of a decision support model, their results indicated that the horizontal collaboration among suppliers allows reducing the aggregate total cost by 17% and the aggregate total CO 2 emissions by 29% in the case studied.

Methodology
The problem addressed in this article deals with the allocation of vehicles for the collection, transport and distribution of fruit from producers to sale points, passing through several intermediate points that are necessary for collection and better distribution.In particular, it is considered that there are two types of collection points, that is, the producers are divided into two groups.Some located on routes (a route is a group of producers with low fruit production) whose initial collection is done in a first location and from there they are transported to a storage center of greater capacity.The producers of the other group directly transport the harvested fruit to these collection centers, which are considered individually in the model because they produce much more fruit than those located on the routes.In addition, the collection on routes and from individual producers can be subcontracted.
Then, the fruit is transported from these collection centers to customers (wholesale centers) located in a city or urban center.For this case, customers are considered to be on a pre-established route, in such a way that the vehicles distribute the fruit by visiting customers one after the other.Figure 1 shows a diagram of the supply chain considered.Then, the fruit is transported from these collection centers to customers (wholesale centers) located in a city or urban center.For this case, customers are considered to be on a pre-established route, in such a way that the vehicles distribute the fruit by visiting customers one after the other.Figure 1 shows a diagram of the supply chain considered.To solve this problem, we designed a multi-objective mixed-integer programming model, which seeks to minimize both the total cost of management and the level of CO2 emissions of vehicles.The formulation of the problem is presented below.To solve this problem, we designed a multi-objective mixed-integer programming model, which seeks to minimize both the total cost of management and the level of CO 2 emissions of vehicles.The formulation of the problem is presented below.CK l = Variable transport cost ($/(kg km)) in vehicle l ∈ L CKm l = Transportation fixed cost ($/km) in vehicle l ∈ L CY i = Fixed cost ($) if the collection is subcontracted in route i ∈ I CYA mj = Fixed cost ($) if collection is subcontracted from the producer m ∈ M to the collection center j ∈ J CP l = Cost of toll payment ($) on the road that vehicle l ∈ L must pay CC kl = Parking cost ($) of vehicle l ∈ L at the customer's k ∈ K facility h j = Fixed cost ($/week) for using collection center j ∈ J EA i = Supply (kg/week) from route i ∈ I EB m = Supply (kg/week) from producer m ∈ M D k = Demand (kg/week) from customer k ∈ K CAP j = Capacity (kg/week) of collection center j ∈ J G l = Capacity (kg) of vehicle l ∈ L DD i = Length (km) of route i ∈ I DA mj = Distance (km) from producer m ∈ M to collection center j ∈ J DB j = Distance (km) from the first location to collection Center 1 DC jk = Distance (km) from collection center j ∈ J to Customer 1 DE kk = Distance (km) from Customer 1 to Customer 2 DT l = Maximum distance (km/day) vehicle l ∈ L can travel NV l = Number of available vehicles type l ∈ L EM l = Quantity of vehicle l ∈ L emissions (kgCO 2 /(kg km)) Mm = Very large number

Variables
x ilt = Quantity of fruit (kg/day) transported on route i ∈ I to the first location in vehicle l ∈ L on day t ∈ T xv ilt = Number of trips to be made on route i ∈ I to the first location in vehicle l ∈ L on day t ∈ T y it = Binary that indicates whether or not the collection is subcontracted in route i ∈ I on day t ∈ T xs it = Quantity of fruit (kg/day) whose harvest is subcontracted in route i ∈ I on day t ∈ T Q jlt = Amount of fruit (kg/day) transported from the first location to the collection Center 1 in vehicle l ∈ L on day t ∈ T QV jlt = Number of trips to be made from the first location to collection Center 1 in vehicle l ∈ L on day t ∈ TP mjlt = Quantity of fruit (kg/day) transported from producer m ∈ M to storage center j ∈ J in vehicle l ∈ L on day t ∈ T PV mjlt = Number of trips to be made from producer m ∈ M to storage center j ∈ J in vehicle l ∈ L on day t ∈ T U mjt = Binary that indicates whether or not the harvest is subcontracted from producer m ∈ M to collection center j ∈ J on day t ∈ T PS mjt = Amount of fruit (kg/day) whose harvest is subcontracted from producer m ∈ M to storage center j ∈ J on day t ∈ T W jklt = Amount of fruit (kg/day) transported from the storage center j ∈ J to Customer 1 in vehicle l ∈ L on day t ∈ T WV jklt = Number of trips to be made from collection center j ∈ J to Client 1 in vehicle l ∈ L on day t ∈ T V kk lt = Quantity of fruit (kg/day) transported from Customer 1 to Customer 2 in vehicle l ∈ L on day t ∈ T VV kk lt = Number of trips to be made from Customer 1 to Customer 2 in vehicle l ∈ L on day t ∈ T I NV jt = Amount of fruit (kg/day) left in collection center j ∈ J at the end of day t ∈ T I NVA j = Binary that indicates whether the storage center j ∈ J is used or not I NVB jt = Binary that indicates whether the storage center j ∈ J is used or not at the end of day t ∈ T DAUX kt = Amount of fruit (kg/day) that is distributed to Customer 1 on day t ∈ T ZTT = Total amount of CO 2 (kg/week) expelled CT = Total cost of management ($/week)

Model
By using the above notation, a multi-objective mixed-integer programming model to minimize the total cost and the total CO 2 emissions is formulated as follows: Objective 2 : ∀y it , U mjt , I NVA j , I NVB jt ∈ {0, 1} In this model, Equation (1) indicates the total cost of the management to be minimized, which includes transport costs among all links in the chain, costs of subcontracting the collection, costs for road tolls and parking at customers' facilities, and costs for collection center use.Note that the fixed costs associated with those variables that represent the number of trips are multiplied by two, because, in the objective, trips are accounted for both ways.Additionally, the variable that represents the amount of fruit collected on the routes is multiplied by 1.5, as the vehicle gradually collects fruit on the route, and returns with a full load.This indicates that, on the outgoing trip, on average, it is loaded at 0.5 of x ilt , while on the return trip, it is loaded at 100% of x ilt .Based on the methodology for calculating emissions by Li [58], Equation (2) minimizes the total amount of CO 2 emitted in a given week.Constraints (3) and ( 4) ensure that all fruit is collected, either on a route or from the individual producers, and by personal or subcontracted vehicles.Constraints ( 5)-( 9) ensure that the amount of fruit transported does not exceed the capacity of the vehicles.Constraints (10) and (11) ensure that subcontracted transport collects fruit only if the decision is made to subcontract.Constraints ( 12)-( 14) ensure that there is a balance in the flow of product in the first location and in the collection centers, respectively.Constraint (15) guarantees that the storage capacity of each storage center is not exceeded on any given day.
Constraint (16) ensures that a collection center is used on a given day only if it is to be used in the week to be planned.Constraint (17) indicates that it is not possible to store the fruit for more than two days, given its perishability.Constraint (18) ensures that a collection center allows product flow only if it is to be used.Constraint (19) ensures that every day there is a balance in the flow of product that passes through Customer 1. Constraint (20) ensures that each day, each type of vehicle transports fruit from Customer 1 to Customer 2 only in the case that it has been used to transport fruit from collection centers to Customer 1, in accordance with the pre-established route.Constraints ( 21) and (22) ensure that all fruit is delivered to customers.This is possible because all produced fruits are sold.Constraint (23) ensures that on any given day the different types of vehicles cannot exceed the maximum total distance traveled round trip.This is necessary because obviously is not possible to travel an unlimited distance in one day.Constraints ( 24)-( 26) indicate the variables that are positive, integer and binary, respectively.

Solution Approach
The multi-objective model was solved via the ε-constraint method [54].This method consists of optimizing a single objective, while the others are formulated as constraints.For the particular case of the present study, the minimization of Objective 1 was maintained, while Objective 2 was reformulated through Equations ( 27)-( 29), for which the additional parameter ε = maximum permitted quantity of emissions (kg CO 2 /week) and additional variable ZT l = total quantity of CO 2 (kg/week) emitted by vehicle type l (ZT l ≥ 0) were previously defined.
ZTT ≤ ε (29) Equation ( 27) calculates the quantity of CO 2 emissions emitted by each type of vehicle.Equation (28) calculates the quantity of emissions added for the type of vehicle.These two equations may be understood as intermediate steps.They are not necessary for the correct functioning of the ε-constraint method, but are important for identification of the quantity of CO 2 emitted by each vehicle, and later analysis of the form in which the scheduling of each of these affects emissions.Constraint (29) limits emission levels.
The methodology used to restrict emissions made use of the parameters shown below: emmax: maximum CO 2 emissions value.This corresponds to the value of the emissions when they are unrestricted, or when the model is run without Equations ( 27) and ( 28), and Constraint (29).emmin is the minimum CO 2 emissions value.This value is found by trial and error, reducing the value of ε as much as possible in Constraint ( 29), until the model becomes unfeasible.δ: integer parameter which varies within 1-10 to controllably limit emissions between emmin and emmax.In general, this parameter does not have to oscillate within 0-10, but the greater the number of values considered, the greater the number of test runs to be performed.The relationship between these parameters is as shown in Equation (30).
The variation of parameter δ within 0-10 in Equation ( 30) results in 11 values for ε, which were tested in each one of the five supply scenarios, for a total of 55 model test runs.Following said test runs, the relationship between both model objectives were established through the construction of a Pareto frontier.

Case Study
The case used for the test runs considers a real case from a fruit producer Association in Cundinamarca, Colombia, whose name cannot be revealed.In this case, there are four routes, eight individual producers, three collection centers, two final customers, two types of vehicles, and collection and distribution may be performed on any day from Monday to Friday.Table 1 shows the parameters associated with each type of vehicle.This Association was established in the area with the objective of increasing bargaining power with customers, as grouping different producers with similar interests results in a large quantity of product to be negotiated.Additionally, the aggregation of this fruit reduces transport costs, which would be higher if each producer transported their goods individually, as the distances within the area are shorter than those to customers (see Table 2, in which all values are noted in km).
Initially, the model was run for five groups (scenarios) of supply values, minimizing only Objective 1.In other words, initially, CO 2 emission minimization was not considered, although it was indeed measured.As supply variations depend mainly on the time of year, and not on each producer or individual route, it is supposed that the same variations exist for all producers in the region.For this study, we consider a scenario of minimal demand, a scenario of maximum demand, and a scenario of average demand, which are based on real data provided by the producer association.Additionally, two scenarios were considered, whose supply values are located between the minimum and the average, and between the average and the maximum.Concretely, if the average supply is taken as base, the values of the supply are as shown in Table 3.On the other hand, since the product to be distributed is fruit, everything that is produced is sold to the collection centers, and, as such, the total demand will always be equal to the total supply.

Results and Analysis
The general results of the test runs are shown in Table 4. Observe that, when the supply is minimal, the minimum management cost is COP 1,155,671.15/week(COP stands for Colombian Peso, local currency), while for the scenario of maximum supply, the cost almost triples to COP 3,098,655.25/week.Similar behavior is observed for CO 2 emissions, growing from 1138.06 kg/week with the minimum supply to 3117.13 kg/week with the maximum supply.Intermediate values are observed for scenarios two, three, and four, showing a directly proportional relationship between the variables.Additionally, relationships between CO 2 emissions and total supply, and between total cost and total supply are shown.It is observed that the differences between said indicators for the various scenarios are small.As one may observe in Objective 2 (Equation ( 2)), the total emissions calculation not only depends on the amount of fruit transported, but also on the type of vehicle used.As is explained further on, based on the results in Table 5, the usage frequency for each type of vehicle varies in each scenario, and, from there, the slight variations in the relationship calculated in Table 4 are derived.Something similar occurs with Objective 1 (Equation ( 1)), in which the total cost depends on the amount transported, type of vehicle used, and number of trips made.
The concrete scheduling required to achieve these costs indicates that it is not necessary to subcontract collection for any route or producer, and as such, adequate assignment of the association's own vehicles is sufficient for collection and distribution.Said scheduling is shown in Table 5, in which R1, R2, R3, and R4 indicate whether each vehicle is assigned to collection routes, CA indicates whether transport is assigned between the first location and Collection Center 1, M1, M2, . . ., M8 indicate whether collection is assigned to individual producers, and in the row entitled "Customers", J1 and J2 indicate whether the vehicle is assigned to transport from Collection Center 1 or 2 to Customer 1. Finally, C2 indicates whether transport is assigned between Customer 1 and Customer 2. In this case, the model test runs never indicated that Collection Center 3 needed to be used, and, for this reason, it was not added to the table.These results indicate that the minimum costs attained making intensive use of 6000 kg vehicles, which is achieved by assigning them to collection and distribution on several days of the week.Thus, for Scenarios 1 and 2, two days of work would be sufficient; for Scenarios 3 and 4, three days of work; and for Scenario 5, four days of work, from Tuesday to Friday.The use of collection centers for fruit storage from one day to the next facilitates this work.Table 7 shows the quantities to be stored (kg) each day in each scenario, in each collection center.Observe that, in Scenario 1, storage from one day to the next is not necessary (although Collection Center 1 is used to collect and distribute, as shown in Table 5), Collection Center 1 is used for storage in Scenarios 2 and 5, and Collection Center 2 is used in Scenarios 3-5.

CO 2 Emission Results
Table 8 shows emissions by type of vehicle, in the case of both maximum and minimum emissions.Additionally, the costs for the case of minimum emissions are shown.It is evident that, to reduce emissions, more intensive use of vehicles with lesser capacities (1500 kg), or even non-use of Vehicle 1 is required in Scenarios 1-4.The more intensive use of Vehicle 2 implies that the total costs increase to more than double the costs shown in Table 4 for each scenario.Likewise, CO 2 emissions can be reduced, at best, to between 39% and 51% of the maximum level, depending on the scenario.The greatest reduction is achieved for minimum supply Scenarios 1 and 2 while the 51% reduction applies to the scenario with highest supply.to compensate for the supply which was covered before with a 6000 kg vehicle.For all scenarios, collection routes and individual producer collection is subcontracted, and the rest of the links are covered by Vehicle 2, with the exception of transport between the first location and Collection Center 1 (J1) on Tuesday, and transport from Collection Center 2 (J2) to both customers on Thursday, which is performed with Vehicle 1 in Scenario 5.
In contrast to the case in which CO 2 emissions are not restricted, this time Collection Center 3 is used to achieve better coverage in the collection-distribution process.Its use limits only those scenarios of greatest supply (4 and 5) and it only receives fruit from two individual producers: M1 and M2 on Monday in Scenario 4, and M2 and M4 on Friday in Scenario 5. Additionally, there is more intensive use of Collection Center 2 (J2), as in Scenarios 1 and 2, Collection Center 1 (J1) is only used to receive fruit coming off of routes, while in the other scenarios, apart from fruit coming off of routes, it is used exclusively for reception of goods from individual producers M4 and M6.For this case, in Scenarios 1-4, 100% of the fruit is transported, whether with Vehicle 2 or with subcontracted vehicles, depending on the link (Table 10).In said table, the minimal use of Vehicle 1 is demonstrated, as it is only used in Scenario 5 to transport a maximum of 16.6% of the fruit in the link which goes from Collection centers to Customer 1.As extreme cases, the scheduling of vehicles has been shown in detail for the cases of maximum and minimum emissions.However, model test runs were performed for all integer values of δ within 0-10, resulting in a Pareto frontier between model objectives for each one of the supply scenarios considered, or between CO 2 emissions and the cost of chain management incurred to achieve said emission levels (Figure 2).
Independently of the scenario considered, a higher cost implies a lower quantity of CO 2 emitted.Further, each scenario's Pareto frontier behaves similarly.For mean values of δ, the chart is practically linear, but on the ends, a change in the slope is shown, which results in each curve presenting upward concavity.This behavior on the ends for large δ values (left end of the curves), is a result of low or null subcontracting, and as such, a simple rescheduling of the association's vehicles allows them to achieve important changes in the quantity of CO 2 emitted, with little additional cost.However, for small δ values (right end of the curves), to guarantee a low level of emissions, all routes and collections from individual producers are subcontracted, and, as such, there are few options for vehicle rescheduling.This implies a considerable increase in cost, when compared to the minimal reduction in emissions.
As extreme cases, the scheduling of vehicles has been shown in detail for the cases of maximum and minimum emissions.However, model test runs were performed for all integer values of within 0-10, resulting in a Pareto frontier between model objectives for each one of the supply scenarios considered, or between CO2 emissions and the cost of chain management incurred to achieve said emission levels (Figure 2).Independently of the scenario considered, a higher cost implies a lower quantity of CO2 emitted.Further, each scenario's Pareto frontier behaves similarly.For mean values of , the chart is practically linear, but on the ends, a change in the slope is shown, which results in each curve presenting upward concavity.This behavior on the ends for large values (left end of the curves), is a result of low or null subcontracting, and as such, a simple rescheduling of the association's vehicles allows them to achieve important changes in the quantity of CO2 emitted, with little additional cost.However, for small values (right end of the curves), to guarantee a low level of emissions, all routes and collections from individual producers are subcontracted, and, as such, there are few options for vehicle rescheduling.This implies a considerable increase in cost, when compared to the minimal reduction in emissions.
Clearly, the principal variable that affects total cost as well as the quantity of emissions is supply, which displaces the Pareto frontier.Figure 2 shows the way in which, as supply increases from Scenario 1 to Scenario 5, the charts move positively, and, for each curve, a greater difference between the ends thereof is presented, i.e., the curve elongates.On the one hand, the motive for the movement of the frontier is evident, given that costs and emissions are directly proportionate to supply.On the other hand, the curve elongates with the change of scenario, given that greater supply increases vehicle scheduling leeway, thus maintaining the objective of minimum cost achievement.For 1,000,000 2,000,000 3,000,000 4,000,000 5,000,000 6,000,000 7,000,000 8,000,000  Clearly, the principal variable that affects total cost as well as the quantity of emissions is supply, which displaces the Pareto frontier.Figure 2 shows the way in which, as supply increases from Scenario 1 to Scenario 5, the charts move positively, and, for each curve, a greater difference between the ends thereof is presented, i.e., the curve elongates.On the one hand, the motive for the movement of the frontier is evident, given that costs and emissions are directly proportionate to supply.On the other hand, the curve elongates with the change of scenario, given that greater supply increases vehicle scheduling leeway, thus maintaining the objective of minimum cost achievement.For example, greater supply implies that the collection-distribution process may be redistributed to other days of the week without incurring excessive cost increases when transporting near capacity each day.For scenarios of low supply, this would not be possible, as making use of additional days would mean sending vehicles with little cargo, which would go against the objective of minimizing total cost.Said increase in leeway, then, absolutely causes emissions to be further reduced, and, as such, costs reduce also.

Quantity of CO
However, in a relative sense, the behavior is not the same, as observed in Figure 3 Pareto frontiers.This shows that the differences between scenarios are not as pronounced when the cost is measured as an increase of the minimum of each scenario, and the quantity of emissions is measured as a percentage of the maximum of each scenario.For example, Scenarios 1 and 2 present the greatest relative CO 2 emission reductions, but Scenario 2 does so with a lesser cost increase (39% of emissions in both scenarios, with a cost increase of 144% in Scenario 1, and 130% in Scenario 2).
This behavior shows that there is no relationship between the supply of fruit and the increase in total cost, which is corroborated in Figure 3, as it does not present a natural order, as observed in Figure 2.Only for δ = 0 (right end of the curves), with the exception of Scenario 2, it is observed that greater supply (Scenario 5) implies a lower increase in cost (119%) and vice versa (see Scenario 1).Said behavior is inverse to that shown in Figure 2.This means that, if supply is kept constant, a high cost could be reduced a great deal, in terms of COP, but not so much in terms of percentage.Likewise, a low cost can be reduced minimally in terms of COP, but a bit more in terms of percentage.
Similarly, there is no relationship between the supply of fruit and the relative quantity of emissions, with the exception of the case in which δ = 0, in which an inverse supply-quantity of emissions relationship is presented as a percentage of the maximum.For example, for Scenario 5 (that with greatest supply) emissions may be up to 51% of the maximum, while in Scenario 4, they may be 45%, and in Scenario 3, 43%.This means that, in general, important reductions in the levels of CO 2 emissions can be achieved just with vehicle rescheduling, which implies an increase in costs associated purely with the collection-distribution process, but investment in additional technology would be unnecessary.
example, greater supply implies that the collection-distribution process may be redistributed to other days of the week without incurring excessive cost increases when transporting near capacity each day.For scenarios of low supply, this would not be possible, as making use of additional days would mean sending vehicles with little cargo, which would go against the objective of minimizing total cost.Said increase in leeway, then, absolutely causes emissions to be further reduced, and, as such, costs reduce also.
However, in a relative sense, the behavior is not the same, as observed in Figure 3 Pareto frontiers.This shows that the differences between scenarios are not as pronounced when the cost is measured as an increase of the minimum of each scenario, and the quantity of emissions is measured as a percentage of the maximum of each scenario.For example, Scenarios 1 and 2 present the greatest relative CO2 emission reductions, but Scenario 2 does so with a lesser cost increase (39% of emissions in both scenarios, with a cost increase of 144% in Scenario 1, and 130% in Scenario 2).This behavior shows that there is no relationship between the supply of fruit and the increase in total cost, which is corroborated in Figure 3, as it does not present a natural order, as observed in Figure 2.Only for = 0 (right end of the curves), with the exception of Scenario 2, it is observed that greater supply (Scenario 5) implies a lower increase in cost (119%) and vice versa (see Scenario 1).Said behavior is inverse to that shown in Figure 2.This means that, if supply is kept constant, a high cost could be reduced a great deal, in terms of COP, but not so much in terms of percentage.Likewise, a low cost can be reduced minimally in terms of COP, but a bit more in terms of percentage.
Similarly, there is no relationship between the supply of fruit and the relative quantity of emissions, with the exception of the case in which = 0, in which an inverse supply-quantity of emissions relationship is presented as a percentage of the maximum.For example, for Scenario 5 (that with greatest supply) emissions may be up to 51% of the maximum, while in Scenario 4, they may be 45%, and in Scenario 3, 43%.This means that, in general, important reductions in the levels of CO2 emissions can be achieved just with vehicle rescheduling, which implies an increase in costs associated purely with the collection-distribution process, but investment in additional technology would be unnecessary.

Conclusions
The present investigation addressed a tactical problem, which was modeled in a supply chain for the collection and distribution of fruit from producers in the countryside to retail customers.In an academic environment, the main contribution of this study is in the analysis of a supply chain through the design of a multi-objective mixed-integer programming model.This model includes decisions about quantities of perishable product to be transported and stored, the number of trips to be made to different links in the chain for vehicles with diverse capacities, use of existing collection centers, subcontracting decisions, and vehicle CO 2 emission levels.All of these decisions must be made from Monday to Friday, given weekly supplies.In the practical scope, the main contribution of this study is in the creation of a plan for weekly vehicle scheduling for the fruit producers Association, which minimizes the total cost of chain management and CO 2 emission levels.
Five supply scenarios were considered, based on real data in which, initially, the objective of minimizing CO 2 emissions was not considered.The results of these test runs revealed a collection-distribution plan for each scenario, each day of the week, in which the use of vehicles with greater capacity (6000 kg) was much greater than use of 1500 kg vehicles, independently of the scenario or link in the chain.Additionally, it is not necessary to subcontract any point or collection route, as the Association's vehicles provide a sufficient weekly capacity, without the need to incur additional costs.This plan provides the absolute minimum cost for each scenario, as non-restriction of emissions may achieve the minimal cost goal.Thus, if the fruit producers Association decides to implement this plan, they should schedule vehicles as described in this article.It is suggested that the Association previously design a forecast model, which would allow them to determine, with the greatest possible precision, the supply of fruit that may be available, and thus schedule collection-distribution in accordance with that information.
Later, CO 2 emissions were minimized using the ε-constraint method [54], in which emissions are limited through the addition of constraints to the model.Said procedure increased weekly management costs, principally owing to a greater use of 1500 kg vehicles than 6000 kg vehicles.Although the former emits less CO 2 per unit of cargo transported and unit of distance covered, their use is costlier, not only per unit of cargo, but also because a greater number of trips must be made to maintain collection and distribution of 100% of the fruit.Secondly, the cost of management also increases, owing to the subcontracting of collection on routes and for individual producers.The representation of the cost of subcontracting in the total cost depends on the scenario considered, and on the limit established in the emissions.However, in the scenario in which this is greater, it presents a maximum of just 17%.Despite this, it is suggested that future investigations include the quantity of subcontracted vehicle emissions, although their cost impact can be minimal.
Finally, although a Pareto frontier was established which shows a general relationship between the quantity of CO 2 emissions being inversely proportional to management cost, it was observed that supply affects this relationship.For future studies, then, the effect of other variables, such as loading and unloading times, use of less contaminating vehicles (which could include decisions about technological investment within the model), traffic in the urban distribution, and the state of roads in the collection chains should be incorporated.In these rural areas, many roads are not paved, and are affected by rain, heavy vehicle traffic, or during intense drought, by the dust generated, which impedes visibility.A tactical-operative or purely operative study could be carried out to consider all of these variables.Additionally, the model could consider that the real supply of fruit is random, and, as such, a study could be carried out which considers the robustness of the schedule designed.

Figure 1 .
Figure 1.Configuration of the supply chain.

3. 1
. Sets I = Set of routes to the first location M = Set of individual fruit producers not located on a route J = Set of collection centers K = Set of final customers L = Set of types of vehicles T = Set of days of the week 3.2.Parameters

Figure 2 .
Figure 2. Relationship between total cost and quantity of CO2 emissions.

Figure 2 .
Figure 2. Relationship between total cost and quantity of CO 2 emissions.

Figure 3 .
Figure 3. Relationship between total cost increase and the relative quantity of CO2 emissions.

Figure 3 .
Figure 3. Relationship between total cost increase and the relative quantity of CO 2 emissions.

Table 2 .
Distances (km) for different links in the chain.

Table 4 .
Results of the first test runs.

Table 6 .
Quantity of fruit (%) transported by every type of vehicle.

Table 7 .
Quantities (kg/day) of fruit inventory to be stored.

Table 8 .
Maximum and minimum CO 2 emission levels.

Table 10 .
Quantity of fruit (%) transported by every type of vehicle for minimum emissions.