Pavement Maintenance Decision Making Based on Optimization Models

: Pavement maintenance prioritization considering both quality and cost is an important decision-making problem. In this paper, the actual pavement condition index of city roads was calculated using municipal patrol data. A linear optimization model that maximized maintenance quality with limited maintenance costs and a multi-objective optimization model that maximized maintenance quality while minimizing maintenance costs were developed based on the pavement condition index. These models were subsequently employed in making decisions for actual pavement maintenance using sequential quadratic programming and a genetic algorithm. The results showed that the proposed decision-making models could effectively address actual pavement maintenance issues. Additionally, the results of the single-objective linear optimization model veriﬁed that the multiobjective optimization model was accurate. Thus, they could provide optimal pavement maintenance schemes for roads according to actual pavement conditions. The reliability of the models was investigated by analyzing their assumptions and validating their optimization results. Furthermore, their applicability in pavement operation-related decision making and preventive maintenance for roads of different grades was conﬁrmed.


Introduction
Continuous urbanization in China has significantly increased the cumulative mileage of roads, which are considered the arteries of economic and social activities at the national and local levels. As an essential component of road infrastructure, a pavement is a type of hard surface made from durable surface material, which is able to withstand traffic and harmful environments. Because of increasing traffic volumes with heavy loads and the impacts of adverse environments, regular pavement maintenance is necessary to repair damage and mitigate degradation. A sufficient budget should be allocated to maintain the pavement at an appropriate condition, but insufficient budget is the primary obstacle of pavement maintenance [1]. When conducting large-scale pavement maintenance, it is important to consider the maintenance costs while ensuring that maintenance quality meets the requisite standards. The selection of a proper pavement maintenance scheme is a multiobjective optimization task. Therefore, determining the optimal maintenance scheme is crucial for decision making in pavement maintenance [2,3]. Decision making in pavement maintenance is a complicated, nondeterministic polynomial time (NP)-hard problem [4].
Maintenance personnel generally make judgments based on their personal experience when facing pavement maintenance tasks. Typical considerations include determination 2 of 15 of which roads need maintenance, measures to be adopted for pavement maintenance, estimation of the maintenance cost, and the deadline for completing the maintenance. However, decision making based on human experience cannot control pavement maintenance costs or evaluate whether the maintenance quality meets the relevant standards. This can subsequently lead to deterioration of road structures or pose challenges in the assessment of maintenance outcomes. Therefore, assessment and prediction of road deterioration [5] as well as establishment of the relationship between maintenance schemes and corresponding influencing factors [6] are crucial in decision making regarding pavement maintenance and fund allocation.
Recently, various decision-making methods for pavement maintenance have been proposed. These methods generally aim at improving quality and reducing costs under diverse scenarios. The most common methods include the analytic hierarchy process (AHP), decision making based on a clustering algorithm, and decision making based on a genetic algorithm (GA). The AHP method was first used for decision making in pavement maintenance by Saaty [7]. Farhan and Fwa proposed a pavement maintenance prioritization method based on AHP and developed a three-level hierarchy model to determine the priority ranking of different pavement maintenance schemes [8]. Li et al. developed an AHP-based pavement maintenance priority model by assigning weights to the factors influencing pavement conditions [9]. However, the selection of the indices and the determination of the index weights in the AHP method remained subjective. Thus, the method of selecting an optimal solution from all available plans determines the effectiveness of the final decision.
In recent years, with the development of machine learning and corresponding algorithms, some researchers began studying pavement maintenance decisions based on machine learning. Han et al. proposed a decision-making framework for pavement maintenance that combined a clustering algorithm with the PageRank algorithm [10]. In this framework, the pavements are first grouped into clusters by the clustering algorithm. Then, the maintenance ranking of the clustered road samples is determined according to the road conditions. Hafez et al. proposed a decision-making algorithm for pavement maintenance based on pattern recognition [11]. This method helps determine the optimal maintenance and repair plans for low-volume paved roads. GAs and traditional mathematical programming are often used to solve road M&R planning problems based on multiobjective optimization (MOO), but they also have their limitations. Elhadidy et al. [12] proposed a two-objective optimization model balancing minimum action costs and maximum conditions for used road networks and developed a GA-based procedure for solving the MOO problem. Hadiwardoyo et al. [13] described the development of a genetic algorithm based on multiobjective programming for pavement and investigated the optimal maintenance strategy options applied as function of road surface distress conditions with a lack of information related to monitoring data and evaluation. Sindi and Agbelie [14] explored the expected accuracy rates of network treatment options through a multiobjective optimization methodology that utilized a GA and mixed-integer programming; their method was capable of effectively assigning pavement maintenance tasks under certain conditions. Hafez et al. [15] performed large-scale optimization to compare the existing maintenance policy with an alternative strategy. Specifically, they incorporated a GA into the optimization model to address the issues associated with specific optimization constraints and the limitations related to low-volume roads during the optimization process. Alqaili et al. offered a new multiobjective stochastic algorithm called the integer search algorithm (ISA) [1]. The ISA and GA were applied to improve the performance condition rating (PCR) of pavement in developing countries and achieved this by maximizing the condition of the pavement with minimum costs at specified constraints. Santos et al. [3] proposed a novel adaptive hybrid GA that incorporated local search techniques to improve the overall efficiency and effectiveness of the search. Hosseininasab et al. [16] proposed two multiobjective evolutionary approaches to solve problems of road construction within a reasonable time. These approaches were combinations of different techniques, such as GA, NSGA-II, the Frank-Wolfe algorithm, and the ordered logit model. The NSGA-II and its improved algorithms are popular for solving MOO problems [17,18]. To address the high computation complexity associated with pavement maintenance at the network level, Hankach et al. [19] developed a model to reduce the search space and formulated the original problem as a generalized assignment problem, which was a well-known problem in mathematical optimization. Ahmed et al. [20] proposed a chaotic particle swarm optimization algorithm to find the optimal solution for pavement maintenance. This algorithm could effectively resolve the maintenance and rehabilitation issues associated with flexible pavements. Ameri and Jarrahi [21] used condition indicators in the form of normalized values and developed technical constraints in a linear integer programming model to improve network-level pavement maintenance and rehabilitation planning. Chen et al. [22] proposed a network-level pavement maintenance and rehabilitation optimization model considering the costs of user travel time and vehicle fuel consumption. The model optimized the asphalt pavement performance evaluation method, including 11 different combinations, which could be easily extended to the study of more complex road networks considering other factors. Mataei et al. [23] proposed a model based on the cloud decision tree (CDT) theory, which included a general decision-making model and various decision trees for every province of the country.
As evident from the foregoing discussion, considerable research has been conducted to improve GAs and to utilize them for decision making in pavement maintenance because of their effectiveness in resolving multiobjective optimization problems. In this study, the pavement maintenance problem was introduced and formulated as a multiobjective decision-making problem. Subsequently, two optimization models were developed based on actual road conditions and applied to a real-world case. Sequential quadratic programming and the NSGA-II algorithm were used to solve the two models. Finally, the characteristics and the optimization space of the proposed models were analyzed based on the results obtained. Future plans for extending the proposed model were discussed in depth.

Problem Formulation
Pavement maintenance usually involves decision making at the macro level for selecting roads to be maintained under objective constraints and decision makers' priorities. In China, pavement maintenance is usually subject to territorial management within each administrative region. For example, each district and county has its own department responsible for decision making during pavement maintenance at the macro level. Typical decision-making tasks involved in pavement maintenance include the formulation of an annual pavement maintenance scheme, budget application, and quality control. When making pavement maintenance decisions in real scenarios, the degree of deterioration of the road network in a certain region is first evaluated based on daily patrol data. Pavement deterioration can be quantified as the pavement condition index (PCI), which is a numerical rating of pavement condition based on the type and severity of distresses observed on the pavement surface. PCI is represented by a numerical value between 0-100, where 0 is the worst condition and 100 is the best.
Subsequently, under the premise of considering the total maintenance cost and time, the annual maintenance scheme is formulated to improve the condition of roads with low PCIs. Thus, the maintenance of regional road networks must take multiple objectives into account, such as minimum maintenance cost and optimum maintenance quality. Therefore, determining a suitable approach to formulating an efficient maintenance scheme is a decision-making challenge.
The multiobjective maintenance problem of a road network in a specific region can be described mathematically as follows. Given the number of roads (N), the area of each road (A) that requires maintenance, the PCI of each road, and the total budget available for pavement maintenance (C) in a certain region, m roads (m ≤ N) are prioritized for maintenance among the total N roads, such that (1) the total maintenance cost is less than or equal to C, (2) the total maintenance cost is minimized, and (3) the maintenance quality is maximized. Selecting the appropriate roads based on human experience results in 2 N combination schemes, because the decision variable is a Boolean representing whether a road is included in the maintenance schedule. Therefore, N cannot be excessively large, for when N is large, it is infeasible to select the final roads manually. This problem is NP-hard. To resolve it, heuristic algorithms were researched and developed aiming at an approximate optimal solution that could be a particularly sensible choice. In this study, a decision-making optimization model for pavement maintenance was solved by using NSGA-II. Optimum maintenance schemes were obtained that could ensure the maximal maintenance quality with minimal maintenance cost. Figure 1 shows the overall procedure of decision-making optimization for pavement maintenance. The foundation is to evaluate road condition by pavement condition index, and the key is to establish optimization models of road maintenance. The models can be applied to provide optimal schemes for multiobjective decision making in regional road network maintenance. Many optimization models have been widely used, as found from engineering management and research literature, but these models of optimizing pavement maintenance have not been fully examined. In this paper, we examined how these models could be extended for our research goals and demonstrated how they could be used in pavement maintenance for an application area. for pavement maintenance (C) in a certain region, m roads (m ≤ N) are prioritized for maintenance among the total N roads, such that (1) the total maintenance cost is less than or equal to C, (2) the total maintenance cost is minimized, and (3) the maintenance quality is maximized. Selecting the appropriate roads based on human experience results in 2 N combination schemes, because the decision variable is a Boolean representing whether a road is included in the maintenance schedule. Therefore, N cannot be excessively large, for when N is large, it is infeasible to select the final roads manually. This problem is NPhard. To resolve it, heuristic algorithms were researched and developed aiming at an approximate optimal solution that could be a particularly sensible choice. In this study, a decision-making optimization model for pavement maintenance was solved by using NSGA-II. Optimum maintenance schemes were obtained that could ensure the maximal maintenance quality with minimal maintenance cost. Figure 1 shows the overall procedure of decision-making optimization for pavement maintenance. The foundation is to evaluate road condition by pavement condition index, and the key is to establish optimization models of road maintenance. The models can be applied to provide optimal schemes for multiobjective decision making in regional road network maintenance. Many optimization models have been widely used, as found from engineering management and research literature, but these models of optimizing pavement maintenance have not been fully examined. In this paper, we examined how these models could be extended for our research goals and demonstrated how they could be used in pavement maintenance for an application area.

Model Construction
The construction of the optimization models involved multiple steps, such as selecting the variables for decision making, setting up the model parameters, determining the objective function, and establishing the constraint conditions. In this section, two optimization models are established to resolve the road network maintenance problem.

Model Construction
The construction of the optimization models involved multiple steps, such as selecting the variables for decision making, setting up the model parameters, determining the objective function, and establishing the constraint conditions. In this section, two optimization models are established to resolve the road network maintenance problem.

Maximum Maintenance Quality-Limited Budget Model
To achieve optimal maintenance quality under a limited maintenance budget, a maximum maintenance quality-limited budget (MMQLB) model is built. The MMQLB model is a single-objective linear optimization model and is expressed as follows.
In these equations, x i denotes the decision-making variable of the model, which represents whether the ith road requires maintenance. For each road, there exist two states in terms of whether maintenance is required. Setting x i to 1 or 0 indicates that the road does or does not require maintenance at present, respectively.
Equation (1) represents the objective function of this model, where n is the total number of roads and PCI i − PCI i represents the level of improvement in PCI (i.e., the maintenance quality) for a specific road. PCI i indicates the PCI of the road after maintenance and is set to 100 in this model; in other words, the road is assumed to be in perfect condition after maintenance. In contrast, PCI i represents the original PCI of the ith road before maintenance. Equation (2) represents the constraint condition of the model, where Area i represents the area (m 2 ) of the ith road that requires maintenance. P is a constant that represents the maintenance cost per square meter of the road and is set as 200 CNY/m 2 . Budget represents the total budget available for pavement maintenance, and D is a constant that represents the maintenance threshold of the road. When the PCI of a specific road drops below this threshold value, the road is judged to be of poor quality and require maintenance. Here, the default value of D is set to 70.
PCI is selected as the parameter in the objective function because it is an important index for evaluating whether the road is in a satisfactory condition according to the Chinese "technical code of maintenance for urban road" [24]. In addition, PCI plays an important role in finalizing the maintenance scheme [25]. In this study, asphalt pavement is selected as the maintenance object. The equation for calculating the PCI of asphalt pavement is as follows: In this equation, PCI ranges from 0 to 100; n is the total number of main deterioration types (here set as 4 for asphalt pavement exhibiting cracks, deformation, loosening, and other forms of deterioration); m is the total number of deterioration subtypes included in each main deterioration type; and DP ij is the deduction caused by the jth subtype of deterioration in the ith main type. ω ij is the weight of the jth subtype of deterioration in the ith main type, which can be calculated by Equations (4) and (5) as follows: When computing the DP ij value, first, the actual deterioration density of a pavement, the ratio of the total area of the subtype j deterioration to that of the pavement, is calculated. The actual deterioration density is compared with the deterioration density of each deterioration subtype in Table 1, and DP ij is found to be proportional to two numbers in two adjacent columns in Table 1. The u ij is calculated by Equation (5) and w ij by Equation (4) based on the DP ij in Table 2. For example, the PCI of AnQing Rd. in the research area is calculated as shown in Table 2. The paved area of AnQing road is about 26,204.1 m 2 , with 1746.94 m of length and 4 lanes. The deterioration density is the ratio of the deteriorated area of each subtype and the paved area of AnQing road. The final PCI of AnQing road is 68.63, subtracting the sum of multiple of DP ij and ω ij from 100.  As mentioned previously, in this model, the road is assumed to be in perfect condition after maintenance. Therefore, the maintenance quality can be represented by 100 − PCI i . In other words, the larger the value of 100 − PCI i , the higher the maintenance quality.

Minimum Budget-Maximum Maintenance Quality Model
Optimal maintenance of multiple roads should seek to minimize the maintenance cost while maximizing the maintenance quality. To this end, a multiobjective optimization model named the minimum budget-maximum maintenance quality (MBMMQ) model is constructed and expressed by the following equations: In Equations (6) and (7), the function Z 1 represents the maintenance quality, which equals the summation of improved PCI for all involved roads, whereas Z 2 represents the total maintenance cost. In this model, PCI i is set as 100 to indicate that roads are assumed to be in perfect condition after maintenance. PCI i denotes the original PCI of the ith road, and x i is the state variable, which indicates whether the ith road requires maintenance. The value of x i can be set only as 1 or 0, which signifies that the road does or does not require maintenance, respectively. Finally, Area i and P have the same definitions as those in Equation (2). In Equation (8), D is the maintenance threshold; when the PCI of a specific road drops below threshold value, it is compulsory to maintain that road. Here, D is set to 70.

Model Application
This section elaborates on the application of the proposed models in decision making for a real-world pavement maintenance problem. Using actual data, the global optimal results were obtained from the model.

Dataset
In total, 149 roads under the jurisdiction of the Shushan District, a county-level district in Hefei, Anhui (China), were taken as the research objects. The deterioration data of these roads in 2019 were collected, and the maintenance area and premaintenance PCI were calculated by Equation (3). The roads were divided into three grades based on the National Standard of Road Classification: 9 expressways, 39 main roads, and 101 branch roads. For brevity, only a few roads in the three grades and their premaintenance PCI are listed in Table 3. Among the 149 roads, the average, minimum, and maximum PCI were 88.43, 67.89, and 100, respectively. There were 5 roads with PCI less than 70 and 26 roads with PCI equaling 100. As can be observed from the data, some roads were in poor condition, while others were in satisfactory condition. This reflects the level of diversity of the road dataset chosen for the study.

Solving the Models
Prior to solving the models, the premaintenance PCI, maintenance area, maintenance cost per unit area, and maintenance threshold for each road were first fed into the model as input parameters. Subsequently, the optimization models were implemented, and the values of the decision-making variables and the objective functions were calculated by using sequential quadratic programming and the NSGA-II algorithm. Finally, the optimal pavement maintenance scheme was determined based on the results obtained from the model. The first model was solved using the sequential quadratic programming method.
Because the multiobjective model could be solved by a GA, the second model was solved by using the NSGA-II algorithm in our research. The calculation steps involved in the NSGA-II algorithm are as follows. First, multiple sets of optimization solutions for pavement maintenance are randomly generated and are referred to as the solution population. Subsequently, these solution sets are ranked based on the maximum pavement maintenance quality and the minimum pavement maintenance cost calculated for each set of solutions. This process is also known as rapid nondominated sorting of the solution population. Next, a new generation of solutions is obtained through selection, crossover, and mutation of the solution population. The new solutions are the offspring of the first-generation solutions. Starting from the second generation, the parent and offspring populations are combined to perform rapid nondominated sorting. Simultaneously, the crowding degrees are calculated for each individual in the nondominated layer. Suitable individuals are then selected according to their crowdedness and the nondominated relationship to form a new parent population. Finally, a new offspring population is generated through selection, crossover, and mutation of the new parent population.
The final solution set of the objective function can be obtained by repeating the aforementioned process continuously until the maximum number of iterations is reached. NSGA-II can select multiple nondominated individuals that form an optimal set of trade-off solutions called the Pareto set. The pavement maintenance scheme can then be formulated according to the solution results. Some researchers have proposed simpler coding methods to represent complex evolutionary phenomena. These methods can realize heuristic searches in the complex search space and can determine the global optimal solution of the objective function with high probability through a simplified genetic process. Table 4 shows the optimal solutions obtained by the MMQLB model with different budgets. The best pavement maintenance quality provided by the model was 1724.29, which equaled the sum of the improved PCI for all roads; meanwhile, the maintenance cost was 10.199 million CNY. Because of actual budget limitations, some roads had to be eliminated from the full maintenance scheme. The maintenance quality was 154.51 with the cost = 1,567,200 CNY if only the worst roads under 70 PCI were to be maintained. The optimal maintenance quality was 324.861 with a cos t = 1,999,600 CNY, and the maximum maintenance was 1025.2 under a 5 million CNY budget. The maximum improvement to PCI was 1479.52 under an 8 million CNY budget. The optimal maintenance was 1702.64, and the cost was 9,991,800 CNY under a budget of 10 million CNY. Table 5 shows the PCI improvements to 149 roads in the optimal scheme obtained by MMQLB model under the 5 million CNY budget limitation. According to the optimal scheme, the PCI improvements of each road equaled the difference between 100 and the original PCI if a road was selected to be maintained; otherwise, its PCI improvement was zero.  The optimal maintenance scheme of 1702.64 PCI improvement under a 10 million CNY budget is shown in Figure 2 using a thematic map of 149 roads in the Shushan District. All roads with a PCI smaller than 70 required maintenance according to the constraint set in the MMQLB model. Each road was labelled by its FID and name in the map. For the roads in red, the corresponding decision-making variables were calculated to be 1 under the optimization model, implying that these roads had to be maintained. For the roads in blue, the decision-making variable was calculated to be 0 under the optimization model, signifying that these roads did not require maintenance.

Results of the Models
Next, Figure 3 shows the resulting Pareto front obtained by the MBMMQ model. The objective function was solved by the NSGA-II algorithm and converged to stable values at approximately the 3000th generation. The first objective function was maximizing maintenance quality in the x-axis direction, and the second was minimizing maintenance cost in the y-axis direction. The solution set included many optimal solutions, represented by star points, which made up the Pareto front. Three solutions are highlighted along with their objective function values in Figure 3.
The two optimal models were applied to the integrated management system for municipal facilities in Shushan Dist., Hefei, Anhui (China). Since 2018, the Municipal Engineering Management Office in Shushan Dist. used the system for pavement decision making. The total improved PCI of road maintenance increased by 15% in 2019, while the budget for maintenance remained almost the same. Optimization scheme by maximum maintenance-quality limit budget model . 0 2,000 4,000 1,000 meters

Discussion
In this section, we discuss the suitability of the developed models, the validation of the obtained results, and extensions of the multiobjective model.

Analysis of Model Suitability
The theoretical suitability of the proposed models was analyzed as follows. A linear optimization model must satisfy the following implicit assumptions.
(1) Assumption of linearization: a function is considered linear [26] when the following equation is satisfied: In the equation, and are constants. The objective functions of the two proposed models can be converted to the form shown in Equation (9). Therefore, these objective functions can be considered as linear functions.
(2) Assumptions of proportionality: A change in the decision-making variable should cause the objective function to change proportionally. Because the relationship between the decision-making variable and objective function can be expressed in the form given in Equation (13) for both models, the assumption of proportionality is satisfied for the decision-making variable and the objective function.
(3) Assumption of additivity: The value of the objective function should equal the sum of the individual contributions from each decision-making variable to the objective function. This assumption is satisfied for all models, as is evident from their equations.
(4) Assumption of the fixed parameter: All the parameters are fixed constants, and no random factor is introduced into the analysis. The PCI used in all models is also a fixed parameter. Therefore, this assumption is also satisfied.

Discussion
In this section, we discuss the suitability of the developed models, the validation of the obtained results, and extensions of the multiobjective model.

Analysis of Model Suitability
The theoretical suitability of the proposed models was analyzed as follows. A linear optimization model must satisfy the following implicit assumptions.
(1) Assumption of linearization: a function is considered linear [26] when the following equation is satisfied: In the equation, K and c i are constants. The objective functions of the two proposed models can be converted to the form shown in Equation (9). Therefore, these objective functions can be considered as linear functions.
(2) Assumptions of proportionality: A change in the decision-making variable should cause the objective function to change proportionally. Because the relationship between the decision-making variable and objective function can be expressed in the form given in Equation (13) for both models, the assumption of proportionality is satisfied for the decision-making variable and the objective function.
(3) Assumption of additivity: The value of the objective function should equal the sum of the individual contributions from each decision-making variable to the objective function. This assumption is satisfied for all models, as is evident from their equations.
(4) Assumption of the fixed parameter: All the parameters are fixed constants, and no random factor is introduced into the analysis. The PCI used in all models is also a fixed parameter. Therefore, this assumption is also satisfied.

Validation of Results Obtained from the Models
The proposed MMQLB model is a single-objective linear optimization model that can provide an exact result. Its result is a maintenance scheme that maximizes the maintenance quality under a limited maintenance budget, such as some key schemes shown in Table 4. Meanwhile, the MBMMQ model, which is a biobjective optimization model, was solved by NSGA-II, a heuristic algorithm. In Figure 3, each point represents an optimal maintenance scheme solved by the biobjective optimization model for 149 roads. Among the optimal maintenance schemes suggested by the MBMMQ model, three schemes were highlighted along with their objective function values, and the costs exactly matched the results obtained by the MMQLB model, satisfying the corresponding cost constraints. The best PCI improvement was 1724.29, with a cost of 10,199,000 CNY; the maximum PCI improvement was 1025.2, with a cost of 4,967,400 CNY; and the necessary PCI improvement was 154.51, with a cost of 1,567,200 CNY. Some solutions of MMQLB were included in the solution set of MBMMQ. This suggests that the results of the MMQLB model verified the part results of the MBMMQ model in this case. For this reason, we believe that the solutions found by the NSGA-II were indeed optimal and that they formed the Pareto front for our biobjective maintenance decision-making problem. The overall shape of the Pareto front indicated the same trend of increasing sum of minimal maintenance cost with the total improved quality. The Pareto front for our problem was discrete, but more importantly, it had a nonconvex shape. This was clearly shown when sum of PCI improvement reached around 1000, where the front took an obvious turn. This suggests that it is often difficult to use a heuristic algorithm to obtain an exact convex Pareto front for a multiobjective decision-making problem.

Extension of the Multiobjective Model
This section focuses on extension of the multiobjective model. Considering different road grades, expressways have a higher maintenance priority than main roads. Moreover, main roads usually have a higher maintenance priority than branch roads. The MBMMQ model was extended accordingly, as follows: Here, w i denotes the grade of road i, which is set to 1.5, 1, or 0.5 when road i is an expressway, a main road, or a branch road, respectively. Based on this revised model, the effects of the different road grades and maintenance schemes were further determined, as shown in Figure 4. Compared with the original solution, the maintenance costs and quality are slightly different, which was primarily due to different road PCI values and maintenance areas between expressways and branch roads. These results indicated that the revised model is suitable for multiple road grades.
Because of constraints on pavement maintenance, a road may be only partially maintained, and its PCI may be lower 100 after maintenance. PCI i could replace x i as a decision-making variable, which represents the PCI of the road after maintenance. Based on the MBMMQ model, the third MOO model is proposed as follows.
In this model, PCI i is selected as the variable, which can vary from PCI i to 100. The constants PCI i , Area i , and P are the same as those in Equation (2). Therefore, when calculating the maintenance cost for a specific road, PCI i /PCI i should be incorporated into Equation (13) to consider the effect of PCI improvement on the maintenance area. This model is suitable to maintain quantitatively each selected road and can be solved by NSGA-II.
In this model, ′ is selected as the variable, which can vary from PCIi to 100. The constants PCIi, Areai, and P are the same as those in Equation (2). Therefore, when calculating the maintenance cost for a specific road, ′ / should be incorporated into Equation (13) to consider the effect of PCI improvement on the maintenance area. This model is suitable to maintain quantitatively each selected road and can be solved by NSGA-II.

Conclusions
In multiobjective decision making of pavement M&R, some main research fields have emerged, such as sustainable pavement maintenance management, optimization of decision models with high-dimensional objectives, nonlinear planning, and high-dimensional variables, which have received increasing attention [27]. Here, actual decision-making problems of pavement maintenance were investigated, and an overall procedure for decision-making optimization was proposed. First, PCI was calculated based on the deterioration type and area to evaluate road conditions using the AHP method. Then, two models, named MMQLB and MBMMQ, were developed. The first model was a single-objective linear optimization model, whereas the other was a MOO model. These models were employed for a real-world case involving the maintenance of 149 roads in Shushan District, Hefei, China. First, the pavement condition data of these roads were recorded and imported into the models. The models were then solved using sequential quadratic programming and a GA. Finally, appropriate pavement maintenance schemes were established based on the solutions of the optimization models. Among the optimal maintenance schemes of the MBMMQ model, some optimal values in the pair of maintenance quality and cost exactly matched the results obtained by the MMQLB model with the corresponding cost constraints. This indicated that the MMQLB model could achieve similar decision making to the MBMMQ model by setting different maintenance constraints. According to one nonlinear constraint of MMQLB, it could be improved to substitute the nonlinear formulation and reduce the difficulty of solving the problem in future research. In practical scenarios, the pavement maintenance cost is also dependent on the time taken and priority levels of the road grades. Therefore, the parameters used in the proposed

Conclusions
In multiobjective decision making of pavement M&R, some main research fields have emerged, such as sustainable pavement maintenance management, optimization of decision models with high-dimensional objectives, nonlinear planning, and high-dimensional variables, which have received increasing attention [27]. Here, actual decision-making problems of pavement maintenance were investigated, and an overall procedure for decisionmaking optimization was proposed. First, PCI was calculated based on the deterioration type and area to evaluate road conditions using the AHP method. Then, two models, named MMQLB and MBMMQ, were developed. The first model was a single-objective linear optimization model, whereas the other was a MOO model. These models were employed for a real-world case involving the maintenance of 149 roads in Shushan District, Hefei, China. First, the pavement condition data of these roads were recorded and imported into the models. The models were then solved using sequential quadratic programming and a GA. Finally, appropriate pavement maintenance schemes were established based on the solutions of the optimization models. Among the optimal maintenance schemes of the MBMMQ model, some optimal values in the pair of maintenance quality and cost exactly matched the results obtained by the MMQLB model with the corresponding cost constraints. This indicated that the MMQLB model could achieve similar decision making to the MBMMQ model by setting different maintenance constraints. According to one nonlinear constraint of MMQLB, it could be improved to substitute the nonlinear formulation and reduce the difficulty of solving the problem in future research. In practical scenarios, the pavement maintenance cost is also dependent on the time taken and priority levels of the road grades. Therefore, the parameters used in the proposed models will be further extended in future studies to reflect real-world scenarios better and to broaden the applicability of the models. Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The raw data of roads and their maintenance cannot be shared at this time, as it belongs to the Hefei Municipal Engineering Management Office. The processed data for solving the multiobjective models in this study are available upon request by contacting the corresponding author.