Multi-Criteria Analysis of Electric Vans for City Logistics

Atmospheric emissions of anthropogenic origin are one of the most important problems in cities. A particularly discrete ecological footprint in urban environment is made by urban freight transport. This problem has become the key challenge for all stakeholder groups involved in freight transport in urban areas. Over the recent years, there has been a growing interest in using alternative fuel vehicles in urban logistics, including those equipped with electric drive systems. This paper presents a multi-criteria analysis (MCA) of selected electric vehicles (EVs) in the context of their application for the purposes of deliveries in cities. In methodological terms, we present a unique MCA-based approach for evaluation of Electric Freight Vehicles (EFVs). A successful attempt is made to build a multistep MCA procedure based on two carefully selected MCDA methods (PROMETHEE II and fuzzy TOPSIS) to handle both certain and uncertain data sets in a single decision process. In practical terms, we successfully demonstrate the usefulness of the proposed approach by creating a set of decision maker’s preference models (based on certain and uncertain data) of carefully selected EFVs and provide sensitivity and robustness analyses of the obtained solutions.


Introduction
Atmospheric emissions arising from human activities are a common problem in medium and large conurbations and urban transport is often a major source of emissions [1,2].Apart from carbon dioxide emissions, the increase in the number of vehicles leads to a rise in emissions of nitrogen oxides as well as air pollution with fine dust resulting from abrasive wear of brake pads [3], tires [4,5] and road surfaces [6].In addition, the consumption of non-renewable fuels translates directly into environmental pollution [7], comprising emissions of: carbon monoxide (CO), nitrogen oxides (NO x ), sulphur dioxide (SO 2 ), hydrocarbons (CH x , including polycyclic aromatic hydrocarbons (PAH)), dioxins and benzene, as well as particulate matter (PM, mainly carbon, cadmium, zinc, nickel, platinum and chromium) [8,9].
The particularly discrete ecological footprint in urban environment is made by urban freight transport [10,11].Urban logistics involve various entities whose interests are often conflicting, which makes it more difficult to implement ecological solutions in this regard [12].Crainic et al. [13] assert that this is mainly due to the fact that the urban sector of freight transport is private.Private shippers often pay less heed to generating negative effects that mainly include noise, environmental pollution, compromised traffic safety and deterioration of urban infrastructure [14].A good example might be the problem with sustaining consolidation centres that often discontinue their activity as soon as the supporting project aimed at their establishing and start-up has expired [15].To a significant extent, this is caused by disparity of interests within public-private partnerships.
In view of the above, the problem of urban logistics operations in the context of its impact on the environment has become the key challenge for all stakeholder groups involved in freight transport in urban areas [16][17][18].Numerous initiatives have emerged in direct response to these problems, including international projects, such as BestUFS (BEST Urban Freight Solutions,), TRAILBLAZER (Transport and Innovation Logistics by Local Authorities with a Zest for Efficiency and Realisation), SUGAR (Sustainable Urban Goods logistics Achieved by Regional and local policies), C-LIEGE (Clean Last mile transport and logistics management for smart and Efficient local Governments in Europe), STRAIGHTSOL (Strategies and measures for smarter urban freight solutions), COFRET (Carbon Footprint of Freight Transport) and GRASS (Green And Sustainable freight transport Systems in cities).While local authorities and inhabitants are the ones most affected by the negative effects of urban freight transport, it is the private companies that are to a large extent able to increase the stability of solutions that reduce them.Dablanc describes it as a badly organized urban logistics system-when delivering goods, shippers adapt to the imposed restrictions, but, from the cities' point of view, the deliveries could be much more optimized in terms of their destinations, which would help reduce the number of vehicles entering the cities [19].Nevertheless, over the recent years there has been a growing interest in using alternative fuel vehicles in urban logistics, including those equipped with electric drive systems [20].
Since the 1970s, there has been a noticeable increase in interest in alternative fuels for cars [21,22].Actions aimed at implementing new technologies were forced by the oil crises [23]: first in 1973, then shortly afterwards in 1979-1982 and in 1990.The result of those actions led to the development of alternative car fuels.Currently, the alternative propulsion systems for motor vehicles include mainly [24][25][26]: gas-powered drives-LPG (Liquefied Petroleum Gas) and CNG (Compressed Natural Gas) [27], biofuels, hydrogen-based technologies [28] and electric drives [29].Gas-powered drives (LPG I CNG), compared to traditional fuels, produce 18% less greenhouse gas and generate ca. 3 dB less noise [26].Their major drawbacks include greater fuel consumption and limited possibilities of filling up the tank [27].Currently, it is the most popular alternative for vehicle powering; moreover, studies on gas fuels are at the most advanced level.Biofuels-in this case IC engines-are modified in such a way that instead of petrol or diesel they use ethanol, being the result of biomass fermentation [1].Currently, the US is the leader in biofuels use [30], as half of the maize grown there is intended for bioethanol production.There are two major ways of obtaining energy in hydrogen drives, where the fuel is hydrogen [31]: hydrogen combustion taking place in a typical piston engine combustion chamber and using fuel cells that produce energy as a result of oxidizing the fuel constantly supplied from outside.Hybrid drives involve a traditional engine combined with an electric one [32].Currently, the basic types of solutions used are, serial drives (IC engine operates all the time, driving the generator which in turn supplies power to the electric motor that propels the vehicle wheels, while any excess power is stored in the battery) and parallel drives (where both the IC engine and the electric one propel the wheels of the vehicle).Electric drives can use only battery powered electric engines [25].They produce less pollution as well as operate more quietly [33].Considering the restrictions that are more and more frequently imposed on freight vehicle traffic, as well as taxes and fees charged by local self-governments, application of electric vehicles may contribute to a reduction in delivery costs in cities, which may account for up to 40% of total transport costs [34].
One of the first research projects regarding the use of electric vehicles in urban freight transport was the ELCIDIS (Electric Vehicle City Distribution) project [35].However, the costs of purchasing electric vehicles are still perceived to be a substantial barrier to their widespread use.Additionally, a substantial difficulty lies in selecting vehicles with operation parameters that fulfil the needs of the logistic processes they are to serve.Therefore, the key challenge is the optimization of the transport fleet while taking into account a multi-criteria evaluation of benefits.City logistics involves implementing schemes in cities for reducing the impacts and costs of freight [36].There are a number of city logistic solutions that involve modifying freight vehicles including alternative fuels such as electric vehicles that can be implemented [37].More and more effective measures have started to be implemented in cities in recent years [38,39].
Although it has been more than 100 years since an electric vehicle first exceeded a speed of 100 km/h [40], the fundamental difficulty still to be overcome is the problem of storing large amounts of energy in batteries with relatively small weight, as well as the problem of how to charge them rapidly [41].The "Jamais Contente" electric vehicle developed by the designer Camille Jenatzy was equipped with 40 kM DC motor and a powerful acid battery bank that ensured the travel range of several kilometres.The problems encountered by Jenatzy have not yet been fully solved.Research on new electric motor designs is still underway, aiming at obtaining high operating parameters as well as advanced systems of drive control, supported by telematic solutions [42].The EU transport policy stipulates-according to the intentions presented in the White Book-to reduce the number of IC engine vehicles to 50% by 2030, and phasing them out by 2060 [43].Although EVs (Electric Vehicles) do not yet have a significant share in the European vehicle market, there is a noticeable, fast and continuous increase in the number of electric vehicles sold [34].
The application of electric vehicles in urban logistics is related to fundamental advantages of using this type of drives.However, it is important to underline some impediments and barriers that must be overcome in order to enable a significant increase in the interest in utilization of them in urban freight transport systems, especially considering the business stakeholders' perspective.One of them is the proper choice of the vehicles.The technical parameters of electric vehicles in city logistics are strongly related to the specificity of the tasks, which are realized by the stakeholder, including type of deliveries, number of deliveries, distance, etc. [44,45].Nowadays, no general method to help properly choose electric vehicles considers the needs of urban freight transport business' stakeholders.Due to the many technical parameters related to the city logistics tasks, these decisions are based on a multi-criteria approach [46].
The general aim of the paper is to provide a multi-criteria analysis approach as a support for this kind of decision processes.The authors' methodological contribution is an attempt to create a unique multistep approach based on MCA methods.This approach accommodates the analysis of certain data, as well as it provides recommendations in cases of both lacking or uncertain data in an integrated multistep decision model.Thus, it produces recommendations based on certain data, consistent with the decision maker's (DM's) preferences, as well as, concurrently, it allows including recommendations for the cases where the data are partially missing or uncertain.In practical terms, we apply the proposed framework in a future-proof problem of the choice of EFVs, for the purpose of deliveries in cities.In addition, generalized guidelines to facilitate effective selection of such solutions in urban logistics are provided.The article also presents the current state of knowledge regarding the extent to which electric vehicles are used in urban logistics and it describes the key parameters for determining the usability of electric vehicles in urban deliveries.Finally, based on the formulated assumptions, a multi-criteria model is presented to enable the selection, out of the defined vehicles catalogue, solutions that are optimal in terms of potential effectiveness.
The article does not discuss the use of passenger vehicles for the goods and passengers delivery at once, focusing exclusively on vehicles dedicated just to the cargo delivery.It is worth noting, however, that passenger vehicles are often used to carry loads and people at the same time.Considering this fact, it may be an interesting extension of our research but it would require expanding the database of electric vehicles with additional models and it requires the introduction of new decision criteria, such as a different payload than for freight vehicles, the number of people, etc.

Overview of the Current State of Knowledge Regarding Application of Electric Drives in Urban Logistics-Selected Examples
According to a report prepared by IHS Automotive, currently, the share of electric and hybrid vehicles registered in most countries worldwide does not exceed 1%.However, their number is growing rapidly, which is being confirmed by statistical data every year.Navigant Research predicts that by 2023 the share of electric vehicles in the global market will have risen to 2.4%.
According to the analyses prepared by Frost&Sullivan, the fastest growing sale markets for electric vehicles are Europe and China, and the number of such vehicles to be used worldwide is expected to approach 10 million within the next five years [47].As for non-European countries, apart from China, the greatest interest in electric vehicles is observed in, inter alia, the USA and Japan.These countries have become leaders in electric automotive industry development.In Europe, a significant increase in electric vehicles sales was not seen until a few years ago.In 2009, as few as 700 electric vehicles were sold, while in 2017 more than 100,000 EVs were registered; this does not include all types of hybrids with the possibility of external charging.In terms of the number of electric vehicles purchased in Europe, the leaders are: Norway, France, Germany and Great Britain.
One of the vital arguments for e-transport is not only ecology, but also economic considerations.The electric power needed to charge traction batteries costs six times less than fuel combusted by IC engine vehicles, assuming similar distances covered by both types of vehicles.Besides, the insurance and servicing costs are also lower.An electric drive is not only more ecological, but also simpler and less susceptible to technical failures compared to IC engines.Despite the higher price of electric vehicles, the maintenance costs of large company fleets based on electric drives may be substantially lower compared to conventionally powered vehicles.However, the prerequisite for introducing e-mobility is a well-developed and well-managed EV charging infrastructure [47].

Applications of Electric Vehicles in Urban Logistics
Traffic generated by freight transport accounts for 20-30% of kilometres covered by vehicles in urban areas and for 16-50% of air pollution derived by all kinds of transport [19,48].It would seem that the optimal solution should be a switchover to electric propulsion and replacing the conventional fleet with electric vehicles.Alas, shippers and drivers have bad experience related to initial implementation of electric vehicles.They primarily point out the lesser performance of electric vehicles (in comparison to the conventional ones) in terms of travel range, speed, acceleration and unreliability [49].Many research studies in that respect date back at least 20 years, e.g., under projects financed by the EU, such as EVD-POST (Electric Vehicle Deliveries in Postal Services) [50].The main conclusions highlighted the high costs of goods delivery, limited choice of vehicle models, lack of support with regard to vehicle sales and long lead time for spare parts, low capacity of batteries, limited travel range, relatively low speed and limited carrying capacity.In short, the authors of those studies indicated that the early versions of EFVs were far from perfect and they were not a serious alternative to conventional vehicles for the purposes of logistic operations in cities.
One of the first successes in terms of EV implementation is described in [51].It pertained to using EFVs within a cooperative framework for the purposes of urban freight transport, wherein its main idea was to establish an EFVs depot in the centre of Osaka and enable their use by multiple distribution companies.The project involved 28 electric freight vehicles, and 79 companies volunteered to participate in the undertaking.As a result of the system implementation, the number of kilometres covered was reduced, which led to a decrease in congestion.Another implementation involved a three-year project described in [52] that was mainly aimed at long-term measurements of fuel and energy consumption, travel range when using an electric motor, drivability, lithium-ion battery capacity and durability.A van prototype was proposed, equipped with a parallel architecture of hybrid drive based on five-cylinder 2.7 L diesel engine.The electric drive system was based on a permanent magnet 90 kW motor and 15 kWh lithium-ion (Li-Ion) battery.Another attempt at developing a modular structure of electric light trucks or vans (ELTVs) was the OPTIBODY project [53].This time, the new architecture was to contribute to improve the passive vehicle safety, in order to reduce the number of fatal accidents and injuries.
The authors of [54] provide a comparison of electric and fuel-powered urban freight vehicles in terms of life cycle, energy consumption, greenhouse gases (GHG) emissions and total cost of ownership (TCO).The authors engaged in studies using the New York City driving cycle that included frequent stopping and a low average speed, electric vehicles emitted 42-61% less greenhouse gas and consumed 32-54% less energy than freight vehicles equipped with diesel engines, depending on performance of individual vehicles.Table 1 presents selected EU projects implemented in 1998-2016, the main goals of those projects and involvement of twin cities.The analysis was based on [34] as well as final reports and web pages of the projects presented in Table 1.The project's central objective was to understand the technical and practical user requirements for using electric vehicles in commercial transport and to develop a set of methods for the fleet management of electric and mixed fleets.
Research studies focused on comparing three-wheel electric freight vehicles of small carrying capacity in relation to freight vehicles equipped with diesel engine were performed by Tipagornwong and Figliozzi [55].The studies were undertaken in urban areas, and the vehicles in question were the so-called trikes, low-emission vehicles using a combination of human effort and an electric motor for propulsion.
Other studies, described in [48], pertained to Amsterdam, a city that for many years has been striving to improve the air quality, road traffic safety and capacity, as well as to reduce noise levels by means of applying intelligent delivery methods and ecological transport.The research found that the current generation of electric vehicles is capable of operating urban freight transport quite efficiently, at the same reducing the distance covered by the vehicles by 19% (when applying consolidation), concurrently improving the air quality and reducing CO 2 emissions by 90%.
Taefi et al. [56] features an extensive analysis regarding the existing EFV initiatives in Denmark, Germany, Holland, Sweden and Great Britain.It concluded that there is a need for such research, as urban freight transport is thought to be one of the most promising areas of electric vehicles application.

Impediments to Application of Alternative Drive Vehicles in Urban Logistics Systems
Application of alternative drive vehicles-in particular electric vehicles-in urban logistics is related to fundamental advantages of using this type of drives.These involve the possibility to produce energy from any source, lack of gaseous or solid pollutant emissions into the atmosphere, lack of noise emissions, higher energy efficiency compared to traditional drives, cheaper production of drives, their maintenance and servicing, providing energy independence, low operation costs depending on the vehicle velocity and the price of 1 kWh.The existing electric power distribution network is the best developed part of the infrastructure.However, it is also important to stress the limitations and impediments which must be overcome in order to enable a significant increase in the interest in alternative drive vehicles.These may be broken down into three basic groups [57]: (a) economic difficulties including costs of purchasing such vehicles, costs of producing electric power for the vehicles, and costs of utilizing used-up batteries; (b) safety-related problems regarding the aspect of hardly audible operation of electric motor compared to conventional engines, which may consequently lead to collisions or accidents; but also potential risk of battery self-ignition; and (c) operational barriers, mainly connected with the fact that the full potential of alternative drives efficiency has not been fully achieved yet, long battery charging time, and limited cargo space in the case of electric vehicles.
What is important in terms of electric vehicles efficiency is the climate of the geographic area where such vehicles are used.It predominantly affects the vehicle's travel range: in winter months, more energy is needed for heating the interior, whereas, in summer months, air-conditioning increases the power demand.An important factor to limit the efficiency of urban freight transport using EFVs is city topography.As indicated in [58], energy consumption was high when driving on motorways and uphill.On the other hand, EFVs did well in urban traffic.Undoubtedly, the greatest barrier is the still high price of purchasing electric vehicles.They are usually more expensive than their traditionally powered counterparts.On the other hand, we can expect that businesses will be able to use different forms of support in purchasing vehicles of this type.This may mean tax relief as well as subsidies from local governments.A good example in this respect is Norway experiencing a substantial increase, the highest in Europe, in the number of electric vehicles.In this context, it is worth noting the conclusion drawn from the studies published in [59] that EFVs will not become competitive until savings derived from the operating costs reduction are sufficient to cover the substantially higher initial costs of purchase.
The most costly component of EFVs is usually batteries, the cost of which accounts for circa half of EFV retail price.Manufacturers focus on the technical capacity of batteries and increasing the operation range while decreasing their weight and charging time.Electric vehicle batteries must feature high Ah capacity; besides, they must be characterized by high power-to-weight and energy-to-weight ratios, as well as high energy density.Research and development in this regard is still underway [60].Currently, economic efficiency of EFVs is characterized by high purchasing costs of batteries, and lower costs of maintenance and fuel purchase.Over a short and medium term, EV costs are considerable, however, their advantages make them more competitive in a long run.The latest research has shown that TCO may be lower for electric vehicles compared to the conventional ones [54,61].It should be noted that EV purchases are currently most often financed with public funds.Private operators will be inclined to replace their fleets if: they notice benefits for their companies, an appropriate number of alternative fuel stations is provided, there are marketing benefits for the company, the given company is in any way connected with environmental protection, appropriate vehicles are available [62].Often, operators are forced to apply solutions of this kind by way of administrative decisions which prohibit in historic city centres or other specified areas vehicles other than those with alternative drives.The limited availability of information in model specifications and its operation is perceived as another factor to inhibit purchasing EFVs and using them in urban logistics.Source London reports that, "at the moment of purchase, there was not much information on basic elements i.e., screws and nuts or if the vehicle required water or oil".Additionally, the impediments that inhibit the use of electric and hybrid vehicles include [62]: higher operating costs of applying electric vehicles, small capacity of batteries, lack of battery charging infrastructure, low reliability and a large number of defects during complicated operation of the vehicles.
Over the recent years, considerable progress has been made with regard to improving the performance of electric vehicles, and many problems identified in earlier attempts at EFVs implementation in urban logistics have been resolved.Unfortunately, many problems entered a new dimension, and solutions have not yet been developed, even though the new EFVs feature longer travel ranges, more carrying capacity, better batteries and it takes less time to charge them completely or partially.It should also be noted that it is not plausible to perform all freight operations by means of EFVs.This relates to particular situations calling for vehicles with a large carrying capacity or long time of effective operation, e.g., when the vehicle must be in continuous operation for more than 12 h.

Key Factors Determining the Usability of Electric Drives in City Logistics
Although today's EFVs demonstrate better performance parameters: longer travel range and more capacious batteries and more carrying capacity, their actual usability in urban logistics still remains limited.The authors of the report [34] specified the challenges and factors that influence successful implementation of EFVs in everyday logistic operations in a city [57]: (a) Technical performance: travel range of EFVs usually does not exceed 100-150 km, even though the values specified by manufacturers tend to be higher.Potential stakeholders declare a greater interest in EFVs if they see improvement of operational parameters such as travel range, operation time with a battery, increasing the number of charging points and charging stations.(b) Operational performance: EFVs show both positive and negative features in comparison to IC engine vehicles.The positive features include the limited impact on the natural environment and decrease in noise level, which translates into their usability in city centres and time windows.At the same time, the charging, carrying capacity, maintenance and the need to adjust the logistic concepts to EFV application are perceived by operators to be the major operational challenges.(c) Economics: Purchase price and TCO still exceed corresponding values related to conventional vehicles.This results mainly from high costs of batteries and EFV small production limits.In the long run, these costs are expected to diminish, which will be connected with improving the operational parameters and efficiency, cutting the purchase prices, e.g., as a result of mass production.Resale value is highly unknown, which stops investors from making purchase on the primary market.(d) Environmental performance: EFVs cause less CO 2 emissions compared to their conventional counterparts, however, their demand for energy as well as energy prices in a long run and energy capacities of individual countries require analyses and research studies.(e) Social and attitudinal impact: EFVs are less noisy and more environmentally friendly than conventional ones, therefore most of the public show a positive view on this direction of development.(f) Impact of local policy and governance structure: Governments of many European countries take up new directives aimed at increasing the use of EFVs while decreasing the use of conventional vehicles.
Technical performance of EFVs (within the scope discussed further on in this article, i.e., carrying capacity, maximum velocity, travel range, engine power, engine torque, battery charging time, battery capacity) is improving with every new design and model.In addition, the economic factor, i.e., price, is changing.This considerably contributes to the growing interest from potential investors who also pay attention to additional aspects that are already implemented in conventional vehicles, e.g., telematics solutions and applications related to safety of drivers, passengers and cargo.
The assumptions for the model introduced below are one of the results of the cooperation with city logistics stakeholders under the work of the Freight Quality Partnership in Szczecin.The partnership has been established as a major activity of international projects C-LIEGE.The activities have been continued under the next international project GRASS.The FQP in Szczecin was focused on identification of major problems and establishing of the proper actions to improve the functioning of urban freight transport in city.The utilization of electric vehicles was one of the most important potential measures.The meetings with UFT stakeholders helped to establish among others the general expectations regarding implementation of electric vehicles in Szczecin city logistics system.

Proposed Conceptual Framework
The advantages and disadvantages enumerated in Section 2.3 with regard to the use of electric vehicles in urban deliveries are not decisive for their possible usability or effectiveness in everyday operation.Despite their undisputable ecological values and strong support with respect to sustainable growth, a separate analysis should be devoted to evaluation of technical performance of individual solutions in view of purchase costs.
The task of constructing a model for evaluating electric freight vehicles, as almost any actual decision-making problem, requires considering many decision alternatives, wherein each alternative should be investigated in view of many criteria.The proposed approach uses the classic MCA decision-making process assumptions.The conceptual framework presented below (see Figure 1) is in accordance with the classical guidelines formulated by Roy [63] or Guitonni [64].The problem of EFVs evaluation was broken down into four successive steps: Structuring Process, Determining and Modelling Preferences, MCA Aggregation Procedure and, Recommendation.who also pay attention to additional aspects that are already implemented in conventional vehicles, e.g., telematics solutions and applications related to safety of drivers, passengers and cargo.
The assumptions for the model introduced below are one of the results of the cooperation with city logistics stakeholders under the work of the Freight Quality Partnership in Szczecin.The partnership has been established as a major activity of international projects C-LIEGE.The activities have been continued under the next international project GRASS.The FQP in Szczecin was focused on identification of major problems and establishing of the proper actions to improve the functioning of urban freight transport in city.The utilization of electric vehicles was one of the most important potential measures.The meetings with UFT stakeholders helped to establish among others the general expectations regarding implementation of electric vehicles in Szczecin city logistics system.

Proposed Conceptual Framework
The advantages and disadvantages enumerated in Section 2.3 with regard to the use of electric vehicles in urban deliveries are not decisive for their possible usability or effectiveness in everyday operation.Despite their undisputable ecological values and strong support with respect to sustainable growth, a separate analysis should be devoted to evaluation of technical performance of individual solutions in view of purchase costs.
The task of constructing a model for evaluating electric freight vehicles, as almost any actual decision-making problem, requires considering many decision alternatives, wherein each alternative should be investigated in view of many criteria.The proposed approach uses the classic MCA decision-making process assumptions.The conceptual framework presented below (see Figure 1) is in accordance with the classical guidelines formulated by Roy [63] or Guitonni [64].The problem of EFVs evaluation was broken down into four successive steps: Structuring Process, Determining and Modelling Preferences, MCA Aggregation Procedure and, Recommendation.In the first step (Structuring Process) of the procedure, the sets of criteria as well as the alternatives available for evaluation (here EFVs) are identified.For this purpose, scientific literature has been used as a reliable source of information and expertise.The construction of the hierarchical family of evaluation criteria was based on a detailed analysis of references that problematically and methodologically cover the decision problem studied in this research.Their detailed analysis was carried out in Section 2.2.
The next two steps (Determining and Modelling Preferences, and MCDA Aggregation Procedure), are related not only to the properties of the studied EFVs' evaluation decision problem, but also to the utilised data aggregation techniques-MCA methods.A detailed analysis of the problem points to the need to analyse two types of data: certain and uncertain.In the above framework, two complementary data aggregation processes were proposed in a single decision support model.It was decided to extract the subset of options (alternatives) containing the full set of crisp data.For the remaining data and the data gaps, their fuzzy versions were developed.To provide a high certainty of the obtained solution, the analysis was performed in two steps.In the first step, the subset In the first step (Structuring Process) of the procedure, the sets of criteria as well as the alternatives available for evaluation (here EFVs) are identified.For this purpose, scientific literature has been used as a reliable source of information and expertise.The construction of the hierarchical family of evaluation criteria was based on a detailed analysis of references that problematically and methodologically cover the decision problem studied in this research.Their detailed analysis was carried out in Section 2.2.
The next two steps (Determining and Modelling Preferences, and MCDA Aggregation Procedure), are related not only to the properties of the studied EFVs' evaluation decision problem, but also to the utilised data aggregation techniques-MCA methods.A detailed analysis of the problem points to the need to analyse two types of data: certain and uncertain.In the above framework, two complementary data aggregation processes were proposed in a single decision support model.It was decided to extract the subset of options (alternatives) containing the full set of crisp data.For the remaining data and the data gaps, their fuzzy versions were developed.To provide a high certainty of the obtained solution, the analysis was performed in two steps.In the first step, the subset of the data with no gaps (crisp data) was analysed with the use of PROMETHEE II method (see Table A1).In the second step, a fuzzy variant of the TOPSIS method was used to evaluate the alternatives represented by fuzzy data (for two input subsets: with data gaps and full data set, respectively).It should be clarified that, while preparing the input data for the second step (fuzzy TOPSIS), it was decided that, for all values of criteria, their fuzzy representations would be developed in a trapezoidal form.It was also decided that the gaps would be filled with trapezoidal fuzzy numbers built based on the extreme domain values for the considered set.Technically, the minimum, maximum, mean and standard deviation values were used during the conversion process.For the remaining crisp values, their fuzzy developments were also prepared.
In both data aggregation scenarios, robustness and sensitivity analyses of the obtained set of decision solutions were performed.The analyses carried out together with their interpretation, as well as the comparison of the rankings obtained in both data aggregation procedures, are the basis for generating the final recommendations.This stage concludes the proposed authors' approach.In the authors' opinion, the individual data aggregation procedure, as well as the MCA numeric algorithms used in them (here PROMETHEE II and fuzzy TOPSIS) require a detailed explanation.
As already indicated, in the first of the two data aggregation runs in the model, PROMETHEE II method was used for the crisp data subset.The PROMETHEE methods are one of the best known and most popular outranking approaches used in sustainable energy planning [65].They have a number of features that are not available in other MCA methods such as the AHP, ELECTRE or TOPSIS [66].With the PROMETHEE method, the decision makers are assumed not to be fully aware of their preferences of the utility of different criteria.This is in contrast to approaches such as the AHP that assume that decision makers are aware of the utility of different criteria values and are able to express the relative importance of different criteria clearly.In addition, the GAIA plane provides a useful analytical tool for visualizing the decision problem that can be used for identifying conflicts between criteria and alternatives, allowing the most effective criteria to be determined.PROMETHEE II also has several other characteristics that make it suitable for evaluating electric freight vehicles such as providing a complete ranking of alternatives and allowing uncertainties in criteria performance values to be modelled.PROMETHEE also requires less input and little interaction with the decision makers, as well as being more transparent and easy to use and understand.The PROMETHEE II method enables a complete order of the resulting ranking of alternatives to be obtained, while in the case of the earlier version of this method, the resulting ranking was rather partial [67].
The PROMETHEE methods are used in determining a synthetic ranking of options.Depending on the implementation, they operate based on real or pseudo-criteria [68].The decision-maker may choose from six preference functions applying: simple criterion, quasi-criterion, criterion with preference level, criterion with linear preference, criterion with linear preference and indifference area, and Gaussian criterion [67].Subsequently, upon determining the compliance coefficients for each pair of options, dominance flows are determined for each of the options [69]:

•
Output dominance flow describing how much option a i outranks the other options: • Input dominance flow informing how much option a i is dominated by the other options: One by one, the decision-maker may establish the complete ranking of options.In the PROMETHEE II method, to establish a complete ranking of options, it is necessary to compute the net dominance flow according to Formula (3) [69]: In this method, equivalence and preference relations are distinguished in a broad sense [67]: The second run of the data aggregation was realized in the proposed authors' approach with the use of a fuzzy development of the TOPSIS method.The TOPSIS (Technique for Order Preference by Similarity to an Ideal Solution) method is based on the concept that the most preferable alternative should have the shortest distance from the ideal solution while, at the same time, be as far as possible from the negative-ideal solution [70].The algorithm of the TOPSIS method is composed of six stages [71].Initially, the decision maker chooses m alternatives and n criteria for use in the decision problem analysis.A decision matrix is created, where the rows represent alternatives and the columns represent the criteria.The data describing each criterion can be expressed in various units and scales; therefore, in the second step of the procedure, the decision matrix is normalized.In the third step, the normalized values are multiplied by the weight of each criterion, thus creating a weighted normalized decision matrix.In the fourth step, the Positive Ideal Solution (PIS) and the Negative Ideal Solution (NIS) are computed, by choosing respectively the best and the worst values from all alternatives describing each criterion.In the fifth step, Euclidean distances are calculated between the positive and negative ideal solutions and the alternative.Eventually, in the last step, the closeness coefficient CCi is calculated, which expresses the relative closeness of the alternative to the ideal solution.In the fuzzy version of the TOPSIS method, the numerical values are substituted by linguistic assessment.Instead of crisp numbers, trapezoidal fuzzy numbers (TFN) n = (n 1 , n 2 , n 3 , n 4 ) are used to express both the attributes of the alternatives, as well as the weights.The TFN's membership function µ n (x) is calculated as follows [72]: The distance between the TFNs m = (m 1 , m 2 , m 3 , m 4 ) and n = (n 1 , n 2 , n 3 , n 4 ) can be obtained with the vertex method [73]: It was also decided that, during the second run of data aggregation (the part containing fuzzy data), the weights of criteria would be expressed by seven linguistic variables: very low, low, medium low, medium, medium high, high, and very high.The TFN representations of the weights are presented in Figure A1.On the other hand, the values of the decision attributes for each alternative and each criterion would be described by a trapezoidal fuzzy number n = (n 1 , n 2 , n 3 , n 4 ), based either on the precise data, if available, or on the aggregated values from all alternatives for each criterion.For the former, n 1 = n 2 = n 3 = n 4 and were equal to the precise value.Therefore, the TFN obtained a form of a singleton.For the latter, n 1 was equal to the minimum value of the aggregated data; n 4 was equal to the maximum value; and n 3 and n 4 were equal, respectively, to the sum and subtraction of the half of the standard deviation value and the mean value.

The Set of Evaluation Criteria
MCA is a very popular methodology widely used in solving many real-life decision problems.It is one of the elemental tools allowing to obtain quantitative-based rankings and to generate decision recommendations on their basis, taking into account a multitude of often conflicting criteria.The literature review allows to identify a number of works using various MCA methods in solving the problems of evaluation or choice of electricity-powered vehicles for both logistic and passenger purposes.
For example, Barfod et al. developed a MCA/SWOT-based decision support system for analysing the challenges, opportunities and policy initiatives.The authors verified their approach in electric commercial vehicles market in Denmark [74].Domingues et al. successfully applied Electre Tri method to the Life-Cycle Assessment of six different powertrains vehicles [75].Oliveira et al. compared a conventional, hybrid and electric vehicles, addressing the Portuguese market and proposed an MCA-based approach during the evaluation process [76].Schmale et al. [77] proposed an Integrated MCA-based Assessment Method for Sustainable Transport System Planning in a Middle Sized German City.The authors of [78] [81].In the study by Pérez et al. [82], a broad state of the art analysis was performed on MCA approach in choice and evaluation of urban passenger problems.
As a result of a profound literature study, the authors chose nine evaluation criteria grouped into four clusters.Table 2 presents the selected criteria along with the abbreviations that were used in the following sections and references to other works justifying their choice.The criteria pertained to the vehicle performance in its broad sense: motor, batteries, and price.The structure of the criteria broken down into groups is presented in Figure A2.

The Set of EFVs Options
The market of vehicles is very diverse.One can enumerate the production of passenger cars, off-road vehicles, trucks, buses or special vehicles.These, in turn, can be broken down into the ones fuelled by internal combustion, electric and hybrid engines.The differentiation of vehicles is very large and is strongly associated to their intended purpose.For this reason, the authors focused on a single group of vehicles.The choice of the selected group of evaluated vehicles was determined by the purpose of the goods' transport for urban logistics, i.e., vehicles with the ability to navigate inside cities and city centres.The determination of such group of vehicles allowed the selection of models with similar technical parameters which resulted in their description with the criteria presented in Figure A2.Based on the above facts, this study considered 36 decision options evaluated from the perspective of nine criteria.
Unfortunately, with regard to many of the considered options, it was not possible to obtain all information on the criteria (car manufacturers do not publish all data, making the selection process more difficult), therefore, there were numerous data gaps and consequently the decision problem was considered under uncertainty.The structure of the decision problem is presented in Table 3.The evaluated vehicles were selected at the time of writing the article.Concept vehicles have been omitted from the list.Currently marketed or ready for sale vehicles have been chosen.

Analyses of Crisp Data with the Use of PROMETHEE II Method
Initially, the set of preferences was identified and modelled.The weights of individual criteria were defined along with their preference direction, preference functions and thresholds values.In this study, equal criteria weights were assumed so that none of them had a greater impact on the decision problem solution.For criteria such as "profit", the preference direction was "maximum", whereas in the case of "cost"-"minimum".The applied preference function was V-shape, so even small differences between the criteria efficiencies of options affected the determination of outranking relation between them.Concurrently, application of the p preference threshold made it possible to change the outranking value fluently.The p threshold was established as the two-fold value of the standard deviation for the option efficiency in relation to a given criterion.The full preference model is presented in Table 4.In the first step of the analysis, a subset of ten alternatives, which contained no gaps in data, was extracted for evaluation from the original set, to assure high certainty of the obtained ranking.The new set was comprised of the alternatives A6, A7, A8, A9, A15, A16, A17, A23, A27 and A28.The solution obtained for this set of alternatives is presented in Table 5.The data collected for this research did not consider how the battery capacity changes as its operating conditions, such as temperature, charging or discharging current or state of charge, as well as its service time vary [88][89][90].Therefore, in the subsequent step of the analysis, three additional rankings were generated with the battery capacity reduced by 10%, 20% and 30% for all the considered alternatives to verify how such change would affect the final ranking.It should be noted that, in these cases, the p preference remained unchanged.The obtained rankings, along with the Phi net outranking values are presented in Table 6.The sequence of options in the presented rankings correspond to the orders of options in the ranking with the original value of the battery capacity criterion.Only minor changes of the Phi net values are observed between the rankings.Thus, it can be assumed that the changes in the battery capacity caused by its operating conditions changes have a minor effect on the final solution of the decision problem.
In the third step of the analysis, the "Battery charging time 80%" criterion was excluded from the set of criteria.The obtained ranking is presented in Table 7.The comparison of the ranking with the one produced in the first step of the analysis is depicted in Figure 2.  The sequence of options in the presented rankings correspond to the orders of options in the ranking with the original value of the battery capacity criterion.Only minor changes of the Phi net values are observed between the rankings.Thus, it can be assumed that the changes in the battery capacity caused by its operating conditions changes have a minor effect on the final solution of the decision problem.
In the third step of the analysis, the "Battery charging time 80%" criterion was excluded from the set of criteria.The obtained ranking is presented in Table 7.The comparison of the ranking with the one produced in the first step of the analysis is depicted in Figure 2. The analysis of Figure 2 allows observing a considerable difference between both rankings.Only four of the ten studied alternatives remain on unchanged positions-A27 as the leading option, A23 on Position 6 and A16 and A28 as the worst alternatives.There was a swap of positions between the alternatives A15 and A17.A considerable change of ranks can be observed on Positions 2-5 among the alternatives A6, A7, A8 and A9.The aforementioned differences between rankings show that the "battery charging time 80%" criterion has a considerable impact on the final rankings.
Another step was the GAIA analysis (Geometrical Analysis for Interactive Assistance).The GAIA diagrams provide a graphical representation of the decision problem, which allows visually The analysis of Figure 2 allows observing a considerable difference between both rankings.Only four of the ten studied alternatives remain on unchanged positions-A27 as the leading option, A23 on Position 6 and A16 and A28 as the worst alternatives.There was a swap of positions between the alternatives A15 and A17.A considerable change of ranks can be observed on Positions 2-5 among the alternatives A6, A7, A8 and A9.The aforementioned differences between rankings show that the "battery charging time 80%" criterion has a considerable impact on the final rankings.
Another step was the GAIA analysis (Geometrical Analysis for Interactive Assistance).The GAIA diagrams provide a graphical representation of the decision problem, which allows visually evaluating how good each alternative is, which criteria support each alternative, and in which areas the alternatives could be improved.The information regarding k criteria and their influence on the alternatives' evaluation in the GAIA methodology, originally represented in k dimensions, undergoes a projection on the GAIA plane, as a result of which, part of the information is lost.The criteria are represented on the GAIA diagram by vectors, and the alternatives by points.The longer a criterion's axis, the more discriminating it is.Criteria expressing similar preferences are oriented in a similar direction, whereas opposite vectors' directions represent conflicting criteria.Orthogonal axes of the criteria symbolize criteria not related to each other in terms of preference.
Figure 3a shows the GAIA analysis of the studied decision problem, broken down into the individual criteria.As a result of projecting the decision problem solution onto the GAIA plane, the quality of 89.2% was obtained, therefore the interdependencies between the criteria and between the alternatives are depicted very reliably.
Sustainability 2017, 9, 1453 16 of 34 the alternatives could be improved.The information regarding k criteria and their influence on the alternatives' evaluation in the GAIA methodology, originally represented in k dimensions, undergoes a projection on the GAIA plane, as a result of which, part of the information is lost.The criteria are represented on the GAIA diagram by vectors, and the alternatives by points.The longer a criterion's axis, the more discriminating it is.Criteria expressing similar preferences are oriented in a similar direction, whereas opposite vectors' directions represent conflicting criteria.Orthogonal axes of the criteria symbolize criteria not related to each other in terms of preference.
Figure 3a shows the GAIA analysis of the studied decision problem, broken down into the individual criteria.As a result of projecting the decision problem solution onto the GAIA plane, the quality of 89.2% was obtained, therefore the interdependencies between the criteria and between the alternatives are depicted very reliably.As a result of the GAIA analysis, it can be recognized that the Battery Capacity criterion is in conflict with the Battery Charging Time criteria.This means that the greater the Battery Capacity, the longer its charging time.On the contrary, the pairs of Carrying Capacity and Battery Capacity, as well as Engine torque and Engine power express similar preferences, i.e., the value of Engine Torque increases along with the Engine Power value and the Carrying Capacity grows together with the Battery Capacity.The Travel Range criterion is unrelated to the Battery Capacity and Carrying Capacity in terms of preferences.The leading alternatives A29, A9 and A8 are most supported by the Max Velocity and Travel Range criteria.On the other hand, the alternatives A15-A17 are supported by the Price and Battery Charging Time criteria.The alternatives A6 and A7 are supported by the Battery Capacity and Carrying Capacity criteria.An analogous analysis was carried out for criteria groups (Figure 3b).
In this case, as a result of projecting the solution onto the plane, only 5% of the information was lost.This analysis indicates that the Price criteria group is in slight conflict with the Engine group.On the other hand, the Engine group is not related to the Battery group, as well as the Performance group values grow independently from the Price group in terms of preference.The Engine group supports the alternatives A6 and A7, whereas the Price group has a positive effect on the ranking position of options A8, A15-A17 and A23.As a result of the GAIA analysis, it can be recognized that the Battery Capacity criterion is in conflict with the Battery Charging Time criteria.This means that the greater the Battery Capacity, the longer its charging time.On the contrary, the pairs of Carrying Capacity and Battery Capacity, as well as Engine torque and Engine power express similar preferences, i.e., the value of Engine Torque increases along with the Engine Power value and the Carrying Capacity grows together with the Battery Capacity.The Travel Range criterion is unrelated to the Battery Capacity and Carrying Capacity in terms of preferences.The leading alternatives A29, A9 and A8 are most supported by the Max Velocity and Travel Range criteria.On the other hand, the alternatives A15-A17 are supported by the Price and Battery Charging Time criteria.The alternatives A6 and A7 are supported by the Battery Capacity and Carrying Capacity criteria.An analogous analysis was carried out for criteria groups (Figure 3b).
In this case, as a result of projecting the solution onto the plane, only 5% of the information was lost.This analysis indicates that the Price criteria group is in slight conflict with the Engine group.On the other hand, the Engine group is not related to the Battery group, as well as the Performance group values grow independently from the Price group in terms of preference.The Engine group supports the alternatives A6 and A7, whereas the Price group has a positive effect on the ranking position of options A8, A15-A17 and A23.
Apart from the GAIA analysis, the sensitivity analysis of the rankings was performed.When the problem is observed from various perspectives (overall performance, engine, battery life, and economical aspects), the relevance of the different sets of criteria can be seen as different, therefore, the weights of the clusters of criteria were examined in this analysis.The results of the analysis are presented in Table 8 and Figure 4.It can be observed that the ranges of stability of the Performance and Battery clusters are much wider than the ranges of stability of the Engine and Price clusters.If the criteria from the Performance cluster were assigned much more important weight (over 65.73%), the A27 alternative would yield to A9, which would become the leading alternative.Furthermore, if the weight of the Engine cluster was above 43.09%, the A6 alternative would become the leading one.On the other hand, if the Engine cluster was considered less important and its weight would drop below 5.46%, the alternative A8 would be positioned on the top of the ranking.the weights of the clusters of criteria were examined in this analysis.The results of the analysis are presented in Table 8 and Figure 4.It can be observed that the ranges of stability of the Performance and Battery clusters are much wider than the ranges of stability of the Engine and Price clusters.If the criteria from the Performance cluster were assigned much more important weight (over 65.73%), the A27 alternative would yield to A9, which would become the leading alternative.Furthermore, if the weight of the Engine cluster was above 43.09%, the A6 alternative would become the leading one.
On the other hand, if the Engine cluster was considered less important and its weight would drop below 5.46%, the alternative A8 would be positioned on the top of the ranking.The sensitivity analyses for distinct criteria are presented in Figure A3.It can be observed that, for the Battery Charging Time 100% (Bat6) criterion, the ranking remains stable only if the weight assigned to this criterion is in the range of 0-25.31%.On the other hand, for the Battery Charging Time 80% (Bat7), the alternative A27 remains the leading one regardless of the weight assigned to this criterion.The sensitivity analyses for distinct criteria are presented in Figure A3.It can be observed that, for the Battery Charging Time 100% (Bat6) criterion, the ranking remains stable only if the weight assigned to this criterion is in the range of 0-25.31%.On the other hand, for the Battery Charging Time 80% (Bat7), the alternative A27 remains the leading one regardless of the weight assigned to this criterion.

Analyses of Fuzzy Data with the Use of Fuzzy TOPSIS Method
Fuzzy TOPSIS Analysis of the Data Gaps Options Subset The PROMETHEE II analysis of the data without uncertainty was followed by the analysis of the alternatives which had some gaps in the data.The analysed subset is presented in Table A2.The analysis was performed with the usage of a fuzzy variant the TOPSIS method.The values of all criteria were represented by trapezoidal fuzzy numbers n = (n 1 , n 2 , n 3 , n 4 ).The gaps in the data were filled with trapezoidal fuzzy numbers built based on the minimum, maximum, mean and standard deviation values of all known values for each criterion.The values that were given were translated into trapezoidal fuzzy numbers where all the points had the same value, i.e., they were represented by singletons.All criteria were assigned a medium trapezoidal fuzzy number as their weight.A fragment of the obtained decision matrix is presented in Table A3.The ranking produced by the TOPSIS method, along with the CCi closeness coefficient values is presented in Table 9 and on the chart on Figure 5.The analysis of the obtained ranking allows observing that the best alternatives, A25, A2 and A20, are approximately two-fold better than the worst alternative, A24.Some groups of alternatives with very similar values of CCi can be observed, i.e., A32 and A12, A26 and A21, and A31 and A33 in the middle of the ranking, as well as A19 and A4, and A35 and A10.The similarities of the CCi value in these pairs of alternatives mean that they can be used interchangeably.
A sensitivity analysis of the ranking was performed by changing the weight of a single criterion between very low and very high, while all the remaining criteria had the medium weight assigned.The results of the sensitivity analysis are charted on Figures 6 and A4.It can be observed that, in the case of modification of the Per3, Eng4, Eng5, Bat7 and Bat8 criteria's weights, the alternative A25 remains unchangeably in the lead position, and alternative A2 receives the second best rank.If criterion Per1 would have a slightly lower weight, i.e., medium low instead of medium, alternative A2 would be better.However, if the weight of this criterion were reduced to very low, A20 would advance to the first position.Similarly, if the weight of Per2 was high or greater, the A20 alternative would outrun the alternative A25.It should be noted, that if medium high to very high weight is assigned to the Pr9 criterion, a significant of the A2 alternative can be observed, from Rank 5 to Rank 1.Other significant rank changes include A13's increase from Rank 9 to Rank 2 with the growth of Per1 criterion, A19's drop from Rank 9 to Rank 19 with Per2 growth, A22's advance from Rank 12 to Rank 4 with Per3 growth or A36's drop from Rank 3 to Rank 13 when Eng4 grows.The rankings remain most stable for criteria Eng5 and Bat7, where, respectively, the first six and four positions remain unchanged.The analysis of the obtained ranking allows observing that the best alternatives, A25, A2 and A20, are approximately two-fold better than the worst alternative, A24.Some groups of alternatives with very similar values of CCi can be observed, i.e., A32 and A12, A26 and A21, and A31 and A33 in the middle of the ranking, as well as A19 and A4, and A35 and A10.The similarities of the CCi value in these pairs of alternatives mean that they can be used interchangeably.
A sensitivity analysis of the ranking was performed by changing the weight of a single criterion between very low and very high, while all the remaining criteria had the medium weight assigned.The results of the sensitivity analysis are charted on Figures 6 and A4.It can be observed that, in the case of modification of the Per3, Eng4, Eng5, Bat7 and Bat8 criteria's weights, the alternative A25 remains unchangeably in the lead position, and alternative A2 receives the second best rank.If criterion Per1 would have a slightly lower weight, i.e., medium low instead of medium, alternative A2 would be better.However, if the weight of this criterion were reduced to very low, A20 would advance to the first position.Similarly, if the weight of Per2 was high or greater, the A20 alternative would outrun the alternative A25.It should be noted, that if medium high to very high weight is assigned to the Pr9 criterion, a significant of the A2 alternative can be observed, from Rank 5 to Rank 1.Other significant rank changes include A13's increase from Rank 9 to Rank 2 with the growth of Per1 criterion, A19's drop from Rank 9 to Rank 19 with Per2 growth, A22's advance from Rank 12 to Rank 4 with Per3 growth or A36's drop from Rank 3 to Rank 13 when Eng4 grows.The rankings remain most stable for criteria Eng5 and Bat7, where, respectively, the first six and four positions remain unchanged.The analysis of the obtained ranking allows observing that the best alternatives, A25, A2 and A20, are approximately two-fold better than the worst alternative, A24.Some groups of alternatives with very similar values of CCi can be observed, i.e., A32 and A12, A26 and A21, and A31 and A33 in the middle of the ranking, as well as A19 and A4, and A35 and A10.The similarities of the CCi value in these pairs of alternatives mean that they can be used interchangeably.
A sensitivity analysis of the ranking was performed by changing the weight of a single criterion between very low and very high, while all the remaining criteria had the medium weight assigned.The results of the sensitivity analysis are charted on Figures 6 and A4.It can be observed that, in the case of modification of the Per3, Eng4, Eng5, Bat7 and Bat8 criteria's weights, the alternative A25 remains unchangeably in the lead position, and alternative A2 receives the second best rank.If criterion Per1 would have a slightly lower weight, i.e., medium low instead of medium, alternative A2 would be better.However, if the weight of this criterion were reduced to very low, A20 would advance to the first position.Similarly, if the weight of Per2 was high or greater, the A20 alternative would outrun the alternative A25.It should be noted, that if medium high to very high weight is assigned to the Pr9 criterion, a significant of the A2 alternative can be observed, from Rank 5 to Rank 1.Other significant rank changes include A13's increase from Rank 9 to Rank 2 with the growth of Per1 criterion, A19's drop from Rank 9 to Rank 19 with Per2 growth, A22's advance from Rank 12 to Rank 4 with Per3 growth or A36's drop from Rank 3 to Rank 13 when Eng4 grows.The rankings remain most stable for criteria Eng5 and Bat7, where, respectively, the first six and four positions remain unchanged.

Fuzzy TOPSIS Analysis of the Complete Options Set
In the third step of the analysis, fuzzy TOPSIS method was used to perform the analysis of the complete alternatives set.Again, the existing data were converted to singletons, and the data gaps were replaced with trapezoidal fuzzy numbers generated on the basis of the minimum, maximum, mean and standard deviation values of all known values for each criterion.A fragment of the decision matrix for solving the decision problem between all 36 alternatives is presented in Table A4.The solution ranking is presented in Table 9 and Figure 7.
In the third step of the analysis, fuzzy TOPSIS method was used to perform the analysis of the complete alternatives set.Again, the existing data were converted to singletons, and the data gaps were replaced with trapezoidal fuzzy numbers generated on the basis of the minimum, maximum, mean and standard deviation values of all known values for each criterion.A fragment of the decision matrix for solving the decision problem between all 36 alternatives is presented in Table A4.The solution ranking is presented in Table 9 and Figure 7.
From the analysis of the ranks and CCi values to which all alternatives were evaluated, it can be observed that the leading alternative A6 has a significant dominance over the second-to-best alternative A2.In addition, the first seven alternatives are at least three times more preferred than the worst alternative A16.The alternatives on Positions 4 and 5, as well as those on Positions 12-15 and 24-26 obtained very close CCi values, which means the difference between them in terms of preference is very minute.Subsequently, the various weights possibly assigned to each of the criteria were analysed in a sensitivity analysis (see Figures 8 and A5).The alternative A6 remains the leading one in regard to all criteria but Per2, Bat6 and Pr9.The changes of the weight assigned to Bat7 criterion appear to have the least influence on the position of the top alternatives in the ranking.If the weight of the Pr9 criterion grows from very low to very high, the alternative A2 advances from the Rank 8 to Rank 1.Other considerable changes of positions can be observed for A5 when the Eng5 criterion weight is changed (growth from Rank 27 to Rank 14), A7 when the Bat6 criterion changes from very low to very high (increase from Position 20 to 4) and A27 when the Bat7 criterion changes (change from Rank 21 to Rank 6).It should be noted, that though the alternative A7 undergoes a significant jump of 16 slots on the sensitivity analysis of the rankings for criterion Bat6, it changes only by two slots in the analysis of the Bat7 criterion.Similarly, although the alternative A27 rank changes by 15 steps for Bat7 stability analysis, in the case of Bat6 analysis it changes only by two slots.Therefore, although the criteria Bat6 and Bat7 express similar property of the analysed vehicles, their influence on the final rankings of preference differs.From the analysis of the ranks and CCi values to which all alternatives were evaluated, it can be observed that the leading alternative A6 has a significant dominance over the second-to-best alternative A2.In addition, the first seven alternatives are at least three times more preferred than the worst alternative A16.The alternatives on Positions 4 and 5, as well as those on Positions 12-15 and 24-26 obtained very close CCi values, which means the difference between them in terms of preference is very minute.
Subsequently, the various weights possibly assigned to each of the criteria were analysed in a sensitivity analysis (see Figures 8 and A5).The alternative A6 remains the leading one in regard to all criteria but Per2, Bat6 and Pr9.The changes of the weight assigned to Bat7 criterion appear to have the least influence on the position of the top alternatives in the ranking.If the weight of the Pr9 criterion grows from very low to very high, the alternative A2 advances from the Rank 8 to Rank 1.Other considerable changes of positions can be observed for A5 when the Eng5 criterion weight is changed (growth from Rank 27 to Rank 14), A7 when the Bat6 criterion changes from very low to very high (increase from Position 20 to 4) and A27 when the Bat7 criterion changes (change from Rank 21 to Rank 6).It should be noted, that though the alternative A7 undergoes a significant jump of 16 slots on the sensitivity analysis of the rankings for criterion Bat6, it changes only by two slots in the analysis of the Bat7 criterion.Similarly, although the alternative A27 rank changes by 15 steps for Bat7 stability analysis, in the case of Bat6 analysis it changes only by two slots.Therefore, although the criteria Bat6 and Bat7 express similar property of the analysed vehicles, their influence on the final rankings of preference differs.

Comparison of the Obtained Solutions
The analysis process was concluded by performing a comparison of the obtained solutions.All alternatives were once again divided into two categories: those which had complete data and those which had some gaps in the data.For the alternatives from the former group, the ranking from Section 3.3.1 was compared with the ranks of selected alternatives from Section 3.3.2.For the alternatives from the latter group, the ranking from the first subsection of Section 3.3.2was compared to the ranks of selected alternatives from the second subsection of Section 3.3.2.Because the latter ranking contained all alternatives, two subsets of the results needed to be extracted in order to perform the comparisons.Table 10 demonstrates the alternatives and ranks used in each comparison.
The Pearson correlation coefficient for the rankings from the first group is equal to 0.5879, which implies some correlation between the rankings.This can be visually verified on the chart in Figure 9a.The horizontal axis of the chart represents the ranks obtained in the PROMETHEE II ranking.The vertical axis represents the ranks of the same alternatives, extracted from the fuzzy TOPSIS ranking containing all alternatives, regardless of the data completion.The closer to the line on the chart each alternative is, the more similarly it was evaluated in both rankings.In this case, only alternative A17 has the unchanged eighth position.However, although their exact position in the rankings differs, alternatives A6, A7, A8 and A27 remain the leaders of the both rankings.The differences in position are caused by the fact that different method was used to aggregate information in both rankings.
In the case of the alternatives where complete data were not available, the Pearson correlation coefficient is equal to 0.9876, which implies that the relationship between both rankings is almost linear.This can be visually confirmed on the chart in Figure 9b.
The analysis of Figure 9b allows observing that almost 30% of the alternatives have an unchanged position and that most of the alternatives are either directly on, or very close to the line.Some of the alternatives, such as A2 and A25, A13 and A20 or A26 and A36, have just swapped their very close positions.This is because the CCi values for the alternatives from each pair were very close, regardless of the ranking analysed.Alternatives A2, A13, A20 and A25 are the leading ones, and alternatives A3, A14, A24 and A34 are the worst ones in both rankings.
The results presented above clearly demonstrate that the combination of the certain and uncertain data in the integrated multistep decision approach provides a way to facilitate the decision problem solution process by expanding the set of alternatives with options where the comprehensive set of data is not always available.

Comparison of the Obtained Solutions
The analysis process was concluded by performing a comparison of the obtained solutions.All alternatives were once again divided into two categories: those which had complete data and those which had some gaps in the data.For the alternatives from the former group, the ranking from Section 3.3.1 was compared with the ranks of selected alternatives from Section 3.3.2.For the alternatives from the latter group, the ranking from the first subsection of Section 3.3.2was compared to the ranks of selected alternatives from the second subsection of Section 3.3.2.Because the latter ranking contained all alternatives, two subsets of the results needed to be extracted in order to perform the comparisons.Table 10 demonstrates the alternatives and ranks used in each comparison.
The Pearson correlation coefficient for the rankings from the first group is equal to 0.5879, which implies some correlation between the rankings.This can be visually verified on the chart in Figure 9a.The horizontal axis of the chart represents the ranks obtained in the PROMETHEE II ranking.The vertical axis represents the ranks of the same alternatives, extracted from the fuzzy TOPSIS ranking containing all alternatives, regardless of the data completion.The closer to the line on the chart each alternative is, the more similarly it was evaluated in both rankings.In this case, only alternative A17 has the unchanged eighth position.However, although their exact position in the rankings differs, alternatives A6, A7, A8 and A27 remain the leaders of the both rankings.The differences in position are caused by the fact that different method was used to aggregate information in both rankings.
In the case of the alternatives where complete data were not available, the Pearson correlation coefficient is equal to 0.9876, which implies that the relationship between both rankings is almost linear.This can be visually confirmed on the chart in Figure 9b.
The analysis of Figure 9b allows observing that almost 30% of the alternatives have an unchanged position and that most of the alternatives are either directly on, or very close to the line.Some of the alternatives, such as A2 and A25, A13 and A20 or A26 and A36, have just swapped their very close positions.This is because the CCi values for the alternatives from each pair were very close, regardless of the ranking analysed.Alternatives A2, A13, A20 and A25 are the leading ones, and alternatives A3, A14, A24 and A34 are the worst ones in both rankings.
The results presented above clearly demonstrate that the combination of the certain and uncertain data in the integrated multistep decision approach provides a way to facilitate the decision problem solution process by expanding the set of alternatives with options where the comprehensive set of data is not always available.

Conclusions
This paper has presented a unique MCA-based approach to the evaluation of electric freight vehicles for use in urban areas.The key factor determining the increase in the number of alternative drive vehicles is adapting the infrastructure to meet their specific needs.This relates in particular to electric vehicles for which it is necessary to develop battery charging systems.Due to the discussed impediments, for alternative drive vehicles to be used in urban deliveries, it is first necessary to take specific actions and implement adjustment processes.However, on the other hand, benefits embodied by a reduction of negative impacts on urban environment make alternative drive vehicles increasingly attractive, all the more that their operational parameters make it already possible to use them effectively.
Undoubtedly, implementation of electric vehicles is the most problematic.For urban freight transport to use electric vehicles and function properly, it is necessary to provide integrated networks of charging terminals in cities, so that all places with limited access to power supply sources are eliminated.The report by Frost and Sullivan, entitled "The strategic analysis of EV charging stations infrastructure in Europe", says that there has recently been a significant progress in that regard, and, in accordance with the forecasts, the number of public chargers provided on the initiative of the governments will exceed two million by 2017, while seven West European countries will spend ca.EUR 700 million over the next seven years [47].The aforementioned report additionally emphasizes the possibilities of developing inductive systems of vehicle battery charging, which may play a key role for freight transport, providing possibilities of rapid charging during short stops connected with deliveries.Therefore, it seems that, as a result of the requirements of the European Union regarding sustainable transport development and the worldwide trends connected with popularization of alternative sources of energy, in the near future the dynamics of urban logistics development based on using alternatively powered vehicles, in particular using electric motors, will be growing and playing an increasingly more important role.
The results introduced in the paper are the first part of the work, realized in Poland.The second stage will focus on utilization of the model in the real business environment.However, nowadays the major obstacle for that is the very low level of development of electric vehicles in Poland, especially taking to the account the UFT systems.Nevertheless, a distinct and steady growth can be observed in the number of electric vehicles and plug-in hybrids being sold in Poland.Most car manufacturers offer electric or hybrid vehicles.The current state of e-mobility in Poland resembles the one observed in Europe 5-6 years ago.Therefore, it is expected that, in the years to come, if sufficiently favourable conditions exist, the electric vehicles sales market in Poland will improve significantly.Moreover, Szczecin Municipality is strongly focused on development of electro-mobility at the city area.Due to that, the model will be introduced under present and future projects and activities realized or supported by the Municipality.
The authors' contribution in this paper was to present a framework to solve MCA problems based on a set of alternatives with both certain and uncertain data, by creating an integrated multistep model reflecting decision maker's preferences.In practical terms, an evaluation of 36 carefully selected Electric Freight Vehicles (EFVs) was performed based on nine indicators.From the set of alternatives, a complete set of measurements data was available only for 10 alternatives.For the remaining set of alternatives, the methods of fuzzy set theory were used to provide uncertain values for the data gaps.A comparative analytical study of the results from rankings based on certain data, uncertain data and mixed data was performed.A verification of correlation between the solutions obtained with the usage of certain and uncertain data was performed.
The presented above methodological contribution included the following highlights:

•
We propose a unique approach that aggregates certain and uncertain data in an integrated multistep model.

•
Following the multi-criteria approach, we incorporate the robustness and sensitivity analyses to the evaluation model in order to achieve an improvement of the practical application of the presented approach.

•
We perform a verification of successful usage of MCA-based methods in the problem of EFVs evaluation.
During the research, some possible areas of improvements of proposed approach and future work directions were identified.The presented model used catalogue information for the battery charging time.However, in "real life" situations, the actual battery charging time can vary as the operating conditions (temperature, charging or discharging current, state of charge, etc.) change, therefore it is important to address this problem in future research.One of the possible approaches of solving this problem would be to use the battery models of each car's battery as a parameter of the MCA model.Moreover, the study was oriented on freight vehicles, while the model could be extended to accommodate research on passenger electric vehicles.Additionally, the utilization of the TOPSIS method provides the researchers with a model of a positive ideal solution.This could be used to form a set of guidelines for creating urban EFVs.Finally, the MCA methods used to create the proposed framework are vulnerable to the rank reversal phenomenon.It would be beneficial to construct a similar model based on methods that are free of this phenomenon, such as the COMET method [91].This is important because the rank reversal problem remains an essential challenge in MCA [92], and even more important in the considered problem, where it should be assumed that, along with the lapse of time, the order of the alternatives would vary because of the changes in the evaluated alternatives' set (commercially available EFVs).

Table 7 .Figure 2 .
Figure 2. Comparison of the rankings based on all nine criteria and on a subset of criteria with the "battery charging time 80%" criterion excluded.

Figure 2 .
Figure 2. Comparison of the rankings based on all nine criteria and on a subset of criteria with the "battery charging time 80%" criterion excluded.

Figure 3 .
Figure 3. GAIA (Geometrical Analysis for Interactive Aid) analysis, broken down into: (a) criteria; and (b) criteria groups.

Figure 3 .
Figure 3. GAIA (Geometrical Analysis for Interactive Aid) analysis, broken down into: (a) criteria; and (b) criteria groups.

Figure 5 .
Figure 5. Values of the closeness coefficients (CCi) for the ranking based on the subset of alternatives with data gaps.

Figure 6 .
Figure 6.Sensitivity analysis of the rankings based on the subset of alternatives with data gaps.Fuzzy TOPSIS Analysis of the Complete Options Set

Figure 5 .
Figure 5. Values of the closeness coefficients (CCi) for the ranking based on the subset of alternatives with data gaps.

Sustainability 2017, 9 , 1453 19 of 34 Figure 5 .
Figure 5. Values of the closeness coefficients (CCi) for the ranking based on the subset of alternatives with data gaps.

Figure 6 .
Figure 6.Sensitivity analysis of the rankings based on the subset of alternatives with data gaps.

Figure 6 .
Figure 6.Sensitivity analysis of the rankings based on the subset of alternatives with data gaps.

Figure 7 .
Figure 7. Values of the closeness coefficients (CCi) for the ranking based on the full set of alternatives.

Figure 7 .
Figure 7. Values of the closeness coefficients (CCi) for the ranking based on the full set of alternatives.

Figure 8 .
Figure 8. Sensitivity analysis of the rankings based on the complete set of alternatives.

Figure 8 .
Figure 8. Sensitivity analysis of the rankings based on the complete set of alternatives.

Figure 9 .
Figure 9. Visual comparison of the rankings: (a) only alternatives without gaps; and (b) only alternatives with gaps.

Figure 9 .
Figure 9. Visual comparison of the rankings: (a) only alternatives without gaps; and (b) only alternatives with gaps.

Figure A1 .
Figure A1.Trapezoidal Fuzzy Number (TFN) representation of the linguistic weights of the criteria.

Figure A2 .
Figure A2.The structure of criteria broken down into groups.

Figure A4 .
Figure A4.Sensitivity analysis of the rankings based on the subset of alternatives with data gaps.Figure A4.Sensitivity analysis of the rankings based on the subset of alternatives with data gaps.

Figure A4 .
Figure A4.Sensitivity analysis of the rankings based on the subset of alternatives with data gaps.Figure A4.Sensitivity analysis of the rankings based on the subset of alternatives with data gaps.

Figure A5 .
Figure A5.Sensitivity analysis of the rankings based on the complete set of alternatives.

Figure A5 .
Figure A5.Sensitivity analysis of the rankings based on the complete set of alternatives.

Table 1 .
Overview of selected EU projects.
analysed Optimal Siting of Charging Stations for Electric Vehicles Based on Fuzzy Delphi and Hybrid Multi-Criteria Decision Making Approaches (fuzzy Grey relation analysis and fuzzy VIKOR methods) from an Extended Sustainability Perspective.Ki et al. used AHP method in the problem of Spatial evaluation of Carsharing Stations [79].MA et al. built a Multi-Objective Optimization Model Considering Battery Degradation for a Multi-Mode Power-Split Electric Vehicle [80].Cai et al. proposed hybrid MCA with life-cycle analysis framework to identify and evaluate sustainable strategies of taxi fleet in Beijing in terms of economic, policy, and environmental implications

Table 2 .
The structure of criteria broken down into groups.

Table 3 .
Table of criteria efficiencies for options with data gaps.

Table 5 .
The ranking and options' efficiencies obtained for the alternatives where complete data were available.

Table 6 .
Rankings and Phi net values obtained from the efficiency data with the battery capacity reduced by 10%, 20% and 30%.

Table 7 .
Ranking and Phi net values obtained from the efficiency data with the "battery charging time 80%" criterion excluded.

Table 6 .
Rankings and Phi net values obtained from the efficiency data with the battery capacity reduced by 10%, 20% and 30%.

Table 8 .
Sensitivity analysis.The ranges of stability for the criteria clusters.

Table 8 .
Sensitivity analysis.The ranges of stability for the criteria clusters.

Table 9 .
Rankings based on: (a) exclusively alternatives with gaps; and (b) all alternatives.Alternatives providing full data are in bold.

Table 10 .
Comparison of the rankings: (a) only alternatives without gaps; and (b) only alternatives with gaps.Alternatives for which the rank in both rankings has not changed are marked in bold.

Table 10 .
Comparison of the rankings: (a) only alternatives without gaps; and (b) only alternatives with gaps.Alternatives for which the rank in both rankings has not changed are marked in bold.
(a) Data without Gaps (b)

Table A1 .
Table of criteria efficiencies for options with no data gaps.

Table A2 .
Table of criteria efficiencies for options with data gaps.