2. Materials and Methods
This section presents the methods used for data collection, the roadway segmentation process, the procedures for the characterization of road conditions, the deterioration models applied in this study and the algorithm proposed to identify the best combination of projects for road preservation. It is important to mention that the present model does not consider political factors but rather purely the life cycle costs related to unpaved and paved roads.
2.1. Road Segmentation and Classification
Four different types of data were collected to develop this research: roadway segmentation data, traffic counts, road width values and road conditions information. Roadway segmentation was conducted by defining segments that start and end at points of new construction, intersections with other roads or other alterations along the road.
Data from unpaved roads were directly collected in the field, which imposes some challenges regarding accessibility and terrain conditions. For example, the average width of the sections was determined using tape measurements. Rutting was obtained using a standardized ruler and straightedge. The average daily traffic volume was determined using the methodology proposed in the DNIT manual for traffic studies. Data regarding paved roads were obtained from the “Rodovias” website managed by the government of the state of MG.
Based on the segmentation process and using the ArcGIS software, a total of 18 unpaved road segments were delineated in the ZMC of MG, alongside 88 paved roads distributed throughout MG.
Figure 1 shows a paved road in MG awaiting maintenance to prolong its service life, while
Figure 2 illustrates an unpaved road in poor condition in MG.
The unpaved road conditions were characterized by seven distinct parameters: cross-section/crown, roadside drainage, rutting, potholes, loose aggregate, dust and corrugation. The Pavement Surface Evaluation and Rating (PASER) manual describes that the surface of unpaved roads can be rated on a five-point scale, where five represents an excellent condition and one indicates a failed condition for unpaved roads. For this paper focusing on MG roads, we proposed transforming this five-point scale into a nine-point scale, following the Gravel Roads Rating System (GRRS) guidelines [
15]. This approach aimed to reduce errors when human raters need to assess the infrastructure conditions. Then, a Pavement Serviceability Index (PSI) was calculated, following the recommendations of the RQRG created by the Wyoming Technology Transfer Center (WYT2/LTAP) [
43]. This nine-points scale was applied for unpaved roads, based on the Overall Road Condition Index (ORCI). This methodology resulted in the ride quality, and the extreme values nine and one represent the best and worst conditions, respectively.
On the other hand, paved roads were assessed using the AASHTO methodology, which ranks conditions from 0 to 100, with 100 representing a newly constructed road and 0 indicating a failed road. According to this rating system, the primary variables are the distresses rutting, fatigue cracking, longitudinal cracking, transverse cracking, patches, potholes, raveling, bleeding and roughness. The evaluation results were reported using the AASHTO guide to road classification, categorizing roads as “Excellent”, “Good”, “Fair”, “Poor”, or “Failed” [
44]. Therefore, a five-point scale was used for the paved roads investigated in this research.
2.2. Deterioration Models
Huntington and Ksaibati [
15] developed a deterioration model for roads in Wyoming (USA), which served as a foundation in the current study for elaborating the deterioration model for the Brazilian roads.
An average deterioration was estimated in the present research, based on overall road conditions. At first, the obtained pavement condition database was appropriately organized according to the defined condition states. Then, the first aspect of the data quantification process involved categorizing the pavements. The subsequent step encompassed generating data matrices for each pavement classification.
The deterioration evolution of unpaved and paved roads was predicted through the Markovian chains. After defining the number of conditions (
k) to be considered in the statistical evaluation, the transition matrix with dimensions
k ×
k was used to evaluate the transition probability of states [
45], as indicated in Equation (1).
The unit value presented in this matrix (
PΔ
t) corresponds to the maximum input value for each event. The probability value of
pkk was defined considering that there was no further condition state after the k state, so it is a limited value [
46]. The transition probability from a state
i to a state
j, for a generic time interval (Δ
t), is represented by
pij, with
i =
j = 1, …,
k, as indicated in Equation (2) [
47].
Then, Markovian chain techniques and regression analyses were conducted using the R statistical software (version 3.5). Using the probabilistic Markovian process and beginning with the data matrix, the transition matrix was calculated to derive the probability distribution matrix. Subsequently, this matrix was used to determine the expected age of the roads based on their conditions, resulting in the development of a statistical model for road deterioration.
2.3. Algorithm Proposed to Identify the Best Combination of Projects for Road Preservation
Five treatment options were proposed for the unpaved roads: general maintenance, chemical and mechanical stabilization, preventive rehabilitation, medium rehabilitation and paving procedures (GM, 1-R, 2-R, 3-R and 4-R, respectively). In addition, six treatment options were proposed for paved roads: general maintenance, preventive rehabilitation, minor rehabilitation, medium rehabilitation, major rehabilitation and reconstruction (GM, 1-R, 2-R, 3-R, 4-R and 5-R, respectively). The characteristics of these treatment options were discussed in
Section 3. The decision trees for the unpaved and paved roads are illustrated in
Figure 3 and
Figure 4, respectively.
Hence, it is evident that the proposed URMS and PMS take into account not only the cost factor but also local conditions, including the ORCI and the average daily traffic (ADT) data. The objective of the developed model was to maximize a function that considers three independent parameters: the road condition, ADT, and cost-factor, as indicated in Equations (3) and (4).
where ADT
i expresses the average daily traffic for road
i, Cost-Factor
i is the function of the treatment type,
xi is an integer equal to one if the project is selected and zero if it is not selected and ORCI represents the minimum index among the distress indices. This is a combinatorial optimization problem where one selects a collection of projects of maximum value while satisfying some weight constraints.
3. Results and Discussion
The collected information for unpaved and paved roads is presented in the comprehensive database summarized in
Table 1 and
Table 2, respectively. The information provided in these tables is also graphically represented in
Figure 5 and
Figure 6, respectively. A detailed analysis of the traffic counts reveals significant patterns in vehicle usage. The unpaved road segments experienced traffic volumes of less than 400 vehicles per day, indicating lighter usage. In contrast, a notable portion of the paved road segments, approximately 20%, handle a heavier traffic load, with daily vehicle counts exceeding 400. The road widths of the segments varied significantly, ranging from 3 to 7 m. The collective length of the road segments was substantial, with the 18 unpaved segments and 88 paved segments together spanning a total distance of 15,236 km. This extensive network underscores the importance of understanding traffic patterns and road characteristics for the effective management of Brazilian roads.
The Markovian process enabled the determination of the following matrices: the transition probability data, transition probability matrix and probability distribution matrix. The results presented in
Table 3 and
Table 4 represent the transition probability data obtained from the unpaved and paved roads datasets, respectively.
Table 5 shows the transition probability matrix for the unpaved roads. A transition matrix of a Markovian process contains information on how the deterioration process moves from one state to another in terms of probabilities. According to the obtained data presented in
Table 5, it was possible to observe that after 1 year, the percentage of roads that stay at level 9 was represented by the element 0.333. Similarly, the 0.481 element represents the percentage of roads that deteriorate from level 9 to level 8, and 0.185 represents the percentage of roads that deteriorate from level 9 to level 7.
Table 6 represents the transition probability matrix for paved roads of MG. These probabilities are assumed to be constant, i.e., the transition matrix is stationary. According to the data reported in
Table 6, it can be observed that the number 0.595 represents the percentage of roads that stay at level five after 1 year. Similarly, the element 0.405 represents the percentage of roads that deteriorate from level five to level four.
The probability distribution derived from the application of the Markovian probability process for unpaved roads is presented in
Table 7. In this table, states of maximum probability were highlighted in red. For example, comparing the obtained data from year 1 to year 5, there was a probability of 0.004 for the unpaved road segments to stay at level nine and a probability of 0.037 for these segments to be at level two after 5 years. In year 5, the highest probability value (0.327) indicated the unpaved road segments of MG at level five.
Table 8 represents the probability distribution derived from the application of the Markovian probability process to the paved road dataset. Again, states of maximum probability are shown in red. Comparing the distribution from year 1 to year 5, there was a probability of 0.075 for the paved road segments to stay at level five and a probability of 0.024 for the segments to be at level one. However, the highest probability value in year 5 was 0.415, which indicated paved roads at level three.
The deterioration models were derived from regression analyses applied to the mean data obtained from the Markovian distribution. Using these models, it was possible to determine the time taken to transition from an excellent condition to a failed condition. Equations (5) and (6) show the deterioration models obtained for the unpaved and paved roads, respectively. In these equations, URCI indicates the unpaved roads condition index, PCI indicates the pavement condition index and AGE is the age (in years).
These results suggest that unpaved roads tended to deteriorate more rapidly over time. The good alignment between the experimental findings derived from the methodology outlined in this paper and the results of previous studies underscores the robustness of the proposed methodology. According to previous literature, the faster deterioration of unpaved roads is affected by permeability, materials, construction, traffic, environment and drainage factors. Among them, the most important factor is the permeability and weather [
33]. As a result of these factors, tropical roads experience faster deterioration, and unpaved road surfaces are more permeable compared to paved surfaces.
Details related to the five treatment options proposed in this methodology for the unpaved roads (GM, 1-R, 2-R, 3-R and 4-R) are summarized in
Table 9, while
Table 10 shows the six possible treatment options (GM, 1-R, 2-R, 3-R, 4-R and 5-R) proposed for the paved roads. The estimated costs were provided in BRL because it is the currency used in MG, thus reflecting the local economic context.
In
Table 11 and
Figure 7, the results provided by the optimization algorithm for unpaved roads are organized into four distinct scenarios, as follows: (i) conditions under current circumstances (referred to as the “current condition”), (ii) conditions following four cycles of 3 months without any maintenance (termed “do nothing”), (iii) conditions after applying treatments without budgetary constraints (designated as “no budget constraints”) and (iv) conditions after applying treatments necessary for maintaining the road’s condition within a limited budget (referred to as “limited budget”). The rationale behind the selection of a 3-month period stems from the observed tendency for unpaved roads to require the next level of maintenance every 3 months [
24].
Table 12 and
Figure 8 present similar results for four different scenarios considered in the optimization algorithm for the paved roads. In this case, the “do nothing” scenario reflects conditions after 12 months without any maintenance. The reason for this 12-month period is that the typical lifespan of paved roads without maintenance is longer than that of unpaved roads [
24].
Comparisons between the maintenance costs of unpaved and paved roads indicated that the full rehabilitation costs of county paved roads were approximately 10 times higher per kilometer compared to those for unpaved roads, considering the difference in maintenance treatment frequency, which is approximately four times higher for unpaved roads (every 3 months) compared to that for paved roads (every 12 months). The models accurately depicted the pattern where major treatments result in a minimal need for further treatments over the following several years in both situations.
It is important to highlight that other countries may have different optimization priorities. For instance, engineers might aim to maximize the overall network road conditions given an unrestricted budget. In the present study, the options for maximizing network road conditions and minimizing costs were thoroughly examined to showcase the potential of the proposed methodology.