Modeling Uncertain Travel Times in Distribution Logistics

: Uncertainty quantiﬁcation is a critical aspect of distribution logistics, particularly unpredictable travel times caused by trafﬁc congestion and varying transportation conditions. This paper explores the modeling of uncertainty in dealing with travel times in the context of distribution logistics using the collocation method. First, we employ Monte Carlo simulations to assess the efﬁcacy of the collocation method in modeling the variability and uncertainty associated with travel times. Second, we implement the collocation method in Casablanca, Morocco, a city renowned for its extensive distribution logistics operations and its dynamic trafﬁc. Four distinct scenarios are considered: morning peak, inter-peak, evening peak, and off-peak periods. Our study explores two scenarios: one with recurrent congestion, representing typical daily conditions, and the other with unpredictable uncertainties in travel times, accounting for unexpected events that may occur during a distribution day. Our research ﬁndings enhance our understanding of the probabilistic nature of travel times in distribution logistics. This knowledge provides valuable insights applicable to both routine situations with recurrent congestion and non-recurrent congestion. The results’ ﬁndings contribute to a better understanding of the probabilistic nature of travel times in distribution logistics, offering valuable insights for optimizing route planning and scheduling.


Introduction
In recent years, the problem of modeling uncertain travel times in distribution logistics has gained increasing attention as a means to improve logistics operations and enhance overall efficiency.The efficient management of distribution logistics is of crucial importance for businesses operating in distribution logistics, especially in cities with high traffic congestion.The ability to accurately predict and account for uncertain travel times is essential for effective logistics planning and optimization.The unpredictable nature of travel times poses significant challenges for logistics professionals.Traditional logistics models that assume fixed travel times often fail to capture the dynamic and uncertain nature of transportation.This can lead to inefficient resource allocation, missed delivery deadlines, increased costs, and customer dissatisfaction.To address these challenges, researchers and practitioners have turned their attention to the modeling of uncertain travel times as a key factor in distribution logistics.The variability and duration of travel time are influenced by various factors, particularly congestion.Congestion can be categorized as recurrent or non-recurrent [1], each with its distinct characteristics and causes.Recurrent congestion arises from daily peak-hour delays caused by the imbalance between traffic demand and the capacity of the transportation infrastructure.On the other hand, non-recurrent congestion results from atypical events, including incidents [2,3], disabled vehicles, road construction works, accidents, adverse weather conditions, and special events.Both categories of congestion introduce variability and uncertainty in travel times, leaving drivers uncertain about their precise arrival times at their destinations.Non-recurrent congestion can have dual effects on traffic conditions.It can generate new congestion during off-peak periods or increase delays during recurring congestion episodes.However, this type of congestion has often been overlooked in traffic engineering and modeling practices due to its irregular and unpredictable nature [4].The modeling of the uncertain time variability has attracted the interest of several researchers.In a study by Ruimin [5], the examination of travel time variability involved the assessment of various factors, including time of day, day of the week, weather conditions, and traffic accidents.To quantify the influence of these factors on travel time parameters, the author employed multiple linear regression models with two-way interactions.The researchers in [6] investigated the prediction of uncertainty in route travel times within urban road networks, utilizing an extensive dataset from floating taxis; in [7], the authors delved into the importance of precise travel time modeling for road-based public transport and provided an overview of the current advancements in statistically modeling bus travel time distributions.They emphasized the significance of aggregating spatial and temporal data and proposed potential avenues for future research in this domain.In [8], researchers investigated travel time variability within the context of Indian traffic conditions and assessed the effectiveness of travel time distributions when considering various temporal and spatial aggregations.Utilizing AVL data gathered from four transit routes in Mysore City, located in Karnataka, India, this research examined travel time distributions across peak and off-peak periods, as well as within different time intervals.The study employed the Anderson-Darling test to evaluate the goodness-of-fit of the distributions, taking into account segments featuring signalized intersections and a variety of land-use types.The findings underscored the effectiveness of the generalized extreme value (GEV) distribution in precisely characterizing travel time variations within public transit systems.Therefore, the quantification of uncertainty remains a pivotal component when modeling and estimating travel time variability.Precisely capturing and forecasting travel time variations is essential in multiple domains, including transportation planning, traffic management, and logistics.Many research endeavors have focused on tackling uncertainties and developing approaches to improve the precision and effectiveness of uncertainty quantification (UQ).UQ tools have been applied in numerous studies to evaluate and mitigate risks effectively.In [9,10], the authors proposed a model that utilizes uncertainty quantification to address risks associated with the customs supply chain.This model employs the moment matching method and takes into account the seasonality of illicit traffics at five different sites in Morocco.In the study [11], uncertainty quantification methods for differential equations were utilized to predict and simulate the transmission of influenza within the setting of a boarding school.In [12], the researchers developed an uncertainty quantification model to assess and manage the risks related to atmospheric dispersion.They utilized five different techniques to depict and analyze uncertainties, encompassing the probability theory, the interval analysis, the fuzzy approach, the mixed probabilistic-fuzzy approach, and the evidence theory.Stochastic uncertainties, which result from the inherent unpredictability of natural events, can originate from two primary sources: sampling uncertainty or measurement uncertainty [13].Epistemic uncertainties, on the other hand, originate from gaps in information or knowledge about these events and may have various causes, such as data limitations, uncertain parameters, algorithmic ambiguities, or methodological uncertainties [13].In [14], the authors introduced a model to predict and model epidemic risks, employing uncertainty quantification (UQ) techniques.They explored two specific methods, namely, the collocation method and the moment matching method.These methods were applied in the context of studying the epidemic risk associated with SARS-CoV-2 in Morocco, serving as an illustrative example.In the humanitarian logistics context, uncertainty quantification (UQ) holds a significant importance.This is primarily because disasters typically exhibit a substantial degree of uncertainty.For example, for flooding, UQ helps in the making of informed decisions, particularly when it comes to predicting disasters [15].One of the fundamental approaches for uncertainty quantification (UQ) is the Monte Carlo simulation (MCS), which has been widely employed in various fields.The basic principle of the MCS involves repeatedly sampling input parameters from their respective probability distributions and propagating them through the model to obtain a distribution of output responses.By generating a large number of samples, the MCS can capture the full range of possible outcomes and provide statistical measures.The computational cost associated with the Monte Carlo simulation (MCS) can be a challenge due to the need for a large number of samples.This requirement can make the MCS computationally prohibitive [16].The Latin hypercube sampling (LHS) technique offers potential improvements in uncertainty quantification.This sampling approach has been shown to provide better accuracy [17].However, it is important to note that the LHS does have its limitations, as pointed out by [18].In this context, a promising approach known as the polynomial chaos expansion (PCE) emerged.The PCE allows for the representation of random variables as multivariate series of Gaussian variables [19][20][21][22].The representation of random variables follows the concept of a Hilbert basis, with approximation achieved by truncating the representation at a specific level.This approach has demonstrated significant advantages over traditional methods such as the Monte Carlo simulation (MCS) and the Latin hypercube sampling (LHS), offering lower computational costs and an increased efficiency.In this paper, we have employed the collocation method, specifically based on approximation interpolation, to model the variability of travel time.The collocation technique involves selecting specific points, known as collocation points, within the domain of interest and determining the corresponding coefficients to approximate the underlying function.To validate the effectiveness of the collocation method in modeling travel time variability, we conducted a comparison with the Monte Carlo simulation (MCS) approach.The MCS method is widely used as a benchmark for uncertainty quantification and serves as a reference for evaluating the accuracy and efficiency of alternative techniques [23].The utilization of collocation methods for modeling the variability of travel time has gained little significant attention in the literature.In [24], the authors utilized the stochastic response surface method (SRSM) to analyze travel time variability under uncertain factors, including traffic volume and rainfall intensity.Applying SRSM to a section of the Kobe Route, Japan, the results showed an improved accuracy in modeling travel time distribution compared to the regression analysis; the collocation method was used for estimating the travel time.In this work, we employed the collocation method to effectively model the travel time variability during the four peak hours, considering both recurrent congestion and non-recurrent congestion scenarios (Figure 1).This approach allowed us to delve deeper into understanding how travel time fluctuations evolve over the course of the day, shedding light on the intricate dynamics of traffic patterns during peak hours.
The collocation method has proven to be a valuable tool for uncertainty quantification in travel time estimation.It enables us to address the challenges posed by the dynamic and complex traffic conditions such as those in Casablanca.By incorporating the variability of traffic volume, incidents like road accidents, and other external factors, we gained insight into how uncertainties impact travel time reliability.The insights derived from our analysis will aid logistics practitioners in making informed decisions, optimizing route planning, and mitigating the potential impacts of travel time uncertainties on delivery schedules.Figure 2 explains the framework methodology of this study: the aim is to model uncertain travel times in distribution logistics using the uncertainty quantification approach with the collocation method.To ensure the reliability and effectiveness of the collocation method, we conducted a thorough validation process using Monte Carlo simulations (MCS).Subsequently, we conducted modeling using the cumulative distribution function (CDF) in both the recurrent and non-recurrent congestion scenarios.Then, a probability density function (PDF) was derived from the CDF using Dirac's delta approximation.To comprehensively assess the variation of travel time uncertainties in distribution logistics, we calculated probabilities across all four peak time scenarios: 'morning peak', 'interpeak', 'evening peak', and 'off-peak'.This encompassing approach helped us gain a holistic understanding of how different factors, including traffic congestion, influence the variability of travel times.The collocation method has proven to be a valuable tool for uncertainty quantification in travel time estimation.It enables us to address the challenges posed by the dynamic and complex traffic conditions such as those in Casablanca.By incorporating the variability of traffic volume, incidents like road accidents, and other external factors, we gained insight into how uncertainties impact travel time reliability.The insights derived from our analysis will aid logistics practitioners in making informed decisions, optimizing route planning, and mitigating the potential impacts of travel time uncertainties on delivery schedules.Figure 2 explains the framework methodology of this study: the aim is to model uncertain travel times in distribution logistics using the uncertainty quantification approach with the collocation method.To ensure the reliability and effectiveness of the collocation method, we conducted a thorough validation process using Monte Carlo simulations (MCS).Subsequently, we conducted modeling using the cumulative distribution function (CDF) in both the recurrent and non-recurrent congestion scenarios.Then, a probability density function (PDF) was derived from the CDF using Dirac's delta approximation.To comprehensively assess the variation of travel time uncertainties in distribution logistics, we calculated probabilities across all four peak time scenarios: 'morning peak', 'inter-peak', 'evening peak', and 'off-peak'.This encompassing approach helped us gain a holistic understanding of how different factors, including traffic congestion, influence the variability of travel times.This paper is organized as follows: Section 3 describes the uncertainty quantification model and contains the collocation model validated by the Monte Carlo simulation (MCS) using an example; Section 4 presents the computational results and discussion; and Section 5 presents the conclusion and discussion.
• Establishing the context involves delving into distribution logistics, specifically the modeling of uncertain travel times.

Defining the context: uncertain travel times in distribution logistics
• A specific example is used to validate our approach, involving an ordinary differential equation (ODE) system.

Validate and compare the collocation method with the Monte Carlo simulation (MCS)
• The collocation method is applied to address the critical aspect of uncertainty in travel time to a specific delivery route within the city of Casablanca.

Uncertainty quantification of travel time using the collocation method
• The objective is to create models for the cumulative distribution function (CDF) of travel time under uncertain conditions.These models will capture variations in travel time within both recurrent and non-recurrent congestion scenarios at different peak times.

Modeling the cumulative distribution dunction (CDF) of uncertain travel time in both recurrent an non-recurrent congestion in various peak hours
• The objective is to provide a detailed view of how the values of random travel time are distributed within the predefined range of travel time data.

Modeling the probability density function (PDF) ) of uncertain travel time in both recurrent an nonrecurrent congestion at various peak times by derivating it from the CDF using Dirac's delta approximation
• The objective is to calculate probabilities associated with different travel time categories, encompassing all four peak times scenarios: "morning peak," "inter-peak," "evening peak," and "off-peak," considering both recurrent and nonrecurrent congestion.This paper is organized as follows: Section 3 describes the uncertainty quantification model and contains the collocation model validated by the Monte Carlo simulation (MCS) using an example; Section 4 presents the computational results and discussion; and Section 5 presents the conclusion and discussion.

Literature Review
Quantifying uncertainty in uncertain travel times is crucial and has wide-ranging applications, especially in logistics and route planning.The variability in travel time within distribution logistics is a consequence of the stochastic nature of numerous operational factors.This variability poses challenges for distribution and logistics firms, introducing uncertainty and associated costs.In [25], the authors conducted a comprehensive systematic review of the existing research on travel time reliability, with a specific focus on the assessment of the value associated with travel time reliability.Additionally, they undertook a meta-analysis to uncover the underlying factors contributing to variations in reliability estimates.Travel time uncertainty varies depending on congestion levels [6].Congestion travel time levels have been considered by researchers.In [26], the authors examined how changes in congestion levels influenced the choice of departure time and they conducted an analysis of costs associated with uncertain travel times.In the study [27], the authors emphasized the relationship between travel time variability and congestion levels.To assess the value of travel time reliability, they developed a comprehensive modeling framework that incorporated trip scheduling, endogenous traffic congestion, travel time uncertainty, and pricing strategies.Hence, to quantify the value of travel time reliability, they integrated trip scheduling, endogenous traffic congestion, travel time uncertainty, and pricing strategies into one modeling framework.The common thread among the models developed in these research studies is that normal traffic flow conditions become unstable when uncertainties arise, especially during peak periods or when other incidents begin to deteriorate the typical traffic conditions.Hence, travel time uncertainty is frequently a result of traffic congestion, and, when congestion occurs, it amplifies the frustration experienced by road users.For the studies [28][29][30], it was important to model and understand the implications of travel time uncertainty on transportation systems.Additionally, studies of transportation networks have shown that travel time uncertainty is an important factor that influences traffic networks [30][31][32].For example, in [31], it was observed that travelers place a greater importance on decreasing travel time variability than on reducing the total travel time for a specific journey.In [32], it was noted that deterministic traffic network models are inadequate for assessing traffic network performance.These models tend to magnify minor discrepancies in travel times while failing to adequately account for the substantial impacts of travel time uncertainty.In [30], it was determined that the loss in utility resulting from uncertainty is comparable in magnitude to the overall travel costs.The majority of the previous research concentrated on either the causes or the effects of travel time uncertainty.Considering the impacts of travel time uncertainty within traffic networks necessitates an explicit approach for quantifying and modeling this uncertainty.Hence, in comparison to previous studies in the literature, our objective is to model travel time uncertainty within the domain of distribution logistics.We aim to demonstrate the variability of travel time at different stages of congestion and during various peak times.To accomplish this, we employ an effective method for uncertainty quantification, namely, the collocation method.In addition, modeling uncertain travel times is crucial for real-world applications, especially in the context of routing optimization, where stochastic scenarios are commonly utilized to account for uncertainties between nodes [33].

Uncertainty Quantification Model
This section develops a model based on uncertainty quantification (UQ) and presents the collocation method, an interpolation-based approximation, which is then validated using the Monte Carlo simulation (MCS).

Model Formulation
A model based on uncertainty quantification (UQ) employs a discrete method that involves approximations by analyzing expressions that depend on the joint distribution of the pair X, Q.Here, Q represents the uncertainties within the system, modeled as a random vector, while X represents the system's response.We assume that X is a function dependent of Q, as X = X(Q).Statistical information can be obtained from a sample of the pair (X, Q).As X is a function of Q, denoted as X = X(Q), we can adopt a Hilbertian approach [16,34]: We assume that = X(Q) ∈ H, where H is a separable Hilbert space; we choose an adequate Hilbert basis {ϕ i } i N * : We consider X and Q as two random variables on the probability space (Ψ, P).Let H be a subspace of functions of X. Find an element PX from H such that PX provides the best approximation of X within W. It is anticipated that PX → X for n → +∞ .The primary objective of a representation is to employ PX instead of X when assessing the probability of significant events.This ensures that PX serves as a practical approximation of X, such that the probabilities obtained with PX are close to the actual values when compared to those obtained with X.In general, achieving this objective entirely is not always feasible, so we seek an approximation where PX can be used as a replacement for X in practical scenarios.This implies that the difference denoted as X − PX, should be small enough, aligning with our specific goals and requirements.If {ϕ i } i N * is not orthonormal, it is necessary to obtain a sequence of orthogonal projections.
where x ∈ R k denotes the unknown coefficients of the expression of PX which describes the evolution of the response of the system and {ϕ i } i N * is a Hilbert basis.

Numerical Method Based on Collocation
For using a sample, {(Q i , X i ) : i = 1, . . . . . .ns } consists of determining In this case, X = (x 1 , x 2 , . . .x k ) t is the solution of a linear system:

Probability Density Function Approximation
The probability density function (PDF) is a mathematical function used to depict the likelihood of a random variable assuming a specific value.It is derived from the cumulative distribution function (CDF) and is obtained through a technique called "Particle derivatives", which involves employing Dirac's delta approximation [35].This approach involves reviewing the derivation using Dirac's delta approximation to identify the PDF.The Dirac delta function δ(y 0 ) can be naturally approximated by a sequence of probabilities that converge to a distribution concentrated at y 0 .
Let F(y) be a function defined by the following: We have We may explore the possibility of employing the approximation technique, which is further elaborated in [35]: so that In the realm of particle derivatives, we use p J to denote individual particles, while N J signifies the associated volume, and f represents the kernel employed in the approximation process.Drawing an analogy with conventional physical parameters, we establish the equation N J = m j /ρ j , where m j represents the mass of the particle and ρ j refers to its volumetric mass density, commonly referred to as "density".Consequently, Equation ( 9) can be derived as follows: To compute the partial derivatives, we follow a similar approach, introducing the approximations F J = Fp J which yields Then, It is important to highlight that a range of kernels can be found in the literature; the most commonly used are defined as follows: where µ(n) is a normalizing factor such that the integral on R n f (x, y, n)dx = 1.In our particular study, we will adopt the quintic kernel, such that µ(n) = 1 120 ; 1 478π ; 1 120π (dimension 1, 2, or 3, respectively)

Model Validation by Monte Carlo Simulation (MCS)
In this study, we employed interpolation-based approximations, specifically the collocation method, to estimate travel time variability.To validate the effectiveness of this approach, we conducted a Monte Carlo simulation (MCS) on a specific example.The Monte Carlo simulation (MCS) is a computational technique used to estimate and analyze complex systems or processes by performing a large number of random simulations.It involves generating random samples from input probability distributions to model uncertain parameters, and then evaluating the system's behavior multiple times to observe the outcomes.In our example, we generated 1000 samples using the Monte Carlo simulation (MCS) and 30 samples using the collocation method.The aim was to compare the results of the two methods.Let us consider a system of ordinary differential equations (ODE) represented as follows: where y(t) = ae bt ; such that y = by y(0) = a over the interval [0, T]; a and b represent two variables generated using triangular distribu- tion to obtain 1000 samples for the Monte Carlo simulation (MCS) and 30 samples for the collocation method; the initial values used in the triangular distribution are represented in Tables 1 and 2. The parameters presented in Tables 1 and 2 represent simulation values of the ordinary differential equation (ODE) system that we used to confirm the reliability of our model.These simulations were conducted while considering various conditions and model parameters, and the resulting values align with the outcomes we obtained in our numerical experiments.In this paper, [0, T] is the time interval over which the solution was computed.
The ODE system was defined on Matlab as follows Algorithms 1: % s is the integrated value of y over the interval [0, T] sol.int = s; % sol is the solution of the ODE system obtained using the ode45 method As a result, Figures 3 and 4 represent two different approximations of the system under consideration using the Monte Carlo simulation (MCS) and collocation methods, respectively.The Monte Carlo simulation allowed us to generate 1000 samples for the variables "a" and "b" using a triangular distribution.On the other hand, the collocation method enabled us to generate 30 samples for the same variables with the same distribution.According to the figures, both the MCS's and the collocation's approximations show very close results to each other.This indicated that even with a small sample size of 30, as used in the collocation method, we achieved high-performance results comparable to those obtained with the MCS method.Additionally, in terms of execution time, the collocation method outperformed the MCS, making it a more efficient choice.Hence, our results revealed that the collocation method provided a robust and precise representation of travel time variability, closely aligning with the MCS approach.Notably, the collocation method achieved comparable results with a smaller sample size, making it a more efficient choice for uncertainty quantification in distribution logistics.The collocation method's concurrence with the Monte Carlo simulation approach reaffirmed its accuracy and dependability.Consequently, for our specific case study aimed at modeling travel time variability in distribution logistics, we opted for the collocation method as it combined accuracy, efficiency, and reliability in order to effectively capture and quantify travel time uncertainties.The decision to utilize the collocation method as our primary modeling technique was grounded in several compelling factors.First and foremost, the collocation method had demonstrated its effectiveness in addressing uncertain parameters, a critical aspect when dealing with the variability of travel times.This effectiveness lies in its ability to produce precise estimates of travel time variability, even when working with relatively smaller sample sizes.One of the significant advantages of the collocation method was its efficiency in capturing the essence of travel time uncertainty.While the Monte Carlo simulation (MCS) is a valuable technique, it often requires a substantial number of random samples to achieve accurate results.In contrast, the collocation method excels in providing reliable estimates without the need for an extensive sample pool.This efficiency translates  According to the figures, both the MCS's and the collocation's approximations show very close results to each other.This indicated that even with a small sample size of 30, as used in the collocation method, we achieved high-performance results comparable to those obtained with the MCS method.Additionally, in terms of execution time, the collocation method outperformed the MCS, making it a more efficient choice.Hence, our results revealed that the collocation method provided a robust and precise representation of travel time variability, closely aligning with the MCS approach.Notably, the collocation method achieved comparable results with a smaller sample size, making it a more efficient choice for uncertainty quantification in distribution logistics.The collocation method's concurrence with the Monte Carlo simulation approach reaffirmed its accuracy and dependability.Consequently, for our specific case study aimed at modeling travel time variability in distribution logistics, we opted for the collocation method as it combined accuracy, efficiency, and reliability in order to effectively capture and quantify travel time uncertainties.The decision to utilize the collocation method as our primary modeling technique was grounded in several compelling factors.First and foremost, the collocation method had demonstrated its effectiveness in addressing uncertain parameters, a critical aspect when dealing with the variability of travel times.This effectiveness lies in its ability to produce precise estimates of travel time variability, even when working with relatively smaller sample sizes.One of the significant advantages of the collocation method was its efficiency in capturing the essence of travel time uncertainty.While the Monte Carlo simulation (MCS) is a valuable technique, it often requires a substantial number of random samples to achieve accurate results.In contrast, the collocation method excels in providing reliable estimates without the need for an extensive sample pool.This efficiency translates According to the figures, both the MCS's and the collocation's approximations show very close results to each other.This indicated that even with a small sample size of 30, as used in the collocation method, we achieved high-performance results comparable to those obtained with the MCS method.Additionally, in terms of execution time, the collocation method outperformed the MCS, making it a more efficient choice.Hence, our results revealed that the collocation method provided a robust and precise representation of travel time variability, closely aligning with the MCS approach.Notably, the collocation method achieved comparable results with a smaller sample size, making it a more efficient choice for uncertainty quantification in distribution logistics.The collocation method's concurrence with the Monte Carlo simulation approach reaffirmed its accuracy and dependability.Consequently, for our specific case study aimed at modeling travel time variability in distribution logistics, we opted for the collocation method as it combined accuracy, efficiency, and reliability in order to effectively capture and quantify travel time uncertainties.The decision to utilize the collocation method as our primary modeling technique was grounded in several compelling factors.First and foremost, the collocation method had demonstrated its effectiveness in addressing uncertain parameters, a critical aspect when dealing with the variability of travel times.This effectiveness lies in its ability to produce precise estimates of travel time variability, even when working with relatively smaller sample sizes.One of the significant advantages of the collocation method was its efficiency in capturing the essence of travel time uncertainty.While the Monte Carlo simulation (MCS) is a valuable technique, it often requires a substantial number of random samples to achieve accurate results.In contrast, the collocation method excels in providing reliable estimates without the need for an extensive sample pool.This efficiency translates into resource savings, making it an attractive choice for practical applications.Moreover, in the context of distribution logistics, where decisions must often be made in real-time, computational efficiency is paramount.The collocation method's computational speed outpaces that of the MCS, allowing for quicker analysis and decision-making.This attribute aligned perfectly with the demands of our study, where timely logistics decisions could significantly impact the efficiency and effectiveness of supply chain operations.

Numerical Analysis for Uncertain Travel Time
This section presents the data analysis related to the travel and computational results of the variability of travel time at different times of the journey, including morning peak, interpeak, evening peak, and off-peak periods.The aim is to explore the details of how travel times behave during these specific phases.This helps us to better understand the patterns and changes that affect distribution logistics operations.The selection of the collocation method stemmed from its proven effectiveness in handling uncertain parameters.We recognized its capability to provide accurate estimates of travel time variability with relatively smaller sample sizes, making it an efficient choice.Additionally, its computational efficiency compared to the MCS is essential for real-time logistics decision-making, a key consideration for our study.

Data Analysis
In today's fast-paced world, the volume of traffic on roadways has become a major concern.As cities continue to expand, the surge in vehicular movement has resulted in significant congestion during peak hours and other periods of high demand.In light of this pressing issue, our study focuses on a specific case in the city of Casablanca, which is renowned for its heavy traffic flow.Casablanca, often hailed as the economic capital of Morocco, is renowned for its significant volume of goods deliveries.Hence, Casablanca, being a bustling economic center, experiences a complex interplay of factors contributing to its traffic congestion.The bustling city serves as a major hub for commerce and trade, attracting numerous businesses engaged in transportation and logistics activities.As a result, the city's intricate transportation network experiences a constant flow of goods deliveries, making it a critical focal point for distribution logistics.The increased traffic not only affects travel time but also the distribution logistics' performance.Currently, the transportation sector in Morocco thrives, with approximately 76,000 companies engaged in road transport, boasting a fleet of around 74,000 vehicles.In Table 3, we show a comprehensive depiction of the number of goods delivery companies operating in various regions across Morocco; the table highlights the distribution of these enterprises across the different regions.This information offers a clear visualization of the regional engagement of goods delivery companies, contributing to a better understanding of the logistics landscape across the country.According to the statistics, it is evident that Casablanca is home to a significant number of goods delivery companies, resulting in a large fleet of delivery vehicles.This makes Casablanca the city with the highest traffic in terms of distribution logistics as well as urban mobility.As a major economic center and the largest city in Morocco, Casablanca has become a focal point for distribution and transportation activities.The concentration of goods delivery companies in the city has led to a high density of delivery vehicles navigating through its busy streets.Our study focused on a delivery route within the city of Casablanca, from the depot (Coopérative Pharmaceutique Marocaine) to the customer (Pharmacie Joyaux de Bouskoura), as shown in Figure 5.Our work delved into the quantification of uncertainties in travel times to address these challenges by employing the collocation method.The aim was to model the variability of travel times in the delivery route.
The collected data comprised travel times between these two points under four scenarios: "Morning peak hour", "Inter-peak", "Evening peak", and "Off-peak".The study considered two traffic states: the first being normal traffic (recurrent congestion) and the second being non-recurrent congestion caused by incidents such as accidents or unexpected disruptions.The aim was to study these variations to better understand the impact of different traffic conditions on travel times and make informed decisions to enhance the efficiency and reliability of distribution logistics operations.The data presented in Tables 4 and 5 represent various values of travel time throughout the day, including during peak hours.These data were collected using the Google Maps routing service, which calculates route-specific travel times based on real-time traffic conditions.goods delivery companies in the city has led to a high density of delivery vehicles navi-gating through its busy streets.Our study focused on a delivery route within the city of Casablanca, from the depot (Coopérative Pharmaceutique Marocaine) to the customer (Pharmacie Joyaux de Bouskoura), as shown in Figure 5.Our work delved into the quantification of uncertainties in travel times to address these challenges by employing the collocation method.The aim was to model the variability of travel times in the delivery route.The collected data comprised travel times between these two points under four scenarios: "Morning peak hour", "Inter-peak", "Evening peak", and "Off-peak".The study considered two traffic states: the first being normal traffic (recurrent congestion) and the second being non-recurrent congestion caused by incidents such as accidents or unexpected disruptions.The aim was to study these variations to better understand the impact of different traffic conditions on travel times and make informed decisions to enhance the efficiency and reliability of distribution logistics operations.The data presented in Tables 4 and 5 represent various values of travel time throughout the day, including during peak hours.These data were collected using the Google Maps routing service, which calculates route-specific travel times based on real-time traffic conditions.For our study, we utilized the extensive and up-to-date travel time data available on Google Maps.This allowed us to capture the dynamic nature of traffic patterns in Casablanca.Considering real-time traffic conditions, we were able to analyze travel times between the depot and the customer in four different scenarios, namely, "Morning peak hour", "Inter-peak", "Evening peak", and "Off-peak".Additionally, we accounted for two distinct states of traffic: "recurrent congestion", representing normal traffic conditions, and "non-recurrent congestion", which corresponded to unpredictable traffic fluctuations.

Computational Results and Discussion
The numerical computations were conducted using Matlab R2018a on a computer equipped with the following specifications: Intel(R) Core (TM) i5-6300U CPU @ 2.40 GHz 2.50 GHz, and 8.00 GB of RAM, running on the Windows operating system.Additionally, to model the variability of travel time, we utilized a polynomial degree of 11.
The study produced travel time variability results for both the recurrent congestion and non-recurrent congestion scenarios.We modeled this variability by employing probability density functions (PDF) and cumulative distribution functions (CDF) for each of the four scenarios: morning peak, inter-peak, evening peak, and off-peak.

•
Scenario 1: Recurrent congestion In this scenario, our focus is on understanding the variability of travel time during recurrent congestion, specifically in four distinct peak time scenarios: morning peak, interpeak, evening peak, and off-peak.Recurrent congestion refers to the regular and predictable traffic jams that occur during specific times of the day when there is a significant increase in traffic volume, such as rush hours.During peak hours, the travel time variability can significantly impact the efficiency of distribution logistics operations, leading to delays, increased costs, and potential customer dissatisfaction.We present in Figures 6 and 7    The cumulative distribution function (CDF) and the probability density function (PDF) of travel time in the case of recurrent congestion are represented in Figures 6 and 7, respectively.The figures show the four scenarios: during the off-peak and inter-peak periods traffic congestion is relatively low compared to the "morning peak" and "evening peak" scenarios, when the traffic volume increases significantly.As a result, it is crucial for logistics and transportation planners to implement effective strategies to reduce delays and ensure on-time deliveries during peak traffic hours.This includes optimizing routes, planning for contingencies, and monitoring real-time traffic conditions.During the four peak time scenarios, we performed probability calculations based on the classification of travel times into three categories: "low" for travel times below 28 min, "medium" for travel times between 28 and 39 min, and "high" for travel times exceeding 42 min.The results of these probability calculations are presented in Table 6.It provides a comprehensive overview of the probabilities associated with each travel time category for all four peak time scenarios: "morning peak", "inter-peak", "evening peak", and "off-peak".The table helps us understand how likely travel times are to fall into each category at different times of the day.The cumulative distribution function (CDF) and the probability density function (PDF) of travel time in the case of recurrent congestion are represented in Figures 6 and 7, respectively.The figures show the four scenarios: during the off-peak and inter-peak periods traffic congestion is relatively low compared to the "morning peak" and "evening peak" scenarios, when the traffic volume increases significantly.As a result, it is crucial for logistics and transportation planners to implement effective strategies to reduce delays and ensure on-time deliveries during peak traffic hours.This includes optimizing routes, planning for contingencies, and monitoring real-time traffic conditions.During the four peak time scenarios, we performed probability calculations based on the classification of travel times into three categories: "low" for travel times below 28 min, "medium" for travel times between 28 and 39 min, and "high" for travel times exceeding 42 min.The results of these probability calculations are presented in Table 6.It provides a  Based on the results, it is observed that the probabilities for "high" traffic levels predominantly occur during the morning peak and evening peak times.On the other hand, during the off-peak and inter-peak periods, the travel times are generally shorter.

•
Scenario 2: Non-recurrent congestion During the non-recurrent congestion, the travel time variability becomes more unpredictable.Incidents such as road accidents can occur at any time of the day, disrupting the traffic flow and leading to unexpected delays.These unpredictable events can significantly impact travel times, causing fluctuations and introducing uncertainty in distribution logistics operations.We will assume that accidents can occur during each of the four peak times, potentially disrupting the traffic flow.These accidents can lead to unexpected road closures, detours, and increased congestion, resulting in longer travel times for distribution logistics vehicles.As a result, the variability in travel times during peak hours becomes even more uncertain and challenging to predict.Figures 8 and 9 present the variability of travel time illustrated in the cumulative distribution function (CDF) and the probability density function (PDF).These graphical representations show the distribution of travel times during different peak hours, taking into account the uncertainties caused by accidents and non-recurrent congestion.
road closures, detours, and increased congestion, resulting in longer travel times for distribution logistics vehicles.As a result, the variability in travel times during peak hours becomes even more uncertain and challenging to predict.Figures 8 and 9 present the variability of travel time illustrated in the cumulative distribution function (CDF) and the probability density function (PDF).These graphical representations show the distribution of travel times during different peak hours, taking into account the uncertainties caused by accidents and non-recurrent congestion.According to Figures 8 and 9, we observe multiple peaks in the PDF, indicating the presence of several disruptions, uncertainties, and random events that affect the travel time.This multimodal distribution suggests that the distribution of travel time is not governed by a single dominant factor but is influenced by various sources of variability.One of the significant factors contributing to these peaks could be road accidents, which can cause unexpected delays and fluctuations in travel time.Road accidents can occur at any time of day and can have a significant impact on the travel time during non-recurrent congestion.When accidents happen, they can lead to unexpected road closures, diversions, or reduced speed limits, causing traffic disruptions and delays.The variability in travel time during non-recurrent congestion becomes more pronounced as these unforeseen events can create fluctuations in the flow of vehicles and introduce additional uncertainties in the travel time estimation.Then, the probabilities obtained for these scenarios were calculated and results are presented below.
The results presented in Table 7 reveal some interesting insights regarding the impact of non-recurrent congestion, such as road accidents, in contrast to recurrent congestion.In non-recurrent congestion scenarios, there is a high probability of travel time disruptions during inter-peak periods in addition to peak hours (morning peak and evening According to Figures 8 and 9, we observe multiple peaks in the PDF, indicating the presence of several disruptions, uncertainties, and random events that affect the travel time.This multimodal distribution suggests that the distribution of travel time is not governed by a single dominant factor but is influenced by various sources of variability.One of the significant factors contributing to these peaks could be road accidents, which can cause unexpected delays and fluctuations in travel time.Road accidents can occur at any time of day and can have a significant impact on the travel time during non-recurrent congestion.When accidents happen, they can lead to unexpected road closures, diversions, or reduced speed limits, causing traffic disruptions and delays.The variability in travel time during non-recurrent congestion becomes more pronounced as these unforeseen events can create fluctuations in the flow of vehicles and introduce additional uncertainties in the travel time estimation.Then, the probabilities obtained for these scenarios were calculated and results are presented below. The results presented in Table 7 reveal some interesting insights regarding the impact of non-recurrent congestion, such as road accidents, in contrast to recurrent congestion.In non-recurrent congestion scenarios, there is a high probability of travel time disruptions during inter-peak periods in addition to peak hours (morning peak and evening peak).This finding highlights the significance of unexpected events, like accidents, in causing traffic disturbances during times that are typically less congested.On the other hand, in recurrent congestion scenarios, we observe high probabilities of travel time disruptions during peak hours, which aligns with the expected patterns of heavy traffic.These results emphasize the importance of considering both types of congestion in distribution logistics planning to effectively manage the variability of travel time and ensure timely deliveries.Travel time variability introduces uncertainties that can significantly impact the performance of distribution logistics.Delays, particularly those caused by road accidents, directly influence the cost of deliveries, resulting in increased expenses.Accurate quantification of travel time variability is crucial for optimizing logistics operations.It helps logistics practitioners allocate resources efficiently, plan for contingencies, and ensure timely deliveries.Managing uncertainties effectively is essential for maintaining high customer satisfaction and building trust in the logistics services.Hence, by acknowledging the influence of non-recurrent congestion events, we can implement better strategies to mitigate the impact of unexpected disruptions and enhance overall logistics efficiency.

Conclusions and Discussion
This study focuses on quantifying uncertainties in modeling travel time variations in two distinct traffic scenarios-recurrent congestion and non-recurrent congestion-observed during four peak delivery periods in Casablanca's distribution logistics.The research highlights the critical importance of addressing uncertainties in travel times within Casablanca's dynamic traffic environment, necessitating reliable methods like collocation for effective management.The probability density functions (PDFs) and the cumulative density functions (CDFs) were utilized to analyze travel time distributions, offering valuable insights into logistics uncertainties.In cases of recurrent congestion during morning peak, interpeak, evening peak, and off-peak times, fluctuations in travel time were observed.The PDFs showed some variability, but, overall, they indicated a relatively stable traffic pattern.Contrastingly, non-recurrent congestion scenarios, often triggered by incidents like accidents, exhibited more complex and multi-modal PDF distributions due to unpredictable traffic disruptions.The presence of multiple peaks in the PDFs during non-recurrent congestion scenarios underscored the diverse and distinct travel time variations caused by such incidents.Accidents and other unforeseen events have a profound impact on distribution logistics, causing delays and uncertainties in delivery times.
The study emphasizes the importance of considering different traffic scenarios when modeling travel time variability for distribution logistics.By quantifying uncertainties and understanding the impact of both recurrent and non-recurrent congestion, logistics managers and planners can develop robust strategies to optimize routes, improve efficiency, and minimize disruptions' impact on delivery schedules.As for the limitations of this study, while it addresses the effects of recurrent and non-recurrent congestion on travel time variability, other factors may also play a role, such as weather conditions, unforeseeable incidents, and driver behaviors.These elements can introduce additional uncertainties that have not been fully explored within the scope of this research.Hence, future improvements

19 Figure 2 .
Figure 2. The methodology of the study framework.

Figure 2 .
Figure 2. The methodology of the study framework.

Figure 4 .
Figure 4. PDF using the collocation method.

Figure 4 .
Figure 4. PDF using the collocation method.

Figure 4 .
Figure 4. PDF using the collocation method.

Figure 5 .
Figure 5.The delivery process arc.

Figure 5 .
Figure 5.The delivery process arc.

19 Figure 6 .
Figure 6.CDF of travel time in recurrent congestion.

Figure 7 .
Figure 7. PDF of travel time in recurrent congestion.

Figure 6 .
Figure 6.CDF of travel time in recurrent congestion.

Figure 6 .
Figure 6.CDF of travel time in recurrent congestion.

Figure 7 .
Figure 7. PDF of travel time in recurrent congestion.

Figure 7 .
Figure 7. PDF of travel time in recurrent congestion.

Figure 8 .
Figure 8. CDF of travel time in non-recurrent congestion.Figure 8. CDF of travel time in non-recurrent congestion.

Figure 8 . 19 Figure 9 .
Figure 8. CDF of travel time in non-recurrent congestion.Figure 8. CDF of travel time in non-recurrent congestion.Appl.Sci.2023, 13, x FOR PEER REVIEW 16 of 19

Figure 9 .
Figure 9. PDF of travel time in non-recurrent congestion.

Table 1 .
Values of variable "a" used to generate samples with triangular distribution.

Table 2 .
Values of variable "b" used to generate samples with triangular distribution.

Table 3 .
Number of goods transportation companies, statistics from December 2022.

Table 4 .
Travel time during recurrent congestion.

Table 5 .
Travel time during non-recurrent congestion.

Table 6 .
Travel time probability in recurrent congestion.

Table 7 .
Travel time probability in non-recurrent congestion.