The Impact of Lead Time Uncertainty on Supply Chain Performance Considering Carbon Cost

: In supply chain operation practices, lead time uncertainty is a common management issue. Uncertain lead time can lead to increased inventory costs and unstable service levels, which will directly affect the overall operation performance of the supply chain. While considering environmental performance in supply chain, it is important to understand how an uncertain lead time will affect sustainable performance. In this paper, we constructed a supply chain model with stochastic lead time and explored the relationship between uncertain lead time and supply chain performance. We considered carbon cost, inventory cost, and service level as a supply chain performance. System dynamics methodology was employed to observe and explore the dynamic change trend of the overall performance in the complicated supply chain model. This was done under both different levels of lead time standard deviation and different order policies. The results demonstrate how stochastic lead times can significantly increase inventory costs and carbon costs. Therefore, we propose appropriate ordering policies which mitigate the negative impacts of stochastic lead times. validation, Z.L., W.F. and E.Z.; formal analysis, Z.L., W.F. and E.Z.; investigation, W.F. and X.C.; resources, Z.L., Y.G. and E.Z.; data curation, W.F. and X.C.; writing—original draft preparation, Z.L. and W.F; writing—review and editing, E.Z., Y.G. and X.C.; visualization, E.Z. and X.C.; supervision, Y.G.; project administration, Y.G.; funding acquisition, Z.L.


Introduction
The concept of a low carbon economy has come into existence to cope with numerous environmental challenges. Due to government regulations and public awareness, enterprises have realized the importance of introducing low-carbon operations and sustainability into their management objectives. Agi et al.
[1] conducted a study on Australian firms and found two main motives that caused enterprises to shift towards green development. The first motive is to comply with environmental regulations, and the second motive is to achieve a competitive edge. Diabat et al. [2] indicated that the sustainable supply chain management (SSCM) system could ensure environmentally friendly practices in traditional supply chains with essential enablers. On the one hand, government policies and regulations, such as carbon caps, are putting pressure on enterprises to incorporate environmental factors into their business models. On the other hand, a number of enterprises have made sustainable green operations as their competitive advancement. In supply chains, several activities produce carbon emissions, for example, transportation and warehousing. Hence, it is critical to manage and control carbon emissions in supply chain operations.
Global trade boosts have shifted the short-distance supply chain within a country to long distance supply chain between countries [3]. This change has resulted in an increase in the lead time and has accelerated the fluctuation. The lead time is also called replenishment delay, which is the time gap between sending an order and receiving it. Existing studies demonstrate that lead time is a key factor that affects supply chain performance [4]. If the lead time is uncertain, the situation is complex from both the cost and service perspectives. Our research is motivated by the growing concern of business, which is forcing business enterprises to reduce carbon emission cost in their supply chain network. Therefore, this paper includes carbon emission costs as one of the supply chain performance indicators. Some studies [5,6] have demonstrated that the uncertain lead time is directly related to supply chain performance. Therefore, this paper examines the effect of uncertain lead time on carbon emission costs. Further, the paper investigates whether the effect of lead time increase on the behavior of carbon cost emissions is the same as the behavior of inventory costs and service levels or not.
When both market demand and lead time are uncertain, the supply chain system becomes more complicated and dynamic with time [7]. The traditional operational method makes it difficult to explore the complicated supply chain system. System dynamics provides a new way to study complex systems. The discipline of system dynamics was founded in 1956 by Professor Forrester of the Massachusetts Institute of Technology. The system dynamics discipline integrates systems theory, cybernetics, information theory, and organization theory with the help of computer technology. The system dynamics model can also show the dynamic behavior of the object system. The system dynamics can be used to better analyze non-linear systems which are caused by many real situations in the supply chain system, such as stochastic lead time, uncertain demand, non-negative order, etc. In this study, we applied a system dynamics method to build a supply chain model and explore the dynamic change trend of supply chain performance indicators, including inventory cost, service level, and carbon cost, under uncertain lead times.
The reminder of this paper is organized as follows: Section 2 presents an overview of the relevant literature. Section 3 presents the assumptions and models and introduces the performance indicators; Section 4 conducts model analysis and numerical simulations. Section 5 focuses on the analysis of experimental results. Finally, Section 6 summarizes the main research conclusions and proposes several future research directions.

Literature Review
Research was conducted on how to reduce the negative impact on the environment during supply chain operations [8,9]. Quite a number of studies focus on transportation studies [10][11][12]. This paper mainly discusses the dynamic change of carbon emission cost in an uncertain supply chain environment with stochastic lead time, and explores the inventory management approach to buffer the fluctuation.
In the research on carbon inventory management, the focus has been on newsvendor models and EOQ (economic order quantity) models [13][14][15]. Song et al. [16] analyzed a newsvendor model and demonstrated that it is possible to achieve multiple targets, including carbon emission reduction, cost saving, and profit improvement at the same time by applying a suitable order policy. Chen et al. [17] also found that there are some conditions under which the reduction in carbon emissions is greater than the increase in cost. Hovelaque et al. [18] presented a novel model which considered the relationship among inventory policy, carbon emissions, and price account. They found that certain environmental policies can significantly decrease total carbon emissions.
Most recent research focuses on different operational factors in the supply chain while simultaneously considering carbon emissions. Benjaafar et al. [19] integrated carbon emissions into supply chain operational process with inventory management, production, and procurement. They also modified traditional models to account for both general cost and carbon emissions. Diabat et al. [20] proposed a genetic algorithm approach for supply chain integration while considering carbon emissions. They considered multiple scenarios such as different demand models and multiple distribution centers. Park et al. [21] investigated punishment mechanism in supply chain structure and found that supply chain members are affected differently by carbon emission reduction. Konur et al. [22] incorporated environmental factors into a supply chain control model with uncertain demand to study the effect of logistics policies and supplier selection on environmental performance. Xu et al. [23] considered demand preferences in low carbon supply chain management to analyze the effect of consumer preference on inventory cost. Hoen et al. [24] studied the effects of transportation facilities on carbon emissions under uncertain demand conditions. They showed that the lead time directly affects carbon emissions on a transport facility. Sarkar et al. [25] reached a similar conclusion. Finally, Arikan et al. [26] studied environmental effects related to lead time uncertainty for transportation facilities. Recent studies on low carbon inventory management have begun to account for an uncertain lead time effect, but the research in this field is not sufficient. Many existing studies were based on order-up-to-level inventory models. For example, Chatfield et al. [27] studied stochastic lead time and established order-up-to-level inventory models of periodic replenishment. They conducted computer experiments to explore the impact of stochastic lead time on the inventory system. The results showed that the uncertainty of lead time aggravates the fluctuation of ordering quantity. Yet, efficient sharing of quality information is beneficial to reduce the fluctuation. Du et al. [28] extended the demand model to an autoregressive moving average model (ARMA) with independent and identically distributed lead times. However, these two studies both used the order fluctuation as the only evaluation indicator to study the uncertain lead time effect. Many studies added the bullwhip effect as the evaluation index. The bullwhip effect is the magnified effect of the demand information distortion while moving from downstream to upstream in the supply chain. The bullwhip effect states that a small fluctuation in demand at the downstream level (i.e., at a retailer level) creates larger fluctuations in demand at the upstream level (i.e., wholesale, distributor, manufacturer, and raw material supplier levels). In supply chain management, each player has limited information over demand and, thus, the fluctuating demand pattern at any supply chain link influences the entire supply chain with its forecasting inaccuracy (i.e., ordering too much or ordering too less). It is considered to be one of the main factors that weakens supply chain performance, because it creates either excess inventory or shortages in the whole supply chain network. A detailed explanation about the bullwhip effect can be found in the literature [27]. Boute et al. [29] considered both order fluctuation and inventory fluctuation. By conducting numerical simulation experiments, they found that the bullwhip effect and inventory costs can be reduced at the same time. However, this study assumed that the supply chain is linear and that there are no order crossovers. In fact, an in-depth study by Saldanha et al. [30] affirmed that the order crossover is a common phenomenon. Chaharsooghi et al. [31] investigated the bullwhip effect in multi-echelon supply chain models and concluded that random lead time has a greater impact on bullwhip effect over fixed lead time.
Jian et al. [32] considered the influence of the randomness of lead time on the inventory model for perishable drugs. Through empirical and sensitivity analysis, they determined that the weightdistribution values of the shelf life of perishable drugs and the service level are reasonable. Spiegler et al. [33] discussed the influence of lead time uncertainty on order quantity, inventory level, and work-in-process quantity. Assuming demand is stable, they followed the small disturbance principle to probe lead time disturbance related to system output, and they found that the order quantity, inventory level, and work-in-process quantity increase with the lead time. Moreover, the reduction in lead time can result in a dynamic change of damping. This result implies that, although reducing the lead time can improve system performance, a flexible management plan is required to avoid unnecessary system turbulence. The above study provides a theoretical basis for setting up the parameter in this paper. Dejonckheere et al. [34] proved that an order-up-to-level order policy creates a bullwhip effect. However, if two feedback loops are introduced to the order policy, the bullwhip effect can be weakened and, sometimes, even eliminated. This order policy is called APVIOBPCS (automatic pipeline, variable inventory, and order based on production control system). The APVIOBPCS model has two feedback loops to control the order quantity by discrepancies in inventory and work-in-process. In addition, by changing the adjustment parameters, APVIOBPCS policy clusters are generated as various ordering policies. When the adjustment parameters of inventory and work-in-process inventory are both 1, it becomes an order-up-to-level policy. The APVIOBPCS model has been applied in practice and in real time [35][36][37].
Sustainability 2018, 10, x; doi: FOR PEER REVIEW www.mdpi.com/journal/sustainability This paper differs from the existing studies in the following aspects. First, the model in this paper considers both uncertain demand and uncertain lead time, which is near reality situation. Correspondingly, there is an order crossover phenomenon in the supply chain model which only emerges when lead time is stochastic. These conditions make the model more complex and nonlinear. Second, this paper introduces carbon cost into the supply chain model as a new systematic performance indicator. It compares average cost and average service level change to carbon cost change under varying supply chain environments. This paper also adopts the system dynamics method to discuss the impact of stochastic lead time on supply chain inventory system, and also find more suitable order policies than the traditional order method, which can dampen the negative effect of uncertain lead time.

Assumptions
In order to explore the impact of stochastic lead times on supply chain inventory systems, this paper constructed a three-stage supply chain with a distributor as the focal point in a periodic continuous review system based on the APVIOBPCS model. The proposed model makes the following assumptions: Assumption 1. At the beginning of a single period, the distributor receives the products from the manufacturer, and the distributor's ordering plan is submitted to the manufacturer at the end of the period; at the same time, the system state is updated. Assumption 2. The upstream manufacturer has abundant production capacity and can fully meet the orders from the distributor without shortages. the demand of each period may not be fully met due to variable market demand. Assumption 4. The demand forecasting for the next period is based on the exponential smoothing method, with the smoothing factor θ .

Assumption 5.
The distributor supplies products to the retailer with the SFS model (supply from stock). That is, the minimum of demand and current inventory is supplied. The shortage will not be replenished in the current period, and it needs to be submitted by downstream retailer in the next period in the form of new demand orders. Hence, inventory is allowed to take negative value. It is noted that this is consistent with business practice. Assumption 6. The distributor's orders are subject to non-negative constraints. In business reality, it can be interpreted as not allowing returns.

Stochastic lead times
This paper investigates the influence of stochastic lead times on the performance of supply chain inventory level. While introducing stochastic lead times into the model, it is essential to consider the allowance of order crossover. In order crossover, the sequence of order arrival can be different from the sequence of order placement. Order crossover is due to stochastic lead time. Due to randomness, multiple order batches from several previous periods could arrive in one period at the same time, or the distributor could receive those orders in advance, which are placed in a later period. According to a global logistics survey [38], crossover orders accounted for more than 40% of total orders. Thus, this paper facilitates the order crossover in the supply chain system.
The logistics data survey in study [38] shows that the lead time distribution is random. Some of the lead time distributions are normal, and some of them are close to normal with long tail. This paper simulates stochastic lead time distributions based on the observations found in the study [38]. We supply chain event process of this model, a single period in the system starts with the distributor receiving the products from the manufacturer and ends at the distributor issuing the current order plan to the supplier. Submission of the order plan is the last step of each period. Thus, there is a natural delay in the periodic cycle when the distributor receives the homologous products of his order plan. Therefore, the earliest and the latest arrival times of the products in the order plan submitted in period t are recorded as 1 t + and 1 t L + + + , respectively. In simulation experiments, the lead time uncertainty can be simulated by the following steps: Step 1: The distribution of random lead time is assumed in the model. MATLAB is used to generate n random numbers representing the lead time in the period; these are stored in an array L . There are five presumed types of lead time distributions corresponding to five experiments.
Step 2: Next, the array L is introduced into the basic inventory simulation model, and the delay, corresponding to the order quantity t OQ and determined at the end of the simulation period t, is denoted as t L .
Step 3: Since each t L is generated in step 1, the arrival time of t OQ can be calculated at the beginning of each period and recorded in the arrival time array as: ( ) 1 t arrivetime t t L = + + .
Step 4: In each period, the arrivetime array from period i to period i L + − is checked to determine whether ( ) arrivetime t equals i . If so, the corresponding arrival quantity can be included in the current inventory. Otherwise, it is counted as the work-in-process inventory.

Models
The basic inventory model of this paper is based on the APVIOBPCS (Automatic Pipeline Inventory and Order based Production Control System) model family. The APVIOBPCS model is a collection of inventory policies, which can express different inventory policies with different parameter settings. The common order-up-to-level inventory order policy, which has been used in practice, is also a form of APVIOBPCS.
The core enterprise of the supply chain model in this paper is the distributor. The upstream of the distributor is the manufacturer and downstream is the retailer. According to the principle of the APVIOBPCS production control, the model is assumed to be a periodic review system. The state of each system variable would be updated as the following process.
At the beginning of the period t , the distributor receives the previously placed orders. Due to the order crossover as mentioned above, in period t , products purchased in multiple historical periods might arrive at the same time, or the distributor may receive the recently ordered products in advance. Accordingly, the arrival quantity t AQ is equal to the sum of orders the manufacturer is delivering t SQ before several periods based on the t L , as shown in Equation (1). Then, the distributor can count the inventory t BINV at the beginning of period t using Equation (2), which is the sum of the inventory at the end of the last period and the current arrivals. Arrivals: Inventory at the beginning: The distributor delivers the products to the retailer from stock to meet market demand t D ，and delivery volume t DQ is the minimum of t D and the initial inventory t BINV at period t , as expressed in Equation (3). If the value of t BINV is less than t D , stock out t SOQ occurs, which can be calculated using Equation (4). In addition, according to the assumptions in Section 3.3, the shortage t SOQ in the current period would not be fulfilled by the distributor in the next period. When there is a shortage t SOQ in the current period, products are replenished in the form of new demand orders submitted by downstream retailers in the next periods. This policy is consistent with the practice used in many manufacturing industries. Therefore, the shortage t SOQ is not in Equation (3).
Delivery quantity from distributor to market: Shortage: Equation (5)  with the smoothing factor θ , as in Equation (7).
Work-in-process inventory: Demand forecast: Equation (8) reveals the mathematical model of the order policy, which is the key point of the APVIOBPCS model. In order to reduce errors caused by the forecast, two feedback loops were introduced into the model: Adjustment for stock and adjustment for work-in-process. The adjustments for stock and work-in-process inventory discrepancies are represented into the order equation Order: The target inventory in the model is denoted by  Table 1. The correlation among the variables is shown in an inventory system dynamics model diagram as in Figure 2.

CC
The carbon cost in period t AC The average inventory cost ASL The average service level    [41], the carbon emission per inventory product can be calculated from the conversion relationship between electricity consumption and carbon emissions. In addition, the sensitivity analysis shows that the carbon emissions per product only affect the total cost of inventory carbon emissions, and would not change its fluctuation trend, which is the focus of this paper. In this paper, Ip is defined as the inventory carbon cost per product. The inventory carbon cost t IP of the period t can be obtained using Equation (12). Equation (13)

Average Inventory Cost
Most scholars define the system inventory cost as the combination of inventory holding cost per unit h and shortage cost per unit b . This paper adopts the same method for the calculation of inventory cost, as defined in Equation (15).

Average Service Level
Macchion et al. [42] assumed the service level t SL of each period to be 0 or 1. In the case of shortages, the service level is considered 0. Otherwise, it is considered 1. In this paper, shortage is introduced to further quantify the service level. As shown in Equation (16), the ratio of shortage t SOQ and demand t D of period t is considered when calculating the service level t SL . Furthermore, the average service level ASL is calculated using Equation (17):

Model Verification
In order to better emulate business practice, a nonlinear model is constructed in this paper. However, it greatly complicates the model, and makes it more difficult to find the optimal analytic solution. Hence, simulation is employed in this paper. In addition, the verification of a simulation model is a vital aspect of using the simulation model. In this paper, the theoretical derivation of the APVIOBPCS model under fixed lead time is used to ensure the correctness of the computer simulation.
Sustainability 2018, 10, x; doi: FOR PEER REVIEW www.mdpi.com/journal/sustainability Disney et al. [43] showed that the bullwhip effect value BWV changes with the adjustment parameters of in-stock and work-in-process stock, when the lead time is fixed, as shown below: T T  T T  T  T  BWV  T  TT  T  T T T  Fixed lead time can be regarded as a special case of stochastic lead time. In order to compare with the theoretical value of Equation (18), the standard deviation of lead time in the simulation verification model was set to 0. The simulation experiment was carried out by MATLAB 2016a to calculate the order fluctuation value under different control parameters. Figure 3a shows a contour diagram generated by Equation (18), which displays the bullwhip effect value with different parameters under fixed lead time and under no order-quantity constraint. Figure 3b was generated by simulation experiments, while Figure 3a was generated using theoretical values. Comparing  Figures 3a and 3b shows that the two graphs are basically the same. This means that the simulation results are consistent with the theoretical values, which validates our simulation model.
It is worth mentioning that the top left and bottom right in Figure 3a are different from those in Figure 3b. The white parts of the upper triangle and lower triangle areas of the coordinate planes in Figure 3a, b represent the large values or negative bullwhip effect values, which implies poor policies. We drew Figure 3a using MATLAB 2016a according to Equation (18), where we set the program judgement conditions about white area. In our program, the white area represents the places where the bullwhip effect value is larger than 18 or the value is negative. In fact, in Equation (18), if both its numerator and denominator are negative, the bullwhip value is incorrect. However, we omitted this condition intentionally to distinguish the two figures. When we consider the situation where both the numerator and the denominator are negative, the top left and bottom right in Figure 3a show as a white area, and the two figures remain almost same.

Model Parameter Analysis
In this paper, Equations (1)- (8) relate to the decisions of the supply chain. Equations (9)-(17) are mainly used to measure the supply chain performance. There are many parameters in the model, but the most critical parameters can be determined by the transfer function derived from the ztransforms of Equations (1)- (8). According to the deduction process shown in the Appendix A, the decision variables are all transformed to the z -threshold. Consequently, the system transfer equation between the order OQ and demand D can be obtained as follows: Equation (21) shows the stability of the system, which is mainly affected by S α , SL α and t L , t L , where ( 1 ) z θ − + represents the stable pole of the system. This deduction further verifies and expands the results reported in [44], which concludes that the parameter combination is the key factor for determining the boundaries among the stability, the stable oscillation, and the exponential oscillation of the system.

Experiment Design
The model verification and the parameter analysis so far attest that the parameters S α and SL α affect the dynamics of the system, and the various parameter values correspond to different ordering policies. Many researchers have explored the practical impact of the parameter setting on system stability performance with the POUT policy, that is, S These researchers argue that POUT policy can improve system stability in a continuous review system [45]. Therefore, this paper adopts the POUT policy, = S S L β α α = , to explore the effect of stochastic lead times on the inventory system in the range of [0.1,1.9] for parameter β .
In order to reflect the different levels of lead time uncertainty, this paper assumes that lead time follows stochastic distribution with the average of 2 L μ = and different standard deviations. After a number of trials, five arrays of lead time with standard deviations of 0, 0.5, 1, 1.8, and 2 were screened. The standard deviation of 0 L σ = represents the fixed lead time. This paper takes the system carbon cost C C , average inventory cost A C , and average service level ASL as the performance evaluation indicators. Under each lead time distribution, the three performance indicators were observed with the help of MATLAB 2016a. We obtained 15 experiments, as depicted in Table 2, and the experimental results are shown with different order policies in Appendix B. Table 2. Experiment design.

Performance indicators
Lead time variance For minimizing the impact of contingency on the test results, the simulation experiment was conducted for 1000 periods with = 0 ω . According to the actual survey data, demand obeys the independent random distribution with the mean value of 5,021 and the standard deviation of 956.48. Oliveira et al. [46] proved that the fluctuation of the system is relatively small with a smaller smoothing factor θ . Hence, the demand forecast parameter is set as = 0.4 θ .
Meanwhile, the values of the constants in the model were set according to the real data obtained from the freight market to ensure the validity of the simulation experiments, as listed in Table 3.     We now discuss the rate of change for each performance indicator caused by the stochastic lead times. The change rate is expressed in the form of growth rate as IR , and AIR represents the average growth rate of the three indicators in Table 4. Overall, when the standard deviation of lead time increases, the IR of the carbon cost C C and the average inventory cost A C remain positive, while the average service level ASL remains negative. In other words, higher lead time uncertainty deteriorates the performance of the inventory system. The lead time uncertainty negatively affects A C , whose average growth rate AIR reaches 17.25%, as shown in Table 4. In addition, the effect of random lead time on ASL is slightly greater than that on the C C , and the average growth rates are -1.62% and 1.29%, respectively.  Table 4 shows that the greater the lead time uncertainty, the greater the change rate of the performance indicators. For the carbon cost C C , when the standard deviation of lead time L σ increases to 0.5, the growth rate IR remains 0.62%. However, when L σ increases from 0.5 to 2, the change rate of C C for all stages remains above 1%. The IR decreases initially, and then rises. This implies that lead time uncertainty will lead to an increase in C C , but the rate of cost increase will fluctuate. When the L σ changes from 0 to 1, average inventory cost A C and average service level ASL rises to 39.32% and 3.38%, respectively. Then, as the fluctuation of lead time intensifies, the rates keep declining. In the transition range of 1.8-2, they fall to 2.82% and 0.23%, respectively. The results summarized in Table 4 are consistent with what Figure 4 shows.

Effect of Stochastic Lead Times on Order Decision
In this subsection, we focus on the relationship between stochastic lead times and inventory systems. It can be observed from Figure 4 that when the standard deviation of lead time L σ increases from 0 to 0.5, the curves of carbon cost C C , average inventory cost A C , and average service level ASL remain relatively smooth, and they are least affected by the change of parameter β . In contrast, when L σ increases from 1 to 2, these three performance indicators change significantly. In particular, the average inventory cost A C and average service level ASL continue to rise and fall with the increase of parameter β . However, the carbon cost curve shows a completely different trend. Figure 4a shows that when L σ is small, the curve remains roughly smooth, which indicates that when the lead time is short, the parameter has little influence on carbon cost. In contrast, with an increase in L σ , the minimum point of the curve gradually decreases, and the decisionmaking area of ordering is reduced. Hence, there exists an optimal order policy to minimize carbon costs. This observation implies that lead time uncertainty influences the adjustment of ordering decisions on carbon cost.
Appendix B is a data table of a computer simulation experiment, which reveals inventory system performance indicators for different standard deviations of lead time and different β values. In comparison, the variance of A C is the largest, followed by the variance of C C , which is then followed by the variance of ASL -and this ranking is not impacted by different values of the standard deviations. Regardless of the uncertainty degree of lead time, the average inventory cost Sustainability 2018, 10, x; doi: FOR PEER REVIEW www.mdpi.com/journal/sustainability A C is most affected by the parameter β , followed by C C and ASL . However, the uncertainty level of lead time changes the strength of the adjusting function of parameter β . See Table 5.  Table 5 depicts the values and the relative change rates of three indicators under = 0, = 0.5， = 1, = 1.8, and = 2 with the parameter β being 0.2, 0.8, and 1.8, respectively. Different values of β represent different order policies. IR1 indicates the growth rate of indicators when β changes from 0.2 to 0.8, while IR2 indicates the growth rate of indicators when β increases from 0.8 to 1.8. Again, C C represents the carbon cost, A C represents the average inventory cost, and ASL represents the average service level. Table 5 shows that when the β increases from 0.2 to 0.8 under fixed lead time, C C decreases by 0.15%, A C increases by 0.95%, while ASL decreases by 0.07%. However, when β changes from 0.8 to 1.8, C C increases by 1.25%, A C increases by 8.35%, and ASL decreases by 0.24%. Hence, it can be concluded that under fixed lead times, the policies of = 0.2 β and = 0.8 β have no significantly different impacts on the inventory system, but they are both better than the policy of = 1.8 β . When the standard deviation of lead time is 0.5, the conclusion is similar.
When the standard deviation of lead time increases to 1, 1.8, and 2, the impacts of the different ordering policies on the system change. When β increases from 0.2 to 0.8, C C decreases more than 3.2%, which is much larger than when = 0 or = 0.5. A C increases by more than 4.3%, except when lead time standard deviation is 1. 8

Summary of Experimental Results
Through the above experimental data analysis, we can see that: 1) Stochastic lead times mainly affect the strength of the inventory system performance. The experimental results show that random lead time has an insignificant impact on the changing trend of each indicator. However, it greatly alters the values of the performance indicators. The sharper the lead time fluctuates, the higher the carbon cost and the average inventory cost of the system and the lower the service level. Thus, the overall performance of the inventory system decreases. The influence of lead time fluctuation on the three system performance indicators is different. Stochastic lead time significantly changes the average inventory cost, while it moderately changes the average service level. In addition, the average service level influenced by the lead time is slightly higher than the carbon cost. Furthermore, with the lead time turning to be more volatile, the inventory system performance continues to be weak, but at a slower pace.
2) Optimizing order policy can mitigate the negative impact of stochastic lead times on the inventory system. Our simulation results show that under the same uncertainty of lead time, the impact of different ordering policies on the system's performance is different, especially when the lead time fluctuates greatly. This phenomenon demonstrates that the inventory system is more sensitive to the order policy when lead time varies. In other words, when the lead time is volatile, the supply chain system is sensitive, and the reverse adjustment function of ordering policies is strong. Hence, there exists appropriate order decisions to weaken the negative impact of lead time fluctuation on the supply chain.
3) The fluctuation of lead time affects the regulation effect of order decisions on an inventory system. Ordering decisions have a moderate effect on the inventory system. However, when the lead time is random, the efficiency of adjustment changes. When the lead time fluctuation is mild, the inventory system is relatively stable, and an order decision does not cause the system to oscillate dramatically. The increase in lead time randomness brings a series of management problems, such as inventory control, procurement planning, order delivery, etc. Furthermore, it puts the supply chain system in a volatile state, which lowers the anti-interference ability of the system. Then, the order decision is the decisive factor in determining the system's performance. Moreover, the average inventory cost is mostly affected by order decisions, followed by carbon cost, and then average service level.
Sustainability 2018, 10, x; doi: FOR PEER REVIEW www.mdpi.com/journal/sustainability 4) Impacts on carbon cost are different for average inventory cost and average service level. In terms of performance indicator fluctuation with lead time variance, the value span of carbon cost is weaker than that of inventory cost and service level, but the change trend is more complex. This brings management challenges to balance the system's performance indicators in a low-carbon economy. Furthermore, when ordering decision-making is constrained by carbon policy, the carbon cost change can be taken as the reference point, so that inventory cost and service level can be balanced to maintain a minimum carbon cost.

Conclusions
Based on the APVIOBPCS model, this paper constructs a three-stage supply chain with the distributors as the focal firm in a periodic continuous review system. The traditional APVIOBPCS model primarily focuses on the impact of parameters of adjustment for inventory and work-inprocess inventory discrepancy on supply chains under fixed lead time. However, in business practice, it is difficult to achieve a fixed lead time. Therefore, this paper makes an important contribution by introducing stochastic lead times into the basic APVIOBPCS model.
Our results show that stochastic lead time only affects the values, but not the trends of the inventory system performance indicators. Furthermore, the three performance indicators have different trends under different order policies. The negative impact of stochastic lead times can be dampened by the appropriate order policy. Moreover, the carbon cost curve is concave with regard to the order policy, while the inventory cost curve is continuously increasing, and service level curve is decreasing. Hence, when the managers make order decisions, they can take carbon cost change as reference point, and then balance inventory cost and service level to maintain a minimum carbon cost.
Carbon emission policies in many countries are still in the exploratory stage. In future studies, the optimal carbon emission policy may be identified. Future study can also focus on the interactions between carbon regulation policies and supply chain conditions, so that appropriate policies can be formulated to achieve green sustainable development.