Integrated Harvest and Farm-to-Door Distribution Scheduling with Postharvest Quality Deterioration for Vegetable Online Retailing

With the rise of vegetable online retailing in recent years, the fulfillment of vegetable online orders has been receiving more and more attention. This paper addresses an integrated optimization model for harvest and farm-to-door distribution scheduling for vegetable online retailing. Firstly, we capture the perishable property of vegetables, and model it as a quadratic postharvest quality deterioration function. Then, we incorporate the postharvest quality deterioration function into the integrated harvest and farm-to-door distribution scheduling and formulate it as a quadratic vehicle routing programming model with time windows. Next, we propose a genetic algorithm with adaptive operators (GAAO) to solve the model. Finally, we carry out numerical experiments to verify the performance of the proposed model and algorithm, and report the results of numerical experiments and sensitivity analyses.


Introduction
With the development of mobile payment and e-commerce platforms, vegetable online retailing has become a trend for urban areas due to its convenience and diversity. According to reports, the online vegetable trading volume in 2017 was 131.3 billion CNY, accounting for 57.7% of the entire fresh e-commerce market [1]. Vegetable online retailing could provide a much time-saving shopping channel for urban citizens, especially official workers in central business districts, and brings in an emerging logistics model called farm-to-door distribution within the context of vegetable online retailing [2]. In fact, the vegetable farms are usually located in some rural areas that relatively closing to the urban districts so as to shorten the distribution distance and keep the freshness, which makes the farm-to-door distribution geographically possible. Compared with the traditional markets, vegetable online retailing and farm-to-door distribution are superior in reducing logistics process and time [3]. However, it brings much more challenges to order fulfillment due to the difficulties in both time and space collaboration. The critical issues in such a context are how to make joint decisions and deliver the vegetables rapidly in order to minimize the cost and quality deterioration. In this study, we aim to investigate an integrated optimization solution for the fulfillment problem of vegetable online retailing, which makes vegetable harvest decisions on the farm side and incorporates them into farm-to-door distribution scheduling with considerations of postharvest quality deterioration.
Specifically, the complication of integrated vegetable harvest and farm-to-door distribution scheduling comes from time constraints, vehicle routing, and vegetable postharvest quality deterioration. Firstly, a critical time collaboration problem decides the harvest time, distribution time, and time windows required in consumer orders. Secondly, the scheduling needs to arrange the number of vehicles to deliver orders as well as the optimal vehicle routings. Finally, since postharvest vegetable quality deteriorates over time, describing quality deterioration in the harvest and farm-to-door distribution process is complicated.
The most relevant literature is that regarding perishable products harvest or distribution problems [4], vegetable harvest planning [5][6][7], the inventory and distribution of agricultural products [8,9], delivery fleet size and truck routes [10], and the spatial and temporal model of the agricultural supply chain [11]. Generally, these previous studies have offered various decision tools to reduce the distribution cost and time of perishable products from different viewpoints, while the literature has seldom addressed the integrated optimization of vegetable harvest and farm-to-door distribution to minimize postharvest quality deterioration within an online retailing context.
Our contributions in this study contain three aspects: (1) model the postharvest quality deterioration of vegetables as a quadratic function; (2) formulate the integrated harvest and farm-to-door distribution scheduling as a quadratic vehicle routing programming model with time windows; (3) propose a genetic algorithm with adaptive operators (GAAO) to solve the model.
Our paper is organized as follows. In Section 2, we review the literature related to our study. In Section 3, we state the problem and formulate the model. In Section 4, we present the GAAO to solve the proposed model. In Section 5, we report the numerical results, sensitivity analyses, and the comparison of computational performance between the genetic algorithm (GA) and GAAO. Finally, we discuss the conclusions and suggest future works in Section 6.

Literature Review
Our work is related to two research streams: the quantitative method of agri-product postharvest quality deterioration, and the harvest and distribution process of agri-products.
The quantitative method can be grouped as a decay function and freshness function. The decay function refers to establishing an accurate model to evaluate the loss of agri-products. Ghare [12], Piramuthu and Zhou [13], and Rahdar [14] proposed a negative exponential function to describe the quality decay of perishable products. Rong et al. [15] presented the first-order reaction function of quality decay in the food supply chain. Wang and Li [16] defined quality decay as a key factor for perishable food and introduced it into a dynamic pricing strategy to control the product quality in the supply chain. Mohapatra and Roy [17], Keizer et al. [18], and Ferrer et al. [19] captured the average decay of heterogeneous agri-products by means of statistical methods. The freshness function refers to describing the relationship between the freshness of agricultural products and time from a positive perspective. Herbon et al. [20] depicted the freshness as a linear function. Cai et al. [21], Díaz-Pérez et al. [22], and Nekvapil et al. [23] introduced a continuous variable to measure the freshness level or deterioration and assumed the fresh index as a differentiable function. Amorim et al. [24] and Wang et al. [25] defined the freshness as the shelf life and incorporated it into the objectives to maximize the freshness of the perishable products.
The harvest and distribution process of agri-products contains harvest optimization and distribution optimization. For harvest optimization, Arnaout and Maatouk [26] and González-Araya et al. [27] studied the harvest operation of agri-products to enhance the utilization of resources and minimize the operational costs. An and Ouyang [28] proposed the relationship between harvest time and quality decay in an uncertain environment to reduce the postharvest loss and transportation cost. Ferreira et al. [29] and Takeda et al. [30] studied the harvest scheduling of fresh fruits and vegetables through machines to increase harvest efficiencies, such as the semi-mechanical system and the mobile platform machine. Ceccarelli et al. [31] discussed the relationship between the harvest maturity stage and cold storage time. For distribution optimization, Lamsal et al. [8] studied the agri-product optimal distribution scheduling from farm to warehouse. Ahumada and Villalobos [32] established an operational model for the distribution of fresh agri-products to maximize the profits of farmers. Amorim and Almada-Lobo [33] and Moons et al. [34] addressed a multi-objective vehicle routing planning problem for perishable foods to minimize the distribution cost and maximize the satisfaction. Many previous studies proposed Agronomy 2019, 9, 724 3 of 13 optimization models to minimize total costs, such as the cold chain distribution model [35], agri-product stochastic model [36], joint production and distribution model [37], and integrated production and distribution model [38].
Generally, in this study, there is a large amount of literature regarding the harvest planning and distribution scheduling of agri-products separately to reduce cost or quality deterioration. Few studies have addressed the integrated optimization of vegetable online retailing with postharvest quality deterioration. In practice, harvest planning and the distribution scheduling of vegetable online retailing are two highly interrelated issues that deserve further comprehensive research. Therefore, our work models a quadratic postharvest quality deterioration function, and incorporates it into the integrated optimization model that combines harvest planning and the farm-to-door distribution scheduling of vegetable online retailing, which can address more challenges of the agri-product supply chain.

Problem Description
The harvest and farm-to-door distribution process of vegetable online retailing contains the orders, harvest, vehicle assignment, and distribution. In this study, we shed light on the harvest, vehicle assignment and distribution to address an integrated harvest and distribution optimization for vegetable online retailing, as shown in Figure 1. In the figure, given some orders from consumers C = {c 1 , c 2 , . . . , c i }, they are firstly sent to a specific farm. For simplicity without loss of generality, assume that there are two kinds of vegetables on the farm (denoted as vegetable p 1 and p 2 ). When the harvest of vegetable p 1 is completed, the farm managers continue to harvest vegetable p 2 ; that is, the ending harvest time of vegetable p 1 (denoted as t kp 1 f ) is equal to the starting harvest time of vegetable p 2 (denoted as t kp 2 s ). After that, the farm managers select the capacity q k and quantity of vehicle K = k 1 , k 2 , . . . , k |K| according to the demand of consumers. Finally, the optimal vehicle routings are programmed considering the geographic location and time window [e i ,l i ] of consumers c i . costs, such as the cold chain distribution model [35], agri-product stochastic model [36], joint production and distribution model [37], and integrated production and distribution model [38]. Generally, in this study, there is a large amount of literature regarding the harvest planning and distribution scheduling of agri-products separately to reduce cost or quality deterioration. Few studies have addressed the integrated optimization of vegetable online retailing with postharvest quality deterioration. In practice, harvest planning and the distribution scheduling of vegetable online retailing are two highly interrelated issues that deserve further comprehensive research. Therefore, our work models a quadratic postharvest quality deterioration function, and incorporates it into the integrated optimization model that combines harvest planning and the farm-to-door distribution scheduling of vegetable online retailing, which can address more challenges of the agri-product supply chain.

Problem Description
The harvest and farm-to-door distribution process of vegetable online retailing contains the orders, harvest, vehicle assignment, and distribution. In this study, we shed light on the harvest, vehicle assignment and distribution to address an integrated harvest and distribution optimization for vegetable online retailing, as shown in Figure 1. In the figure, given some orders from consumers = , , … , , they are firstly sent to a specific farm. For simplicity without loss of generality, assume that there are two kinds of vegetables on the farm (denoted as vegetable and ). When the harvest of vegetable is completed, the farm managers continue to harvest vegetable ; that is, the ending harvest time of vegetable (denoted as ) is equal to the starting harvest time of vegetable (denoted as ). After that, the farm managers select the capacity and quantity of vehicle = , , … , | | according to the demand of consumers. Finally, the optimal vehicle routings are programmed considering the geographic location and time window [ , ] of consumers . Generally, the integrated harvest and farm-to-door distribution involves three decisions, i.e., (1) when to harvest, (2) which vehicles to use, and (3) how to arrange the optimal distribution routings. Generally, the integrated harvest and farm-to-door distribution involves three decisions, i.e., (1) when to harvest, (2) which vehicles to use, and (3) how to arrange the optimal distribution routings.

Objective Function
In this study, we consider two objectives, i.e., postharvest quality deterioration cost and distribution cost.

Postharvest Quality Deterioration Cost
Let γ p denote the speed of vegetable p quality deterioration with time t after harvest, and let β p denote a constant quality deterioration of vegetable p after harvest. Therefore, the rate of postharvest quality deterioration can be modeled as a linear function v(t), as given in Equation (1): For vegetable online retailing, the manager needs to control the postharvest quality deterioration of the vegetables from the starting harvest time t kps to the time t kpi at which they are received by consumer c i . In other words, the postharvest quality deterioration equals the integral of v p (t) in [t kps , t kpi ]. Let V p (t) denote the postharvest quality deterioration, which can be formulated as a quadratic function in Equation (2): Therefore, the objective function of postharvest quality deterioration cost (denoted as ob j 1 ) can be calculated in Equation (3): where d ip is the demand of consumer c i for vegetable p, and ω p is the unit postharvest quality deterioration cost of vegetable p. In the objective, we would minimize ob j 1 .

Distribution Cost
Let α k and λ k denote the fixed cost and distribution cost of vehicle k in unit time, and let τ ij denote the distribution time from node i to node j. The distribution cost objective function can be represented in Equation (4): where N = C ∪ {0} is a set of node numbers, x k ij ∈ {0, 1} is a decision variable, and the arc (i, j) belongs to the route of vehicle k if x k ij = 1. In the objective, we would minimize ob j 2 .

Constraints
To satisfy the time windows of consumers, we formulate the model as a vehicle routing problem with time windows (VRPTW). Let b p be the unit harvest time of vegetable p. x k i ∈ {0, 1} is a decision variable and the consumer c i is served by vehicle k if x k i = 1. Constraints are given as follows: i∈N k∈k j∈N k∈k Agronomy 2019, 9, 724 j∈C k∈K Constraints (5) and (6) ensure that each consumer can be served by one vehicle only. Constraint (7) represents the network flow balance constraint; that is, the vehicle arrives at a certain point and then leaves this point. Constraints (8) and (9) represent that each vehicle can only leave and return to the harvest location at most once. Constraint (10) indicates that vehicle k provides service to consumer c i if it is a node assigned to this route. Constraint (11) is a vehicle capacity limit. Constraint (12) is a vehicle quantity limit. Constraint (13) indicates the relationship between the starting harvest time and ending harvest time for each vehicle k and vegetable p. Constraint (14) means that the ending harvest time for vegetable p is equal to the starting harvest time for vegetable p + 1. Constraint (15) denotes the relationship between the ending harvest time and starting delivery time from the harvest location. Constraint (16) states that the distribution time of vehicle k cannot be later than the arriving time of consumer c i . Constraint (17) is a constraint eliminating the sub-tours problem. (17) is equivalent to t ki + τ ij − t k j ≤ M. Constraints (18)- (20) limit the time windows on harvest location and consumers. Constraint (21) defines x k ij , x k i as binary decision variables.

Model Formulation
According to the discussions in Sections 3.2 and 3.3, our formulation represents a constrained multi-objective optimization model. The integrated harvest and farm-to-door distribution scheduling model can be expressed as: This is a mixed integer nonlinear programming model that has proven to be an NP-hard problem [39]. The model can be described as a general problem including the traveling salesman problem (TSP) and vehicle routing problem (VRP). If Constraints (18)- (20) are removed, the model is equivalent to the VRP problem. If only one vehicle is provided for delivery, the model is equivalent to the TSP problem.

Genetic Algorithm with Adaptive Operators
With decision variables changing, the solution spatial structure of the mixed integer nonlinear programming model will change accordingly, which increases the difficulty of solving the model. We adopt GAAO to solve the integrated harvest and farm-to-door distribution model. As an improved version of GA, GAAO has a better performance in finding global optimal solutions and dealing with the large-scale decision variables problem. More details about the GAAO are described in Sections 4.1 and 4.2. Section 4.3 shows an implementation and flow chart of GAAO.

Controlling Parameter of Population Prematurity
In this section, we propose an assessment function to evaluate the degree of population prematurity. F tmax is the optimal individual fitness, and F tmax is the average fitness of optimal individuals; then: ω presents the premature degree of the population. The value of ω is calculated in order to avoid the negative effects caused by poor individuals, and to clearly show the convergence degree of individuals with greater fitness. It could describe the prematurity degree of individuals more accurately.

Self-Adaptive Crossover Probability and Mutation Probability
We use a nonlinear curve to adjust the adaptive crossover probability and mutation probability. According to this improvement, the excellent individuals are effectively preserved, and the ability of the poor individuals is enhanced. GAAO could leave out the local optimal solution and overcome the shortcomings of prematurity.
Based on ω described in Section 4.1, we propose the new strategy for adaptive adjustment coefficients. The crossover probability p c and the mutation probability p m could change with ω, as described as Equations (25) and (26): where l 1 , l 2 > 0, and ω ≥ 0, the range of p c is between [0.5, 1], and the range of p m is between [0, 0.5].
When the individual population tends to be discrete, the ability of developing better individuals for the population is enhanced. When the individual population tends to be convergent, the ability of producing new individuals for the population is increased.

Implementation of GAAO
Based on the discussion in Sections 4.1 and 4.2, the detailed steps of GAAO are as follows: Step 0: Initialize the population size N, evolution iteration T, crossover probability p c , and mutation probability p m .
Step 2: Calculate the fitness value for each individual according to Equation (22).
Step 3: Based on the fitness value, select the parents (denoted as P 1 ) from the current population through roulette selection.
Step 4: Generate the offspring (denoted as P 2 ) from parents P 1 according to the following steps.
Step 4.1: Based on Equations (24)- (26), implement the crossover and mutation operators adaptively to generate the offspring P 2 .
Step 4.2: Check each individual chromosome in the offspring P 2 to judge whether it satisfies Constraints (18)- (20). If it satisfies conditions, go to next step; otherwise, go to Step 4.1.
Step 5: Calculate the fitness value for each individual again according to Equation (22).
Step 6: Combine the initial population P and the offspring P 2 (denoted as P ∪ P 2 ).
Step 7: Replace the individual chromosome in the current population P, and select the first N chromosomes.
Step 8: Check the termination condition. If Gen < T, then go to Step 3; otherwise, output the optimal solution.
The flow chart of GAAO is shown as Figure 2:

Numerical Experiments
We examine the efficiency of the proposed model and algorithm through numerical experiments. In Subsection 5.1, we calculate the basic numerical results of GAAO by 15 pieces of consumer information and give the optimal vehicle routings. In Subsection 5.2, we test the sensitivity by varying the postharvest quality deterioration factors and give managerial insights. In

Numerical Experiments
We examine the efficiency of the proposed model and algorithm through numerical experiments. In Section 5.1, we calculate the basic numerical results of GAAO by 15 pieces of consumer information and give the optimal vehicle routings. In Section 5.2, we test the sensitivity by varying the postharvest quality deterioration factors and give managerial insights. In Section 5.3, we compare the results of medium-large scale numerical experiments between GA and GAAO from aspects of CPU time and iteration. All experiments are performed on an Intel Core i3-8100 3.60 GHz processor with 8GB RAM and the MATLAB 2017a. The details of the results are explained below.

Basic Numerical Experiments
In this case, the logistics network includes a vegetable farm (0) (which is viewed as the distribution center) and 15 consumers ( 1 ∼ 15). The consumer locations and online order demands, as well as the time windows of consumer orders, are all known.

Basic Numerical Experiments
In this case, the logistics network includes a vegetable farm (0) (which is viewed as the distribution center) and 15 consumers (1~15). The consumer locations and online order demands, as well as the time windows of consumer orders, are all known.
The distribution center has two kinds of vegetables, and  Figure 3 shows the spatial coordinate of each consumer (cited by Solomon's problem [40]). Table 1 gives the consumer information.    We test the basic numerical experiments 20 times, and the optimal results of GAAO are shown in Table 2 and Figure 4. Four vehicles are provided in this distribution service, and the objective value is 2687.3 CNY. The iteration time of GAAO is 186.3 s. The consumer with an earlier time window has a service priority, which indicates that optimal routes are not according to the shortest paths but are selected by the time window constraints in the model. We test the basic numerical experiments 20 times, and the optimal results of GAAO are shown in Table 2 and Figure 4. Four vehicles are provided in this distribution service, and the objective value is 2687.3 CNY. The iteration time of GAAO is 186.3 s. The consumer with an earlier time window has a service priority, which indicates that optimal routes are not according to the shortest paths but are selected by the time window constraints in the model.

Sensitivity Analyses on Postharvest Quality Deterioration Factors
In order to research the relationship between the postharvest quality deterioration factor γ and the objective value, we vary γ and γ from 0.1 to 0.9, respectively. As shown in Figure 5a, the objective value increases as the postharvest quality deterioration factors increase. Figure 5b demonstrates the magnitude of variation in objective value with γ . In the case of γ and γ both being small, due to the slow decay of vegetables, the farm managers have more time to fulfill orders, which makes little gaps between the objective values. In other words, the farm managers can ignore the postharvest quality deterioration cost and select the optimal route with the minimum distribution cost. In the case where γ and γ are both large, there is a significant jump in the objective value, since the time sensitivity of the vegetables increases substantially. That is to say, the farm managers need to arrange more vehicles in this case to fulfill the vegetable distribution in a shorter time, which leads to higher costs inevitably.

Sensitivity Analyses on Postharvest Quality Deterioration Factors
In order to research the relationship between the postharvest quality deterioration factor γ p and the objective value, we vary γ 1 and γ 2 from 0.1 to 0.9, respectively. As shown in Figure 5a, the objective value increases as the postharvest quality deterioration factors increase. Figure 5b demonstrates the magnitude of variation in objective value with γ p . In the case of γ 1 and γ 2 both being small, due to the slow decay of vegetables, the farm managers have more time to fulfill orders, which makes little gaps between the objective values. In other words, the farm managers can ignore the postharvest quality deterioration cost and select the optimal route with the minimum distribution cost. In the case where γ 1 and γ 2 are both large, there is a significant jump in the objective value, since the time sensitivity of the vegetables increases substantially. That is to say, the farm managers need to arrange more vehicles in this case to fulfill the vegetable distribution in a shorter time, which leads to higher costs inevitably. In general, we can draw two managerial insights. On one hand, controlling the postharvest quality deterioration factor of vegetables is an effective method to reduce the total cost, such as cold-chain transportation or food preservation techniques. On the other hand, when the perishability of vegetable is high, the farm managers could shorten the distribution time by increasing the vehicle number to satisfy the fresh quality requirement of consumers.

Medium to Large-Scale Numerical Experiments and Comparisons Result between GA and GAAO
In this section, we use 10 different medium-large scale numerical experiments to demonstrate the advantage and efficiency of GAAO. The number of consumers increases from 20 to 110.
We report the computational results of GA and GAAO, and compare the time gap and iteration gap in Table 3. The maximum iteration number of GAAO is 432, while the minimum iteration number of GA is 425. The iteration number of GAAO is significantly less than that of GA, which indicates that GAAO is more efficient. We draw a conclusion that GAAO is superior to GA because GAAO saves at least 10% CPU time and has at least a 20% iteration gap. In other words, compared with GA, GAAO can reduce the iteration and converge the optimal solution in a shorter amount of time.  In general, we can draw two managerial insights. On one hand, controlling the postharvest quality deterioration factor of vegetables is an effective method to reduce the total cost, such as cold-chain transportation or food preservation techniques. On the other hand, when the perishability of vegetable is high, the farm managers could shorten the distribution time by increasing the vehicle number to satisfy the fresh quality requirement of consumers.

Medium to Large-Scale Numerical Experiments and Comparisons Result between GA and GAAO
In this section, we use 10 different medium-large scale numerical experiments to demonstrate the advantage and efficiency of GAAO. The number of consumers increases from 20 to 110.
We report the computational results of GA and GAAO, and compare the time gap and iteration gap in Table 3. The maximum iteration number of GAAO is 432, while the minimum iteration number of GA is 425. The iteration number of GAAO is significantly less than that of GA, which indicates that GAAO is more efficient. We draw a conclusion that GAAO is superior to GA because GAAO saves at least 10% CPU time and has at least a 20% iteration gap. In other words, compared with GA, GAAO can reduce the iteration and converge the optimal solution in a shorter amount of time.  Figure 6a shows that the trend of GAAO iterations changes from 100 to 300 gradually with the increasing experiment scales, and that of GA iterations ranges from 420 to 900. From Case 1 to Case 6, the iteration gaps between GAAO and GA are all above 50%, but in Case 7 to Case 10, the iteration gaps are smaller. This difference indicates that the GAAO has a distinct solution advantage when the consumer number is less than 70. Figure 6b demonstrates the comparison of CPU time between GAAO and GA. In these cases, GAAO is superior to GA significantly in CPU time. Generally, the performance of GAAO proposed in this paper is consistently better than that of GA in solving the integrated optimization of harvest and farm-to-door distribution, and could obtain an optimal solution more quickly.  Figure 6a shows that the trend of GAAO iterations changes from 100 to 300 gradually with the increasing experiment scales, and that of GA iterations ranges from 420 to 900. From Case 1 to Case 6, the iteration gaps between GAAO and GA are all above 50%, but in Case 7 to Case 10, the iteration gaps are smaller. This difference indicates that the GAAO has a distinct solution advantage when the consumer number is less than 70. Figure 6b demonstrates the comparison of CPU time between GAAO and GA. In these cases, GAAO is superior to GA significantly in CPU time. Generally, the performance of GAAO proposed in this paper is consistently better than that of GA in solving the integrated optimization of harvest and farm-to-door distribution, and could obtain an optimal solution more quickly.

Conclusions
In this paper, we propose an integrated optimization model that simultaneously considers vegetable harvest and farm-to-door distribution. To capture the perishable feature of vegetables, we formulate a quadratic postharvest quality deterioration function. The optimization model integrates the decisions regarding the harvest time, delivery vehicle assignment, and farm-to-door distribution planning in terms of consumers' demand requirements and time windows. Moreover, we design a GAAO to solve the proposed model. Finally, numerical results are reported to show the efficiency and effectiveness of the proposed model and solution algorithm, and managerial insights are offered for vegetable online retailing in practice. The numerical results show that the GAAO can obtain the optimal solution with an average 24% CPU time and an iteration number that is 48% less than that of GA. Generally, compared with the traditional vegetable supply chain, our study suggests that integrating harvest and farm-to-door distribution scheduling can reduce the postharvest quality deterioration of vegetables and subsequent distribution costs dramatically.
Our work serves as an initial effort on integrated harvest and farm-to-door distribution scheduling for vegetable online retailing. Several promising directions for future research remain. Firstly, all order information for a certain period in this model is assumed to be determined. In practice, the consumer orders arrive randomly, and such information could be forecasted in advance, which will lead to a scenario-based stochastic programming model. Secondly, we only consider common heterogeneous vehicles in this paper. Considering different types of refrigerated vehicles with variable capacity and multiple compartments for perishable products is also a future

Conclusions
In this paper, we propose an integrated optimization model that simultaneously considers vegetable harvest and farm-to-door distribution. To capture the perishable feature of vegetables, we formulate a quadratic postharvest quality deterioration function. The optimization model integrates the decisions regarding the harvest time, delivery vehicle assignment, and farm-to-door distribution planning in terms of consumers' demand requirements and time windows. Moreover, we design a GAAO to solve the proposed model. Finally, numerical results are reported to show the efficiency and effectiveness of the proposed model and solution algorithm, and managerial insights are offered for vegetable online retailing in practice. The numerical results show that the GAAO can obtain the optimal solution with an average 24% CPU time and an iteration number that is 48% less than that of GA. Generally, compared with the traditional vegetable supply chain, our study suggests that integrating harvest and farm-to-door distribution scheduling can reduce the postharvest quality deterioration of vegetables and subsequent distribution costs dramatically.
Our work serves as an initial effort on integrated harvest and farm-to-door distribution scheduling for vegetable online retailing. Several promising directions for future research remain. Firstly, all order information for a certain period in this model is assumed to be determined. In practice, the consumer orders arrive randomly, and such information could be forecasted in advance, which will lead to a scenario-based stochastic programming model. Secondly, we only consider common heterogeneous vehicles in this paper. Considering different types of refrigerated vehicles with variable capacity and multiple compartments for perishable products is also a future research direction. Finally, our model sets most factors to certainty. Incorporating some uncertain factors into harvest and distribution scheduling is an important future research direction, such as the uncertainty of temperature or road conditions.