Towards Road Sustainability—Part I: Principles and Holistic Assessment Method for Pavement Maintenance Policies

: Assessing the holistic sustainability of public policies remains a challenge rarely taken up due to a lack of adequate assessing methods. Frequently, only environmental and/or ﬁnancial aspects are addressed, rather than the three pillars, including macro- and micro-economic as well as social performance. This paper presents an assessment method to fully compare the performance of pavement resurfacing policies for all its stakeholders and considering pavement–vehicle interactions. First, an analytical and then systemic approach to road maintenance highlights all its stakeholders, and a complete set of sustainability indicators is proposed to quantify the various impacts of maintenance programs: tax revenues, road operator’s and users’ savings, domestic production and employment, net present value, users’ time savings and noise reduction health beneﬁts, as well as protection of natural resources, biodiversity and human health. Second, speciﬁc physical models of road condition (International Roughness Index) and its role in pavement–vehicle interaction in terms of vehicle consumption and wear as well as trafﬁc noise are introduced. Then, equations to calculate these indicators are presented based on a comparison of existing assessment methods. The ﬁnal transdisciplinary method pulls from road engineering, industrial ecology, acoustics and economics. It especially combines environmental and economic life cycle assessments and economic input– output analysis, as well as ﬁnancial and socioeconomic appraisals. Finally, this article takes up the interdisciplinary challenge of building a fully holistic assessment method to help decision makers properly address sustainability, and its general algorithm can be adapted to assess a variety of transportation policies.


Introduction
The sustainability of road transportation strongly depends on infrastructure maintenance policies, but the impacts of these policies have never been quantified with a comprehensive triple bottom line framework. As an example, in France, road transportation accounts for 8% of jobs and 13% of the gross domestic product, and it generates annual tax revenues of EUR 45 billion [1]. The average household spends more than 10% of its budget on cars [2]), and the average individual spends a cumulative period of 4 years on the road in their life (calculated from the French households transportation survey [3]). These socioeconomic figures have concomitant environmental implications, as 29% of primary energy consumption [4] and 28% of greenhouse gas (GHG) emissions [5] are attributable to road transportation in the country. Road vehicles are also responsible for local pollution, particularly particulate matter emissions, including PM2.5, that alone would cause almost 10,000 annual deaths in France (see Supplementary Material). Moreover, the impact of road noise affects 16 million people, around 1 million of them at critical health thresholds [6].
Road conditions-and thus pavement maintenance policies-influence this performance. For instance, road surface deterioration increases consumption and wear for the vehicles traveling over them due to pavement-vehicle interactions (PVI) [7], thereby affecting the overall efficiency of the road and especially its environmental performance [8]. Nevertheless, the holistic sustainability of pavement maintenance policies has never been assessed comprehensively. While some methods were developed to make triple bottom line assessments of specific activities or products, such as buildings in Turkey [9], small and medium manufacturing companies [10] or urban food system governance [11], methods are still needed for transportation [12] and its infrastructure maintenance plans. Numerous methods claim a sustainability approach while being in fact mono-pillar, i.e., considering only one of the three dimensions: social, environmental or economic.
The main objective of this article is to build a comprehensive method to assess the holistic sustainability of road maintenance policies, with a focus on pavement resurfacing and PVI. Indeed, besides major refurbishments, pavement resurfacing is the most common and important kind of road maintenance operation, consisting of (re)building the top layer of a road surface, by spreading road materials over a thickness of 1 to less than 10 cm. This paper's subobjectives are as follows: first, to define road maintenance stakeholders and the categories of holistic sustainable impacts affecting them; second, to introduce specific physical models of road condition (international roughness index (IRI)) and its role on the pavement-vehicle interaction in terms of vehicle consumption and wear as well as traffic noise and third, to put forward a set of equations to fully specify the assessment model and make it computable. A companion article "Part II" will present a case study on maintenance policies for an interurban highway in France and show the practicability of the method and its concrete usefulness.

Mono-Pillar Sustainability Methods and Hybridization
Fairly extensive literature exists on the environmental evaluation of roads through life cycle assessment (LCA) [13], sometimes used specifically for the question of resurfacing [14][15][16][17][18][19][20][21]. Life cycle cost analysis (LCCA) can also be used to financially optimize roads [22][23][24][25]. Next, the input-output analysis (IOA) approach ( [26] is an operational method for grasping the macroeconomic impacts of an economic policy on domestic production and employment [27]. It has been used in the transportation sector to assess the impact of a new road in the U.S. [28] or in France [29], of rolling out electric vehicles in France [30] or of the Paris-Bordeaux high-speed rail line in France [31], but never for road resurfacing. Finally, hybrid methods are used to assess the impact of roads excluding road maintenance. For example, a socioeconomic appraisal method is used by the French Ministry of Transportation to assess the public utility of a new transportation infrastructure, which combines macroeconomic, financial, social (time saved) and environmental dimensions [29]. This method has been hybridized with LCA in the case of an urban road rehabilitation operation [32]. However, it is not suitable for analyzing maintenance and there are deficiencies in its scope. Moreover, whole life cost analysis-hybridizing LCA and LCCA-was recently used in a pavement management study [33].

PVI Models and Road Maintenance Sustainability
Many models of vehicle operating costs have been developed since the 1970s [7]: incremental models linking the IRI with vehicle consumption and wear were developed for the World Bank between 1971 and 2005 by the Massachusetts Institute of Technology (MIT) [34], the Transportation Research Lab [35][36][37] and Birmingham University [38]. Originally developed for roads and vehicles in developing countries, these models were first calibrated for a developed country in the U.S. by Michigan State University [7]. This calibration reveals consumption differences of around 3%, 8% and 130% in fuel, tires and suspension systems, respectively, between vehicles traveling on a very smooth surface (IRI = 1 m/km) and vehicles running on a road in fairly poor condition (IRI = 6 m/km), these percentages depending on the vehicle type and its roll speed.
Finally, several studies have assessed some sustainability consequences of vehicle consumption due to road surface conditions. We summarize the characteristics of these studies in Table 1 and describe them further. The models of vehicle consumption sensitivities to the IRI established by Chatti and Zaabar [7] have been used several times to assess the effectiveness of road maintenance practices. Wang et al. [8] compared the improvements in GHG emissions and energy consumption obtained from a resurfacing program leading to a lower than usual IRI, including in the scope of their assessment both the impact of roadwork and the fuel consumption of the vehicles traveling on the roads. Chong and Wang [39] conducted a fairly similar study, extending the latter approach with an optimization program considering the initial dimensions of the road. Wang, Harvey and Kendall [40] assessed the optimum potential gain in GHG emissions in relation to road maintenance policy in the State of California and calculated the cost of mitigation in dollars per ton of CO 2 equivalent saved. Bryce et al. [41] presented a five-year optimization model that minimizes total energy consumed in both road maintenance and fuel consumption, constrained by the road operator's costs and road conditions. Yu et al. [42] developed an optimization program of the same type, extending the environmental scope to GHG emissions, acidification potential and respiratory effects, using a weighting method. Yang et al. [43] compared the impact of resurfacing programs with differences in the proportion of reclaimed asphalt pavement included in the new resurfacing material from an environmental perspective-GHG emissions and primary energy consumption-and from the perspective of the financial impact on the road operator. The only environmental effect from PVI considered was fuel consumption. Similarly, Santos et al. [44] compared the environmental performance of several resurfacing techniques but also left out the impact of differences in tire and suspension wear arising from PVI. This performance was assessed using the environmentally extended input-output method but not LCA in a second comparative study considering tires and suspension parts [33]. The latter study also includes an economic performance indicator, in the form of net present value (NPV), including the economic flows associated with road maintenance as well as vehicle consumption, calculated with a single discount rate of 2.3%. A decision support method for road surface management, developed by Pellecuer et al. [45], calculates the monetized health impacts of road noise and atmospheric pollution, as well as the cost of carbon dioxide emissions associated with traffic, under different road surface conditions. It also includes the variation in fuel consumption but not wear on tires and suspension. This study considers the monetary cost of noise pollution.
Finally, the marginal effect of vehicle consumptions according to road condition on the time spent by the users to operate their vehicle is considered by Guevara [46]. Nevertheless, no studies have compared the excess time spent because of roadwork with the time saved on the logistics of vehicle consumptions (journeys to gas stations to fill up and to service stations for tire and suspension maintenance).

Article's Scientific Contributions
As shown in the literature review, the methods developed to assess the sustainability of road maintenance policies are numerous but fragmentary with respect to the range of impacts on sustainable development and on stakeholders, as well as in consideration of the consequences of PVI. We thus propose a method that completes the existing approaches by adding the social, macroeconomic and financial aspects missing in the literature and by considering all the consequences of PVI-in particular the previously neglected environmental contributions of tire and suspension system wear. A novel method is also integrated to quantify the impact of road pavement resurfacing on traffic noise over time and ultimately damage to road local residents' health.

Identifying Road Maintenance Stakeholders and a Comprehensive Classification of Impacts
We built the method based on a reflection on which fundamental goals a public policy should follow towards a sustainable pathway, considering all four thinking schemes of socioeconomic ethics [47], main findings of the economics of welfare [48] and happiness [49] and the analysis of the original French texts theoretically driving the national social contract [50][51][52]. After a thorough analysis, the socioeconomic factors of happiness which are related to the strategies of road maintenance are: health, safety and security, employment, leisure and free time, resources and natural environment (details are available in the Supplementary Material and a Ph.D. thesis manuscript [53]). These factors affect the following stakeholders: the asset manager (i.e., road operator), the users, the residents, the State (or Government) and the environment.
Then, from a road maintenance stakeholder analysis and the multiple interests underlined in the literature review, we selected a set of indicators to holistically assess the multiattribute performance of a road resurfacing policy. The number of indicators was arbitrarily restricted to around 10 to avoid decision making cognitive saturation [54]. Nonquantitative impacts such as user comfort were not considered, and road safety was not accounted for due to its inherent multifactorial aspect, only partially depending on the infrastructure itself [55], thus difficult to separate from the other factors. Finally, the set of eleven indicators is composed of: users' time savings, health protection related to road noise, protection of resources, biodiversity and human health, road operator's and users' savings, domestic production and employment, tax revenue and national economic savings.
The stakeholders, impact causality chains and selected indicators are illustrated in the Figure 1. It shows how PVI are at the center of the stakeholders' impacts.

Overview of the Method's Architecture
The architecture of the integrated method that we developed to assess the different indicators selected in Figure 1 is represented in Figure 2. As an input, it uses IRI data of the road studied at the time t0. Depending on the pavement resurfacing scheme, it simulates the evolution of the IRI over time and calculates the associated inventories of consumptions and emissions on the assessment period: road works and vehicle consumptions, as well as noise emissions. It then combines assessment methods covering all aspects relating to road maintenance sustainability for the stakeholders: whole cost analysis, financial and socioeconomic assessments, LCA and Leontief's IOA. Combining these methods allows for calculating the main impacts caused by a maintenance policy as a result of the speed and intensity of road surface deterioration throughout its lifespan and on the different stakeholders highlighted previously.

Overview of the Method's Architecture
The architecture of the integrated method that we developed to assess the different indicators selected in Figure 1 is represented in Figure 2. As an input, it uses IRI data of the road studied at the time t 0 . Depending on the pavement resurfacing scheme, it simulates the evolution of the IRI over time and calculates the associated inventories of consumptions and emissions on the assessment period: road works and vehicle consumptions, as well as noise emissions. It then combines assessment methods covering all aspects relating to road maintenance sustainability for the stakeholders: whole cost analysis, financial and socioeconomic assessments, LCA and Leontief's IOA. Combining these methods allows for calculating the main impacts caused by a maintenance policy as a result of the speed and intensity of road surface deterioration throughout its lifespan and on the different stakeholders highlighted previously.

Overview of the Method's Architecture
The architecture of the integrated method that we developed to assess the different indicators selected in Figure 1 is represented in Figure 2. As an input, it uses IRI data of the road studied at the time t0. Depending on the pavement resurfacing scheme, it simulates the evolution of the IRI over time and calculates the associated inventories of consumptions and emissions on the assessment period: road works and vehicle consumptions, as well as noise emissions. It then combines assessment methods covering all aspects relating to road maintenance sustainability for the stakeholders: whole cost analysis, financial and socioeconomic assessments, LCA and Leontief's IOA. Combining these methods allows for calculating the main impacts caused by a maintenance policy as a result of the speed and intensity of road surface deterioration throughout its lifespan and on the different stakeholders highlighted previously.

Practical Use of the Method
The method was developed to be applied at the scale of the road section. The section must be considered over an observation period that is long enough to contain several resurfacing cycles and where the physical condition of the section is identical at the end of the analysis period for all the maintenance plans considered. A maintenance plan consists of a temporal sequence of resurfacing operations. Several resurfacing techniques exist, differing in terms of the thickness of the materials used, the quality of these materials and the construction technique of the surface layer. Each technique presents a specific financial, macroeconomic and environmental performance profile from the point of view of the maintenance operation. On the other hand, each maintenance technique has a potentially different impact on road surface condition and therefore the consumption and wear of the vehicles it supports. The vehicle population carried by the road section is separated into different categories of vehicles. The different consequences of PVI on these categories in terms of consumption and wear are accounted for. The impacts of this consumption and wear, as well as those of road works, are considered over their entire life cycles, as macroeconomic, financial and environmental consequences happen over their entire supply chain. The environmental aspects quantified are therefore calculated on a global scale. On the other hand, the health impact of road noise for residents relates to a local and not a global indicator. Finally, macroeconomic indicators focus on a national approach.

Capturing the Evolution of Road Surface Condition
Considering the literature, the surface condition is expected to be a key parameter of road maintenance sustainability. From the IRI data for the road section to be assessed from time t 0 , we calculate the evolution of the IRI over the entire lifespan of the road surface, between two resurfacing operations and before/after resurfacing operations.

IRI Evolution between Two Resurfacing Operations
The studies published on changes in the IRI between two resurfacings [8,23,43,[56][57][58][59][60] propose different formulas in the mathematical form described in Equation (1)-where a, b and d are calibration parameters, and t is time.
In previous studies, this ultimately corresponds to deterioration speeds ranging from 0.019 m/km·year [60] to 0.23 m/km·year [59]. We recommend calibrating the equation for every road section or network assessed using field measures.

Effect of Resurfacing Operations on the IRI
The effect of resurfacing works on the IRI can be estimated by calibrating the American equation established by Wang et al. [59] using the IRI data for the road assessed as recalled in Equation (2), where e and g are calibration coefficients depending on the resurfacing technique, t R+ the time just after resurfacing and t R− the time just before.
By default, for a resurfacing technique with a thickness between 3 and 7.5 cm, one can take the parameter values of the study referred to, recalculated here in international system units: e = 0.40 and g = 0.67 m/km [59].

Consumptions and Emissions over the Pavement Surface's Lifetime
Vehicle consumptions and emissions partly depend on the IRI, which partly depends on the resurfacing program and thus on maintenance operations. To compare the sustainability of different road maintenance programs, we assess for each option the resurfacing works demand, the vehicle consumptions-in fuel, tires and shock absorbers-as well as the road noise emissions. PVI intensity and consequences depend on the kind of vehicle considered. To model the traffic, we choose a typology with 4 categories of vehicles like Chatti and Zaabar's method [7]: passenger car (PC), light commercial vehicle (LCV), small heavy vehicle (2 axles) (SHV) and large heavy vehicles (3 axles and more) (LHV).

Resurfacing Works Demand
Based on the resurfacing programs to be tested, with Equation (3), we calculate the resurfacing works demand RWD k (in square meters) using a resurfacing technique k depending on the width ω of the pavement to be resurfaced in meters, its length λ in meters and N k , the number of resurfacing operations of type k over the assessment period.

Excess Vehicle Consumptions
Based on the data for pavement condition and dynamic traffic over time, we calculate the excess consumptions due to pavement surface irregularities over the pavement's lifetime for each vehicle category and consumable type-fuel, tires and suspensions. To do this, we consider the French average consumptions and select, after a model comparison, the IRI consumption sensitivity factor from HDM-4 [38] calibrated for the U.S. [61], which we adapt to the speeds actually observed on the different types of French road networks by statistical regression (see Supplementary Material). French speed limits are very similar to those in other European countries, but the equation adaptations can be conducted using the same method for each country's real conditions.

•
Excess fuel consumption The equations to quantify excess fuel consumption (EFC), depending on IRI, vehicle type and road type, are calculated in the Supplementary Material and presented in Table 2 for each kind of intercity network in France considering the average speeds on these networks in km/h for light vehicles (LV)-PC and LCV-and heavy vehicles (HV)-SHV and LHV. Again, regressions can be calculated using the same method for each country's real conditions. Fuel consumption due to the deterioration of the surface of the road over time can be calculated using Equation (4): with FC j,k the fuel volume of type k overconsumed on a road of type i over a period from t 0 to t f in hundredths of liters, λ the length of the road section evaluated in kilometers, EFC j,i the EFC of a vehicle of type j on a road of type i in percentage, AFC i,j the dynamic average fuel consumption of a vehicle of type j on a road of type i in liters per 100 km and Q j the traffic flow of type j vehicles on the road section considered in number of vehicles.
• Excess tire wear The equations to calculate excess tire wear (ETW), depending on IRI, vehicle type and road type, are calculated in the Supplementary Material and presented in Table 3. for each kind of intercity network in France and LV and HV average speeds, respectively. Specific national regressions can be calculated using the same method for any country. Tire consumption due to the deterioration of the surface of the road over time can be calculated using Equation (5): with TC i,k the number of tires of type k overconsumed on the road of type i, λ the length of the road section evaluated (in kilometers), ETW i,j the ETW for a vehicle of type j on a road of type i, TWR i,j the kilometric tire wear ratio for a vehicle of type j on a road of type i and Q j the traffic flow of type i vehicles.

• Excess suspension wear
According to Chatti and Zaabar's model, vehicles age faster on roads only with IRIs over 3 m/km. In our model, we will consider that only the suspension systems age faster under high IRIs. Over the 3 m/km threshold, the equations to calculate the shock absorber lifespan, depending on IRI, vehicle type and road type, are calculated in the Supplementary Material and presented in Table 4, for each kind of intercity network in France and LV and HV average speeds, respectively. These equations can be recalculated for any network in any country using the same process presented in the Supplementary Material. Shock absorber consumption for IRIs above 3 m/km can be calculated using Equation (6): with SAC i,k the number of shock absorbers of type k consumed on the road type i, λ the length of the road section evaluated (in kilometers), SAL i,j the SAL for a vehicle of type j on a road of type i and Q j the traffic flow of type j vehicles. Note that this time, we directly calculate a consumption rather than an overconsumption: this does not change the final calculation, as it will be conducted on a differential approach between resurfacing programs.

•
Unitary noise linear power level over time Noise characterization can be calculated through the unitary noise linear power level (Lw/m) corresponding to the noise power emitted by one vehicle on one meter. To understand its evolution over time, the statistical models of tire-road noise evolution from the European benchmark method [62] are updated using the most recent version of the French road noise database established by the CETE of East, a public technical center that reports to the French Ministry of Transportation (see Supplementary Material). These new models, different for heavy vehicles (HVs) and light vehicles (LVs), are presented in Table 5. Rx corresponds to acoustic categories of pavement rolling course, these categories being detailed in the Supplementary Material. Although the model is updated using a French database, physical acoustic phenomena is the same everywhere; thus, it must be usable in other countries with similar climate conditions, where pavement ageing is equivalent.

Rx
Light Vehicles Heavy Vehicles The noise of one motorized vehicle rolling on pavement is calculated by adding two noise components: the noise emitted by the motor and the noise emitted by the tireroad interaction. To calculate the unitary noise linear power level at t = 2 years, Lw/m (t = 2 years) per component and type of vehicle, we use formulas from Sétra [62] detailed in the Supplementary Material in the conditions of French interurban roads. We obtain Table 6. Table 6. Average noise power level-motor and tire-road components-per meter of source line for a two-year-old rolling course.

Lw/m-Motor Component (dB(A))
Lw/m- We can then calculate Lw/m (t = 2 years) per type of vehicle (LV or HV) by adding the noise source u-motor and tire-road noise-using the generic noise addition Equation (7).  (8) is the formula to calculate the linear density ρ j in vehicles of type j per meter, with q j the traffic in number of vehicles j per hour and V j the average speed for vehicles of type j.
Then, Equation (9) is the formula to calculate L, the noise level of a road per meter depending on the linear density of heavy and light traffic d LV and d HV during the day (6 a.m. to 10 p.m.) or the night (10 p.m. to 6 a.m.) and the noise level of one HV and LV L LV and L HV per meter. • Additivity of the line sources' noises and temporal weighting In the case of multiple road lanes, the noise of each line source must be added to obtain the linear power level of the road. Several methods can be considered to add the line sources' noise levels, depending on the road and noise receptor configuration. The method must be selected depending on the configuration studied.

•
Calculation of the sound power With Equation (10), we then calculate the sound power W /m emitted by the traffic from the linear power levels of the road, using Equation (9), with W 0,linear = 1 pW/m = 1 × 10 −12 W/m [63].
The sound level is then integrated over the length λ of the road section studied to obtain the average sound powers during the day (Wj) or the night (Wn) according to the following Equation (11): This result is then converted into acoustic energy emitted by multiplying by the time period considered. Our sound energy varying with the aging of the road surface and the latter being calculated over a time step of one year, we therefore multiply W tot by 3600 s × 365 days × the number of hours of the period considered (16 h for the day and 8 h for the night) to obtain Ej and En, respectively, the daytime and nighttime acoustic energies emitted in J(A).

Three Pillar Indicators' Algorithm
The classification of performance indicators between the three pillars of sustainable development is a matter of perspective when an indicator can be attached to two or three pillars. We propose the following classification, recognizing this potential intrinsic overlap. The method aims to assess the benefits of an alternative maintenance scheme (Alt) compared to a business-as-usual (BAU) program. Thus, performance indicators will always be calculated as the difference between the alternative scheme's performance and the BAU scheme, using the generic Equation (12):

LCA: Characterization Methods and Indicators
We use LCA to evaluate the environmental impacts, following ISO standards 14040 and 14044 [64,65]. This method quantifies the environmental impact of a system-being a product or an activity-by inventorying input and output flows crossing the system over its life cycle and calculating their impact on the environment using characterization methods. The inputs and outputs are product, energy or material flows. Input flows come from the natural environment or the technosphere. Based on this inventory, characterization methods relate the flows to their potential impact on the environment, potentially considering their fate and the exposure of specific ecosystems. These methods make it possible to calculate many different environmental impacts, and the same type of indicator-e.g., impact on climate change, acidification or eutrophication-can be calculated in multiple ways depending on the characterization method chosen. The quality of a characterization method depends on the good representativeness of the physical realities modeled: the method must be scientifically up-to-date and adapted to the area of assessment. Two kinds of LCA indicators exist: midpoint and endpoint indicators. Midpoint indicators focus on single environmental problem-such as climate change or acidification-while endpoint indicators aggregate midpoint indicators in the three areas of environmental protection: natural resources, ecosystems and human health. A complete midpoint indicator set has a dozen indicators. As we want to restrict our indicator set to around ten, we choose to select endpoint indicators, which allow covering the total environmental damage concisely. Selecting LCA endpoint indicators rather than midpoint indicators is doubly beneficial: focusing on meaningful environmental indicators for the preservation of our planet and restricting the number of performance indicators to help decision makers while encompassing all the midpoint impact categories, including the most popular one: climate changecontribution. Further adaptations of this method do not exclude calculating a complete midpoint set to assess the environmental impacts of resurfacing policies and, for instance, a climate change contribution indicator.

Selection of Endpoint Indicators
We select the two operational environmental endpoint indicators from the method IMPACT World+ [66]: damage to ecosystems (renamed "biodiversity") and to human health. IMPACT World+ is the update of IMPACT 2002+, LUCAS and EDIP, and currently the most scientifically advanced characterization method in LCA [67]. The IMPACT World+ biodiversity damage indicator covers the effects on biodiversity of (short-and long-term) marine, freshwater and soil acidification; freshwater and marine eutrophication; land use; ecotoxicity (short-and long-term); climate change (short-and long-term); ionizing radiation and heat pollution affecting water quality and water availability (for terrestrial and aquatic freshwater ecosystems). The biodiversity loss calculated is expressed in PDF.m 2 .yr: it corresponds to the potentially disappeared fraction of species over a surface area of one squared meter over one year. The IMPACT World+ indicator for damage to human health includes the following effects: climate change (short-and long-term), carcinogenic and non-carcinogenic toxicities (short-and long-term), ionizing radiation, ozone layer depletion, formation of fine particulates and photochemical oxidants and availability of water. This indicator is expressed in DALY, meaning disability-adjusted life years. One DALY is equivalent either to one year of potential life lost through premature death or one year of productivity lost as a result of incapacity. Note that the human health damage indicator includes the impact of air pollution over the entire supply chain of road maintenance consumption and road usage. Air pollution due to roads is a major public health issue, even more when highly trafficked and localized in densely populated areas. However, manufacturing is also a major source of air pollution: we chose an indicator that does not discriminate the damage to human health due to air pollution around roads and around facilities worldwide (where, for instance, spare parts or fuel are produced), whereas socioeconomic appraisals only consider local pollutions.
A third endpoint indicator focusing on damage to non-renewable resources and calculated with the ReCiPe method [68] comes to complete this environmental set of indicators. The latter, expressed in dollars, counts the economic impacts relating to the consumption of fossil and mineral resources, using a marginal cost approach. This approach is underpinned by the notion of resource scarcity: the scarcer a resource, the more its consumption will generate additional future extraction costs. ReCiPe offers the most advanced resource damage indicator, as it is an update of the Eco-Indicator and CML methods [67].

Calculation of the Metrics
Equation (13) is the formula to calculate the impact of a unitary consumption I j,k,o for each type of consumption k-a type of tire, suspension system or fuel needed to travel one kilometer-per type of vehicle j, for each type o of the three damage indicators. It sums the products of each flow of consumption or emissions f low c listed in the corresponding consumption life cycle inventories by characterization factor CF c,o obtained from the methods IMPACT World+ or ReCiPe, we.
The environmental impact of type o related to switching from a BAU resurfacing scenario to an alternative scenario is then calculated by multiplying the unitary impact of each type of consumption k by the number of units consumed over the assessment period for the total traffic. Consumption in one resurfacing scenario can be calculated with Equations (3)-(6).

Social Metrics
Two social indicators are considered in the method: users' time savings and the impact of road noise on residents' health. The latter could be considered an environmental indicator. However, as it involves the local population exposed to the road's traffic, it is considered a social indicator.

Road Noise Health Impact Indicator
The method includes an indicator of the impact of road noise on local residents' health DALY ∆t . It combines the road noise model developed in the noise inventory section of the article with a road noise LCA model developed by Meyer (2017). More specifically, we will use Meyer's characterization factors CF ∆t that relate the noise energy E ∆t emitted by the road to human health damage according to Equation (14), with CF day = 6.61 × 10 −7 DALY/J(A) and CF night = 1.25 × 10 −5 DALY/J(A).
The residents' health damage due to traffic noise will be calculated over the assessment period according to Equation (15).
Contrary to air pollution impacts on human health that are considered over the entire supply chain, for this indicator, we calculate the impact of road noise on human health of people living around the road, due to a lack of database on noise emissions of different activities.

Users' Time Saving Indicator
This indicator considers, based on the resurfacing program as well as fuel consumption, tire and suspension wear, the time spent by users in roadwork zones but also in journeys to gas stations or garages for vehicle maintenance and operation.
We use Equation (16) to calculate the time loss for a vehicle of type j due to a work zone, with θ alternative,j the time to cross the work zone during construction for the type of vehicle j and θ standard,j the standard time to cross the section. θ alternative,j can relate to time lost due to congestion, construction traffic lights, reduced speed limitations or detours.
TimeLoss works,j = θ works,j − θ standard,j Then, the time lost on the road section over the assessment horizon is calculated by summing the time lost for the four types of vehicles j and all the traffic during the n resurfacing operations, with q j (t) the hourly traffic of type j vehicles, ∆t_m the duration of the resurfacing operation k and N k the number of resurfacings of type k over the assessment period, as shown in Equation (17).
On the specific road network of type i, the total time lost due to vehicle consumption of type c (tires, shock absorbers or fuel) is calculated with Equation (18) from N c , the vehicle spare part consumption in number of total replacements of tires and shock absorbers or the volumes of 100 L of fuel for the vehicle of type j, and the estimated duration of the various operating and maintenance activities OD c,j , with c = 1 relating to tire consumption, c = 2 to shock absorber consumption and c = 3 to fuel consumption.
The users' time savings generated by an alternative resurfacing scheme compared to a reference are calculated using Equation (12) based on Equation (19): TimeLoss users = TimeLoss vehicle consumptions + TimeLoss works (19)

Economic Metrics
Based on resurfacing operations and vehicle consumption, we calculate the economic impacts of resurfacing programs, i.e., their macroeconomic effects and their financial consequences for the stakeholders concerned: the road operator, users and government. This requires cost models for vehicles and resurfacing consumption.

Users' Costs
The excess cost to use their vehicles for the road users VehicleExpenses users , discounted over time at a rate r, can be calculated over the assessment period by multiplying each type c of excess consumption EE k -namely, FC k the excess fuel consumption of type k, TC k the number of tires of type k overconsumed and SAC k the number of shock absorbers of type k consumed-by their costs C k (t) in constant currency, following Equation (20).
By default, we propose to consider a discount rate a equal to the government's 10-year borrowing rate, as it represents the household financial market reality.

Road Operator's Costs
The cost to maintain the road for the operator MaintenanceExpenses operator , discounted over time at a rate r, can be calculated over the assessment period by multiplying the surfaces of the road maintained using the resurfacing technique k, RWD k (t), by C k (t), the cost of this technique per square meter over time in constant currency, following Equation (21).
By default, we propose to consider a discount rate r equal to the rate of return of the operator if it is a private company or to the government's 10-year borrowing rate if the road is publicly operated.

Domestic Production and Employment
To calculate domestic production and employment content indicators, we use Leontief's IOA method [26,69]. This method is used by governments to analyze the national accounts and produce their macroeconomic projections. It uses input-output (I/O) tables representing production by economic sectors and their mutual dependencies to calculate the macroeconomic impact of an economic shock, i.e., a change in demand on certain economic sectors.
First, we convert the previously calculated physical consumption items from Equations (3)-(6) into monetary flows employing cost models. These models are specific to a country and a period. The calculation can be conducted on an annual basis to potentially consider dynamic cost models. Then we apply the Leontief method using the I/O tables-generally supplied by the national ministries for the economy-to calculate an indicator of domestic production. It includes the annual sum of direct and indirect domestic production resulting from all demands for roadwork, fuel and service station changes to tires and suspension. Two formulas are central. First is that of the technical coefficients in Equation (22) to calculate the share of intermediate consumption of a sector per unit of production of this sector, with α ij the coefficient corresponding to the share of expenditure on the product i in production of branch j, ter ij the matrix of intermediate consumption and Y = [y j ] the matrix of production.
Then, the inverse Leontief relation, noted in Equation (23), makes it possible, from the matrix A = [a ij ] of technical coefficients and the vector f of demand in a sector, to calculate the production sectorial effects P, with I the identity matrix.
We can calculate an annual vector f following Equation (24), where each element of the vector f corresponds to the differential consumption (in currency) of a maintenance scenario compared to the reference. Hence, only the elements corresponding to pavement construction, vehicle maintenance, tires, shock absorbers and fuel consumption (or larger activities depending on the country economic nomenclature), will be nonempty.
The direct and indirect domestic production related to the maintenance scheme assessed can be calculated by summing the elements of vector P.
Then, based on the production vector, Equation (25) uses a national job content vector for the different economic sectors [JobContent], which is, most of the time, supplied by national ministries for one country economy to calculate domestic employment in Full-Time Equivalent (FTE) jobs, relating to these demands, on the supply chain within the country where road maintenance is performed.

Tax Revenues
The tax revenues collected on consumption for roadwork, fuel, tires and suspensions, discounted over time at a rate r, are calculated using Equation (26), with f t the demand vector for the year T in constant currency and %tax t the tax rate vector (or a matrix in the case of tax rate variations over time), indicating the tax rate applicable for each type of consumption for the year T.
By default, we propose to consider a discount rate a relative to the government's 10-year borrowing rate, as it represents the State financial market reality.

Integrated National Economic Indicator
A government is often interested in the total cost of a policy such as a road maintenance program, combining the financial flows for the different stakeholders. The NPV indicator can be used to estimate the financial interest for an actor or group of actors in an operation: the larger it is, the more financially interesting the operation. It is used in particular by decision makers and road managers to justify their decisions. We use this concept to calculate the multistakeholder financial impact of a maintenance scheme-for the operator, users and tax authority. We recall in Equation (27) the generic formula for calculating an NPV, i.e., the sum of the financial flows-incomes and expenses that we detailed earlier for all the stakeholders-in constant currency, which is actualized each year over the entire duration of the assessment period, with E T the expenses, I T the incomes planned for the year t, and a the discount rate.
We suggest setting r at 2.5%, which is the risk-free discount rate recommended for the valuation of major investments in France. A risk-free rate choice allows us not to override the impact of future events [29] and, thus, to limit the burden of negative impacts for future generations.

Discussion
The geographical scope of the calculated impacts varies amongst the indicators: this scope goes from the close area around the road studied with the indicator of damage to human health due to road noise to a worldwide perimeter with the three LCA indicators, passing through a national scope with macroeconomic indicators. This heterogeneity of perimeters is partly justified by the scope of the impacts or the contribution of the various phenomena involved. For example, it is preferable to assess the environmental impacts at a global scale to limit the burden's geographic shifts. In addition, climate change, an important component of damage indicators, is based on global dynamics. On the contrary, it is likely that most of the noise changes linked to road maintenance policies are located around roads rather than around the manufacturing sites of road maintenance consumables because the impact of road maintenance policies on the production of these sites is marginal. Then, economic solidarity mostly occurs at a national level. However, the impacts of air pollution around roads and, more generally, in the country generate national socioeconomic impacts due to the cost of healthcare. Thus, comparing local or even national impacts to impacts in the rest of the world would advocate for or against further discretizing the two scales and potentially aggregate health damage due to noise and local air pollution into a resident human health damage indicator.
From a road management point of view, this method could benefit from future developments. It is for now restricted to the evaluation of resurfacing, but combined optimization of structural design and maintenance treatment scheduling over the lifespan of a road would allow for capturing the feedback effects between the mechanical behavior of the road structure and its surface layers. Such an approach will require advances in the field of road aging modeling. Next, the models linking road surface condition to vehicle fuel consumption and wear on tires and suspensions are central to the assessment mechanics in the method: these need to be developed and validated in regional conditions to enhance reliability of results. Uncertainty calculations could also be conducted to assess this reliability. Finally, a few indicators could be added. First is a user comfort indicator: riding comfort partly depends on the rolling course condition; nevertheless, we did not find any conclusive approach to propose an adequate quantitative indicator. Ideally, at the level of the road network, this user comfort indicator could be exploited in a more systemic way to calculate the impact of a resurfacing policy by considering the modification of the users' routes that may result from the evolution of the surface conditions. The selection of the route according to comfort can also imply differences in time, speed and consumption of travel, which are challenging to capture. Second is a pavement safety and reliability indicator quantifying the consequence of the road condition and, for example, cumulating the following costs: noninjury incidents, such as vehicle damage, goods breakage and damage for the freight, and monetized impact of injury accidents and fatalities. However, in practice, road safety depends on many factors, and, to our knowledge, the impact of the surface condition has not been decoupled from other factors [55]. Third is the tax revenue indicator, which could be refined to consider the financial impact of maintenance policies on social allowances from public authorities, e.g., on the unemployment allowances budget, depending on the employment performance of the policy.

Conclusions
In this article, we present the first method of holistic sustainability assessment developed to our knowledge in the field of transportation. Notably, it integrates the three pillars of sustainability and considers important macroeconomic and other state-of-the-art indicators compared to existing sustainability rating systems of pavement and transportation projects [70]. It especially completes the existing approaches by adding the social, macroeconomic and financial aspects missing in the literature and by considering all the consequences of PVI-in particular the previously neglected environmental contributions of tire and suspension system wear. A novel method is also integrated to quantify the impact of road pavement resurfacing on traffic noise over time and, ultimately, damage to local residents' health. A set of indicators that is both comprehensive and concise is needed to support decision makers for whom triple bottom line performance remains a theoretical concept that they struggle to grasp in its entirety. Indeed, holistic sustainability can hardly be defined "above ground", i.e., in a generic way without looking at a specific system or object. Therefore, we have developed this method on the specific example of road maintenance policies, considering sustainability under the French social contract. While the calculation of the inventories (Section 4) is specific to road maintenance, and the PVI equations are calibrated in French conditions as an example, the performance calculation algorithm (Section 5) can be used to assess all kinds of transportation policies in any country. The selection of the relevant quantitative indicators is based on an integrated vision of sustainability by including all stakeholders highlighted by the literature review. The benefits of an alternative road maintenance policy compared to standard practices for these stakeholders are assessed on all three pillars of sustainable development. Thus, this report offers a positive or upbeat consequential approach to triple bottom line decision making, in the manner of the environmental handprint concept [71], that reverses the vision of the classic environmental footprint calculated in attributional LCA. Finally, this method can be used by any road manager or road owner willing to fully understand and tackle its sustainability responsibility.
Supplementary Materials: Detailed information referred to in the main text as accessible in the supplementary material can be found in the supplementary material document following this article. The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3 390/su14031513/s1. The code for the method is available on the GitHub repository accessible at the following link: https://github.com/Anne2B/PhD.git (accessed on 15 January 2022). Funding: This research was funded by the industrial chair ParisTech-VINCI in "eco-design of buildings and infrastructure".