Road safety assessment and risks prioritization using an integrated SWARA and MARCOS approach under spherical fuzzy environment

There are a lot of elements that make road safety assessment situations unpredictable and hard to understand. This could put people's lives in danger, hurt the mental health of a society, and cause permanent financial and human losses. Due to the ambiguity and uncertainty of the risk assessment process, a multi-criteria decision-making technique for dealing with complex systems that involves choosing one of many options is an important strategy of assessing road safety. In this study, an integrated stepwise weight assessment ratio analysis (SWARA) with measurement of alternatives and ranking according to compromise solution (MARCOS) approach under a spherical fuzzy (SF) set was considered. Then, the proposed methodology was applied to develop the approach of failure mode and effect analysis (FMEA) for rural roads in Cosenza, southern Italy. Also, the results of modified FMEA by SF-SWARA-MARCOS were compared with the results of conventional FMEA. The risk score results demonstrated that the source of risk (human) plays a significant role in crashes compared to other sources of risk. The two risks, including landslides and floods, had the lowest values among the factors affecting rural road safety in Calabria, respectively. The correlation between scenario outcomes and main ranking orders in weight values was also investigated. This study was done in line with the goals of sustainable development and the goal of sustainable mobility, which was to find risks and lower the number of accidents on the road. As a result, it is thus essential to reconsider laws and measures necessary to reduce human risks on the regional road network of Calabria to improve road safety.


Introduction
Failure mode and effect analysis (FMEA) is a common risk analysis method used to find and fix possible failures, problems, and mistakes in a process, system, design, or service in order to make it safer and more reliable [1]. System safety is enhanced when FMEA aims to recognize all potential risk factors in a system, prioritize the potential risk factors, and take corrective actions to eliminate or decrease the high-risk ones. FMEA unlike other critical analysis techniques is a preventive tool that discover a solution before failure occurs [2]. Therefore, it can prepare valuable information for risk managers, guide them to modify existing plans, take estimates to decrease the possibility of failure, and avoid dangerous accidents. A crossfunctional team should be formed first before analyzing a particular system or product. This will allow an FMEA to be carried out. Then, FMEA team members should determine all potential product or system failure modes [3].
The risk factors prioritization is evaluated with Risk Priority Number (RPN) by employing Severity (S), Occurrence (O), and Detection (D) [4], which is defined as Eq. (1). Each parameter takes values as 1 lowest and 10 highest.
Risk factors with higher RPN are more significant and will require incite activities to avoid or relieve probable risks [5]. Though a systematic way to rank risk assessments, this conventional FMEA logic has some disadvantages. The following shortcomings can be attributed to conventional FMEA: (A) Apart from triple risk factors (S, O, and D), additional parameters like Cost (C) that have an impact on risk assessment have not been considered [6,7]. (B) Weights of S, O, and D are not considered when computing RPN in conventional FMEA [8,9]. (C) It is challenging to examine precisely the S, O, and D parameters due to their subjective evaluation on a 1-10 scale. Using linguistic terms in fuzzy numbers, the FMEA can be better driven [6,9,10]. (D) The RPN elements have duplicate numbers. Sometimes two or more risk factors may have same RPN scores. Conventional FMEA does not thoroughly prioritize risk assessments and confuses risk management decision-makers [11].
In this study, an adjusted FMEA model based on the spherical fuzzy extension of Measurement of Alternatives and Ranking according to COmpromise Solution (SF-MARCOS) is proposed to overcome the conventional RPN assessment limitations style.
The fuzzy set theory was introduced by Lotfi [12] to overcome ambiguity and uncertainty in difficult situations to determine whether the crisp state is true or false, which is not a clear frontier of information. Since then, fuzzy sets have been used in many research fields [13,14], and the theory has been extended to have more reliable results. SFS has recently been utilized in several applications, such as landfill site selection for medical waste [15], waste disposal location selection [16], and design evaluation and technology of a linear delta robot [17]. Ashraf et al. [18] developed and applied spherical fuzzy operators to Dombi aggregation to solve group decision-making problems and used spherical fuzzy to detect the cov-19. Also, spherical fuzzy sets are applied in multi-criteria issues. Kutlu Gündogdu and Kahraman [16] presented the spherical fuzzy VIKOR (SF-VIKOR) method to show spherical fuzzy applicability. They also extended the conventional WAS-PAS method to the spherical fuzzy WASPAS (SF-WAS-PAS) method [19]. Jafarzadeh Ghoushchi et al. [11] applied spherical fuzzy MOORA to prioritize circular economy implementation barriers in the designing of sustainable medical waste management systems. Boltürk [20] applied spherical fuzzy TOPSIS and compared the results with the neutrosophic TOPSIS method. Gündogdu and Kahraman [17] proposed a spherical fuzzy QFD (SF-QFD) method that includes linguistic assessment under certainty and uncertainty.
Nevertheless, the MARCOS technique has not been merged with the SWARA method based on the FMEA technique within the context of SFSs, though SFSs have been proven as one of the valuable tools to handle the uncertainty and vagueness that occur in real-life concerns. Consequently, the present study focuses on SFSs. This type of fuzzy set eliminates some aspects of neutrosophic sets and Pythagorean fuzzy sets by requiring the sum of membership, non-membership, and hesitancy degrees not to exceed one. It does not disregard independent hesitancy, in contrast to Pythagorean fuzzy sets. Even though many researchers have focused on the SFS context, no one has looked at the serious risks that could be present in rural road prioritization problems in a spherical fuzzy environment. Existing literature shows that there is a need to evaluate the critical potential risk factors in rural roads because the evaluation of critical potential risk is not only a vital issue but also an uncertain subject [21][22][23][24]. The rural roads' critical potential risk factors selection process includes many objective and subjective attributes with conflicting goals. Thus, the above problem requires a systematic and suitable approach to evaluate rural roads' critical potential risk factors. To address this concern, an integrated SF-SWARA-MARCOS method is developed that allows DMs to specify a membership function on a spherical area to generalize other fuzzy set components and assign the membership performance parameters independently of the larger domain. Hence, the developed approach can provide more relevant and accurate results by utilizing the advantage of the SFS set, which reflects uncertainty more appropriately and is equivalent to judgments made by decision-makers in assessing potential risks on rural roads. The main contributions of this study are as follows: (A) Provide an FMEA based SF-SWARA-MARCOS integration framework (B) Consider Cost (C) as an essential management indicator with SOD factors (C) Using the SFS-based SWARA to evaluate criteria weights (D) Presenting and evaluating a real problem in prioritizing the existing potential severe risk factors in rural roads of the Calabria region (Italy) using the proposed method SF-SWARA-MARCOS and showing the application of the introduced method (E) Sensitivity and comparison analysis to validate and reveal the usefulness of the suggested method The rest of this research is structured as follows: In Sect. 2, several researches are investigated in three subsections as FMEA model and methods based on it, applications of SWARA and MARCOS methods separately. In Sect. 3, additional descriptions of Spherical fuzzy sets and SF-SWARA are provided. In Sect. 4, extended Spherical fuzzy MARCOS (SF-MARCOS) is introduced. Then, the suggested approach to this study is presented. Finally, a case study is introduced, and the results from implementing the suggested approach are analyzed in Sect. 4. In the end, the conclusions and evolution suggestions of this research are presented in Sect. 5.
2 Literature review 2.1 Hybrid FMEA approach Failure mode and effect analysis (FMEA) is a method widely utilized in different industries to detect potential failures. Reliability analysis is increasing in today's complex systems to find and prevent probable failures to improve the reliability of products, processes, safety, and systems. Apart from different reliability management tools, the main goal of the failure mode and effect analysis (FMEA) is to identify and prioritize the potential failure modes before a failure arises [25]. Failure modes with a higher RPN are more important and need to be fixed quickly to avoid or reduce possible risks [5]. FMEA has been used in various fields due to its application. The tunnel construction process is a very costly and risky project built in uncertain conditions, and there may be risks in the tunnel construction process. Therefore, Amini and Mojtaba [26] have used the FMEA technique to identify, classify, and analyze potential hazards in the tunnel construction process. The vehicle, the environment, and the car's driver are the three elements associated with driving, and the failure of each leads to an accident. By examining potential errors, effects, and causes, solutions can be offered to prevent them. For this reason, the FMAE technique has been used to analyze driving errors [27]. The traditional RPN index used in this method suffers from significant drawbacks, such as allocating the same weight to different risk factors and an absence of complete ranking in the risk evaluation procedure [9]. Several multiple criteria decision-making (MCDM)-based FMEA have been proposed [28]. Moreover, the literature has made many efforts to improve FMEA using fuzzy inference. FMEA has been used to identify significant human errors, causes, and effects during Computer numerical control (CNC) machining operations. The Analytic Hierarchy Process (AHP) method is used to calculate the weight of risk factors, and the Measurement of Alternatives and Ranking according to the compromise Solution (MARCOS) method is used to rank the risks. In addition, a Fuzzy environment has been used to reduce the uncertainty in the data [29].
Also, due to the shortcomings of the FMEA method, Alvand, Mirhosseini et al. [30] have used the integrated Stepwise Weight Assessment Ratio Analysis (SWARA) and Weighted Aggregated Sum Product Assessment (WASPAS) approaches in a fuzzy environment in order to identify and rank the risks in construction projects. Finally, the results show that the combination of MCDM and FMEA methods is more capable and accurate than conventional FMEA. Also, a combination of MCDM and FMEA methods in a fuzzy environment has been used to identify risks in the automotive spare parts industry. In this study, the Best-Worst Method (BWM) method is used to weight the criteria, and the Multi-Objective Optimization by Ratio Analysis (MOORA) method is used to prioritize the risks [31].

SWARA applications
The stepwise weight assessment ratio analysis (SWARA) approach was introduced by Keršuliene et al. [32], to appointing suitable weights to the criteria. In this method, decision-makers (DMs) have an important role in assessing the criteria weights. At first, each DM prioritizes the criteria according to their significance [33]. The major attribute of the SWARA method is that DMs' views about the importance ratio of the criteria are approximated through the weights find out process. The capability to evaluate experts' opinion about significance ratio of the criteria in the process of their weights calculation is the main element of the SWARA method. Furthermore, SWARA is helpful for gathering data from experts. Moreover, SWARA method is uncomplicated and the experts can comfortably work together [34]. SWARA can be applied in many complicated MCDM problems, such as selection of supplier [35], packaging design [36], selection of machine tool [37], logistic problems [38], and selection of personnel [39,40]. Moreover, SWARA is a beneficial method where issues ranks are already clarified and criteria assessment is not needed anymore, while other criteria weighting methods, like the Analytic Hierarchy Process (AHP), are highly related to the criteria assessments. A SWARA-COPRAS method was recommended by Hashemkhani Zolfani and Bahrami [41] for prioritizing high-tech industries investments. Vojinović et al. [42] applied a novel approach include IMF SWRA method to assess healthcare system. A SWARA-WASPAS approach has been proposed to solve a problem of location selection in turkey [43].
Most of the time, DMs' judgments are not express because of hesitancy and vagueness of their views or lack of information; thus, to manage this issue, researchers have extended MCDM methods in uncertain environments [44,45]. A hybrid fuzzy SWARA-COPRAS approach was proposed by Zarbakhshnia et al. [46], to assess and select the sustainable provider. Ghoushchi et al. [47] applied developed SWARA and CoCoSo methods to evaluate wind turbine failure modes. A Pythagorean fuzzy SWARA and VIKOR framework was recommended in another study to choose the solar panel performance assessment [48]. Ghoushchi et al. [49] proposed a hybrid Z-SWARA and Z-MOORA approach with a developed FMEA method to cope some shortcomings of traditional FMEA. An extended interval-valued PF-SWARA-MULTIMOORA was presented by He et al. [50] to examine the status of tourism community in India. Ghoushchi et al. [15] proposed an integrated SF-SWARA-WASPAS approach to select the suitable location for medical waste landfill site. Rahmati et al. [51] discussed a new methodology combining Z-SWARA and Z-WASPAS techniques with the FMEA method to evaluate the priority of risk factors of financial measurement for production companies' management control systems.

MARCOS applications
MARCOS is a new MCDM method introduced by Stević et al. [52]. MARCOS accelerates decision-making in the literature by solving a wide range of different problems. Defining the relationship between the options and the ideal and counter-ideal degrees as reference points determines the utility functions of the options and obtains an adaptive ranking of the options. Higher efficiency, ease in structuring and optimizing the decision process, more accurate determination of the degree of desirability concerning the reference point, more excellent stability and robustness of the results in terms of changing measurement scales, and no problem of ranking inversion [52] are some of the advantages of MARCOS method compared to other MCDM methods such as WASPAS, SAW, and TOPSIS. The MARCOS method has been used in various fields. In this regard, Stević and Brković [53] used the MARCOS method to evaluate the transportation system of an international transportation company. They used this method to rank 25 drivers based on five criteria. Finally, a sensitivity analysis was performed between MARCOS methods and other methods, which shows the results of the superiority of MARCOS method and the reliability of the ranking. In addition, the MARCOS method has also been used to evaluate battery-powered electric vehicles. In this work, 10 cars have been selected as alternatives and have been ranked based on technical specifications such as price, battery, energy, and allowable load. Sensitivity analysis has also been performed to show the validity and robustness of the results [54]. Chakraborty et al. [55] applied a modified MARCOS technique with the application of D numbers to select suitable supplier. Also, an integrated F-PIPRECIA and F-MARCOS have been used in regional aircraft selection problem [56].
MCDM methods allow experts to select the best option based on multiple criteria, the relative importance of the criteria, qualitative and quantitative information, and the preferences of DMs [57]. On the other hand, the opinions and preferences of experts are usually vague and inaccurate when evaluating and selecting options based on various criteria, so the fuzzy concept is combined with MCDM methods to consider the uncertainty. For these reasons, Stanković et al. [58] have developed the MARCOS method in a fuzzy environment. MARCOS has been used to get rid of ambiguity and uncertainty in a number of different subjects. For example, to evaluate the quality of electronic services in the aviation industry from the perspective of consumers, the AHP and MARCOS methods have been used. Due to the ambiguous nature of e-services evaluation, the fuzzy concept has been used. First, the three criteria of reliability, security, and comprehensibility were weighed by the Fuzzy AHP method, and then, the ranking was done by the Fuzzy MARCOS method [59]. Ecer and Pamucar [60] have used the MARCOS method in an intuitionistic fuzzy environment to evaluate and rank insurance companies in healthcare services.
To date, there are few studies that have investigated the applications and capabilities of MARCOS method in uncertain environments. As a result, based on the literature, there is no study reported to evaluate rural roads' critical potential risk factors by applying FMEA-based methodology in the context of spherical fuzzy. Hence, it motivated us to focus on this topic and propose a modified SF-SWARA-MARCOS approach base on the FMEA technique as a contribution of this study.

Prelimination of spherical fuzzy sets
The concept of spherical fuzzy set (SFS) is one of the latest fuzzy sets proposed by Kutlu Gündogdu and Kahraman [61]. Zadeh [62] introduced Type-2 fuzzy sets; compared to the membership function of type-1 fuzzy sets, which are two-dimensional and fixed, type-2 fuzzy sets have a fuzzy three-dimensional membership function that provides more freedom to the decision-making process [63]. Intuitionistic fuzzy set (IFS) is a two-dimensional that defines both membership and non-membership degrees. Atanassov [64] evolved intuitionistic fuzzy sets (IFSs) with the aid of introducing indices, membership degree (l(x)), and nonmembership degree (t(x)) such that 0 B l(x) ? t(x) B 1. Yager [65] developed the Pythagorean fuzzy set (PFS) such that 0 B l 2 (x) ? t 2 (x) B 1. It is a generalization of IFSs. Finally, the Spherical fuzzy set (SFS) was proposed by Kutlu Gündogdu and Kahraman [61], which is the generalized form of Neutrosophic sets (NSs), PFSs, and picture fuzzy sets (PiFSs). Spherical sets offer an innovative approach to a selection of problems that have historically proven very difficult to resolve through existing other extensions of fuzzy set theory, such as human suppositions, including answers such as yes, no, abstain, and refusal. According to Ashraf et al. [66], the main difference between PFS and SFS is that in SFS, we study a neutral degree, while in PFS, it does not [67,68]. This section presents some of the properties, arithmetic operations, and principles of SFSs.
Definition 1 Spherical Fuzzy Sets (SFS) S of the universe of discourse X is given by: where l S : X ! 0:1 ½ :v S : X ! 0:1 ½ :p S : X ! 0:1 ½ For each x, the numbers l S , v S , and p S represent the degrees of membership, non-membership, and hesitance for every x 2 X in the SFS S, respectively. Also, be two SFSs and k is greater than 0 as a constant number. Basic mathematical operations are defined as given as follows: These rules under the condition k; k 1 ; k 2 [ 0; are valid: Definition 4 The score value (SV) and accuracy function (AF) of the number S ¼{l s :m s :p s } are calculated as follows: Note that: Definition 6 Spherical Weighted Geometric Mean (SWGM) with respect to, w ¼ ðw 1 :w 2 . . .w n Þ:w i 2 ½0:1; P n ði¼1Þ w i ¼ 1; is computed as follows: Neural Computing and Applications (2023) 35:4549-4567 4553 Definition 7 The distance between two spherical fuzzy numbers A and B is calculated as follows: This expression can also be used to calculate the spherical distance between two spherical fuzzy sets.
The normalized spherical distance between A and B is defined as:

Spherical fuzzy SWARA
Weight assignment to criteria is one of the significant stages in MCDM problems. The SWARA method was first introduced by Keršuliene, Zavadskas [32]. In this method, a group of experts express their opinions about the options freely while a researcher takes notes and determines the relative weights by ranking them based on the evaluation of experts [5]. The SF-SWARA is applied in this research to evaluate the weights of FMEA parameters considering the probable uncertainty in the views of DMs. The SF-SWARA steps are given as follows: Step 1: Sort the evaluation criteria in a descending order according to the DM's opinions and linguistic variables (LVs).
Step 2: Convert the LVs expressed by DMs to SFNs using Table 1 and make a SFN decision matrix.
Step 3: Aggregate the preferences of each DM using the SWAM or SWGM operator, as shown in Eqs. (17 (18).
Step 4: Evaluating the comparative significance of score value: It starts from the criteria in the second place, a score between 0 and 1 allocated by DMs to the factor j concerning the last criterion (j -1). After applying this process to all the criteria, the comparative importance of score value (s j ) is obtained.
Step 7: Estimate the relative weights of the evaluation criteria as follows: where w j represents the relative weight of criterion j and n represents criteria number.

Methodology
This study aims to propose a hybrid integrated SF-based decision-making approach for evaluating and prioritizing potential risk factors. Therefore, the proposed approach is discussed in this section.

Problem formulation (spherical fuzzy MARCOS method)
Step 1: Construction of the decision matrix. The first step in all MCDM methods aimed at prioritizing is to set up a decision matrix. In the MARCOS method, n criteria evaluate m alternatives, and each alternative gets a value based on each criterion.
. . .; b m f g specify our alternatives and C n ¼ c 1 ; c 2 ; . . .; c n f gspecify the criteria. Therefore, the decision Step 2: Conversion of linguistic variables to SF numbers.
The decision matrix formed in the first step is converted to SFNs utilizing Table 1.
Step 3: Determining the ideal and the counter-ideal.
In this step, ideal (A ai ) and anti-ideal (A id ) values are find out based on Eqs. (26) and (27), where B signifies the benefit aspect criteria, and C signifies the criteria with a cost aspect.
Step 4: Weighing decision matrix. Form the weighted decision matrix v ij by multiplying the elements of matrix K ij with corresponding criteria weight coefficients based on Eq. (29).
Step 5: Determine the utility degrees (UG) of alternatives.
In this step, the UGs of alternatives are calculated utilizing Eqs. (30) and (31).
Step 6: Determine the utility functions of each alternative and ranking.
Finally, the utility functions of each alternative are calculated using Eq. (32).
where f K À i À Á and f K þ i À Á represent the utility functions of alternatives, calculated according to: Step 7: Determine the ranks of the alternative. Rank the alternatives according to the parameter F i in descending order.

Proposed method
In this section, a systematic approach of extended SWARA and MARCOS methods is proposed to assess the potential risk factors. The proposed approach is presented in three steps as follows: According to Fig. 1, in phase one, the FMEA team identifies potential risk factors, and a value is assigned to each of the four criteria described in the FMEA process (severity, occurrence, detection, and cost). The FMEA team then uses Table 1 and linguistic variables (LVs) to figure out the values of SODC factors for each identified risk factor. So that the lowest and highest scores of SODC factors for each potential risk factor are determined by LVs of Absolutely Low Importance (ALI) and Absolutely More Importance (AMI), respectively. Also, the significance values are shown using a 9-point rating scale from ALI to AMI.
Then, in the second phase, the importance of the criteria (severity, occurrence, detection, and cost) is determined using spherical fuzzy LVs by the experts' team, which consists of 3 experts with experience and expertise in various fields. Using Table 1, the LVs are turned into spherical fuzzy numbers. Then, the preferences of the experts are combined, and the SFSWARA method is used to get the final weights of the SODC factors as a spherical fuzzy number.
In the third phase, based on the outputs obtained in the first and second phases, an attempt is made to prioritize the FMEA team's identified potential risk factors based on the weights of the criteria. In the SF-MARCOS method, after figuring out the decision matrix (the result of the first step), which is made up of spherical fuzzy LVs, Table 1 is used to turn these values into spherical fuzzy numbers. Then, a weighted spherical fuzzy decision matrix is organized using weights of SODC factors (output of the second phase). The SF-MARCOS method is utilized to prioritize the potential risk factors in the following. In prioritizing the potential risk factors based on the score obtained from the proposed approach, potential risk factors with a higher score will be ranked as the first priority. The flowchart in Fig. 1 shows the proposed approach to rank the potential risk factors.

Results and discussion
In recent years, in line with a global goal to reduce road accident casualties, more extensive research and studies have been conducted to increase the level of road safety [69][70][71][72][73][74][75][76]. This target is one of the United Nations Sustainable Development Goals (SDGs) and is connected to wellness and urban living. Extensive and valuable studies have been conducted on road safety and the parameters involved in causing accidents. However, due to the uncertainty in examining these factors, the need for further studies is fully felt. Therefore, considering the roads of the Calabria region in southern Italy as the case study in this research, the methodology proposed in the previous sections is applied to this case. More discussion regarding the risks of the case study will be given in the following section.

Context definition
The approach introduced in the paper is applied to a study context that coincides with the rural road network of the Calabria region. Based on the Piano Regionale dei Trasporti [77], the Calabrian rural road network has an extension of 9066 km, divided as follows (Fig. 2 From an administrative point of view, except for the provincial roads, the other Calabrian rural roads (1640 km) are managed by ANAS, the concessionaire company for managing the road and motorway network of national interest in Calabria.
The A2 Salerno-Reggio Calabria motorway, currently being modernized, is the main road in Calabria; it ensures the region's connections with Italy (and more generally with Europe), absorbs traffic flows in transit through Sicily, and guarantees long-distance connections within the regional territory.
The state roads represent the main axes of the regional road network of national interest and are required to ensure rapid inter-provincial or intra-provincial connections.
The provincial roads constitute the preponderant component of the Calabrian road infrastructure patrimony; they integrate the network constituted by the motorway and the state roads, ensuring the connection of the municipal areas to them.
According to ISTAT [78], in Calabria, from 2001 to 2010 road fatalities fell by 20.2%, less than the national average (-42.0%). Between 2010 and 2020, there were -55.8% and -41.8% variations. In the same period, the mortality rate decreased from 4.1 to 2.9 deaths per 100 accidents, while the national average recorded a slight increase (1.9 to 2.0 deaths per 100 accidents).
In 2020, the most significant number of accidents occurred on urban roads (1286, 61.8% of the total), causing 26 dead (42.6% of the total) and 1904 injured (58.3%). The most serious accidents occur on other roads (4.8 deaths per 100 accidents) and motorways (2.8 deaths per 100 accidents).
Despite the indexes of injury, mortality, and severity decreased in the 2 years 2019-2020 1 (respectively from Fig. 1 The flow of the proposed approach 1 The COVID-19 pandemic and the measures adopted to contain it has led to a significant decrease in the traffic flow volumes on the roads, with a consequent decrease in exposure to the accidents risk. 164.2 to 157.0, from 3.8 to 2.9, and from 2.2 to 1.8), the accident rate remains high on the A2 Salerno-Reggio Calabria motorway (117 accidents, 4 deaths, and 177 injured in 2020) and along the two main coastal roads, especially for the Ionic SS106 state road, along which the highest number of accidents is recorded (200 accidents, 9 deaths, and 335 injuries), and of the Tyrrhenian SS18 state road. In inland areas, the accident rate remains high along the SS280 state road called ''dei Due Mari,'' the SS283 state road called ''delle Terme Luigiane'' and the SS682 state road called ''del Passo della Limina.'' Under these assumptions, intending to assess the risk of an accident and the related components, following the approach introduced in the previous sections, the team of experts identified a set of serious potential risk factors in the rural roads of Calabria above described (Table 2). In fact, this team of experts, in the form of a research group, by identifying and collecting risks that are generally considered in road safety, considered more than 25 risks in the first stage. Then, in several meetings of these experts, among these collected risks, 16 risks were selected as potential risks that could be taken in this case study. It is important to note that the goals of achieving sustainable mobility were taken into account when these risks were first looked at and chosen.

Results
In this section, the results of the proposed approach in assessing serious potential risk factors in rural roads of Calabria are examined.
This study aims to introduce a new approach using MCDM methods in an environment of uncertainty.
According to the first step of the proposed approach, experts first express the importance of each criterion according to Table 3 in linguistic terms using Table 1.
The LVs assigned by the experts are then transformed to SFNs according to Table 1. Utilizing Eqs. (17) and (18) and the experts' weights, which are 0.2, 0.3, and 0.5, respectively, the aggregated decision matrix is created (see Table 4).
The SVs of SFNs are obtained using Eq. (14), and the matrix S is formed in the next step. Then, according to the score function value, the criteria are sorted in descending order, and k j and p j are calculated utilizing Eqs. (22) and   (23). Finally, w j , which is the final weight of the criteria, is calculated by Eq. (24) ( Table 5). The SF-MARCOS method is then applied to prioritize risk factors. According to the first step of this method, the decision matrix is formed by the FMEA team with SF linguistic terms ( Table 6). The FMEA team have been constructed by selecting some experts in the safety and risk management field and some specialists from traffic department.
Then, decision matrices based on SF variables are transformed to SFNs utilizing Table 1. The weights obtained from the SWARA method are applied (see Table 7).
For example, K þ i and K À i for the first risk factor are calculated using Eqs. (30,31) as follows: In the next step, f K þ 1 À Á and f K À 1 À Á values which present utility functions of first risk factor, calculated using Eqs. (34,33), respectively.  According to Table 8, we find that R.F. 7, with a score of 0.291, has a higher priority than other risk factors. R.F. 2 and R.F. 6, with 0.288 and 0.282 scores, are in the second and third priority, and R.F. 15, with score 0.255, is in the last priority. Therefore, based on this prioritization, specialists can take preventive and corrective measures to avoid the negative effects of these potential risk factors. In fact, this ranking, which is based on the opinions of experts, can give a comprehensive picture of all the possible risks that have been looked into for this case study, so that preventive measures can be taken and the situation can be made better. This ranking was independent of the classification of risks based on controlled or uncontrolled risks and was done only based on the experience and expertise of experts who were fully familiar with the conditions of the roads and the environment of the case study.

Validation of results
This section provides some validation tests of the obtained rankings. For the first validation study, a comparative study between the proposed approach's results (modified FMEA by SF-SWARA-MARCOS) and some existing approaches include SF-TOPSIS, SF-MOORA, and conventional FMEA (RPN score) has been performed. According to Table 9, based on conventional FMEA (RPN score), risk factor 2 (R.F. 2) has a higher priority than other risk factors. In addition, R.F. 6 and R.F. 7 are both in second priority. By examining the priority of potential risk factors based on the RPN score, it can be concluded that the ranking of potential risk factors is done so that the potential risk factors are placed in 14 categories instead of 16 categories. This issue indicates that rank based on the conventional index is not entirely done and perplexes decisionmaker in risk management and planning corrective/preventive estimates. In line with the conventional FMEA approach results shown in Table 9, incomplete rank can be because of the lack of allocation of various weights (according to experts) to SODC factors and lack of uncertainty in the values of these factors. The problem of duplicate prioritizes has been solved by applying the SF-SWARA-TOPSIS method and risk factors are placed in 16   been demonstrated that spherical fuzzy sets serve as the theoretical basis for MCDM's greater capability in dealing with ambiguous problems. Thus, more than classical fuzzy sets, it reflects uncertainty in real-world issues, which leads to the project's critical path being determined correctly, and as a result, project planning is closer to reality. By comparing the results obtained by applied approaches, it can be observed from Fig. 4 that SF-MOORA method results have high variations compared to other investigated methods. It seems that this method is not suitable to solve the problem of this study with correlation coefficient as 0.18334. On the other hand, when we compare the results obtained by SF-SWARA-TOPSIS and proposed approach (SF-SWARA-MARCOS), we observe that there are few rank variations between each of them with conventional FMEA. The correlation coefficient is observed as 0.87656 for proposed approach and 0.91982 for SF-SWARA-TOPSIS. It can be obtained from this results that this proposed approach is valid because the detected results of this approach is similar to the results of SF-TOPSIS approach as a reliable method in literature [61,79]. Although the coefficient of SF-SWARA-TOPSIS is higher than the proposed approach of this study, in the detailed examination, the result of this examination is different. The critical risk factors that obtained from conventional FMEA are more close to critical risk factors that obtained from the proposed approach from SF-SWARA-TOPSIS, where R.F. 7 is in the first and second priority (Table 9) in the results of the conventional FMEA and SF-SWARA-MARCOS, but it is in the seventh priority according to SF-SWARA-TOPSIS approach. Although we do detect extreme rank variations between the proposed approach in this study and some existing approaches, it can be declared that the application of this proposed approach is new in the FMEA domain.
A sensitivity analysis was done as a second validation study. Here, the variation in risk factor ranking is analyzed with regard to changes in criterion weights. We exchange the weight values of the modified FMEA parameters for the sensitivity analysis. As there are four risk parameters in the modified FMEA, in our case study, eight combinations (current scenario and seven different scenarios) are created.
The ranking orders of 16 risk factors with reference to the eight weight values are shown in Table 10. It can be observed from Table 10 that when the weight value varies, there are variations in the ranking orders of risks. As a result, our method is sensitive to the weights of the modified FMEA risk parameters. While R.F. 7 and R.F. 2 are ranked as the most crucial risk factors, the least critical risk factor varies depending on the weight values. Based on the comparative evaluation of the results obtained from our proposed approach with related literature studies [61,80,81], we can state that the ranking result obtained by our approach is credible and valid.
The correlation between the scenario results and the main ranking orders with regard to weight values is also examined. According to Fig. 5, SC.7 and SC.4 has highest correlation with the main current scenario results with 0.994 coefficient. Also, it can observed that SC. 6 has the least correlation with the main current scenario results with 0.950 coefficient. As a result of this analysis, it is clearly visible that ranking results based on the seven different weight values have a very high correlation. Figure 5 displays a correlation coefficient more significant than 95% for all the results that show a high and positive correlation.
As previously stated, the risk assessment process in road safety is an essential step in achieving sustainable mobility. Nevertheless, uncertainty is the most undesirable effect affecting the risk assessment process. So, an integrated approach in a spherical fuzzy environment for accurate assessment of risks can be helpful to minimize the uncertainty effects of the risk assessment process. The present study aims to present an integrated approach for risk assessment and prioritizing rural road safety in Calabria, in southern Italy. For this aim, the developed approach of FMEA with the SF-SWARA-MARCOS method was presented in this study. In this regard, numerous risk factors were identified and then consulted with transportation and road safety experts who were familiar with the study area. Then, sixteen risk factors were selected from six risk sources, including driver conditions, vehicle conditions, environmental conditions, engineering and technological conditions, policies, and unpredictable risks. Finally, the developed approach was applied to assess and rank the sixteen risk factors. The following is an explanation of the study's findings. According to the proposed approach, the first source of risk affecting rural road safety in Calabria is the driver's condition. As many as five risk factors associated with driver conditions are at the top of the ranking (Lack of driving skills, Presence of alcohol, medicinal or   [82][83][84][85]. In order to reduce the risks caused by drivers, effective measures should be taken to improve and control the behavior of drivers themselves. At present, the common measures are mainly financial penalties. These measures can include periodic training, control, and monitoring of drivers' behavior at regular or unexpected times. Also, drivers can get training and rules-related reminders all the time by using modern technologies like smartphones, text messages, and training applications. Furthermore, enacting and enforcing legislation on key risk factors could be advisable. Also, it will be important to come up with ways to measure some key performance indicators (KPIs) set by Agenda 2030. These KPIs include speed, the use of protection systems (helmets, seat belts, and child seats), the use of alcohol and drugs, the safety level of the vehicle fleet and the national road network, distractions while driving, and how well rescue systems work in case of an accident.
The last two risk factors had the lowest level of risk (unpredictable risks). The ranking results of these two risks were in full accordance with the conditions prevailing on the roads.
It is also suggested that for further studies, the methodology used in this study should not only be used for other rural roads in Italy and their results should be compared with each other, but it is also suggested to be used on urban roads as well, so that by examining urban road risks, we can use this method to check the safety of roads in urban areas in fuzzy environments and under uncertainty. Although the results are generally in line with national data, it should be emphasized that they are specific, especially concerning the individual factors analyzed; they should only be utilized on rural roads in the Calabria region and not on other rural roads. A more detailed analysis could be carried out by examining other factors that have a more significant impact on the risk of a road accident in urban areas.

Conclusions
This study presented the developed approach of FMEA with SF-SWARA and SF-MARCOS methods. Each method was utilized to cover several shortcomings of the traditional FMEA method in order that after determining the probable flaw scenarios are based on FMEA, SF-SWARA is used to count the weight of factors, and SF-MARCOS is utilized to prioritize potential risk factors. Actually, an integrated SF context decision framework is provided in the weighting phase and the prioritization phase. It was the first time that integrated SF-SWARA-MARCOS methodology have been developed base on the FMEA technique. Given the fact that decision-making is essential in the real world, providing effective approaches is critical. For this purpose, in this study, in line with the goals of sustainable development and sustainable mobility, the proposed method was applied to rural roads in the Calabria region in southern Italy. The results showed that most of the risk factors in which drivers had a direct role had the highest level of risk. This study has identified that ''lack of driving skills'' has the highest priority with 0.291 score among risk factors. Also, ''presence of alcohol, medicinal or recreational drugs'' and ''lack of attention to the longitudinal safety distance'' are in the second and third priority with the scores of 0.288 and 0.282, respectively. The proposed method was verified in two ways. As a consequence of these verifications, it was evident that the rating results were reliable, with adequate levels of accuracy and robustness. Based on this, a set of necessary measures to increase the level of awareness of drivers was proposed. The suggested method might be a specific approach but an effective tool for evaluating risk factors impacting rural road safety and, therefore, beneficial in decision-making under uncertainty.
The main weakness of this study was not considering the relations between criteria. The criteria probably have linear or non-linear relations. We may suggest a novel approach for this. In addition, the number of experts was limited. So, increasing the number of experts will change the results to be more reliable. Also, the expert had limited 9 scale linguistic variables to present their opinions. Therefore, increasing the scale of linguistic variables may give them more freedom and increase certainty and confidence. For future research, this study provides directions: the proposed approach can be directly applied to other types of roads or any other issues that related to risk assessment. Additionally, we recommend new aggregation operators to be extended by applying SFS. Another recommendation is to extend other MCDM methods to their spherical fuzzy environments such as BWM, VIKOR, SECA, and TRUST.
Funding Open access funding provided by Università della Calabria within the CRUI-CARE Agreement.
Data availability All data generated or analyzed during this study are included in this published article.

Declarations
Conflicts of interest The authors declare that they have no conflict of interest or financial conflicts to disclose.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.