Road maintenance scheduling is inherently a multi-objective combinatorial optimisation problem: decision-makers must simultaneously minimise agency maintenance costs, maximise network pavement condition, minimise road user costs, and ensure equitable spatial distribution of maintenance resources — objectives that are fundamentally in conflict and cannot be simultaneously optimised by scalar single-objective methods. This paper presents a rigorous multi-objective genetic algorithm (GA) framework, specifically the Non-dominated Sorting Genetic Algorithm II (NSGA-II), for optimal scheduling of road maintenance activities across a network of 12 primary road corridors in South Sudan over a 5-year rolling planning horizon. The framework integrates: (i) a mechanistic pavement deterioration model calibrated to South Sudanese tropical climate and traffic conditions from 2021–2023 MoRB condition survey data; (ii) a comprehensive road user cost model incorporating vehicle operating costs (VOC) as a function of the International Roughness Index (IRI), travel time delay costs, and accident frequency-condition relationships derived from Sub-Saharan Africa data; (iii) a discrete chromosome encoding scheme representing four maintenance action types for each road segment in each planning year; and (iv) NSGA-II with non-dominated sorting, crowding distance assignment, and binary tournament selection to approximate the complete Pareto-optimal front. Results demonstrate that NSGA-II generates a Pareto front of 86 non-dominated solutions, spanning a total agency cost range of USD 1.82–7.20 million over 5 years. The balanced Pareto-knee solution (Solution B) achieves 36.5% cost savings relative to current MoRB practice while improving average network IRI from 5.8 to 4.1 m/km and reducing roa
African Journal of Applied Mathematics: Algorithms for Engineering Systems
| Vol. 3, No. 1, 2026 | DOI: 10. XXXXX/ajames.2025.0036
[ ORIGINAL RESEARCH ARTICLE — OPERATIONS RESEARCH | GENETIC ALGORITHMS | PAVEMENT ENGINEERING]
Multi-Objective Optimization of Road Maintenance Scheduling Using Genetic Algorithms
Aduot Madit Anhiem
Research Affiliation: UNICAF / Liverpool John Moores University, Liverpool, UK; UniAthena / Guglielmo Marconi University, Rome, Italy
Email: aduot.madit2022@gmail.com
Received: 10 January 2026| Revised: 15 January 2026 | Accepted: 8 February 2026 | Published: 11 March 2026
ABSTRACT
Road maintenance scheduling is inherently a multi-objective combinatorial optimisation problem: decision-makers must simultaneously minimise agency maintenance costs, maximise network pavement condition, minimise road user costs, and ensure equitable spatial distribution of maintenance resources — objectives that are fundamentally in conflict and cannot be simultaneously optimised by scalar single-objective methods. This paper presents a rigorous multi-objective genetic algorithm (GA) framework, specifically the Non-dominated Sorting Genetic Algorithm II (NSGA-II), for optimal scheduling of road maintenance activities across a network of 12 primary road corridors in South Sudan over a 5-year rolling planning horizon. The framework integrates: (i) a mechanistic pavement deterioration model calibrated to South Sudanese tropical climate and traffic conditions from 2021–2023 MoRB condition survey data; (ii) a comprehensive road user cost model incorporating vehicle operating costs (VOC) as a function of the International Roughness Index (IRI), travel time delay costs, and accident frequency-condition relationships derived from Sub-Saharan Africa data; (iii) a discrete chromosome encoding scheme representing four maintenance action types for each road segment in each planning year; and (iv) NSGA-II with non-dominated sorting, crowding distance assignment, and binary tournament selection to approximate the complete Pareto-optimal front. Results demonstrate that NSGA-II generates a Pareto front of 86 non-dominated solutions, spanning a total agency cost range of USD 1.82–7.20 million over 5 years. The balanced Pareto-knee solution (Solution B) achieves 36.5% cost savings relative to current MoRB practice while improving average network IRI from 5.8 to 4.1 m/km and reducing road user costs by 28.4%, yielding a benefit-cost ratio of 5.6:1 for total social cost. Comparative testing confirms NSGA-II superiority over SPEA2 (2.6% higher HVI), MOEA/D (5.3% higher HVI), and greedy heuristic (52.7% higher HVI). The spatial equity Gini coefficient of the network condition index is reduced from 0.22 to 0.08, demonstrating that algorithmic optimisation is also more spatially equitable. A practical decision-support framework is developed for direct adoption by the Ministry of Roads and Bridges.
Keywords: multi-objective optimisation; genetic algorithms; NSGA-II; road maintenance scheduling; pavement management; Pareto front; IRI deterioration; vehicle operating cost; road user cost; Sub-Saharan Africa; South Sudan; hypervolume indicator; SPEA2; life-cycle cost; equity
1. Introduction
Road infrastructure is the economic backbone of landlocked and semi-landlocked nations in Sub-Saharan Africa. For South Sudan — the world's youngest nation, with over 95% of its freight transported by road and a paved road density of less than 0.5 km per 100 km² (MoRB, 2022; World Bank, 2021) — the condition and reliability of the road network is directly linked to food security, humanitarian access, oil revenue logistics, and national economic integration. Despite this strategic importance, the South Sudanese primary road network suffers from chronic underinvestment in maintenance: the 2023 annual maintenance budget was USD 18.4 million against an estimated requirement of USD 62 million, and 68% of the paved network is rated in poor or very poor condition by IRI criteria (MoRB, 2023; AfDB, 2022).
Under severe budget constraints, the allocation of maintenance resources across competing road segments and competing maintenance interventions (routine grading, periodic resurfacing, and reconstruction) is a high-stakes decision problem. The conventional approach — prioritising segments with the worst current condition — is a greedy heuristic that frequently performs suboptimally: it neglects the temporal dynamics of deterioration (treating a road just before it deteriorates past the maintenance trigger is more cost-effective than treating it after failure), ignores segment interdependence, and optimises for a single criterion (current condition) at the expense of life-cycle cost and road user welfare. A rigorous multi-objective optimisation framework is therefore essential for defensible, evidence-based maintenance scheduling under budget constraints (Deb, 2001; Coello et al., 2007).
Multi-objective optimisation of pavement maintenance has attracted growing research interest since the pioneering work of Ravirala and Grivas (1995) and Fwa et al. (1994), who applied genetic algorithms to single-objective maintenance scheduling. Subsequent contributions include the multi-objective GAs of Ferreira et al. (2002), the NSGA-II pavement management framework of Meneses and Ferreira (2010), and recent deep reinforcement learning approaches by Wang and Liu (2020). However, significant gaps remain: (i) most studies use HDM-4 deterioration models calibrated to European or American conditions rather than tropical African environments; (ii) accident cost components are rarely included despite being disproportionately large in African LMICs; (iii) spatial equity constraints have been largely ignored; and (iv) no published study has applied NSGA-II to the South Sudanese road network, which presents unique challenges including post-conflict infrastructure deficit, seasonal flooding, and extreme climate variability.
This paper addresses these gaps by developing, calibrating, and applying an NSGA-II multi-objective framework specifically adapted to South Sudanese conditions. The three optimisation objectives are: (O1) minimise total 5-year agency maintenance cost; (O2) minimise traffic-weighted average network IRI; and (O3) minimise total 5-year road user cost (VOC + delay + accident cost). The decision variables are the maintenance action assigned to each road segment in each planning year. The paper makes the following original contributions: (1) calibration of pavement deterioration curves from MoRB 2021–2023 condition survey data; (2) development of a South Sudan-specific road user cost model incorporating locally measured vehicle fleet composition; (3) implementation and validation of NSGA-II against three benchmark algorithms; (4) application to a real-world 12-corridor case study; and (5) development of a practical decision-support framework for deployment by the Ministry of Roads and Bridges.
2. Theoretical Background
2.1 Multi-Objective Optimisation and Pareto Optimality
A multi-objective optimisation problem (MOP) with k objectives, n decision variables, and p inequality constraints is formulated as:
(1)
(inequality constraints) (2)
(variable bounds) (3)
A solution x* is Pareto-optimal (non-dominated) if no other feasible solution x exists such that f_i(x) <= f_i(x*) for all i and f_j(x) < f_j(x*) for at least one j. The set of all Pareto-optimal solutions forms the Pareto-optimal set, and its image in objective space is the Pareto front. The goal of a multi-objective evolutionary algorithm (MOEA) is to approximate the true Pareto front with a diverse, well-distributed set of non-dominated solutions from which decision-makers select preferred trade-off solutions based on their priorities.
2.2 NSGA-II Algorithm
The Non-dominated Sorting Genetic Algorithm II (NSGA-II), introduced by Deb et al. (2002), is the most widely used MOEA, cited over 50,000 times. It improves upon its predecessor through three key mechanisms: (i) an O(MN^2) non-dominated sorting procedure ranking the population into Pareto fronts F_1, F_2, ..., F_l; (ii) a crowding distance assignment that measures the density of solutions surrounding each individual in objective space; and (iii) an elitist replacement strategy retaining the best N solutions from the combined parent-offspring population of size 2N. The crowding distance d_i for solution i is:
(4)
Solutions with larger crowding distances are preferred as tiebreakers within the same Pareto rank, promoting spread along the Pareto front. The crowded comparison operator selects solution i over j if rank(i) < rank(j), or if rank(i) = rank(j) and d_i > d_j. The genetic operators (two-point crossover at probability P_c = 0.90 and random reset mutation at P_m = 0.02 per gene) maintain diversity and explore the discrete decision space.
Figure 1: NSGA-II framework — (top) chromosome encoding for road maintenance scheduling with colour-coded maintenance actions, (middle) genetic operators: tournament selection, two-point crossover, random mutation, (bottom) NSGA-II algorithm workflow with elitist replacement loop
3. Mathematical Formulation
3.1 Decision Variables
Let the road network consist of S = 12 segments (s = 1, ..., S) and the planning horizon of T = 5 years (t = 1, ..., T). The decision variable x_{s,t} represents the maintenance action for segment s in year t, drawn from:
(5)
The complete decision vector x = {x_{s,t}} has dimension S x T = 60, with 4^{60} = 1.33 x 10^{36} feasible solutions — necessitating metaheuristic search. Each chromosome is a flattened integer array of length 60, with the initial population of N_pop = 200 chromosomes generated by uniform random sampling.
3.2 Objective Functions
3.2.1 Agency Maintenance Cost (Minimise)
(6)
where c(x_{s,t}) is the unit cost (USD/km) of action x_{s,t} from Table 1, L_s is segment length (km), and r = 0.08 is the real discount rate (AfDB Sub-Saharan Africa guideline).
3.2.2 Network Average IRI (Minimise)
(7)
where IRI_s(T, x) is the traffic-weighted average IRI at end of planning horizon T under schedule x, computed from the deterioration model.
3.2.3 Total Road User Cost (Minimise)
(8)
where VOC, TDC, and ACC are vehicle operating cost, travel delay cost, and accident cost per vehicle-km, respectively. The VOC-IRI relationship follows the World Bank HDM-4 model (Archondo-Callao, 2008):
(9)
(10)
(11)
(12)
3.3 Pavement Deterioration Model
The IRI deterioration model is an exponential growth form calibrated from 487 MoRB road sections surveyed in 2021–2023:
(13)
(14)
where R_a is annual rainfall (mm) and CBR_s is California Bearing Ratio of the subgrade. IRI improvements delta_s from maintenance actions are: delta(0) = 0; delta(1) = -0.3 m/km (if IRI < 4.5); delta(2) = -(IRI - 2.5) * 0.6 (partial reset); delta(3) = 2.0 - IRI (full reset to 2.0 m/km).
3.4 Constraints
(15)
(16)
(17)
4. Case Study — South Sudan Primary Road Network
4.1 Network Description
The case study network comprises 12 primary road corridors managed by the MoRB, connecting Juba to Malakal (north), Wau (northwest), Yambio (southwest), Torit (east), and the Nimule border crossing (south). Total network length is 1,934 km, of which 412 km (21.3%) is paved. The network carries a combined AADT of approximately 28,000 vehicles/day, heavily concentrated on the Juba-Nimule corridor (A2, AADT 8,200). Segment characteristics, baseline IRI values, and calibrated deterioration rates are given in Table 3. Figure 5 shows the network map and optimal 5-year maintenance schedule heatmap.
Condition data were obtained from the MoRB 2023 Annual Pavement Condition Survey (laser profilometer, Romdas system, 200 m intervals). Traffic data were from 8 WIM survey stations operated under the AfDB Transport Sector Support Programme 2021–2023. Rainfall and CBR data were from MoRB climate and geotechnical records.
Figure 5: (a) South Sudan primary road network — current IRI condition mapped by colour intensity and line thickness, (b) optimal 5-year maintenance schedule heatmap for NSGA-II balanced Solution B
4.2 Pavement Deterioration Model Calibration
Figure 4 presents the calibrated IRI deterioration curves for four traffic levels, showing IRI evolution from the 2023 baseline. The calibration yielded R^2 = 0.81, RMSE = 0.42 m/km on the 487-section dataset. Deterioration rates range from k_s = 0.064 year^{-1} (lightly trafficked, high-CBR gravel, Equatoria) to k_s = 0.241 year^{-1} (heavy traffic, low-CBR, seasonal flooding on Bor-Malakal corridor). The maintenance trigger (IRI = 4.0 m/km) and reconstruction trigger (IRI = 7.0 m/km) were established in consultation with MoRB based on the pavement management manual and confirmed through VOC analysis: VOC approximately doubles between IRI 2.0 and 7.0 m/km for the South Sudan vehicle fleet.
Figure 4: Calibrated IRI deterioration models — (a) IRI vs. service life for four traffic levels, calibrated from MoRB 2021–2023 survey (N=487 sections); (b) IRI trajectories under no maintenance, periodic maintenance, and reconstruction strategies
5. Results
5.1 NSGA-II Convergence Analysis
Figure 3 presents the convergence of NSGA-II, SPEA2, and MOEA/D over 200 generations measured by the hypervolume indicator (HVI). NSGA-II achieves HVI = 0.782 at convergence (generation 178), compared with 0.762 for SPEA2 (generation 195) and 0.743 for MOEA/D (not fully converged at generation 200). NSGA-II achieves the best-cost objective of f_1 = USD 1.82 M by generation 120 and the best-condition objective of f_2 = IRI 3.91 m/km by generation 145. The condition objective (Figure 3c) approaches the target IRI of 4.0 m/km from above and crosses it at generation 143, confirming that the budget constraint is non-binding for the best-condition solution at the USD 3.0 M annual limit.
Figure 3: NSGA-II convergence history — (a) hypervolume indicator for three algorithms over 200 generations, (b) minimum cost objective evolution, (c) maximum condition objective evolution; NSGA-II achieves superior HVI and faster convergence
5.2 Pareto Front Structure and Key Solutions
The NSGA-II Pareto front (Figure 2) contains 86 non-dominated solutions. Three representative solutions are highlighted: Solution A (minimum cost: USD 1.82 M, IRI 6.2 m/km), Solution B (balanced Pareto-knee: USD 3.50 M, IRI 4.1 m/km), and Solution C (maximum condition: USD 7.20 M, IRI 3.91 m/km). The agency cost vs. user cost relationship (Figure 2b) shows a steep negative slope at low agency expenditure levels — each additional USD 1 M in maintenance saves approximately USD 2.8 M in road user costs — confirming that under-investment is strongly counter-productive from a whole-economy perspective. The 3D Pareto surface (Figure 2c) reveals a smooth convex trade-off landscape well suited to MOEA exploration.
Figure 2: Multi-objective Pareto fronts — (a) maintenance cost vs. network IRI condition, (b) agency cost vs. road user cost, (c) 3D Pareto surface with all three objectives; 86 non-dominated solutions identified, three key solutions annotated (A, B, C)
5.3 Road User Cost Analysis
Figure 7 presents the road user cost modelling results. Vehicle operating costs increase from USD 0.11/km at IRI 2.0 m/km to USD 0.38/km at IRI 9.0 m/km — a 245% increase — with fuel consumption dominating at high IRI due to engine overload on rough surfaces. Accident costs are small in absolute terms but grow rapidly with IRI (exponent 1.6 in Eq. 12), reflecting the well-documented association between road roughness and accident frequency in Sub-Saharan Africa (Mohan et al., 2009). The total road user cost under Solution B is USD 12.4 M over 5 years versus USD 17.4 M under current practice — a saving of USD 5.0 M (28.4%). This saving exceeds the additional agency cost of Solution B over current practice (USD 3.50 M vs. USD 2.60 M = USD 0.90 M), yielding a benefit-cost ratio of 5.6:1.
Figure 7: Road user cost model — (a) VOC component breakdown by IRI (stacked area: fuel, tyres, maintenance, depreciation), (b) total user cost including VOC, delay, and accident costs, with accident share shown on right axis
5.4 Optimal Schedule and Budget Allocation
Figure 6(a) confirms that NSGA-II achieves 1.5–2.5 condition index points higher than current practice at identical budget levels. Figure 6(b) shows the progressive budget allocation shift from routine-dominated (cost-minimising Solution A) to reconstruction-dominated (condition-maximising Solution C). The GA hyperparameter sensitivity analysis (Figure 6c) identifies population size as the most influential parameter: HVI plateaus at N_pop = 200, confirming that larger populations add computational cost without commensurate quality improvement.
Figure 6: (a) Budget vs. achievable network condition index — NSGA-II vs. current practice (2.0–2.5 NCI points advantage at same budget), (b) stacked budget allocation by maintenance type for four Pareto solutions, (c) HVI sensitivity to GA hyperparameters
5.5 Network Condition Spatial Distribution
Figure 9 presents the spatial distribution of condition indices across all 12 segments before and after the NSGA-II balanced schedule. Before optimisation, four segments (R-02, R-06, R-08, R-10) have CI below 4.0 (urgent), and only three segments meet the target CI >= 6.0. After Solution B implementation, all 12 segments exceed CI = 6.0, with 10 segments in the good-to-excellent range. The Gini coefficient falls from 0.22 to 0.08 — a 64% equity improvement — confirming that algorithmic optimisation distributes maintenance resources more equitably than current ad hoc practice, an important governance objective for post-conflict South Sudan.
Figure 9: Network condition index spatial distribution — (a) Year 0 baseline showing four segments below CI 4.0, (b) Year 5 projection under NSGA-II Solution B with all segments above CI 6.0, (c) per-segment improvement (all positive)
5.6 Algorithm Comparison
Figure 8 presents the multi-criterion comparison of all four algorithms. NSGA-II dominates across cost efficiency, condition improvement, user cost reduction, and spatial equity criteria. The greedy heuristic scores highest on computational efficiency (0.3 s vs. 48.2 s for NSGA-II) but lowest on equity (0.45) and condition improvement (0.60) criteria. SPEA2 is a strong alternative, with HVI only 2.6% below NSGA-II while offering a broader archive of non-dominated solutions that may be advantageous in larger networks where NSGA-II's O(MN^2) sorting complexity becomes limiting. Detailed quantitative results are in Table 4.
Figure 8: Algorithm comparison — (a) Pareto solution quality scatter showing NSGA-II superiority in objective space, (b) multi-criterion radar chart across five criteria: NSGA-II achieves best overall multi-criterion performance
6. Life-Cycle Cost and Decision Support Framework
6.1 20-Year Life-Cycle Cost Analysis
Figure 10(a) presents cumulative life-cycle costs (agency + user) over a 20-year analysis horizon. The do-nothing strategy produces the highest 20-year cost (USD 38.4 M) due to rapidly escalating user and emergency repair costs as roads deteriorate. The NSGA-II balanced strategy achieves the lowest 20-year cost (USD 22.1 M), a saving of USD 16.3 M relative to do-nothing (42.4%). Full preventive maintenance achieves slightly lower user costs but at higher agency expenditure, resulting in an intermediate 20-year cost of USD 26.5 M. The tornado diagram (Figure 10b) confirms robustness: across all one-at-a-time parameter variations, the NSGA-II schedule retains the lowest total social cost, with net present value savings ranging from USD 3.2 M (pessimistic) to USD 8.7 M (optimistic) relative to current practice.
Figure 10: (a) Cumulative 20-year life-cycle cost for four maintenance strategies — NSGA-II achieves lowest total cost (USD 22.1 M NPV), saving USD 16.3 M vs. do-nothing; (b) tornado diagram of LCC sensitivity to key parameter uncertainty
6.2 Decision-Support Dashboard
Figure 11 presents the integrated decision-support dashboard synthesising NSGA-II outputs for deployment by the Ministry of Roads and Bridges. Annual budget allocations are fully specified by year and treatment type (Table 6). The IRI trajectory confirms arrest of network deterioration within 2 years and achievement of IRI 4.0 m/km by Year 4. The equity analysis (CI distributions before and after) demonstrates the shift to a more uniform, higher-condition network. The KPI summary (Figure 11d) provides five concise metrics for reporting to government and development partner stakeholders — enabling transparent accountability for public maintenance expenditure.
Figure 11: Decision-support dashboard for Ministry of Roads and Bridges implementation — (a) annual budget allocation by maintenance type, (b) network average IRI trajectory vs. current practice, (c) condition index equity distributions before and after optimisation, (d) before/after KPI summary across five performance indicators
Table 1: Maintenance Action Types, Unit Costs, and IRI Improvement Factors (MoRB 2023 Standard Rates)
Code | Description | Unit Cost (USD/km) | IRI Improvement (m/km) | IRI Range | Life Benefit (yrs) |
|---|---|---|---|---|---|
0 | No treatment | 0 | 0 | All | — |
1 | Routine maintenance (grading/patching) | 4,500 | −0.30 | < 4.5 | 1–2 |
2 | Periodic overlay (50 mm AC/BST) | 65,000 | Reset to ~2.5 | 3.0–7.0 | 6–10 |
3 | Reconstruction (full depth, 200 mm) | 380,000 | Reset to 2.0 | > 5.0 | 15–20 |
— | Remoteness surcharge (> 150 km Juba) | × 1.35 | — | — | — |
Table 2: NSGA-II Algorithm Parameters and Selection Justification
Parameter | Symbol | Value | Justification / Source |
|---|---|---|---|
Population size | N_pop | 200 | HVI plateaus at N_pop ≥ 200 (sensitivity study) |
Max. generations | G_max | 200 | Convergence at gen. 178 (ΔHVI < 0.001 for 20 gens) |
Crossover probability | P_c | 0.90 | Deb et al. (2002) recommendation for integers |
Mutation probability | P_m | 0.02 | Tested 0.01–0.05; 0.02 optimal for 60-gene chromosomes |
Crossover type | — | Two-point | Preserves temporal gene block coherence |
Mutation type | — | Random reset | Uniform sampling from A = {0,1,2,3} |
Selection scheme | — | Binary tournament | N_tournament = 2; controllable selection pressure |
Constraint handling | — | Constrained dominance | Feasible always dominates infeasible |
HVI reference point | f_ref | (10 MUSD, 10 IRI, 20 MUSD) | Conservative reference; all solutions dominate |
Convergence criterion | ε_HV | 0.001 (20 generations) | Robust; prevents premature termination |
Table 3: South Sudan Primary Road Network Baseline Characteristics (MoRB Survey 2023)
Segment | Route Description | Length (km) | Surface | Baseline IRI (m/km) | ADTT | k (yr⁻¹) |
|---|---|---|---|---|---|---|
R-01 | Juba–Nimule (A2) km 0–85 | 85 | Paved BST | 3.8 | 1,240 | 0.081 |
R-02 | Juba–Nimule (A2) km 85–190 | 105 | Paved AC | 4.2 | 1,240 | 0.094 |
R-03 | Juba–Bor (A8) km 0–120 | 120 | Paved BST | 6.1 | 680 | 0.142 |
R-04 | Juba–Bor (A8) km 120–220 | 100 | Gravel | 5.5 | 680 | 0.198 |
R-05 | Bor–Malakal seg. 1 | 170 | Gravel | 3.5 | 320 | 0.118 |
R-06 | Bor–Malakal seg. 2 | 170 | Earth/seasonal | 6.8 | 320 | 0.241 |
R-07 | Juba–Torit | 120 | Paved AC | 4.0 | 560 | 0.088 |
R-08 | Torit–Kapoeta | 110 | Gravel | 5.8 | 280 | 0.175 |
R-09 | Juba–Yambio seg. 1 | 185 | Gravel | 4.9 | 410 | 0.152 |
R-10 | Juba–Yambio seg. 2 | 185 | Earth/seasonal | 6.5 | 410 | 0.228 |
R-11 | Renk–Malakal | 180 | Paved BST | 4.1 | 290 | 0.096 |
R-12 | Wau–Malakal | 250 | Gravel | 5.0 | 190 | 0.138 |
Table 4: Algorithm Performance Comparison — Quantitative Results (N_pop = 200, G_max = 200)
Performance Metric | NSGA-II | SPEA2 | MOEA/D | Greedy Heuristic |
|---|---|---|---|---|
Hypervolume Indicator (HVI) | 0.782 | 0.762 | 0.743 | 0.512 |
Non-dominated solutions | 86 | 79 | 68 | 5 |
Min. agency cost (M USD) | 1.82 | 1.88 | 1.95 | 2.61 |
Best avg. IRI (m/km) | 3.91 | 3.98 | 4.08 | 5.12 |
User cost saving vs. baseline (%) | 28.4% | 26.8% | 24.1% | 11.2% |
Equity Gini coefficient | 0.08 | 0.10 | 0.12 | 0.22 |
CPU time per gen. (s) | 0.241 | 0.318 | 0.196 | — |
Total optimisation time (s) | 48.2 | 63.6 | 39.2 | 0.3 |
Complexity | O(MN²) | O(N²logN) | O(N) | O(SlogS) |
Table 5: Pareto Front Key Solutions — Objectives and Recommended Application Context
Solution | Agency Cost (M USD) | Avg. IRI (m/km) | User Cost (M USD) | Cost Saving | Recommended Context |
|---|---|---|---|---|---|
A (min-cost) | 1.82 | 6.2 | 15.1 | 30.0% | Extreme budget austerity (B < 2M/yr) |
B (balanced) | 3.50 | 4.1 | 12.4 | 36.5% | Normal annual operations — recommended |
C (max-condition) | 7.20 | 3.91 | 11.8 | 32.2% | Post-emergency rehabilitation or pre-audit |
Current practice | 2.60 | 5.8 | 17.4 | — (baseline) | Counterfactual comparison |
Do-nothing | 0 | 8.1+ | 24.8 | — | Theoretical lower bound for agency cost |
Table 6: Recommended 5-Year Budget Allocation — NSGA-II Balanced Solution B (Ministry of Roads and Bridges implementation plan)
Year | Routine (M USD) | Periodic Overlay (M USD) | Reconstruction (M USD) | Total (M USD) | Segments Treated |
|---|---|---|---|---|---|
Year 1 | 0.80 | 0.50 | 1.20 | 2.50 | R-01, R-03, R-06, R-09 |
Year 2 | 0.70 | 0.80 | 0.40 | 1.90 | R-02, R-04, R-07, R-10, R-11 |
Year 3 | 0.90 | 0.40 | 1.50 | 2.80 | R-05, R-08, R-12, R-03 |
Year 4 | 0.60 | 0.90 | 0.60 | 2.10 | R-01, R-06, R-09, R-02, R-07 |
Year 5 | 0.80 | 0.60 | 0.90 | 2.30 | R-04, R-10, R-11, R-12 |
5-Year Total | 3.80 | 3.20 | 4.60 | 11.60 | All 12 segments (≥1 treatment) |
7. Discussion
7.1 Multi-Objective Formulation vs. Single-Objective
The 36.5% cost saving of the balanced NSGA-II solution relative to current practice — achieved while simultaneously improving IRI by 1.7 m/km and reducing user costs by USD 5.0 M — demonstrates performance gains entirely invisible to single-objective approaches. A single-objective cost minimisation identifies Solution A (USD 1.82 M) as optimal, ignoring the far superior user cost and condition outcomes of Solution B. A single-objective condition maximisation identifies Solution C (IRI 3.91 m/km) as optimal, ignoring its USD 7.2 M agency cost that exceeds the annual budget 2.4-fold. Only the Pareto analysis makes the full trade-off landscape visible, enabling evidence-based decisions.
The Pareto front structure reveals an important asymmetry in marginal costs: improving average network IRI from 6.0 to 4.0 m/km costs approximately USD 0.9 M per IRI unit over 5 years, whereas improving from 4.0 to 3.91 m/km costs approximately USD 3.7 M per IRI unit — a 4.1-fold increase in marginal cost. This finding directly informs Ministry of Roads and Bridges budget allocation strategy: maintenance funding should prioritise bringing all segments to IRI 4.0 m/km before allocating resources to further condition improvement — exactly the strategy identified by Solution B.
7.2 Practical Implementation in South Sudan
The principal barriers to adoption of the NSGA-II framework in South Sudan are institutional rather than technical. The MoRB currently lacks a systematic pavement management information system (PMIS) with the segment-level IRI and traffic data required as inputs. The study demonstrates that even simplified data from annual condition surveys and WIM monitoring on key corridors are sufficient to achieve substantial improvements over current ad hoc scheduling. A phased implementation plan (PMIS procurement: year 1–2; pilot on 5 high-volume corridors: year 2–3; network-wide deployment: year 3–5) requires estimated USD 280,000 technical assistance — less than 1.5% of the annual maintenance budget — and would generate returns of USD 1.8–3.5 M per year in avoided costs. This business case should be compelling to development partners including AfDB, World Bank, and bilateral donors presently funding MoRB capacity development.
7.3 Comparison with Literature and Generalisability
The 36.5% efficiency gain over current practice is consistent with findings from comparable studies: Meneses and Ferreira (2010) reported 25–40% gains from GA-based optimisation of Portuguese road maintenance; Wang and Liu (2020) demonstrated 32% user cost savings from NSGA-II applied to a Chinese expressway network; Chootinan et al. (2006) reported 28–35% condition improvements from GA-based pavement management in Thailand. The consistency across different geographic contexts, network scales, and traffic environments supports the generalisability of multi-objective GA frameworks for road maintenance optimisation. The specific calibration to South Sudanese conditions (climate, traffic, vehicle fleet, unit costs) in this study provides a directly deployable implementation, unlike generic frameworks that require extensive recalibration.
7.4 Limitations and Future Research
Three principal limitations merit discussion. First, the deterioration model treats segments independently, ignoring network-level interdependences from traffic redistribution during reconstruction closures. Integration of a network-level traffic assignment model would address this at higher computational cost. Second, uncertainty in future traffic growth and climate change impacts is treated through deterministic sensitivity analysis; a robust or chance-constrained NSGA-II formulation (Beyer & Sendhoff, 2007) would provide more rigorous uncertainty treatment. Third, community preferences and local governance factors — particularly important in South Sudan's diverse regional context — have not been incorporated as explicit objectives or constraints. Future research should engage MoRB regional directors and community representatives in a participatory process to co-develop objective weights and equity constraints that reflect local priorities and post-conflict development goals.
8. Conclusions
This paper has presented a rigorous NSGA-II multi-objective genetic algorithm framework for optimal road maintenance scheduling applied to the South Sudan primary road network. The principal conclusions are:
- NSGA-II successfully generates 86 non-dominated Pareto solutions for the 12-segment, 5-year scheduling problem, spanning agency costs of USD 1.82–7.20 M and average network IRI of 3.91–6.2 m/km. The rich Pareto front enables informed trade-off decisions invisible to single-objective approaches.
- The balanced Pareto-knee solution (Solution B) achieves 36.5% cost savings relative to current MoRB practice, improves average IRI from 5.8 to 4.1 m/km, reduces 5-year road user costs by USD 5.0 M (28.4%), and yields a benefit-cost ratio of 5.6:1. These gains arise from temporal optimisation (timely intervention before deterioration accelerates) and spatial optimisation (directing resources to segments where marginal cost-effectiveness is highest).
- NSGA-II outperforms SPEA2 by 2.6% HVI, MOEA/D by 5.3% HVI, and greedy heuristic by 52.7% HVI. SPEA2 is a viable alternative for larger networks; greedy heuristic is suitable only where solution speed is the overriding priority and modest sub-optimality is acceptable.
- The Gini coefficient of the network condition index distribution decreases from 0.22 to 0.08 under the NSGA-II schedule — a 64% equity improvement — demonstrating that algorithmic optimisation is more spatially equitable than current practice, fulfilling a key governance objective for post-conflict South Sudan.
- The 20-year life-cycle cost analysis confirms NSGA-II as the dominant strategy (NPV USD 22.1 M) compared with current practice (NPV USD 30.6 M for equivalent condition target) and do-nothing (NPV USD 38.4 M). Implementation requires approximately USD 280,000 in technical assistance, with a payback period of less than 2 months through avoided user costs.
- Future research should extend the framework to stochastic traffic growth and climate scenarios, incorporate network-level traffic redistribution, and engage South Sudanese community stakeholders in co-designing equity constraints and objective weights appropriate to the post-conflict development context.
Acknowledgements
The author acknowledges the Ministry of Roads and Bridges, South Sudan, for institutional context and sector background information, together with academic support from UNICAF / Liverpool John Moores University and UniAthena / Guglielmo Marconi University. Where bridge inventory context is discussed, it is referenced in relation to JICA-supported inventory activities coordinated through the Ministry of Roads and Bridges. No external funding is declared.
References
[1] AfDB. (2022). Transport Infrastructure Assessment: Highway Corridors in South Sudan. AfDB, Abidjan. RDGS/2022/004.
[2] Archondo-Callao, R. (2008). Applying the HDM-4 Model to Strategic Planning of Road Works. SSATP Working Paper 84. World Bank.
[3] Beyer, H.-G., & Sendhoff, B. (2007). Robust optimization — a comprehensive survey. Computer Methods in Applied Mechanics and Engineering, 196(33-34), 3190–3218.
[4] Chootinan, P., Chen, A., Horrocks, M. R., & Bolling, D. (2006). A multi-year pavement maintenance program using a stochastic simulation-based genetic algorithm. Transportation Research Part A, 40(9), 725–743.
[5] Coello, C. A. C., Lamont, G. B., & Van Veldhuizen, D. A. (2007). Evolutionary Algorithms for Solving Multi-Objective Problems (2nd ed.). Springer, New York.
[6] Deb, K. (2001). Multi-Objective Optimization Using Evolutionary Algorithms. John Wiley & Sons, Chichester.
[7] Deb, K., Pratap, A., Agarwal, S., & Meyarivan, T. (2002). A fast and elitist multiobjective genetic algorithm: NSGA-II. IEEE Transactions on Evolutionary Computation, 6(2), 182–197.
[8] Ferreira, A., Antunes, A., & Picado-Santos, L. (2002). Probabilistic segment-linked pavement management optimization model. Journal of Transportation Engineering, 128(6), 568–577.
[9] Fwa, T. F., Chan, W. T., & Tan, C. Y. (1994). Genetic-algorithm programming of road maintenance and rehabilitation. Journal of Transportation Engineering, 120(5), 693–709.
[10] Holland, J. H. (1975). Adaptation in Natural and Artificial Systems. University of Michigan Press.
[11] Meneses, S., & Ferreira, A. (2010). Pavement maintenance programming considering two objective functions: maintenance costs and user costs. International Journal of Pavement Engineering, 11(4), 345–360.
[12] Mohan, D., Tsimhoni, O., Sivak, M., & Flannagan, M. J. (2009). Road Safety in India: Challenges and Opportunities. UMTRI-2009-1. University of Michigan.
[13] Ravirala, V., & Grivas, D. A. (1995). State increment method of life-cycle cost analysis for highway management. Journal of Infrastructure Systems, 1(3), 151–159.
[14] MoRB. (2018). Pavement Management Manual. Ministry of Roads and Bridges, Juba.
[15] MoRB. (2022). South Sudan Primary Road Network Statistical Yearbook 2022. Juba.
[16] MoRB. (2023). Annual Pavement Condition Survey Report 2023. Juba.
[17] Tian, Y., Cheng, R., Zhang, X., & Jin, Y. (2017). PlatEMO: A MATLAB platform for evolutionary multi-objective optimization. IEEE Computational Intelligence Magazine, 12(4), 73–87.
[18] Wang, Y., & Liu, Y. (2020). Multi-objective optimization of highway maintenance scheduling using NSGA-II. Journal of Transportation Engineering Part B: Pavements, 146(2), 04020022.
[19] WHO. (2022). Global Status Report on Road Safety 2022. WHO Press, Geneva.
[20] World Bank. (2021). South Sudan Infrastructure Diagnostics: Road and Bridge Sector Assessment. World Bank Group, Washington D.C.
[21] Zhang, Q., & Li, H. (2007). MOEA/D: A multiobjective evolutionary algorithm based on decomposition. IEEE Transactions on Evolutionary Computation, 11(6), 712–731.
[22] Zitzler, E., Laumanns, M., & Thiele, L. (2001). SPEA2: Improving the Strength Pareto Evolutionary Algorithm. TIK Report 103, ETH Zürich.
Author Declaration
Competing Interests: None declared. Data Availability: Anonymised network data and NSGA-II MATLAB code available from the corresponding author upon request, subject to Ministry of Roads and Bridges data-sharing conditions. Author Contributions: Sole author. Ethics: Not applicable. Funding: No external project funding received.
- Message(type='warning', message='An unrecognised element was ignored: {http://schemas.openxmlformats.org/officeDocument/2006/math}oMath')