Toward a Comprehensive Pavement Reliability Analysis Approach

Reliability has been incorporated in pavement design tools to account for input variability influence on predicted performance. As they are not based on a probabilistic method of uncertainty propagation, the reliability analysis methodologies that are currently implemented in pavement performance tools lack rigor and robustness. This paper investigates the potential of three reliability analysis methodologies for pavement application: the Pavement ME reliability analysis methodology, Monte Carlo simulation (MCS), and the first-order reliability method (FORM). The MCS and FORM involve a response surface method for the generation of a second-order surrogate model. The investigation was performed using inputs and performance data from accelerated pavement testing structures. Inputs that were identified as significant were characterized as random variables and their associated variability was established using measured structural and material properties. Pavement performance with respect to rutting was predicted using the ERAPave performance prediction tool, while MCS was used to generate the actual variability of the distress. The reliability analysis results have shown that a comprehensive reliability analysis methodology is required that effectively captures input variabilities and the error associated with surrogate models.

The pavement design process is probabilistic, as it involves uncertainties that originate from different sources (1). The impact of these uncertainties on predicted performance should be addressed and mitigated as this would allow maintenance and rehabilitation interventions that are timely and cost effective. This will consequently lower the total life cycle cost estimate of highway projects (2). Many pavement design procedures have incorporated reliability into the design process to account for input variability influence on predicted performance. Reliability is computed in these design procedures by utilizing different reliability analysis methodologies and using variabilities that are characterized and quantified on the basis of various assumptions, making it challenging to perform a direct comparison across pavement design procedures (3)(4)(5). A reliability analysis methodology that is comprehensive in nature and that characterizes input variability through statistically sound methods is required for enhancing the implementation of reliability for pavement applications.
The 1993 American Association of State Highway and Transportation Officials (AASHTO) pavement design guide adopted the reliability design factor (Fr) as a positive spacing parameter between the allowable and expected traffic repetitions (3). This way of lumping together all the involved uncertainties into a single factor might not deliver sections of uniform reliability, as design inputs have different levels of impact on estimated performance. To address and overcome this problem, different reliability analysis methods were suggested. These methodologies recognize and consider predicted distress variability influence on reliability. The AASHTO Mechanistic-empirical Pavement Design Guide (MEPDG) employs a simplified reliability analysis methodology that assumes each predicted distress to be a key output of interest, while other methodologies utilize Monte Carlo simulation (MCS) or an analytical approximation method such as the first-order reliability method (FORM) (4,(6)(7)(8). Reliability-based design procedures in the format of the load and resistance factor design (LRFD) format were also suggested (9). The AASHTO MEPDG reliability analysis approach does not utilize a probabilistic method of uncertainty propagation to estimate input variability influence on predicted distress. Instead, it computes reliability by assuming each pavement distress to be an independent random variable. The variability of each individual pavement failure mode is approximated by using a normal distribution function. Furthermore, the associated variability of each distress is characterized with a standard error (SE) value, which is estimated using input variability, uncertainty caused by the construction process, and model error (4). It is primarily the unavailability of a comprehensive reliability analysis methodology that has led many pavement design tools to adopt simplified approaches with respect to the variability of predicted distress (8,10,11).
Accelerated pavement testing (APT) involves a continuous measurement of material and structural properties along with observed performance, providing an excellent opportunity to study the input variability influence on predicted distress (12,13). APT is also very convenient for controlling the level of variability associated with inputs, delivering sections with homogeneous structural and mixture properties. As full-scale pavement structures are constructed using the same kind of procedure as field pavements, APT sections are expected to reflect actual field conditions (14). As testing is performed in a controlled environment, APT allows measured pavement response and performance to reflect accurately the impact of both inputs and their associated variability. Furthermore, the short period of testing allows material properties not to exhibit variations that are temporal in nature.

Objective
The main objective of this research is to investigate and establish the impact of a reliability analysis methodology on estimated reliability level. The reliability analysis approach in Pavement ME, MCS, and the FORM are the three investigated methodologies. In addition, the study characterizes the variability associated with inputs, such as layer modulus and layer thickness. For this purpose, two APT sections that were tested using heavy vehicle simulation (HVS) and that are representative of actual in-service conditions were used. The variability of the inputs that were identified as significant through a sensitivity analysis was assessed measured falling weight deflectometer (FWD) and elevation level measurements. ERAPave performance prediction (PP) was used to determine pavement performance with respect to rutting.

Pavement Reliability
The variability nature of pavement performance has motivated many researchers and practitioners to adopt reliability for pavement performance evaluation (3)(4)(5)15). Reliability estimates the probability that a pavement structure will perform as planned during its design life. Reliability estimate of a given distress and design condition is highly dependent on the manner in which the performance function is formulated, and the way failure is defined. Furthermore, reliability analysis requires the statistical characterization of the performance function, which is the product of the independent random variables. Therefore, the design inputs with significant impact on predicted performance need to be identified and statistically characterized. In cases where the full probability density function (PDF) of the random variables (f x (x)) is available, the probability of failure (P f ) and its complement, reliability (R), can be computed using Equations 1 and 2, respectively. The probability of failure, as can be seen in Equation 1, integrates the joint PDF of the independent random variables over the failure region g(x) \ 0: An exact solution for the integration problem in Equation 1 can be obtained only for special cases and finding a solution through numerical integration will be impractical once the number of variables exceed two or three. Most engineering problems are multi-dimensional and involve high nonlinearity, making it very difficult to obtain an exact solution. For these reasons, pavement reliability is performed using simplified approaches that involve either MCS or the FORM (4,9,10).
Reliability was estimated in this study using MCS, the FORM, and the reliability analysis approach in Pavement ME. For the FORM analysis, a twocomponent reliability analysis methodology was used (5,16). The first component generates using a central composite design (CCD) response surface method (RSM), a simplified surrogate model to represent the performance function. The surrogate model is a second-degree polynomial function and considers as inputs all the significant random variables (i.e., the thickness and modulus of each layer). The second component computes the reliability of the performance function using the Rackwitz-Fiessler FORM (17). Equation 3 presents the performance equation (PE) for the ERAPave PP rutting predictive model: where RD limit and RD predicted are the limit or failure rut depth and predicted rut depth, respectively.

Reliability Methods
Pavement ME Reliability. Pavement ME is a reliabilitybased mechanistic-empirical pavement design approach and is routinely implemented to optimize pavement sections for distresses, such as fatigue cracking and rutting. The reliability approach in Pavement ME does not directly consider input variability influence on predicted performance. It assumes instead each major distress to be a key output of interest and models these outputs as random variables. A normal distribution is assumed to represent the randomness associated with these distresses (4). The mean value of the distress, which represents 50% reliability, is estimated using the nominal input values. The SE, which defines the overall variance of the predicted distress, is estimated using outputs that are obtained during the calibration process and includes errors such as input variability, construction error, and model bias. In cases where a deterministic failure criterion is used, the reliability of the performance equation depends on the first and second moments of the predicted distress. For rutting, the following equation is used to calculate the amount of rut that should be permissible at the specified design period and designated target reliability level: where RD ZR represents rutting at the designated target reliability level, RD m is the rutting predicted using average nominal values, RD SE is the SE of the rutting prediction method, and Z R is the standard normal deviate. The SE for the total predicted rutting is estimated using the following equation: where RD SE,AC , RD SE,GB , and RD SE,SG represent the SE associated with the predictions for the asphalt concrete (AC) layer, the unbound granular base (GB) layers, and the subgrade (SG), respectively. The following three equations are used to estimate the SE for each respective layer (4): where RD m,AC , RD m,GB , and RD m,SG are the rutting values determined for the AC, GBs, and SG using nominal average input values, respectively.
Monte Carlo Simulation. MCS uses randomly sampled input variables to determine the reliability of engineering systems. MCS generates the required sampling points using the PDF of the individual random variables. In the simulation, many randomly generated set of the basic random variables X (i.e., x i=1,2.,n ) are evaluated deterministically through numerical experimentation to determine whether or not each of the realizations fulfils the requirements of the limit state condition. Those outcomes that do not fulfill the requirements of the limit state function, in the case when g(x) ł 0 defines failure, are considered to represent a failure condition. Reliability is estimated by dividing the number of simulation cycles that fulfil the condition g(x) . 0 (N P ) by the total number of simulation cycles (N) as follows: A MCS outcome is highly dependent on the number of realizations that are used to evaluate the performance function. It is expected that the accuracy of this outcome increases as the number of cycles increases and this value would attain its true value as the number of samples approaches infinity. However, this might be very hard to achieve, especially for analysis such as pavement performance evaluation where a single computation requires a considerable amount of processing time. MCS can also be used to generate the actual PDF of distresses, such as fatigue cracking and rutting, which would otherwise be very difficult to determine through numerical analysis (6).
The First-Order Reliability Method. The FORM is an analytical approximation approach and has been implemented for a variety of pavement reliability analysis problems (7,9,11). The FORM utilizes two-level idealization and simplification steps to find a solution for the reliability problem. In the first simplification, all the random variables are transferred from the original random space to the standard normal space where variables are independent and uncorrelated. In the second step, a linear function is used to approximate the limit state function at the most probable point (MPP), forming a hyperplane that divides the random space into safe and failure regions (17,18).
The failure or design point is determined through an optimization and reliability is estimated by taking the minimum distance from the origin to the failure point. This minimum distance is called the reliability index (b) and reliability (R) is estimated as follows: The FORM includes two approaches that are employed depending on the complexity of the performance function and the amount of information available with respect to input variability (17,18). The Rackwitz-Fiessler FORM is well suited to problems where the performance function features high nonlinearity, and the random variables are characterized by a non-normal PDF. At the checking point, the Rackwitz-Fiessler FORM approximates the non-normal random variables using an equivalent normal variable and estimates reliability through an integration procedure that utilizes the partial derivatives of the performance function (17).
Analytical-based reliability analysis methods such as the FORM require the pavement performance equation to be expressed by an explicit closed-form function of the design input variables, which in most cases is not available. RSM can be implemented to overcome this problem and to establish an explicit mathematical expression for the implicit pavement performance equation. RSMs utilize mathematical and statistical techniques and regression analysis to generate a first-or second-order polynomial surrogate model (19). The CCD RSM is well adapted for generating higher order surrogate models that can readily be implemented for pavement reliability analysis (16). These surrogate models are continuously differentiable in the regions of the MPP and capture the interaction between parameters and critical points of the function. Equation 11 presents the mathematical formulation for a second-degree performance equation (y) that takes independent variables x 1 , x 2 ,.,x k as an input: where b 0 is the model constant and e represents the residual associated to the experiments. Here, b i , b ij , b ii are the coefficients for the linear, interaction, and quadratic terms, respectively.

Distress Model and Pavement Sections
ERAPave PP ERAPave PP employs mechanistic-empirical design principles for the prediction of distresses in flexible pavements. For the prevailing environmental and traffic conditions, the tool primarily optimizes flexible pavements for fatigue cracking and rutting distresses. It can also analyze pavement structures for conditions that are typical of cold climates, such as frost heave and studded tire wear. Several studies that have been conducted using laboratory and APT investigations have shown that the pavement response and performance models of ERAPave PP are capable of delivering acceptable results (20)(21)(22). The process that is used in ERAPave PP for the performance evaluation of flexible pavements is presented in Figure 1. As can be seen in Figure 1, the input module provides all the required information for pavement response and pavement performance evaluations. These inputs comprise mixture properties, structural inputs, environmental conditions, and traffic factors. A multi-layered linear elastic theory (MLET)-based analysis is used to determine the field variables, such as strain, stress, and deflection, for the prevailing pavement analysis condition (1,23). Damage is computed by the empirical component of the tool and this damage is accumulated on the basis of time or number of traffic cycles.
ERAPave PP uses the layer strain approach for the prediction of rutting in flexible pavements. This approach obtains the total observed surface rutting by computing and aggregating the contribution of each layer (24). This requires separate models or model coefficients for the AC and unbound granular layers (UGLs). This is mainly because pavement structures are constructed using layers that have different material compositions. ERAPave PP takes into account the stress history dependence behavior of the material while computing the evolution of permanent deformation for each layer of the pavement structure. For this, a time-hardening approach that accounts for the stress path or stress level influence on permanent deformation accumulation is implemented (25,26).

Permanent Deformation Prediction for the Asphalt Concrete
Layer. The permanent deformation behavior of the AC layers is affected by factors such as temperature, binder type and content, aggregate gradation, and structural thickness. Factors related to traffic are also observed to influence this behavior (11). ERAPave PP predicts the development and evolution of permanent deformation in the AC layers using the model that was originally developed for Pavement ME (4,27). The mathematical formulation of this model is presented in Equation 12. The model captures the permanent deformation behavior of the AC layers using factors such as traffic volume (N), temperature (T), and applied load level: where e p and e r are the resilient and accumulated permanent strains, respectively. Here, a = 0.03, b = 1.85, and c = 0.27 are model coefficients

Permanent Deformation Prediction in Unbound Granular
Layers. Permanent deformation or rutting is the primary failure mode for the layers that are comprised of unbound granular materials (UGMs). The permanent deformation behavior of these layers is highly dependent on factors such as structural thickness, traffic volume, stress level, stress history, moisture content, and aggregate gradation (28)(29)(30). ERAPave PP employs a strainbased model for the prediction of permanent deformation in layers such as the base, subbase, and SG (21,31).
As can be seen in Equation 13, the model predicts permanent strain (e p ) using as an input number of applied traffic cycles (N) and resilient strain (e r ): where a = 3.0 and b-which has values of 100, 200, and 250 for the base, subbase, and SG model coefficients, respectively.

Pavement Sections
This study was carried out using two APT sections. The studied sections (i.e., SE-14 and SE-18) represent different design features and are representative of Swedish inservice conditions. The sections were constructed using workmanship and equipment that are employed for the construction of actual field pavements. A schematic representation of the two sections is presented in Figure 2.
As can be seen in Figure 2, the two structures consist of four distinct layers: a bituminous surface, a crushed aggregate base and subbase, and a sand SG, which is constructed over a rigid concrete floor.
As can be seen in Figure 2, the studied sections were embedded with sensors and gauges of different kinds. Instrumentation was provided in the longitudinal direction along the main loading path and a multiple set of sensors and gauges was provided for each variable of interest. The instruments provided a continuous measurement of the field variables for the whole testing period. For the measurement of the field variables, asphalt strain gauges (ASGs), soil pressure cells (SPCs), and vertical strain gauges (eMU coils) were provided. Measurements of interest were the longitudinal and transverse horizontal strains at the bottom of the AC layer, the compressive vertical pressure at different locations within the unbound and SG layers, and the fulldepth vertical strain.
The APT facility is a fully insulated system where environmental factors such as moisture and temperature are fully controlled and continuously measured. Traffic loading is applied during the main accelerated loading stage using mobile HVS (Mark IV). The HVS is capable of applying different load levels and wheel configurations while moving at a speed of 12 km/h (20,22). The lateral application of the traffic in the transverse direction is also carried out during loading.
Surface rutting is measured using laser beams that are installed at five different locations along the longitudinal direction, while the permanent deformation of each layer is measured using the eMU coils. FWD measurements are performed on numerous occasions during construction and testing to assess the structure integrity of the pavement structure. It is also customary to perform surface elevation measurements that characterize the homogeneity of the structural thickness.  Table 1 provides the layer thickness and modulus of the individual layers for each APT section. These values represent average nominal conditions, and the moduli were calculated using average FWD measurements. For pavement response and performance calculations, the layer moduli obtained through the FWD were optimized further using the full-depth vertical strain measurements. This is required so that the material properties can reflect accurately the loading condition of the HVS. Testing was performed at a constant pavement temperature of 10°C while the moisture condition in the UGLs and SG was unaffected. During the main accelerated loading phase, the sections were subjected to a traffic loading that consists of single axles of different magnitudes and volumes. A dual wheel configuration with a tire pressure of 800 kPa was used for the traffic application. SE-14 was subjected in the first phase to 580,000 traffic cycles of 80 kN, while in the second phase 600,000 cycles of 100 kN were applied. For SE-18, the traffic was applied in three phases. Axle load magnitudes of 80, 100, and 120 kN, with traffic volumes of 481,000, 434,000, and 300,000 cycles, were applied during the three loading phases, respectively.

Variability Characterization
The design of pavement structures requires a multitude of inputs. As all of these inputs do not have the same level of impact on predicted performance, it will be uneconomical to consider all of them as random variables. A sensitivity analysis has shown that inputs related to traffic, environmental conditions, material properties, and structure influence significantly the rutting development in flexible pavements (11). As the APT was performed at a constant temperature and using a fixed volume and magnitude of traffic loading, the only inputs that are expected to be variable are the inputs related to the material and structure. Therefore, the thickness and modulus of each layer were modeled as random variables. The full probability approach, which uses a PDF along with the first and second moments of the data (i.e., mean sand standard deviation), was used to characterize the variability. In addition to the most common PDFs, such as normal and lognormal, PDFs such as Weibull, gamma, and general extreme value were used.
Layer Modulus. FWD testing were conducted using three dynamic load levels (30,50 , and 65 kN). The pavement response caused by the applied dynamic loads was received using geophones. As a variability study, multiple sets of testing were conducted along the longitudinal and transverse directions. All testing was conducted following the recommended guidelines and test protocols. Figure 3 provides the surface deflections for the load impact point along the longitudinal direction. As can be seen in Figure 3, a and b, these deflections show a remarkable   variation, and this variation is consistent among the three dynamic loads. A backcalculation analysis that involves MLET-based analysis was employed to determine the modulus of each layer (32). To simplify the optimization and to reduce prediction error, the layers were grouped together to form a three-layered pavement structure: the AC, UGL, and SG. Furthermore, to eliminate the uncertainty associated with pavement thickness, the homogeneity of the structural thicknesses was carefully assessed. Table 2 presents the statistical analysis for each section layer modulus. As can be seen in Table 2, the SG is more homogenous than the asphalt and UGLs. The same pattern of variability was also observed in previously studied APT sections (12). However, for field sections, the AC and UGLs are more homogenous than the SG (6). For this study, coefficients of variation (CVs) of 20%, 15%, and 5% were selected to represent the variability associated with the AC, UGLs, and SG, respectively.
The dispersion associated with the measured modulus was tested with various types of PDFs. Figure 4, a-f, presents the measured moduli of each layer fitted with a lognormal distribution. It is not uncommon for the measured modulus to be described and modeled concurrently by multiple PDFs. In addition, the sample size has a significant influence on the type of PDF that best fits the measurement. The histogram plots have also shown that the measured moduli of each layer exhibit a wide range of values. For this study, a lognormal PDF was selected to characterize the layer modulus variability.
Layer Thickness. Elevation level measurements that characterize the uniformity of the compacted surface were carried out after the placement and compaction of each layer. As a variability study, multiple elevation level measurements were performed in both longitudinal and transversal directions. Figure 5 provides the center line average elevation level measurements for the longitudinal direction. As can be seen in Figure 5, a and b, there is a negative correlation between the layers. This is expected, as the intent during construction is to attain a smooth flat surface.
The thickness profile of each layer was obtained by deducting its elevation level from the layers adjacent to it. Table 3 presents the statistical analysis for each section layer thickness. As can be seen in Table 3, the CV for each layer falls within a narrow range, showing the consistent construction practice at the APT facility. For this study, CVs of 10%, 15%, and 10% were selected to represent the variability associated with the AC, UGL base, and UGL subbase, respectively. Figure 6 presents the measured thickness of each layer fitted with a normal PDF for the two studied sections. As can be seen in Figure 6, a-f, the variability of layer thickness can be described by a normal PDF. The variability of the layer thickness was also observed as exhibiting a wide range of values. For this study, a normal PDF was selected to characterize layer thickness variability.

Reliability Analysis
Predicted Rutting Variability. MCS was carried out to generate the frequency distribution of predicted surface rutting. As ERAPave PP requires a considerable amount of time to carry out many MCSs, the analysis was limited to 500 cycles. Although this number might not be enough to truly reflect the actual variability of predicted rutting, it will provide valuable insights with respect to the expected range of values. Statistical analysis has shown that a lognormal PDF with a CV range of 15%-30% can be used to describe the variability of the MCS-generated rutting. Figure 7 presents the MCS-generated rutting histogram fitted with a lognormal PDF for the two sections. The averaged measured value of the rutting is also included in the figures. As can be seen in Figure 7, a and b, there is an overlap between the measured average values and most frequent values generated by MCS, which shows the capability of the ERAPave PP rutting prediction model.
Surrogate Model Validation. The CCD RSM generates the different level matrices using experiments that involve a full factorial design and an additional design where experimental points are at a certain distance from the Note: PDF = probability density function; AC = asphalt concrete; CV = coefficient of variation; UGL = unbound granular layer; GEV = Generalaized extreme value.
center. This was achieved using the MATLAB inbuilt CCD function. A one-way analysis of variance (ANOVA), which was performed on 100 randomly generated experimental points, has shown that the generated surrogate models can deliver acceptable predictions. A direct comparison was also made between the rutting that was predicted by the actual and surrogate models. Figure 8 presents the comparison between the two values for the two studied sections. As can be seen in Figure 8, a and   Reliability Comparison. The limit state equation in Equation 3 was used to estimate the reliability of the two sections. As experimental sections that represent actual field conditions, the failure criteria for the two sections were established by carefully studying the existing guidelines and specifications. Thus, rut depths of 8 and 15 mm were established as failure criteria for sections SE-14 and SE-18, respectively. Two separate failure criteria were used for the two sections, as the sections represent different design and functional conditions. The reliability of the two sections was assessed using the three approaches. In the first approach, MCS in conjunction with a surrogate model was used to determine the ''exact'' reliability of the two sections. MCS cycles of 10,000 were used, as this number of cycles was observed to be reasonable. The second approach utilizes the reliability methodology incorporated in Pavement ME. Equations 5-8 were used to determine the required SE for each layer and for the total observed surface rutting. The third approach estimates the reliability of the two sections using the FORM. Table 4 presents the estimated reliabilities according to the three approaches for the two sections.
As can be seen in Table 4, the estimated reliabilities of the three approaches differ significantly for both sections. The estimated reliabilities are more precise for SE-14 in comparison with SE-18. For SE-14, when compared with MCS, the reliabilities of Pavement ME and the FORM have normalized errors (NE) of 6.5% and 2.9%, respectively. For SE-18, these values are 12.9% and 16.8%, respectively. As can be seen in Figure 7b, the MCS-generated histogram of SE-18 exhibits a wide dispersion and might be the reason for the wide gap between the Pavement ME and FORM reliabilities.
In addition, for the Pavement ME approach, CV values of 29% and 23% were obtained for SE-14 and SE-18 rutting, respectively. This difference in CV values might be the reason why the Pavement ME approach delivered a higher reliability for SE-18 than SE-14. For the FORM approach, when compared with SE-14, the reliability estimate of SE-18 diverges by a significant margin from its corresponding MCS reliability. This might be attributed to the SE-18 surrogate model, which as can be seen in Figure 8b is less accurate and precise than the surrogate model in SE-14.

Conclusions
Reliability analysis approaches that rely on a probabilistic method of uncertainty propagation have not been  fully utilized in pavement design. This is mainly attributed to factors such as the level of uncertainty associated with the pavement design process, the implicit nature of pavement performance evaluation, and the lack of accurate data on input variability. This has led to the adoption of simplified pavement reliability analysis methods that may not deliver the intended benefits. A comprehensive reliability analysis methodology that captures the combined variance of input variabilities on predicted performance is required. This would deliver designs of uniform performance, providing the right platform for the technical, economic, and environmental assessment of pavement structures. The variability characterization of the APT sections has shown the challenge of obtaining uniform and homogeneous structures even in constructions that are performed in a controlled environment. The layer modulus and thickness of the APT sections exhibit less variation in comparison with field pavement sections. Unlike field sections, the SG of the APT structures is relatively uniform and homogenous than the AC and UGLs. Lognormal and normal PDFs can be used to describe the variability of the layer modulus and thickness, respectively. A lognormal PDF is also observed to be appropriate for describing the dispersion associated with MCS-generated rut depth.
The CCD RSM has managed to generate surrogate models that are relatively accurate, while the evaluated reliability analysis approaches are observed to provide reasonable estimates. The reliability estimates of the Pavement ME approach are highly dependent on the SE of each layer, as it lacks the mechanism to consider directly the variability associated with inputs. It is also   observed that the currently utilized SE equations are positively skewed to higher rut depths, requiring further adjustments. MCS on the basis of a surrogate model can be a good alternative to obtain an ''exact'' estimate of the reliability, minimizing the forbidding computational time. As estimated ''exact'' values are highly dependent on the chosen number of simulated cycles, a further interpretation of these ''exact'' results is required. The FORM coupled with properly validated surrogate models can provide acceptable reliability estimates. In addition to reliability estimates, the FORM provides information on failure points and directional cosines of each random variable, making it preferable for the development of a deterministic reliability-based design procedure for pavement application.

Author Contributions
The authors confirm contribution to the paper as follows: study conception and design: Y. Dinegdae, S. Erlingsson; data collection: Y. Dinegdae; analysis and interpretation of results: Y. Dinegdae, A. Ahmed, S. Erlingsson; draft manuscript preparation: Y. Dinegdae, S. Erlingsson. All authors reviewed the results and approved the final version of the manuscript.

Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research was sponsored by the Swedish Transport Administration (Trafikverket).

Data Accessibility Statement
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.