key: cord-0900486-wfy6685k
authors: Keshavarz-Ghorbani, Fatemeh; Pasandideh, Seyed Hamid Reza
title: Modeling and optimizing an agro-supply chain considering different quality grades and storage systems for fresh products: a Benders decomposition solution approach
date: 2021-10-04
journal: J Comb Optim
DOI: 10.1007/s10878-021-00802-5
sha: 1aba0d471abf399e842c3af51b8595c90db2b013
doc_id: 900486
cord_uid: wfy6685k

This paper proposes a mathematical model in the context of agro-supply chain management, considering specific characteristics of agro-products to assist purchase, storage, and transportation decisions. In addition, a new method for determining the required quality score of different types of products is proposed based on their loss factors and purchasing costs. The model aims to minimize total cost imposed by purchasing fresh products, opening warehouses, holding inventories, operational activities, and transportation. Two sets of examples, including small and medium-sized problems, are implemented by general algebraic modeling language (GAMS) software to evaluate the model. Then, Benders decomposition (BD) algorithm is applied to tackle the complexity of solving large-sized instances. The results of both GAMS and BD are compared in terms of objective function values and computational time to demonstrate the efficiency of the BD algorithm. Finally, the model is applied in a real case study involving an apple supply chain to obtain managerial insights.

packaging, storing, and distributing agro-products, can add significant value to logistics while satisfying customers' demands and maintaining food security. Therefore, research on agricultural logistics can speed up modernization and competition abilities in agricultural processes (Otchia 2014) .

The specific characteristic of agro-products, such as seasonality, perishability, and long procurement time, make it different from other supply chains, therefore, it is necessary to have accurate information about their effects on supply chain management. These characteristics put restrictions on logistic decisions and make deep concerns over resource depletion and the occurrence of agricultural wastes and losses. According to Widodo et al. (2006) , 20-60% of agro-products are wasted along supply chains in any country. Therefore, designing efficient inventory and transportation systems for perishable products plays a vital role in cost-saving and freshness-keeping (Coelho and Laporte 2014) . From this perspective, cold chains with proper temperature adjustment can significantly reduce losses (Yu and Xiao 2017) .

Although efficient agro-supply chain management prevents storage and transportation losses, facility location problem is another issue strongly affecting the security and safety of foods, especially in developing countries (Morton et al. 2008) . Given capacity constraints in warehouses, determining warehouse locations and assigning cold storage rooms under controllable temperature conditions play a pivotal role in loss and waste reduction over a supply chain that may yield a complicated problem (Aung and Chang 2014) . Forasmuch as the current study deals with agricultural products with a short shelf life, optimal facility location and storage condition decisions can reduce losses due to the dispersed distribution centers, long distances transportation, etc.

This study presents an agro-supply chain based on a case study involving a wide range of decisions such as supplier selection, facility location, inventory control, and transportation planning. This research aims to optimize the shelf life of products via assigning proper cold storage systems that can prevent postharvest losses. In addition, different quality grades are considered for each type of product, whereby those farmers who have gotten the average quality score are selected to satisfy warehouses' demands. For these purposes, a mixed-integer nonlinear programming (MINLP) problem is formulated considering several operational constraints. Then, the model is linearized and transformed to mixed-integer linear programming (MILP) problem and solved by GAMS/CPLEX solver in small and medium sizes. Moreover, Benders decomposition (BD) algorithm is proposed to tackle the complexity of the problem caused by increasing the number of binary and integer variables in large sized-instances. Then, the proposed model is implemented for a real case study in Iran covering eight months of a year.

The rest of the paper is organized as follows: Sect. 2 gives an overview in the context of agro-supply chains. In Sect. 3, problem definition, notations, assumptions, and the mathematical model are presented. In Sect. 4, BD algorithm is proposed to tackle the complexity of large-sized instances. In Sect. 5, numerical experiments and a case study are performed to demonstrate the practicality of the proposed model. Section 5 draws conclusions and remarks on some future research directions.

Economic costs in agro-supply chains have been implied by inefficient management in purchasing, storage, and distribution of agro-products (Ahumada and Villalobos 2009) . Although there are copious studies dealing with facility location and transportation problems, they have not been widely discussed for agricultural products. The facility location problem is a critical element of strategic decisions in agro-supply chains, especially in developing countries (Morton et al. 2008) . This is because transportation and storage time will become shorter by properly determining facility locations (Morganti and Gonzalez-Feliu 2015) . Therefore, dealing with the capacitated facility location problems with applications in storing and distributing agro-products is an essential component to adapt to environmental and competitive aspects (Validi et al. 2014) .

Here, related literature with emphasis on three core parts of agro-supply chains, purchasing, facility location, and transportation decisions, has been reviewed to pinpoint the research gap. In this context, Paam et al. (2016) and Ganesh Kumar et al. (2017) have comprehensively reviewed agro-supply chain models. Then, Banasik et al. (2017) analyzed tradeoffs between economic and environmental benefits of an agro-food closed-loop supply chain with location-allocation decisions and production planning. They attempted to valorize agro-wastes by applying recycling concepts. Ghezavati et al. (2017) proposed a mathematical model for a distribution planning system considering the freshness and ripeness of products. They proposed BD algorithm to solve their problem in a reasonable time. Orjuela-Castro et al. (2017) presented a MILP model for a multi-product facility location problem. They aim to minimize the costs of transportation, installation, and fruit losses and show how intermediate logistics platforms affect fruit losses. Soto-Silva et al. (2017) formulated a two-stage agro-supply chain model for a case study in Chile. The first stage contains particular decisions on supplier selection and transportation planning and the second stage associates with establishing different types of cold chambers to store products and transportation planning. Allaoui et al. (2018) proposed a two-stage programming model for a sustainable agro-food supply chain. In the first stage, a hybrid analytic hierarchy process with an ordered weighted averaging method was used for assigning weights to sustainable indicators. In the second stage, a multi-objective mathematical model was proposed with respect to the results in the first stage. Paam et al. (2018) proposed an agro-supply chain model to optimize the shelf lives of apples with different quality grades in cold storage rooms. Liu et al. (2018) studied optimal purchase decisions and inventory retrieval policies for a single-product finite period agro-supply chain to maximize total profit under fluctuating prices. Jonkman et al. (2019) investigated the impacts of location-allocation, inventory, and transportation decisions on an agro-food supply chain under uncertain conditions. Liu et al. (2019) developed a transportation model for fresh agricultural products considering the impacts of ethephon concentration, time, and temperature on the maturity of soluble solids. They used ant colony algorithm to solve their model. Onggo et al. (2019) assessed stochastic demand in a multi-period inventory routing problem and attempted to minimize inventory, transportation, and shortage costs. They used a local metaheuristic algorithm to tackle the complexity of their model. Paam et al. (2019) proposed a new model for inventory management optimization in an apple supply chain. They analyzed the effects of two scenarios in storage rooms on fruit losses. Yu et al. (2019) investigated the impacts of horizontal and vertical integration on the quantity/quality losses and retailer profits.

The prior work on agro-supply chains is summarized in Table 1 to reveal existing gaps in the literature. The majority of the research has investigated multi-level supply chains under several constraints, without sufficient consideration on specific characteristics of agro-products, farmers, and distribution processes. However, a survey of practical experiences illustrates significant impacts of these properties on spoilage rate, in turn, on the amount of loss and system cost. Therefore, this paper takes several actions to fill the gaps, which are vividly explained as follows:

This study is an extension of the work done by Soto-Silva et al. (2017) in a way to be applicable for an apple supply chain in Iran. In their study, unlike many agro-supply chain models, three types of refrigeration technologies for maintaining different types of agro-products have been proposed that can make a privileged contribution in this context. In spite of their primary step towards sustainability of agro-products, it lacks an in-depth discussion on adjustment of storage conditions and its effects on inventory decisions. However, practical accessibility gave us valuable insights to design a mathematical model for adjusting storage conditions according to costs imposed on the system and its effects on the durability of agro-products. In addition, the storage of agro-products for a long-time mainly depends on their initial quality. Therefore, this study aims to manage the quality of agro-products by giving scores to producers and accordingly selecting them. To the best of the authors' knowledge, there is no previous effort on quality management despite its importance in handling the shelf life of agro-products. Overall, the main contributions of the paper are as follows:

• A MINLP model has been presented to jointly optimize facility location, transportation, inventory decisions in an agro-supply chain. • A new method is proposed to assign weights to suppliers based on their loss factor. Accordingly, decisions on supplier selection are made. • The current research seeks to optimize the shelf life of products through adjusting appropriate storage conditions. • The proposed MINLP method is linearized and transformed to a MILP model. • BD algorithm has been utilized to tackle the complexity of the proposed model and solve it in a reasonable time. • The proposed model is implemented in a real case study in Iran to demonstrate its applicability and efficiency.

This study is an extension of the proposed model by Soto-Silva et al. (2017) , who addressed supplier selection, transportation, and facility location decisions in an agro supply chain. However, their study lacks in-depth discussions on quality management and storage and inventory planning. Therefore, we formulate a multi-period MINLP model under some operational constraints to intertwine supplier selection, transportation, storage operation, inventory planning, and facility location decisions in an agro-supply chain. Afterwards, we reformulate it to a MILP model and propose BD algorithm to solve it. Figure 1 is meant to show the process of this paper. Figure 2 is a perspective of the proposed agro-supply chain, including producers, warehouses, and a wholesaler. Fresh products should be purchased from the producers who gain the required score based on their quality grades. This approach satisfies the required quality of products to be kept for a long time in cold chambers, which is more profitable for the system. The products are shipped to the warehouses by the private fleet. The reason is that producers, either small or medium sizes, are often located in rural and sparsely populated areas. The long distances between orchards and facilities and busy harvest seasons require the producers to provide flexible operating by their private fleet. The processing operations should carry before the storage of products, by which products are sorted based on their varieties, sizes, etc. Then, each type of products regardless of their sizes, color, etc., are kept in a cold storage room, and the temperature conditions should be adjusted to optimize the shelf life of products. Afterward, the products are delivered to the wholesaler to satisfy customers' demands. Herein, transportation decisions are made by warehouses. In total, the proposed model aims to minimize the total cost of the agro-supply chain model while satisfying the constraints.

In order to formulate the mathematical model, the following assumptions are considered. In addition, indices, parameters, and decision variables are presented in Table  2 .

• Shortages are not allowed.

• Three different quality grades are considered for products in orchards.

• Each of the quality grades for each variety of products imposes different amounts of losses. • Each variety of products is available for purchase at a specific period. Therefore, they are purchased with sufficient quality to be stored in cold chambers in order to meet customers' demands for a while. • Decisions on the type of temperature conditions should be made for each cold chamber in warehouses.

In this section, a MINLP model is presented as follows:

Min Z i∈I j∈J s∈S t∈T 

Parameters CT jw : Operational costs to transport products from producer j to warehouse w 

The objective function (1) minimizes total cost consisting of six components, purchasing cost, operational costs including packing, laborers, and transportation costs (from producers to warehouses), operational costs in warehouses, fixed cost of opening warehouses, fixed cost of transportation (from warehouses to the wholesaler), inventory holding costs. Constraints (2) ensure that the purchased products are of satisfactory quality with reasonable prices. Constraints (3) ensure that only one type of product should be stored in each cold chamber under a specific temperature condition. Constraints (4) indicate the effects of temperature conditions on the shelf life of products. Constraints (5) ensure that products with the same quality grade are purchased from a producer. Constraints (6) show the relation between variables A i jst and Q F i jst . Constraints (7) ensure that the amount of product in warehouses should not exceed the amount of purchased product. Equations (8) and (9) show the inventory level of products in warehouses where g < t and g t, respectively. Constraints (10) state the limited storage capacity of cold chambers in warehouses. Constraints (11) are the demand response. Constraints (12) represent the limitation in the number of cold chambers in a warehouse. Constraints (13) are capacity constraints of vehicles for transporting products from warehouses to the wholesaler. Equations (14) calculate the losses after purchasing products from producers. Constraints (15) indicate the available capacity of producers to supply product i. Constraints (16) define types of decision variables which are non-negative continuous, binary, or integer variables.

In the proposed model, the upper limit of index g is considered as a decision variable making the model nonlinear. In these cases, if-then decisions should be made to remove spoiled products from warehouses. Indeed, if the storage time exceeds the shelf life of products, the products should be removed from warehouses (If t − g > T E iwvct → I iwvctg 0; QW iwvcktg 0). Conversely, if the storage time is less than the shelf life of products, Eq. (8) are used to calculate the inventory level of products

. In order to handle the nonlinear formulation and transform the MINLP model to MILP, a linearization technique is derived from Chen et al. (2010) .

Consider the following nonlinear equivalents ( f (x) ≤ b → g j (x) ≤ b j ) where b and b' are the right-hand side of inequalities and f (x) and g j (x) are the functions of variable x. In this case, if f (x) is less than b, then g j (x) ≤ b j . The linearization of the if-then constraint is applied as follows:

where BIG and y are a positive big number and binary variable, respectively. Accordingly, the following equations are presented to linearize the if-then constraints in this paper.

The complexity of MILP problems rises by increasing the number of binary and integer variables, consequently, has a strong influence upon the computational time (Keshavarz-Ghorbani and Pasandideh 2021). Therefore, it is necessary to apply efficient solution procedures to overcome the computational burden. BD algorithm is known as one of the powerful algorithms to tackle the complexity of NP-hard problems. BD algorithm uses the concept of "divide and conquers" to reformulate the original problem into a master problem (MP) and dual subproblem ( The primary steps of the BD algorithm work iteratively by outer linearization and relaxation. In other words, the original problem is decomposed into an MP and subproblem (SP). The MP is a pure integer problem and includes integer and binary variables. The MP is solved and puts its solutions among the SP. In this step, the SP turns into DSP to make the feasible solution space independent from integer solutions. Then the DSP is solved, and its optimal bounded solutions/ unbounded solutions are added as feasibility/optimality cutting planes into the MP. This procedure reduces the computational time compared to the original problem. The convergence of BD comes from closing upper and lower bounds obtained from solving the DSP and MP, respectively. Figure 3 illustrates the steps of the iterative BD algorithm.

To implement BD for the proposed model, the MP is called to minimize a part of the objective function associated with binary and integer variables with respect to constraints (3)-(5), and (12). That is a relaxation of the original problem to obtain lower bounds. 

w∈W C W w AW w + k∈K t∈T C FT k y tk i c AC iwvct ≤ 1 ∀w ∈ W ; v ∈ V ; t ∈ TYesT E iwvct ≤ T ic AC iwvct ∀i ∈ I ; w ∈ W ; v ∈ V ; c ∈ C; t ∈ T ; s A i jst ≤ 1 ∀i ∈ I ; j ∈ J ; t ∈ T ; i v c AC iwvct ≤ F W w AW w ∀w ∈ W ; t ∈ T AW w , AC iwvct , A i jst ∈ {0, 1}; y kt , T E iwvct ∈ I nteger +(23)

The SP includes a part of the objective function associated with continuous variables subject to constraints (2), (6)-(11), (13)-(15) that can be seen in Table 3 . 

The DSP makes the feasible solution space independent from integer solutions and reaches the same objective function value as the SP. Therefore, the DSP is used instead of the SP to obtain upper bounds for the original problem. In this paper, the DSP is formulated as follows (Table 4) :

The DSP and its corresponding variables are presented in Table 5 . After solving the DSP, an optimal or feasible cut is added to the MP according to Eq. (24) and an extreme ray of DSP. BD algorithm terminates whenever the stop

where S is the feasible solution space of the MP, ε is a small non-zero number. In addition, UB iter and LB iter are the upper and lower bounds obtained from solving DSP and MP, respectively. The pseudocode of BD is given as follows: Table 4 Dual variables associated with the sub-problem 

This section is presented in three parts. In the first part, numerical test problems, including small, medium, and large-sized instances, are executed for assessing the efficiency and applicability of the proposed model and solution approach. These test problems are randomly generated in Excel and presented in Tables 6 and 7. In the second part, the model is implemented in a case study in Iran to derive valuable insights. Finally, a sensitivity analysis is carried out to help decision-makers determining key parameters which should be carefully tracked in the presented model.

The presented MILP model and BD algorithm are coded in GAMS 27.2 optimization software and solved by a CPLEX solver. Both CPLEX solver and BD algorithm are executed on a computer with Intel @ Core™ i7-CPU 2.20 GHz and RAM 8.00 GB. Different sizes are presented in Table 8 , and numerical results are presented in Table  9 , including optimal values of objective functions, the lower and upper bounds of BD, computational time, and gaps between upper and lower bounds of BD and CPLEX solver. BD algorithm continues once it reaches termination criteria which are as follows:

• The number of iterations reaches 200 • A relative gap between lower and upper bounds become less than 0.3%

In addition, the CPLEX solver terminates once the absolute and relative gaps become less than 0.001.

The results in Table 9 point out that the CPLEX solver reaches optimal solutions in less than one minute, whereas the computational time of BD is more than the CPLEX solver in several small sized-instances. As increasing the number of variables in medium sizes, the lower and upper bounds of BD converge to optimal solutions less than an hour, which is remarkably faster than the CPLEX solver. It should be noted, for 3  3  3  3  3  5  3   Max  25  20  20  3  3  10  3 large sized-instances, the CPLEX solver cannot find optimal solutions and releases an "OUT OF MEMORY" error. In these cases, BD can be utilized to cope with complexity and reach optimal solutions in a reasonable time. Figure 4 demonstrates the changes of computational time for BD and CPLEX solver, in which BD outperforms CPLEX solver, especially in medium and large-sized instances. According to the results, the computational time of the CPLEX solver increases whatever the size of the problem increases. This affirms the efficiency of BD in solving the model in a reasonable time. Besides, the average gap between the upper bounds of BD and the CPLEX solver is about 0.01684. It turns out BD converges to optimal solutions with a near-zero gap. In order to verify the efficiency of BD and its convergence, three instances in small, medium, and large sizes are derived from GAMS/CPLEX solver and presented as follows: Figure 5 shows that BD provides optimality convergence in less than 17 iterations for a small-sized test (no. p5). Figure 6 shows that BD provides optimality convergence in less than 12 iterations for a medium-sized test (no. p8). Figure 7 demonstrates that BD provides optimality convergence in less than 23 iterations for a large-sized test problem (no. p14).

This section devotes a discussion on an apple supply chain belonging to a domestic distributor in Iran to demonstrate the application of the proposed model. The company decides to establish warehouses at four potential locations, namely Tehran, Fars, Isfahan, and Razavi Khorasan, and activate cold chambers under proper temperature conditions. They often purchase different varieties of apples from producers who truck products by private vehicle fleet and located in Tehran, Alborz, Isfahan, West Azerbaijan, East Azerbaijan, and Semnan provinces. The company purchases good storage apples based on quality characters (i.e., texture, cosmetic appearance) to remain for a while in warehouses and satisfy the wholesaler's demands in different months of a year. Tables 10, 11, 12, 13, 14, 15 and 16 represent input data related to cold storage construction projects (www.damoon-co.com) and purchasing, transportation, and storage planning for eight months. The average demand is obtained based on historical data when the company faces shortages, and also demands should be met by any variety of apples. In addition, the purchasing cost of apples with different quality grades are available in fruit and vegetable fields organization in Iran, and information about temperature conditions in cold chambers are available in (https://pdfsara.ir/containerstorage-services-atmospheric-plan/).

The quality score of fruits supplied by different producers (W E P i js ) is calculated based on the minimum and maximum values of loss factor for a product with quality grade s as follows: where M I L is is the minimum loss factor for product i with quality grade s and M AL is is the maximum loss factor for product i with quality grade s. Equation (26) illustrates the quality score of producer j2 supplying product i1 with quality grade s3, which equals 1.

W E P i js 1 − 0.06 − 0.06 0.09 − 0.06 1 (26) 

The model has 184,572 discrete variables, 330,049 continuous variables, and 642,748 equations for the case study and is implemented by both CPLEX solver and BD. Herein, the efficiency of the proposed model and related solutions have been discussed. The lower and upper bounds of BD become 25,289,297 and 25,289,302, respectively. In addition, the optimal solution is obtained by the CPLEX solver, which is 25,289,302.

The results indicate that BD converges to the optimal solution in a reasonable time. Figure 8 shows the convergence of BD for the case study. First, producers who gain good scores are selected to supply products. The results show a 44% of apples are of first quality grade, 45% are of second grade, and 12% are of third quality grade. In this respect, the total cost of purchasing apples with different quality grades is about US$ 996,000.

According to the results, it is beneficial to open warehouses in Tehran, Fars, and Isfahan provinces at a total cost of US$ 6,600,000. In this case, the operational costs in warehouses become US$ 147. It is worth mentioning that the presented model is able to show the temperature condition in each cold chamber. Three temperature conditions are utilized for maintaining apples with different shelf life. Accordingly, 73.554%, 19.008%, and 7.438% of products are kept under the temperature levels c1, c2, and c3, respectively at a total cost of US$ 17,520,000. In addition, operational costs, including laborer, packing, and transportation costs from producers to warehouses are about US$ 67,422. Transportation cost between warehouses to the wholesaler is about Table 15 Loss factors for different types of apples with different quality grades and minimum acceptable quality score based on expert's opinion 

Sensitivity analysis is performed on the case study to determine the behavior of the model by changing the main parameters. Clearly, there are many reasons affecting the selling price of products and may lead to fluctuation in demands. Changes in the selling price of products may be caused by specific economic conditions (e.g., economic impacts of COVID-19 outbreak), irreversible damage during production (e.g., frost, floods), etc. The customers' demand may be triggered by these factors or changes in taste and expectation, thereby becoming lesser or higher. From this point of view, it is vital to assess the model behavior during changes in parameters D it and C Q i js . In addition, this paper gives four vehicle types for fruits transportation and two cold storage rooms for 500 or 1000 Tons of apples that may change. The containers of different capacities can easily fit onto different forms of trucks and provide other options for transportation planning. Besides, redesign of cold rooms of less capacity may affect the feasibility of the problem. Finally, two remaining parameters are inventory holding costs and operational costs, which may be changed by embracing new technologies, identifying inefficiencies and unused services, etc. Therefore, this section provides an overview of the changes and their impacts to revise the plan after the occurrence. The results are summarized in Table 17 .

The results in Table 17 draw up the sensitivity of the model to parameters D it and H W c , where an increase of 25% in parameter D it makes the model infeasible, and a decrease of 50% in parameter D it reduce the total cost by 59%. Indeed, it is Fig. 9 Impacts of the changes on operational costs necessary to adjust capacity to match it with the changes in demands. We suppose that the available capacity of warehouses or producers is not sufficient in response to increased demands, therefore, the model becomes infeasible. In this case, the capacity of either warehouses or producers may step up to meet the demand. In addition, it is expected to purchase fewer products and avoid overstocks when demand decreases. This may cause a decrease in purchasing costs, operational costs, and inventory holding costs. Besides, a 50% increase in H W c leads to a 34.63% increase in the total cost. This indicates the importance of temperature condition adjustment and its impacts on the total cost. For further investigation, the effects of the changes in parameters H W c, C Q i js , D it on operational costs, inventory holding costs, and purchasing costs are presented as follows: Figure 9 demonstrates the effects of the changes on operational costs.

• The changes in parameter C Q i js have insignificant opposite impacts on the operational costs. • The changes in parameter D it have the greatest impacts on the operational costs, and a 25% decrease in demand leads to an 18% decrease in operational costs. Figure 10 shows the impacts of the changes on inventory holding costs.

• A decrease in parameter H W c significantly reduces inventory holding costs. Indeed, a decrease of 25% in H W c leads to a 25% decrease in inventory holding costs. • An increase in parameter D it significantly affects the inventory holding costs, and a 25% decrease in D it leads to a 36% decrease in inventory holding costs. Figure 11 shows the impacts of the changes on purchase costs.

• Changes in parameter D it have the greatest impacts on the total purchase cost. Indeed, a 25% decrease in D it leads to a 27% decrease in total purchase cost. • Total purchase cost increases by increasing the amount of parameter C Q i js . 

Nowadays, an appropriate configuration of agro-supply chains becomes all the more critical issues owing to specific characteristics of agro-products. Hence, in this paper, a multi-product multi-period MINLP model has been developed for an agro-supply chain contributing to purchase, facility location, and transportation decisions. Forasmuch as adjustment of temperature condition affects the quality of products and consequently their shelf lives, this study determines proper temperature conditions based on the varieties of products. This model aims to minimize the total cost of the supply chain considering several operational constraints. The primary step to solve the proposed model is to linearize the MINLP model and reformulate it into a MILP problem. Then, three sets of numerical examples, small, medium, and large sizes, have been presented to evaluate the performance of the mathematical model. The problem was coded in GAMS software and executed for small and medium test problems. Whatever the size of the problem increases, the computational time increases and the problem becomes difficult to solve. Therefore, BD algorithm has been implemented to tackle the complexity of the problem. Then GAMS/CPLEX solver and BD were compared in terms of objective function values and computational time. The results depict the convergence of BD to optimal solutions in a reasonable time.

Afterward, a precise assessment of the model was carried out by implementing a real case study. The results turn out the optimal purchase decisions of required quality score to be kept for a long time and number of cold chambers under proper temperature conditions. In addition, a sensitivity analysis was performed on the main parameters to predict the effects of the variations.

For future research, first, the model can be revisited for conditions under uncertain demands, thereby proposing efficient scenarios to overcome their associated shortage. Another possible extension of the model is to analyze purchase decisions under discount. This paper assumes only the products of the required quality for long time storage should be purchased. However, a new framework is promising to accommodate the distribution of products with short and long shelf lives that need precise inventory management.

Application of planning models in the agri-food supply chain: a review

Benders decomposition for the inventory vehicle routing problem with perishable products and environmental costs

Sustainable agro-food supply chain design using two-stage hybrid multi-objective decision-making approach

Temperature management for the quality assurance of a perishable food supply chain

Stochastic optimization models for joint pricing and inventory replenishment of perishable products

A benders-local branching algorithm for second-generation biodiesel supply chain network design under epistemic uncertainty

Closing loops in agricultural supply chains using multi-objective optimization: a case study of an industrial mushroom supply chain

Optimal joint replenishment, delivery and inventory management policies for perishable products

Agri-food supply chain management: literature review

A Benders' decomposition algorithm for optimizing distribution of perishable products considering postharvest biological behavior in agri-food supply chain: a case study of tomato

Benders decomposition for the distributionally robust optimization of pricing and reverse logistics network design in remanufacturing systems

Integrating harvesting decisions in the design of agro-food supply chains

A multi-echelon multi-product stochastic model to supply chain of small-and-medium enterprises in industrial clusters

A Lagrangian relaxation algorithm for optimizing a biobjective agro-supply chain model considering CO 2 emissions

Optimal purchase and inventory retrieval policies for perishable seasonal agricultural products

Optimization of distribution planning for agricultural products in logistics based on degree of maturity

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

City logistics for perishable products. The case of the Parma's food hub

Accessing food resources: rural and urban patterns of giving and getting food

A Benders decomposition approach for a real case supply chain network design with capacity acquisition and transporter planning: wheat distribution network

Agri-food supply chains with stochastic demands: a multiperiod inventory routing problem with perishable products

Coupling facility location models in the supply chain of perishable fruits

Agricultural modernization, structural change and pro-poor growth: policy options for the Democratic Republic of Congo

Planning models to optimize the agri-fresh food supply chain for loss minimization: a review

An integrated loss-based optimization model for apple supply chain

The impact of inventory management on economic and environmental sustainability in the apple industry

A memetic algorithm for bi-objective integrated forward/reverse logistics network design

Low carbon chance constrained supply chain network design problem: a Benders decomposition based approach

Optimizing fresh food logistics for processing: application for a large Chilean apple supply chain

A case analysis of a sustainable food supply chain distribution system-a multi-objective approach

A periodical flowering-harvesting model for delivering agricultural fresh products

Pricing and cold-chain service level decisions in a fresh agri-products supply chain with logistics outsourcing

Price and cold-chain service decisions versus integration in a fresh agri-product supply chain with competing retailers

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations