Selection of EVs as Tourist and Logistic Means of Transportation in Bosnia and Herzegovina’s Nature Protected Areas Using Z -Number and Rough Set Modeling

Tis research focuses on the use of electric vehicles (EVs) to transport visitors and cargo within Bosnia and Herzegovina’s Kozara National Park. Reduced air pollution and the preservation of natural resources are required to help protect this aerial spa. Together with the expert employees of this NP, the EV that would best suit their needs was chosen. Te process of decision-making combines subjective and objective methods. Employees frst chose the criteria and alternatives and then weighed their importance. On that occasion, Z -numbers were used to include uncertainty in the decision, because it is not always possible to make decisions with complete certainty. Furthermore, the weight of these criteria was determined using the fuzzy PIPRECIA (PIvot Pairwise Relative Criteria Importance Assessment) method. Range ( C 1) became the most important criterion, followed by vehicle cost ( C 2), and the technical specifcations of these EVs were used to compare them. Because these specifcations vary, a rough set was used in which the minimum and maximum EV characteristics were taken based on specifc criteria. To rank the alternatives, the R-CRADIS (Rough Compromise Ranking of Alternatives with Distance to Ideal Solution) method was used. According to the research results, the Mercedes eVito Tourer 90kWh is the highest ranked EV and the validation of the results confrmed these fndings. According to the research results, the Mercedes eVito Tourer 90kWh is the highest ranked EV and the validation of the results confrmed these fndings. Te sensitivity analysis revealed that if criterion C 1 is not as important, the other EVs are ranked higher. Tis research`s methodology has demonstrated fexibility, therefore it is recommended for use in similar research.


Introduction
Protected nature areas and transportation are strongly interconnected in terms of tourism as transportation cannot be separated from the aspect of tourist ecological footprint.In addition to the carbon footprint, transportation has a signifcant infuence on the quality of the tourist's experience [1], but rather than being solely driven by tourist demand, transportation in protected nature areas should be planned to meet clearly-defned management objectives on the matter.Te transportation component ultimately makes up a major part of the calculation of the ecological footprint of tourism, and the energy consumption factored into the transportation calculation will likewise make up a sizable share [2].Te private vehicles cannot always be the principal means of transport after scenic spots and other natural resources are what attract visitors there, not the transportation system.
In terms of the impacts of transport in tourist areas, carbon pollution is one of the most important factors.It is claimed that tourism is directly responsible for 5% of global CO 2 emissions and that mobility accounts for three quarters of that amount, note that while worldwide emissions are beginning to decrease, those related to the tourism sector are increasing [3].As a result, it is important to maintain a sufcient transportation infrastructure to decrease the distance factor with appropriate access quality and reduce journey time [4], as well as to encourage more sustainable solutions.Tere ought to be various systems of transportation, as well as the adoption of more environmentally friendly vehicles.Tis could be achieved with awarenessraising initiatives and the demystifcation of EVs as social status badges of individuals or businesses.
Despite the fact that during the past 40 years, sustainable transportation strategies have become more prevalent in many urban contexts, they have only lately been applied to natural environments [5].Tis move toward ecological means of transportation is critical for nature-based locations as it helps reconcile tourism activities with protected nature area sustainability [6].On the other hand, there are some setbacks.McGinlay et al. [7] reported that managing European protected natural areas in epidemic times was fraught with transportation and trafc challenges such as incidents of negligent parking and trafc volume as a result of the visitor's unwillingness to use organized shuttles or public transportation out of concern that a virus would be spread [5].
Te most common means of transportation, the variables that afect tourist mobility, and an assessment of visitors' perceptions of the issues most signifcant to them when selecting a means of transportation must all be thoroughly examined in order to analyze tourist mobility [3].Sustainable planning and managing must also take into account the "journey to destination" component in addition to minimizing the ecological footprint produced by the visitor's stay at the destination.
Regarding mobility in general and visitors' transfers, the highest ecological footprint values were noticed in transfers that are longer than 100 km, such as from an airport or harbors, especially when numerous pickups are required owing to visitors' individual arrival/departure schedules.On the other hand, a low-mobility ecological footprint was observed where public transportation was promoted over the use of private vehicles.Te common idea is that transportation to and from a location has the highest ecological footprint and that the greatest carbon reduction may be accomplished by transitioning from air to ground travel, notably by train [8], as the least carbon-intensive mode of transportation.
Tis paper is focused on analyzing the mobility of tourists in protected natural areas.By doing so, we intend to determine the best environmental and economic sustainable alternative mode of transportation for tourist purposes in these areas.Decisions concerning the transportation industry are often riddled with disputes since they might result in both benefts and drawbacks for various stakeholders.To achieve sustainable transportation solutions, multiple viewpoints from various stakeholders must be included in the decision-making cycle [9].Stakeholder engagement in decision making is critical as it assists in determining and assessing the objectives among many stakeholders, improving decision acceptance, and boosting the efciency of decisions.
Based on the literature on alternative transportation, available electric vehicles, and protected natural areas policies, this research ofers a selection of EVs as tourist and logistic means of transportation in Bosnia and Herzegovina's nature protected areas.When choosing an EV, the decision maker (DM) has at his disposal several alternatives that are observed through certain criteria.In this sense, selecting an EV falls within the scope of the multicriteria decisionmaking (MCDM) problem and applying the multicriteria analysis (MCDA) approach allows for this problem to be solved.
Te following issue is that while selecting a vehicle, one faces client perception, because each customer has their own vision of how a vehicle should feel [10].In this research, this problem will be addressed in such a way that the technological qualities of the EV will be taken into account.However, these technical parameters of EVs vary depending on charging time and method, vehicle model, and so forth.In order to perform a more comprehensive research, data will be obtained in the form of intervals.Furthermore, group decision-making will be utilized to determine the relevance of the weights of the criteria with uncertainty factored into the decision-making process.Tese problems serve as the foundation for the following research objectives: (1) Applying a novel decision-making model based on Z-numbers and a rough approach.
(2) Performing EV selection for Kozara National Park needs.
(3) Contributing to environmental preservation by application of sustainable tourism.(4) Providing a methodological framework for selecting EV.
To achieve these objectives, a hybrid MCDM model based on the PIvot Pairwise Relative Criteria Importance Assessment (Z-PIPRECIA) and Rough Compromise Ranking of Alternatives with Distance to Ideal Solution (R-CRADIS) methodologies will be utilized.Te Z-PIPRECIA method is an upgrade to the standard PIPRECIA method that uses Z-numbers to help the reduction of uncertainty in decision-making.Te R-CRADIS method is an upgrade to the classical CRADIS method in the rough set.Rough set is used here due to the alternatives' values being expressed as an interval.Tis interval is necessary to include all technical data regarding the observed EVs.
Te paper is organized as follows.Section 2 reviews literature about grand narratives (GNs) for achieving sustainable mobility and ofers a literature overview about diferent EV selection problems with listing various approaches and methods that were used by the authors.Section 3 consists of a research background that closely describes the subject of the conducted research, followed by EV selection methodology developed for NP Kozara's needs.In Section 4, the results of the ranking are presented and then validated by comparison with fve other methods' ranking results.Results of the sensitivity analysis are presented in this section as well.Section 5 discusses the efect of more realistic decisionmaking on Z-number employment, and Section 6, as the conclusion to this research, summarizes the most important fndings, limits, and recommendations for further research.

Literature Review
In comparison to a system based on the use of private vehicles, a sustainable transportation system based on public transportation could be more easily achieved.In order to maintain the ecological and cultural worth of nature tourist spots, public transportation is used to entice visitors to leave their private vehicles outside of these protected nature areas.To create tourist mobility sustainable, the majority of stakeholders in protected natural areas use a number of transportation strategies and policies, including private vehicle restrictions, park and ride initiatives, and public transportation support [9].
It is important to talk about sustainable mobility narratives [11] that include technological and social practices that reduce the need for driving, promote the use of shared transportation instead of private transportation, and support the widespread adoption of low-carbon transportation, especially electric vehicles (EVs) [12] in order to address these issues.
As electromobility, collective transport 2.0, and lowmobility society are proposed as the three GNs for achieving sustainable mobility [11], the transportation industry needs all stakeholders to act immediately and on a scale that has not existed before in order to coordinate and support each other.At all decision-making levels, solutions that recognize the assistance of experts, the developing market, and the support of citizens can function efectively.

GN 1:
Electromobility.Te term "electromobility" refers to the replacement of traditional ICE vehicles by EVs.EVs are available in a variety of forms, including battery (BEV), plug-in hybrid (PHEV), range-extended (REV), and fuel cell (FCEV) [13,14] and including hydrogen hybrid options such as FCHEV or FCPHEV as well [15].Electromobility, as the frst grand narrative, is not just about replacing individual private internal combustion engine (ICE) vehicles, but all current fossil fuel-powered vans, buses, light and heavy-duty vehicles, trains, ships, and short-haul airplanes with equivalent EV drive [11].It is worth mentioning that electromobility neither challenges nor demands that tourists travel less.

GN 2:
Collective Transport 2.0.Public transportation has often been the solution to the problem of sustainable mobility.New categories of collective transportation are proposed by Collective transport 2.0, and this approach involves a change from owning to using a vehicle, which has also been referred to as shared or as "as a service" mobility [11,16,17].For example, in Kruger National Park in South Africa, according to Brett [1], the present infrastructure's capacity is frequently surpassed and an increasing number of tourists require smart and organized transportation alternatives.Te research presented that using shared vehicles, in this case open safari vehicles, had advantages as they could carry 5.9 tourists on average, compared to 2.4 tourists in private vehicles.Te use of autonomous vehicles is also advocated as a transport service innovation in the second grand narrative: Collective Transport 2.0.

GN 3:
Low-Mobility Societies.Te numerous efects on the sustainability of private vehicle use could be observed through increased air and noise pollution, destroyed roadside vegetation, shortage of parking spots, visitor discomfort, growing trafc [9], and overall environmental degradation, and the transportation policy for protected natural areas is required.Te third grand narrative, lowmobility societies, suggests car-free options through gradually reducing, and in the end, eliminating the use of private vehicles.In the process, it is suggested that low-mobility societies should travel on fewer and shorter trips and that the preferred choice of transportation, when needed, should be an EV alternative.
As an overview about diferent EV selection problems, various approaches and methods will be listed hereinafter (Table 1).
Te fact that various approaches were utilized to obtain the selection of diferent EVs demonstrates a weakness to the decision-making process.Tus, this particular research proposes a novel decision-making model based on Z-numbers and a rough approach presented in the following section.

Methodology
3.1.Research Background.Te amount of protected areas in Europe, both on land and at sea, has been steadily expanding since the turn of the 20th century.As of 2020, the 38 European Economic Area member nations have roughly 130,000 protected areas [42].Tese locations span roughly 1.27 million km 2 (22.7%) of the land and approximately 570,000 km 2 (8%) of water, which ranks Europe as the continent with the highest protected area coverage rates.Moreover, the new EU Biodiversity Strategy 2030 [43] mandates EU member states to keep expanding conservation eforts to protect 30% of both land and water coverage by 2030, with 10% of them strictly protected [42][43][44].On a global scale, the most popular NPs attract millions of visitors each year, which has economic advantages but also has environmental consequences [45].
In Bosnia and Herzegovina, 37 nature areas are classifed as protected.Te protected area covers 105,602.18hectares or less than 3% of the entire land area of the country [46].In Bosnia and Herzegovina, two areas have been designated as strict nature reserves (SNRs): SNR Prašuma Janj and SNR Prašuma Lom.Furthermore, the country is home to a total of four national parks, including National Park (NP) Kozara, NP Sutjeska, NP Una, and the most recently established NP Discrete Dynamics in Nature and Society  [20] Advanced battery technologies selection MCDM-WPM Babar et al. [21] Viability of EV in the current market Fuzzy SWOT, fuzzy LP Büyüközkan and Uztürk [22] Sustainable urban logistics SAW, VIKOR Hamurcu and Eren [23] Electric bus selection AHP, TOPSIS Khan et al. [24] HEV selection TOPSIS Ziemba [25] EV selection for government use PROSA-C, Monte Carlo Ecer [26] Comprehensively assessing BEV alternatives SECA, MARCOS, MAIRCA, COCOSO, ARAS, COPRAS Khan and Ali [27] Smart waste management adoption framework Fuzzy SWARA, fuzzy VIKTOR Ren et al. [28] Selection strategies for BEVs LDA, DEMATEL, DANP, VIKTOR Ziemba [29] Selection of city and compact EVs NEAT F-PROMETHEE Ziemba [30] Analysis and recommendation of the EV Monte Carlo, fuzzy TOPSIS, fuzzy SAW, NEAT F-PROMETHEE II Hamurcu and Eren [31] EV selection for transportation in inner city AHP, GP, TOPSIS He [32] Selection of battery electric bus EWM Oztaysi et al. [33] EV selection problem F-SMART Ozdagoglu et al. [34] Bus selection for intercity transportation PIPRECIA, COPRAS-G Stopka et al. [35] Evaluation of selected passenger EVs AHP Štilić et al. [36] Taxi service EV selection SWARA, MSDM, MABAC Wei and Zhou [37] EV supplier selection BWM, fuzzy VIKOR Ziemba and Gago [38] Selection and analysis of e-scooters PROMETHEE GDSS, GAIA Puška et al. [39] EV selection DNMEREC, DNCRADIS Ba ˛czkiewicz and Wa ˛tróbski [40] EV selection Entropy, standard deviation (SD), CRITIC, Gini coefcient-based, MEREC, statistical variance, CILOS, IDOCRIW, VIKOR Dwivedi and Sharma [41] EV selection Entropy, TOPSIS Drina.In this research, the focus will be on creating a sustainable transportation system based on the selection of EVs for the needs of the NP Kozara.According to the 2013 Law on Nature Protection [47] and Te Spatial Plan of the Special Purpose NP Kozara, the protected area of Mount Kozara is located within the municipalities and towns of Prijedor, Gradiška, and Kozarska Dubica [48,49] (Figure 1).Te majority of the karstic Kozara [50] falls under IUCN Category II, classifying it as a National Park [42].
In 1967, Kozara Mountain's center region was designated a national park in order to preserve its historical, cultural, and ecological qualities.Te NP Kozara spans 3.907,54 acres and is situated in the heart of the aforementioned mountain.Te park's territory is extended in a north-south direction, with a length of 7 km and a width that varies from 3 to 6 km.In NP Kozara, 865 plant species have been identifed [51], which include 117 fungus species, 11 lichens, 80 mosses, and 657 species of vascular plants [52].Tere are 19 rare and endangered plant species among the overall number, such as Northern white straw (Galium boreale), Holly (Ilex aquifolium), Cardamine (Cardamine bulbifera), Telekia (Telekia speciosa), and wild cyclamen (Cyclamen purpurascens).[51].Even though detailed research on the fauna has not yet been conducted, we know that the NP Kozara is home to many autochthonous animal species, as well as many species that pass through it and stay for a short time.Roe deer (Capreolus capreolus), wild boar (Sus scrofa), wolf (Canis lupus), fox (Canis vulpes), badger (Meles meles) [51], and other autochthonous animal species may also be found in the park.Kozara also homes native bird species such as owls (Strigidae), magpies (Pica pica), fnches (Fringilla coelebs), crested tits (Parus cristatus), and nightingales (Luscinia luscinia) [51].
In order to reach NP Kozara, visitors need to follow the main Banja Luka-Prijedor road, and then for the last 12 kilometers, the winding asphalt road will take them to the central parking on the Mrakovica plateau.Te Mrakovica plateau, which sits at an elevation of 806 meters above sea level and dominates the center region of NP Kozara, is where many park visitors gather and where, as a result of its urban characteristics, cultural manifestations are frequently held.Every year, Kozara National Park's great biodiversity and ecological worth draw a large number of visitors, who may very well impact particular ecosystems [52].

Research Methodology.
In all countries, environmental protection is becoming increasingly important.Bosnia and Herzegovina recognises this through undertaking the frst steps of environmental protection.To preserve protected areas of nature, it is vital to reduce pollution and use environmentally friendly means of transportation.EVs have become the best alternative, and it is essential to adopt as many of them as possible for tourist and business transport.As a result, Kozara National Park initiated a campaign to deploy EVs for visitors and business transport in order to protect the NP`s environment.Tis research was conducted in collaboration with NP Kozara employees.Te research methodology is presented in Figure 2.
Te research was conducted in such a way that NP Kozara initially identifed fve of its employees as practice experts who dealt with this problem.Tese employees took part in group decision-making.Te selection of the criteria came frst.To prevent certain personal preferences of workers while selecting a vehicle, the technical parameters of the EVs were used to minimize subjectivity.Out of a total of 15 criteria that were initially identifed as potential use cases, seven criteria, in the viewpoint of the employees, were chosen as the most crucial for the deployment of EVs at Kozara NP (Table 2).Tese criteria outline the fundamental criteria for EVs to meet the needs of NP Kozara.
Te frst criterion that was considered crucial is the EV range.Te EV's range should be as long as possible in order to cover as many kilometers with a single battery charge.Te range, however, varies depending on the temperature fuctuations [25].Terefore, the range that a single vehicle can cover with a single battery charge is measured at both its maximum and minimum values.Vehicle cost is another signifcant criterion.Te cost price of individual vehicle models was used to calculate the vehicle cost [31].Te cost of a vehicle varies depending on its size and equipment, and the values of two versions of the same EV are taken here as it is essential to keep expenses as low as possible in order to obtain a high-quality EV at the lowest possible cost.Te following criterion is vehicle consumption or how much electricity is used for every 100 kilometers travelled.Tis consumption rises or decreases depending on the outside temperature and the driving mode used.According to Stopka et al. [35], the consumption rate is reduced as the temperature rises and the velocity decreases, and vice versa.To lower the cost of utilizing EVs, this consumption must be as low as feasible.As a result, the minimum and maximum consumption for each EV were calculated.Cargo volume was used as the fourth criterium.Given that each traveller has a set amount of luggage, the cargo volume of EVs is crucial for their usage [29,30].To be able to transport as many items as feasible with the EV, the cargo volume should be as large as possible.Because the vehicle comes in many confgurations, the minimum and maximum volumes were taken.Te carrying capacity is the sixth criterion.It should be as high as possible to accommodate as many passengers and belongings as necessary.When determining the carrying capacity, we must keep in mind that the weight of the battery is a challenge for all EVs, therefore the total carrying capacity varies [26].Tis criterium depends on the performance of a certain EV, and the minimum and maximum carrying capacities of specifc EVs were defned.Te fnal two criteria address battery charging using home chargers and charging stations.Because no charging station has yet been established in NP Kozara, how long it takes to charge an EV using home chargers is of key importance.Tis charging cycle should be made as short as possible in order for the vehicles to be used.Te charging time is also afected by the charger and its voltage [23], so as a result, charging using 220 v and 400 v chargers is taken into account here.Fast charging time makes up the fnal criterion.Charging stations are used for fast charging [26] and the installation of such stations is required for both the NP Kozara's vehicles and the private visitors' vehicles.Due to the difering voltages of these charging stations, 50 kw and 100 kw, charging rates at these voltages were taken into account.
Once the criteria were defned, their relevance in the selection of EV had to be decided as the weight of the criterion defned as the signifcance of the individual criteria.Te expert employees evaluated the criteria using linguistic values, comparing the signifcance of the criteria to the frst and last criterion.Tey evaluated the importance of individual criteria regarding the frst criterion, and then in relation to the last criterion.Te value scale was used to calculate these parameters (Table 3).Te value scale is determined for values when a specifc criterion is more important than the frst or fnal criterion (scale of 1-2) and when this criterion is less signifcant than the frst and last criteria (scale of 0-1).Te expert employees then decided the degree of certainty in their judgment after considering the value of the criteria and the level of certainty is assessed using a fve-point scale ranging from very small (VS) to very high (VH).Using these scales, the Z-PIPRECIA method will be used to calculate the weight of the criteria.
Te selection of alternatives, i.e., EVs, was initiated after the criteria were decided upon by the employees.Alternatives were chosen based on their capacity to transport a number of visitors and their cargo volume for transporting goods.Te small van type of EVs, which have a minimum passenger capacity of 8, was used for this reason.Moreover, the availability of these vehicles on the Bosnia and Herzegovina`s market as well as the availability of authorized services for them was taken into consideration.Tis way, 13 alternatives were selected by the employees of NP Kozara, namely, Opel Vivaro-e M50 (A1), Citroen e-Jumpy (A2),    4).Tese values are presented in the form of intervals, and a rough set was used to rank these alternatives and select the one that best meets NP Kozar's needs.To rank the alternatives, the CRADIS method will be applied.
After determining the weights and ranking of the alternatives, the results are evaluated in such a way that the ranking of the alternatives is tested using 5 other methods.Tese methods are Rough SAW (simple additive weighting), Rough ARAS (additive ratio assessment), Rough MABAC (multiattributive border approximation area comparison), Rough MARCOS (measurement of alternatives and ranking according to the compromise solution) and Rough TOPSIS (technique for order performance by similarity to ideal solution).Te aim of the evaluation of the results is to compare the results obtained with diferent methods in order to confrm the results obtained with the rough CRADIS method as each of the methods has its own specifcs by which it ranks the alternatives.Te Rough SAW method determines the ranking based on the value of the weighted normalized decision matrix [53].Te Rough ARAS method ranks alternatives based on the degree of usefulness in relation to the optimal alternative and Rough MABAC ranks the alternatives based on the average value of the alternatives.Rough MARCOS ranks alternatives using a utility function in relation to ideal and anti-ideal solutions and Rough TOPSIS ranks alternatives in relation to ideal and "anti-ideal" solutions.Terefore, it is important to evaluate how these method's characteristics infuence the alternatives' ranking and determine whether this ranking difers signifcantly from the ofered CRADIS method ranking.
Te research results will be followed by a sensitivity analysis with the main objective to determine how much each criterion afects the fnal ranking of the alternatives.Te weights of the alternatives are altered based on the scenarios, the ranking of the alternatives is established, and how sensitive the alternatives are to particular criteria will be shown by the sensitivity analysis's results.

Z-number represents an ordered pair of fuzzy numbers Z � (A, B)
, where A is a fuzzy number that represents the limitation of the variable X, while B is a fuzzy number that represents the reliability of the fuzzy number A [57].In general, Z-numbers could be presented as  Z � (a 1 , a 2 , a 3 ;  w A ), (b 1 , b 2 , b 3 ; w B )}. Te transformation of the Z-number into a classical fuzzy number is performed by applying the following steps [58]: Step 1. Converting a B fuzzy number to a crisp number, Step 2. Adding the weight of B fuzzy number to A fuzzy number, Step 3. Converting the Z-number weights to an ordinary fuzzy number, ( By using these steps, the Z-number will be converted into an ordinary fuzzy number and the other operations will be performed as with ordinary fuzzy numbers.
Te PIPRECIA method was developed by Stanujkić et al. [59].Tis method enables the ability to evaluate criteria without frst ranking them in terms of importance; instead, the value is compared to the frst and fnal criteria.Te fuzzy PIPRECIA method is applied as follows: Step 1. Forming a set of criteria for comparison.
Step 2. Decision makers initially evaluate the criteria by comparing them to the frst criterion and determining the relevance of those other criteria in relation to the frst criterion.Tey next compare the other criteria to the last criterion and determine their relevance in relation to the last criterion.Te value scale shown in Table 3 is used on this occasion as follows: where s r j represents the evaluation of the criteria by the r-th decision maker.In order to generate matrix s r j , it is necessary to calculate the average matrix s r j using the geometric mean.Step 3. Determining the coefcient k j .
Step 4. Determining the fuzzy weight q j .
Step 5. Determining the relative weight of the criterion w j .
Tis method is used for both grading regarding the frst criterion and grading in relation to the last criterion.Te fnal weight of the criteria is calculated by averaging the weights determined with respect to the frst and last criteria.

Rough CRADIS Method.
Te CRADIS method was developed by Puška et al. [60], whereas the Rough CRADIS method was developed by Dordevic et al. [61].Te CRADIS method calculates the deviation of the alternatives from the ideal and anti-ideal alternatives, as well as the deviation of the alternatives from the optimal alternative.Te steps of the Rough CRADIS approach are as follows: Step 1. Formation of initial rough decision matrix.
Step 2. Normalization of the initial rough decision matrix.
Step 3. Weighting of normalized decision matrix.
Step 5. Calculation of ideal and anti-ideal solutions' deviations.
Discrete Dynamics in Nature and Society Step 6. Calculation of individual alternatives' deviations from ideal and anti-ideal solutions.
Step 7. Te utility function calculation in relation to deviations from optimal alternatives. where where ϑ − 0 � minϑ − i .
Step 8. Calculation of the values of the CRADIS method.
Highest value represents the highest ranked alternative.

Results
In order to determine which EV best meets the objectives of NP Kozara, EVs were ranked based on their technical characteristics and using selected criteria.Tese criteria were evaluated by expert employees of NP Kozara in such a way that they assessed which criteria were more and which were less important.Tey based their evaluations in accordance with the scale (Table 3).Firstly, criterion C1 was given a neutral value (NO), and the other criteria were given a value in relation to this frst criterion, as well as the values for other criteria.At the same time, they expressed their level of confdence in their assessment.Uncertainty was included in decision-making in this manner.Following that, they gave the last criterion a neutral rating, while the other criteria were rated in relation to criterion C7.In this case, uncertainty is included (Table 5).Te Z-PIPRECIA method's steps are then carried out.To begin, the linguistic values are transformed into fuzzy numbers, with A representing fuzzy number A and B representing fuzzy number B. To accomplish this, linguistic values must be transformed into fuzzy numbers through the membership function (Table 3 and Table 6).
Te fuzzy number B is defuzzifed once the linguistic values are transformed into fuzzy numbers A and B. Te number is then rooted and that value is multiplied by the value of the fuzzy number A [55].Te average value of the fuzzy number A was calculated to account for each employee's opinion when determining the weight of the criteria.Te initial decision matrix for the fuzzy PIPRECIA method is thus formed in this manner.
When multiplying the fuzzy number A by the value of number B, a table is generated for calculating the PIPRECIA method.Initially, the value s j is computed, representing the starting point for calculating the criteria weights.Tis value is then subtracted from two (2) to derive the value k j .Te value q j is obtained by reusing the value of k j , for the frst criterion.For criterion C2, q j is calculated by dividing the previous criterion's (C1) q j value by its corresponding k j value.Similarly, this process is applied to all criteria to obtain q j values, which are subsequently summed up.Te weight value w j is determined by dividing the individual q j values by the total sum of these values.Tis same procedure is also conducted in reverse for the second part of Table 7.Here, criterion C7 serves as the primary criterion, and the remaining criteria are calculated based on it.
Weight values for individual criteria were obtained by transforming fuzzy numbers into crisp numbers using the defuzzifcation process (Table 8).According to NP Kozara employee opinion, the most important criterion for selecting an EV is range (C1), followed by vehicle cost (C2), while the fast charge time (C7) criterion has the least importance.Tese results indicate that, in this case, it is of key importance that the vehicle can travel as many kilometers on a single charge as possible while also being as inexpensive as possible.In terms of charging, the employees believe that vehicles will be charged overnight while visitors are not visiting NP Kozara, thus these criteria are not of great importance.
After determining the weights of the criteria, the alternatives are ranked using the Rough CRADIS method.After determining the weights of the criteria, the alternatives are ranked using the Rough CRADIS method.Te rough CRADIS method begins with the creation of a rough decision matrix (Table 4).Te rough decision matrix is then normalized as the following step (Table 9).When performing normalization, the type of criteria, whether beneft or cost, must be taken into account.
Furthermore, the normalized rough decision matrix is weighted.In this step, the elements of the normalized rough decision matrix are multiplied by the appropriate weights.Te next step is to determine the ideal and anti-ideal solutions.Tis involves selecting the ideal value, which corresponds to the maximum value among all alternatives based on all criteria, and the anti-ideal value, which represents the minimum value among alternatives when considering all criteria (expressions 13 and 14).Following this, the deviation from these solutions is calculated by subtracting the values of the weighted normalized decision matrix from the maximum value, i.e., the anti-ideal solution subtracted from the weighted normalized decision matrix (expressions 15 and 16).Subsequently, values are computed, optimal alternatives 10 Discrete Dynamics in Nature and Society are determined, and these derived values are summed for all alternatives, including the optimal ones (expressions 17 and 18).Tis leads to the calculation of the utility function (expressions 19 and 20), and the average values of these functions are subsequently computed (Table 10).
Te fnal results of the Rough CRADIS method are obtained by calculating the average value of the rough numbers (Table 11).Te results revealed that alternative A12-Mercedes eVito Tourer, a 90 kWh model, best satisfes the specifc objectives of the research, while alternative A6-Mercedes eVito Tourer 60 kWh received the lowest ranking.In order to validate these results, they will be examined using other rough methods.
In order to minimize the impact of data normalization on the ranking order, the same normalization process was used in all methods.Tis is because the original form of the MABAC, ARAS, and TOPSIS methods has a diferent   [62] research, which had a large deviation (Figure 3).Te ranking order for alternatives A1 and A10 was diferent when using the rough ARAS technique, but only the rank order has been changed for these alternatives.Te obtained validation results demonstrated that the rough CRADIS ranking did not difer from those of other methods, thus the ranking obtained by this method was confrmed.Following the validation of the research results, a sensitivity analysis was carried out by developing and executing scenarios.Many recent research studies have used this form of sensitivity analysis, i.e., Durmić et al. [63], Božanić et al. [64], Aytekin [65], Jokić et al. [66], Bakir et al. [67], Alosta et al. [68], Pamučar et al. [69], Tešić et al. [70], Ðukić et al. [71], and Stojanović, et al.'s [72] studies.Te scenarios were designed in such a manner that the individual criterium's value decreased by 15% while the weights of the other criteria grew to compensate for the reduction of these individual criteria.As a result, 42 scenarios were developed, along with scenario S0, which represents the rank order of the rough CRADIS method with specifed initial weights.In this manner, the extent to which a certain criterion infuences the ranking of alternatives was tested.Te sensitivity analysis revealed that none of the alternatives are immune to the conducted analysis (Figure 4).Te sensitivity analysis results revealed that the A12 alternative achieved the best results with the highest number of scenarios where it was ranked frst, and as such, it is the preferred option for EV in NP Kozara`s case.Tis alternative demonstrates how the ranking order changed noticeably when the range criterion's weight was changed.As a result, in the S4-S6 scenarios, A12 was ranked 12th.Tus, it was proven that this alternative would not be selected as the best if the maximum required EV range was between 300 and 400 km.Moreover, worth mentioning is the fact that among all alternatives, alternative A7 has the most consistent ranking.Its poorest rank was sixth, but in two of the scenarios, it was ranked frst.

Discussion
Despite having a wealth of biodiversity and countless natural attractions, Bosnia and Herzegovina treats its environment and natural resources improperly [73].For the beneft of future generations, the environment must be protected [74].Since NP Kozara is also an air spa, it is important to protect the air in the area.Internal combustion engines contribute to air pollution by emitting greenhouse gasses and other harmful gasses [75].Tis must be avoided, hence vehicles  [36].For selection of EV NP, Kozara appointed fve expert employees in order to determine the most important criteria for this selection.Tey decided on 7 criteria and defned limitations for possible EV purchases.Tese limitations were in terms of a minimum carrying capacity of 8 people and that they are suitable for carrying cargo as well.As a consequence, only small van types of vehicles were considered and 13 alternatives that are available on Bosnia and Herzegovina`s market were examined.
To avoid absolute subjectivity in selecting an EV, the technical values of these vehicles were combined with NP Kozara employees`evaluation of the criteria.In decisionmaking, Z-numbers are intended to include uncertainty [56], thus Z-numbers are utilized in the combination of these subjective evaluations of the criteria.Integrating uncertainty into the decision-making process is imperative, given the inherent difculty in achieving absolute certainty in one's decisions.Tis strategy fosters a more pragmatic approach to decision-making.Decision makers may be doubtful about the criterium's evaluation, whether it is accurate, or whether it should be diferent.It is therefore necessary to include additional evaluations that will indicate how certain one is of his/her decision [76].Employees must frst determine the signifcance of each criterion in relation to the frst and last criterion, as well as their level of certainty.For the mutual relationship and decision-making certainty, they provided two gradings in the form of linguistic scales.
To calculate the weight of the criteria, i.e., the importance of individual criteria for decision makers, the Z-PIPRECIA method was used, that is, the classic fuzzy PIPRECIA method was extended with Z-numbers.Te PIPRECIA method was used in this case because, unlike the FUCOM (full consistency method), it is not necessary to sort the criteria by importance beforehand [77].Tis contrasts with the FUCOM method, where it is necessary to rank the criteria in order of importance during the frst step.While employing the PIPRECIA method, employees were instructed to compare all criteria with the frst criterion and the last criterion.Te weight value of the criteria was then calculated by averaging these two comparisons [78].Te fact that every employee was given the same importance in decision-making must be emphasized.
Te results obtained based on employees' opinions showed that the most important criterion for them is range (C1), followed by vehicle cost (C2), while the least important criterion is fast charge time (C7).However, with EVs, these two most important criteria are in confict.To have a longer range, EVs need larger capacity batteries, which cost more and raise the total vehicle's price.In order to select an EV that best meets all criteria, these conficting criteria must be balanced.Furthermore, the EV criteria values also vary, particularly when considering the range criterion.Te outside temperature afects the range of electric vehicles, so when it is colder outside, their range is also shorter, and vice versa [18].Likewise, the vehicle consumption criterion is afected by the outside temperature.Since the minimum and maximum values are provided, these criteria are set in the form of an interval.In addition, it is necessary to develop new innovative methods for energy production [79], while the focus should be on renewable energy sources to power these EVs.
Using the Rough CRADIS method, the results indicated that the EV Mercedes eVito Tourer 90 kWh (A12) ranked frst.Tis was due to the fact that this vehicle has the largest battery capacity and thus the longest range, which was the most important criterion.Te validation of the results, which used other rough methods, also confrmed these research results.Te sensitivity analysis, however, revealed that if a lower value of the rank criterion is assigned, alternative A12 is not the best, but it did fnish second to last in three scenarios.Tis analysis revealed that depending on the weight of the criteria, other alternatives can be selected as the best.Consequently, in order to choose the best alternative as objectively as possible, it is necessary to consider which criteria are used and how important they are.

Conclusion
Creating a sustainable transportation system based on electromobility, collective transport 2.0, and low-mobility society is not a myth, but the limitations lie in the fact that the majority of stakeholders in protected areas, such as NP Kozara, must develop a number of strategies and policies that would restrict private vehicle usage and public vehicles support to create a sustainable tourist mobility environment and to reduce the ecological footprint of tourism.
Following the grand narrative, this research ofered a contribution to environmental preservation by the application of sustainable tourism through a methodological framework for selecting EVs for the needs of NP Kozara.As in the process of vehicle selection, decision makers could have their own vision on how a vehicle should feel, this research addressed the problem by taking into account technical and measurable parameters.Te objectives of the research were achieved by using a hybrid MCDM model based on Z-PIPRECIA and R-CRADIS, with the goal of reduction of uncertainty in decision-making and by using the interval sets to express the alternative values needed to present EVs technical data.
Te results showed that alternative A12, Mercedes eVito Tourer, a 90 kWh model, was the best ranked among the thirteen alternatives, owing primarily to the battery capacity, therefore this alternative's range was the greatest when compared to others.As a result of the outcomes of the analysis carried out using Z-PIPRECIA, the highest weight was given to this criterion.With the notion that the sensitivity analysis revealed that if this criterion is not the most important, then alternative A12 is also not the best in terms of ranking.Te criteria that were used are where this study's limitations lie.However, these criteria were chosen by NP Kozara employees because they are the most important to them.In future research, other criteria should be used to see how they afect the fnal decision.Furthermore, the research was limited to only 13 alternatives as only those met the expert employees' set limitations.Future research must expand the number of these alternatives, especially as the number of EVs on the market grows by day.Additionally, there are limitations in terms of the number of criteria utilized in the study.For future research, it is essential to incorporate subjective evaluations from experts to assess these electric vehicles.However, in order to fully exploit the advantages of EVs in tourist transportation, future research papers should also delve into researching chargers that generate electric energy from renewable sources.Tis approach would serve to mitigate the adverse environmental impact, preserving nature for forthcoming generations.
Tis research presented a novel EV selection methodology.Te methodology demonstrated great fexibility, therefore it is recommended for use in similar decisionmaking problems.Moreover, the obtained results demonstrated that the R-CRADIS method results do not difer from the results obtained by other methods, implying that it can be used in future research.During the validation process, it was also discovered that the rough TOPSIS method's ranking order deviates the most from the other ways utilized in this research; thus, future research could address the cause for this discrepancy and attempt to make the outcomes of this method closer to the results of other methods.

Table 1 :
EV selection problems, an overview of various approaches and methods.

Table 2 :
Criteria for EV selection.

Table 6 :
Fuzzy number B membership function.

Table 7 :
Calculation of the Z-PIPRECIA method.

Table 11 :
Final ranking of the alternatives.

Table 10 :
Ideal and anti-ideal solution and utility function.needs of NP Kozara.Even though it is important to encourage the electric vehicle purchase decision in Bosnia and Herzegovina, EVs are not yet widely available in BiH, thus only about 50 EVs were imported in 2021