A Fuzzy AHP-MAIRCA Model for Overtourism Assessment: The Case of Malaga Province

: The unexpected outbreak of COVID-19 has placed the tourism sector in an unfavourable position. It has altered the reality of destinations worldwide, changing the paradigm of tourism and showing an unusual scenario within the tourism sector, the undertourism. This paper conﬁrms how this global pandemic has put the Spanish tourism sector at a halt, making it essential to establish public–private response strategies that reactivate it and, in turn, avoid overtourism in the post-COVID-19 scenario. This is done to emphasise the necessity of embracing an adapted systematic approach that entails a realistic set of factors to quantitatively measure the impact of overtourism on the local community. In this sense, a fuzzy modiﬁed hybrid multi-criteria decision-making (MCDM) approach is proposed, combining the analytical hierarchy process (AHP) model, the multi-attribute ideal–real comparative analysis (MAIRCA) model, and the sensitivity analysis. The modiﬁed F-AHP model is used to identify the criteria importance. For the ﬁnal ranking of the considered alternative municipalities, the modiﬁed F-MAIRCA is applied. To validate the applicability of the proposed approach, we designated the province of Malaga (Costa del Sol) as an empirical study. The ﬁndings show that the municipality of Nerja was the most affected by overtourism in the pre-pandemic scenario, and today it is found in undertourism. This calls for transformative decisions of fully engaged tourism policymakers and practitioners to thrive in the post-COVID-19 tourism scenarios in terms of sustainability. ‘Consider the Bayesian BWM and VIKOR methods to discuss the importance of sustainable sport tourism and recommend potential sport tourism attractions.’ ‘Propose a modiﬁed DEMATEL approach to examine the cause–effect relationships among the barriers hindering the tourism industry development and assess the critical barriers.’ the ELECTRE method to facilitate the land-use suitability analysis for rural


Introduction
With over 70% of the borders of destinations closed to international tourism and restrictions on 100% of world destinations [1] and the socio-economic unparalleled impact caused by COVID-19, on a global scale, the pandemic has begun to make its effects tangible, reaching 'zero tourism' in most tourist destinations. The tourism sector is being hit hard by the health crisis with global shutdowns, thereby causing significant social and economic impacts. More than 90% of air traffic, more than 300 million jobs, as well as entire tourist destinations have been paralysed [2]. In this scenario, the estimated international tourist arrivals point to a decrease of between 60% and 80% to add losses to 1.20 trillion US dollars [3]. Strielkowski [4] confirmed that "the impact of the corona COVID-19 pandemic on the tourism industry is going to be devastating and similar but with a higher magnitude to what happened in the case of Severe Acute Respiratory Syndrome (SARS) IN 2002IN -2003". Fuentes et al. [5] acknowledge that the pandemic will cause major economic impacts, leading to the steepest decline since the Great Economic Recession of 2008. The COVID-19 pandemic has put 75 million jobs at risk in travel and tourism [6] as the aviation passenger revenues drop USD 252 billion amidst the global restrictions [7] and dragging international travel to levels of 2012/2014 [6]. The aforementioned impact of the COVID-19 outbreak on tourism has not been enclosed to this, but sustainability is also on the verge of stalling as acclaimed by all,

•
Addressing the necessity for identifying the impact of overtourism on the local community employing decision science methodology, mainly MCDM, due to its consistency and reliability in empirical studies; • Considering an adapted two-module model for overtourism assessment that can positively contribute towards a reliable decision within the frame of sustainability of destination; • Filling the gaps by proposing an adapted approach that encompasses a fuzzy hybrid model, AHP-MAIRCA, to assist the tourism policymakers to identify overtourism and make optimal decisions in the future; • Empirically verifying the novel integrated model employing real data taken from experts working in the tourism industry in Malaga province.
With that aim, the paper is structured as follows: a comprehensive literature review of relevant existing studies is given in Section 2. In Section 3, a fuzzy hybrid model, combining AHP and MAIRCA model customised to the overtourism evaluation, is proposed. Section 4 Sustainability 2021, 13, 6394 3 of 32 entails a presentation of the numerical application, whereas Section 5 presents a discussion, sensitivity and comparative analysis. Finally, Section 6 presents conclusions, research boundaries, and directions of future studies.

Review of Prior Methodologies Used in the Tourism Industry
As this paper underpins tourism concerns, diverse methodologies have been employed in the tourism industry towards sustainable tourism management. The domain of Multi-Criteria Decision-Making (MCDM) has been significantly developed, thanks to the vast number of academic and scientific works dealing with the application of the specific decision by adopted models that pertain to the particular field. In addition, MCDM is a convenient methodology for solving complex problems due to its capacity to assess numerous alternatives employing a particular set of criteria. The main reason behind developing MCDM by scholars is to establish mathematical formulas aimed at facilitating the criteria evaluation and the selection of the most adequate alternative. Besides this, MCDM techniques have been applied to smooth the decision-making process in various business areas, having lately gained attention in the domain of tourism as presented in Table 1.
'Propose an integrated case-based reasoning and MCDM system for web-based tourism destination planning.' [25] DEMATEL, NRM, DANP, VIKOR. 'Improve tourism policy implementation by the use of hybrid MCDM model.' [26] TDC, DEMATEL, DANP. 'Adopt strategies for promoting tourism competitiveness using a hybrid MCDM model.' [27] DEMATEL, DANP, VIKOR. 'Improve metro-airport connection service for tourism development by using hybrid MCDM models.' [28] DEMATEL, ANP. 'Propose a hybrid MCDM model for improving the tourism entrepreneurship policy.' [29] Qualitative investigation. 'Analyse the overtourism problem by the use of qualitative investigation among 80 stakeholders in 13 European cities.' [30] Descriptive statistical methods. 'Analyse tourism data in Barcelona to detect the elements behind the popular aversion of tourism.' [31] Fuzzy set theory, DELPHI 'Adopt a fuzzy Delphi method in establishing sustainable ecotourism indicators. The fuzzy set theory for addressing vagueness and uncertainty of judgments.' [32] SAW, TOPSIS. 'Application of MCDM methods to tourism evaluation of cultural sites.' [33] Observation method. 'Identify the weaknesses of tourist accommodation density as a parameter of overtourism in Spain.' [34] Integrative review. 'Examine the overtourism issue between the right of the tourists' and residents' rights.' [35] Descriptive statistical methods.
'Determining the development of local, small-scale, and ideally more sustainable experiences can contribute to economic growth without necessarily increasing numbers of arrivals.' [36] MCDM, UTASTAR. 'Apply a priority-based intuitionistic multiplicative UTASTAR method and its application in low-carbon tourism destination selection.' [37] Bayesian BWM, Rough DEMATEL 'Develop a novel two-stage multi-criteria decision-making (MCDM) model to incorporate the concept of sustainable development into sports tourism.' [38] TODIM, Rough SWARA 'Apply a hybrid MCDM method, which comprises rough SWARA and TODIM, to help policymakers and stakeholders to gain exact data about tourism and travel sector.' [39] SVR-SARIMA; PROMETEE II 'Select the best forecasting model for long-term tourism demand by using hybrid MCDM approach.' Table 1. Cont.

Authors
Methods Abstract [40] Bayesian BWM, VIKOR 'Consider the Bayesian BWM and VIKOR methods to discuss the importance of sustainable sport tourism and recommend potential sport tourism attractions.' [41] 2-Tuple DEMATEL 'Propose a modified DEMATEL approach to examine the cause-effect relationships among the barriers hindering the tourism industry development and assess the critical barriers.' [42] ELECTRE 'Conduct the ELECTRE method to facilitate the land-use suitability analysis for rural tourism.' Heretofore, MCDM techniques have been applied in the tourism area to contribute and facilitate the policymaking process concerning diverse issues. The most common matter that occupied academic attention is the case of the selection of eco-tourism destinations. Hence, other authors investigated the theme of tourists' behaviour by using MCDM methods [43,44]. The implementation of the best tourism planning strategies has been performed by applying MCDM techniques. However, it can be deduced from this review that there is no research that has been conducted on the MCDM approach to assess the overtourism phenomenon and its negative impact on the community. To our knowledge, most of the previous works used traditional approaches to analyse the overtourism issue and did not pay as enough attention to assessing its negative impacts with advanced numerical approaches [45][46][47][48]. In other words, there is a dearth of consistent works, and even developing a systematic model and examining the relative importance of overtourism has been ignored. As a solution to this problem, a need remains for developing a novel methodical model and offering an adapted version of the MCDM approach to suit the examination of overtourism as well as the analysis of its impact on the local community.
The increasing number of MCDM use in the past few decades has pushed scholars to develop and propose a range of various techniques to solve different sorts of complex problems in business. However, we cannot deny the fact that these developed techniques hold pros, cons, and limitations. Thus, it is compulsory to designate an adequate technique for the research problem that is being resolved [49]. The most prominent MCDM techniques are the analytic hierarchy process (AHP) method [50]; the analytic network process (ANP) method [51]; the complex proportional assessment of alternatives (COPRAS) method [52]; the preference ranking organization method for enrichment evaluation (PROMETHEE) method [53]; elimination et choix traduisant la réalité (ELECTRE) method [54]; VlseKriterijumska Optimizcija I Kaompromisno Resenje (VIKOR) method [55]; the technique for order performance by similarity to ideal solutions (TOPSIS) method [56], and so on. Additionally, there is a contemporary set of MCDM techniques that can deal with a broader range of complex problems, such as a new multi-attribute ideal-real comparative analysis (MAIRCA) method [57]; the step-wise weight assessment ratio analysis (SWARA) method [58]; the additive ratio assessment (ARAS) method [59]; the multi-objective optimization on basis of ratio analysis (MOORA) method [60]; the multi-objective optimization by ratio analysis plus the full multiplicative form (MULTIMOORA) method [61], and so forth. It is important to highlight that some of the aforesaid techniques are applied for calculating the weights of the criteria and some of them are used for prioritising the alternatives.
In the actual context of overtourism assessment, the analytic hierarchy process (AHP) is an appropriate technique of 'measurement through pairwise comparisons' and considered as one of the most resilient multiple criteria decision-making techniques [62]. The nature of the AHP method is analytical and allows the decision-makers to turn the abstract form of the problem into a constituent component throughout a 'bottom-up from the alternatives to the criteria' procedure. Due to its practicability, this study adopts the AHP technique to solve the research questions. Some studies highlighted the importance of this technique and described in detail the procedure of its application [63]. On the other hand, the Multi-attribute Ideal-Real Comparative Analysis (MAIRCA) was initially  [64], and it is a novel potential technique that can be softly merged with different MCDM techniques. The main purpose of selecting this method in this study is to determine the gap between ideal/real assessments, sum the gaps for each criterion in the decision matrix, and calculate the total gap for each alternative in the evaluation process. At the end of the evaluation process, the ranking of the alternatives depends on the obtained gap values, and the alternative with the lowest gap value is the best one in the filtration process and vice versa [65]. Advocating the effectiveness of the proposed methodology, the core challenge of the present paper is to give deep insights into assessing the overtourism issue in this exceptional situation of the pandemic and applying corresponding sustainable actions in the future.

Dimensions of Overtourism
Hitherto, overtourism is at a very early stage of research and often allied with the number of tourists' arrival and the destination carrying capacity [66]. In practical terms, overtourism encompasses multidimensional push factors and entails other stakeholders like residents and local businesses. As tourism activity rises in a region, the negative impacts on the local area become evident [67]. This conforms to the definition provided by Mowforth [68]: "the result of carrying capacity measurements will always depend on the context of the situation being measured and that this context will vary not just with the physical and social environments, but also with the values of those asking the questions and establishing the conditions for measurement". This involves that all destinations will know manifold capacities, not just in the context of the balance between different dimensions of carrying capacity but within these dimensions as well. However, at the time of realising this work, the number of international tourists arriving in Spain diminished by 64.3% compared to March 2019. The most touristic places in Spain, the Balearic Islands and Catalonia, registered a decline of 68% while the Canary Islands recorded a drop of over 66% in international arrivals and with a national average of almost 65% [69]. The COVID-19 outbreak has intensively affected one of the vivid sectors for the Spanish economy and had declined the tourism activity by 98% in March 2020 [70]. As a result, the estimated evolution of the output gap of the Spanish economy knows again a negative gap, registering rates close to −10% by 2020, which refers to the sudden change in the cyclical position of the national economy. The evolution of the output gap was gradually closed since 2013, when it reached its minimum, and in 2018 and 2019 reached slightly higher than their potential levels [69]. However, the Spanish government is drafting a budget policy that allows the reconstruction of the country and underpins economic recovery where the policy will lay the basis for economic, social, and environmental reform [69]. From this context, the criteria employed in this study were in the first place determined in line with a detailed comprehensive literature review [25,30,66,[71][72][73][74][75][76][77] and then by 14 tourism experts' approvals. Similarly, Peeters et al.'s [11] study also acknowledges that overtourism evaluation is deemed suitable by considering economic, social, and environmental dimensions. The designation process is categorised into three dimensions and the drivers leading to them are presented in Table 2. Table 2. Description of the set of criteria for overtourism assessment.

Dimension
Criteria Description It shows the number of tourists' arrival that participates more in environmental degradation. High level of carbon dioxide emission and marine pollution. A sustainable tourism strategy should be established to ensure that the region's natural ecological space is not damaged by tourism activities.
It applies the consumption of inorganic waste in touristic events that hinders the environmental and ecological balance of the area. Inorganic waste should be managed and controlled in type and amount.
Natural resources usage C 13 It involves the harsh use of natural resources, which shows some reverse effect. Tourism activities must take preventive measures and be organised in an environmentally friendly mode to ensure the conservation, preservation, and sustainability of the natural resources.

Social C 2
Safety and security C 21 It encompasses the negative impacts of the tourism perceived by the local population like insecurity, noises, alcohol and drug consumption. Thus, the capability of local law to safeguard the safety and security of both citi-zens and tourists is essential.
Accessibility/facilities to touristic sites C 22 This multidimensional aspect entails the improvement of the transportation system, facilitate travellers and residents to arrive at all destinations safely, and permit the local population to access touristic spaces.

Infrastructure degradation C 23
Some existing touristic buildings and facilities often cause damages and degradation to cities' infrastructure. Therefore, sustainable actions of the tourism sector should be implemented and build ecologic constructions to maintain the original look of the region.
Culture identity loss C 24 As the tourism demand inclines, the loss of cultural identity and traditional activities in the local area becomes higher. This would be observed negatively by the local community and may create an intensive feeling of rejection and change the citizens' attitudes towards tourism. Thus, this impact may lead to a loss of community moral values, pride, and a loss of local cultural identity.

Economic C 3
Inflation/low Salaries C 31 Inflation and low salaries are strongly related to tourism activities in the local area. Inflation of prices and decrease of the handiness of some goods, services, and factors of production aimed at citizens. The youth are likely vulnerable to low wages, discrimination in job opportunities, and employment quality.

Infrastructures costs C 32
Potential tourism capital investments are dedicated to infrastructures' restructuration. These expenses create a strong challenge and may cause an imbalance in budget handling. Economic dependence on tourism sector C 33 Rely heavily on the tourism sector to generate incomes more than any other sector and is strongly impacted by seasonality and the degradation of other sectors of employment.

Study Area
Malaga province is a seaside capital in southern Spain's Andalusia autonomous community, also called Costa del Sol (Sun Coast) and it is in 6th place among 50 provinces and 2 autonomous cities by demographic size with a land area of 802.7 km 2 . The province is located in 37 • north latitudes and 5 • west longitudes. This coast has a Mediterranean climate with mild, relatively wet winters and hot and sunny summers. It has a warm climate and hardly any precipitation. The average maximum temperature in Costa del Sol is 28 • C in summer and 14 • C in winter. The water temperature is between 17 • C and 25 • C. Figure 1 describes the cartographic representation of the study area by the use of the software packaging Qgis, version 3.10.6. In 2019 the population of Spain did not reach 46.5 million and the number of international tourists who visited Spain in that year was 83.7 million, generating 12.3% of GDP and 12.7% of employment [10]. According to the Ministry of Industry, Commerce and Tourism [78], the tourism sector became, after the economic crisis (2015-2019), one of the reliable strategic sectors. Andalusía is the first destination in terms of receiving foreigners worldwide. The annual number of visitors exceeds 29.5 million, more than triple the regional population and revenues reach 20.4 billion euros. According to the latest available data, the region has some 546,000 regulated accommodation places, of which around 250,000 are offered in the hotel category [79]. In Malaga, the number of foreign tourists who came to travel exceeded 12 million, which represents an increase of 3.4% compared to 2018 [80]. The capital has some 337,200 regulated accommodation places, of which around 100,395 are offered in the hotel category. Tourism activity has a great weight in Costa del Sol for its contribution to regional GDP, the employment it generates, and its ability to drag on other economic activities [81]. In recent years, the favourite destinations for tourists are the beaches of Malaga province, which represents the maximum exponent totals a third of the hotel beds in the region. The lower pressure of foreign demand favoured the hiring of Spanish tourists due to greater accessibility to the accommodation offer in Spanish destinations on the Spanish coast with the highest presence of tourist packages [82].
In 2020, for the first time after the Second World War, the tourism sector reduced its activity to almost zero overall in a controlled way. Uncertainty has been the protagonist of this unprecedented situation. The movements of people were reduced to their minimum expression, called essential and applied to all scales. The first reports predicted a pessimistic tourist future, with millions in losses, unemployment, destroyed companies, and a reactivation of the sector conditional on the mass production of a valid vaccine. In any case, Navarro et al. [83] declare that the new tourism will go through adaptation not only conceptual and process, but planning as well. The Costa del Sol is divided into two areas: the 'Western Costa del Sol' and the 'Eastern Costa del Sol', the capital and centre of both being the city of Malaga. Based on the tourism activity intensity reported by the statistical registry of Andalusia tourism [84], the municipalities being considered in this study are Fuengirola, Torremolinos, Benalmadena, Malaga capital, Marbella, Estepona, Mijas, and Nerja. Figure 2 describes the cartographic representation by the use of the software packaging Qgis, version 3.10.6 to explain the geographical location of the study area and the densities of the accommodation spaces. pessimistic tourist future, with millions in losses, unemployment, destroyed companies, and a reactivation of the sector conditional on the mass production of a valid vaccine. In any case, Navarro et al. [83] declare that the new tourism will go through adaptation not only conceptual and process, but planning as well. The Costa del Sol is divided into two areas: the 'Western Costa del Sol' and the 'Eastern Costa del Sol', the capital and centre of both being the city of Malaga. Based on the tourism activity intensity reported by the statistical registry of Andalusia tourism [84], the municipalities being considered in this study are Fuengirola, Torremolinos, Benalmadena, Malaga capital, Marbella, Estepona, Mijas, and Nerja. Figure 2 describes the cartographic representation by the use of the software packaging Qgis, version 3.10.6 to explain the geographical location of the study area and the densities of the accommodation spaces.

Methodological Overview
This paper proposes the use of modified F-AHP and F-MAIRCA as an adapted methodology to evaluate overtourism issue of the case of Malaga province and accordingly propose the post-COVID-19 tourism sustainable/potential actions. The hierarchy framework consists of three succeeding steps, which are graphically presented in Figure  3. This framework can help tourism policymakers and manager analysts to identify the proposed set of criteria for overtourism evaluation by the study area, conclude the relative importance of criteria weights, and determine the most impacted alternative. First, the determination of the criteria and alternatives comes in the second place directly after the research problem and study aim identification. This set of criteria and alternatives have been designated through a deep investigation of relevant literature, current empirical studies, and inputs concluded from brainstorming meetings with experts in the tourism industry. As a result, a set of 10 criteria have been arranged and divided into three dimensions regarding environmental, social, and economic aspects. Additionally, eight municipalities of Malaga province have been selected due to their immense touristic activity in the peak seasons. Second, the dimension and criteria weights are calculated separately and then the total weights are generated. The main role of the F-AHP model is to identify relevant weights of each criterion and check the consistency rate of pairwise comparison matrices obtained from expert's judgments. After obtaining the fuzzy global weights, the initial decision matrices are generated with the assistance of the decision-makers and experts in the tourism sector. Furthermore, the obtained weights are integrated with the F-MAIRCA model to assess and select the most impacted municipality by overtourism.

Methodological Overview
This paper proposes the use of modified F-AHP and F-MAIRCA as an adapted methodology to evaluate overtourism issue of the case of Malaga province and accordingly propose the post-COVID-19 tourism sustainable/potential actions. The hierarchy framework consists of three succeeding steps, which are graphically presented in Figure 3. This framework can help tourism policymakers and manager analysts to identify the proposed set of criteria for overtourism evaluation by the study area, conclude the relative importance of criteria weights, and determine the most impacted alternative. First, the determination of the criteria and alternatives comes in the second place directly after the research problem and study aim identification. This set of criteria and alternatives have been designated through a deep investigation of relevant literature, current empirical studies, and inputs concluded from brainstorming meetings with experts in the tourism industry. As a result, a set of 10 criteria have been arranged and divided into three dimensions regarding environmental, social, and economic aspects. Additionally, eight municipalities of Malaga province have been selected due to their immense touristic activity in the peak seasons. Second, the dimension and criteria weights are calculated separately and then the total weights are generated. The main role of the F-AHP model is to identify relevant weights of each criterion and check the consistency rate of pairwise comparison matrices obtained from expert's judgments. After obtaining the fuzzy global weights, the initial decision matrices are generated with the assistance of the decision-makers and experts in the tourism sector. Furthermore, the obtained weights are integrated with the F-MAIRCA model to assess and select the most impacted municipality by overtourism.

Data Acquisition
To evaluate the set of pertinent criteria identified through literature review and prospective alternatives, the background of experienced experts in the tourism sector was considered. A team of twelve highly proficient decision-makers was created. Each decision-maker, industrious in a differently situated yet important Costa del Sol region (Malaga province), was invited to furnish subjective judgments and decisions based on self-knowledge and proficiency. In this research, a combination of literature review and questionnaire survey was employed involving the appraisal of rich judgments obtained from experts as a sample with quantitative data [85]. A pilot survey was performed over semi-structured interviews with the team to eliminate the equivocal viewpoints, rectify errors, and test the compatibility of the identified indicators with the actual concerns in Spain. The questionnaire was handed out among tourism practitioners working in different public administrations including tourism faculty, municipalities, tourism delegations, and national tourism organisations as illustrated in Table A1. These selected experts, at a political level, face daily the challenge generated by the impact of tourism in their destinations, addressing the complex and non-linear interaction between the 'right to the city' and the disruption of 'the tourist city'. The survey was constructed with two passages dedicated to creating the decision matrices to be employed in both AHP and MAIRCA method as explained in Table S1. However, the decision-makers' subjective judgements might entail ambiguity and some linguistic variables can refer to the same importance [86]. As a solution to this problem, triangular fuzzy numbers (TFNs) were utilised to overcome the vagueness and determine the linguistic importance weights. During the interviews, the experts stated the methods used by the tourism planning institutions to measure the overtourism in Malaga province as follows: • An analysis of tourist flows by territories and roads in the city is carried out using Big Data tools, evaluating the carrying capacity of each one of them. Tourism statistics are quantified and scaled by creating micro-statistics associated with the impact of tourism activity and its relationship with residential statistics, tourists, and visitors, by residents; • An analysis of behaviour in terms of spending and visitor profiles is carried out. From other areas of Government of the City Council such as environmental sustainability, urban planning, commerce, mobility, security, and others, different ratios are evaluated and a comparison is made with other European cities; • Evaluation of the reports issued by Tourism and Planning Costa del Sol through meetings with professionals in the sector; • Cooperation with neighbourhoods' service offices that enable them to deal directly with incidents suffered by residents due to excess tourism in the summer months. Consider monthly reports of hotel occupancy, tourist apartments, and analysis of the distribution of the VFT (housing for tourism purposes), and field work-study Profile of the Visitor to the city of Malaga.

Preliminaries
The triangular fuzzy number is used in this study. Therefore, some prime definitions and operations are concisely furnished in this section for the sake of the reader understanding. The Fuzzy set theory was primarily introduced by Zadeh [87] and has been widely employed to handle judgments extracted from real-world issues that contain ambiguity and vagueness. The main role of this theory is to deal with inaccurate data mathematically. A fuzzy set could be expressed arithmetically by giving a value to each discrete in the universe N = x, µ N (x) , x ∈ R , R 1 = −∞ < x < +∞ and where each value is referring to its grade of membership in terms of membership function µ N (x) to the close interval of [0, 1]. This is made to respond in a logical way to the primary imprecise judgments taken from various sources of information that might be ambiguous. The triangular fuzzy numbers (TFN) are the pillar of the fuzzy set and it is expressed with borders rather than crisp numbers to conduct a pairwise comparison matrix of the decision-makers. A triangular fuzzy number (TFN), symbolized by N = (l, m, u), entails the subsequent function [88]: The triplet values (l, m, u) determine the fuzzy situation and border the area of possible assessment to the close interval of [0, 1] where the parameter m shows the most favourable value, l represents the lower possible value, and the parameter u denotes the upper possible value. The basic operational laws for two triangular fuzzy numbers (TFNs) N 1 and N 2 , N 1 = (l 1 , m 1 , u 2 ) and N 2 = (l 1 , m 1 , u 2 ), are defined as follows [79]:

Modified Fuzzy AHP Method
The present study proposes a modified fuzzy AHP model, supports the evaluation process by using triangular fuzzy numbers to translate the judgments of experts, and proposes the geometrical mean method to obtain the value of the total weights. The following basic steps are essential to calculate the criteria weights [89]: Step 1: Construct pairwise comparison matrices for dimensions and criteria of the hierarchy structure using linguistic variables in Table 3. Absolutely Important (AI) (7/2,4,9/2) In Equation (7), kth denotes the number of the decision-makers who take part in the decision-making process. Each element d ij of the pairwise comparison matrix D (k) is a fuzzy number equivalent to the linguistic variables.
A linguistic term is a sentence or a word formed to refer to the decision-maker's preferences, which can be realised through the questionnaire. To translate the experts' subjective pairwise comparisons, the triangular fuzzy number (TFN) for importance weights of each criterion is structured as shown in Figure 4: Sustainability 2021, 13, x FOR PEER REVIEW 12 of 32 subjective pairwise comparisons, the triangular fuzzy number (TFN) for importance weights of each criterion is structured as shown in Figure 4: Step 2: In this step, the aggregated fuzzy pairwise comparison matrix is determined with the help of the decision-makers as represented in Equation (8).
Step 3: This step is necessary to be done to check the consistency of the decision-maker's judgments throughout the evaluation procedure. The consistency index (CI) of the fuzzy aggregated pairwise comparison matrix = , and consistency ratio (CR) are computed by applying Equations (9) and (10), respectively.
where, denotes the largest eigenvalue, the number of dimensions being used in the comparison matrix, and represents the random index. The values of varies correspondingly to the order of the judgments as provided in Table 4. The evaluation process is acceptable (consistent) when the value of the consistency ratio (CR) is equal to or less than 0.1. In case the obtained results are not consistent, the experts are requested to re-evaluate the pairwise comparison matrix.
Step 4: Next, the geometrical mean method is proposed in this study to determine the fuzzy geometric means of the triangular fuzzy numbers and the fuzzy weights for each row in the matrix as explained in Equation (11).
where, denotes the number of the decision-makers. Furthermore, the obtained fuzzy criteria weights are integrated with the proposed fuzzy MAIRCA method in the next sub-section to assess the alternatives. Step 2: In this step, the aggregated fuzzy pairwise comparison matrix is determined with the help of the decision-makers as represented in Equation (8).
Step 3: This step is necessary to be done to check the consistency of the decisionmaker's judgments throughout the evaluation procedure. The consistency index (CI) of the fuzzy aggregated pairwise comparison matrix D = d ij , and consistency ratio (CR) are computed by applying Equations (9) and (10), respectively.
where, λ max denotes the largest eigenvalue, N the number of dimensions being used in the comparison matrix, and RI represents the random index. The values of RI varies correspondingly to the order of the judgments N as provided in Table 4. The evaluation process is acceptable (consistent) when the value of the consistency ratio (CR) is equal to or less than 0.1. In case the obtained results are not consistent, the experts are requested to re-evaluate the pairwise comparison matrix.
Step 4: Next, the geometrical mean method is proposed in this study to determine the fuzzy geometric means of the triangular fuzzy numbers and the fuzzy weights for each row in the matrix as explained in Equation (11). where, k denotes the number of the decision-makers. Furthermore, the obtained fuzzy criteria weights are integrated with the proposed fuzzy MAIRCA method in the next sub-section to assess the alternatives.

Modified Fuzzy MAIRCA Method
The main shortcoming of the traditional MAIRCA method is in responding to the preference of decision-makers by using crisp values and its inability to calculate the total aggregated performance of multiple decision-makers. This study proposes a modified fuzzy MAIRCA method for the alternative's assessment with the use of triangular fuzzy numbers and the aggregated performance method [92]. The selection of the most impacted province by overtourism is conducted based on the gap values. The ranking is accomplished to measure every region with highlighted objectives of the analysis. The following steps are required to apply the proposed fuzzy MAIRCA model: Step 1: In this step, the linguistic variables are translated into triangular fuzzy numbers (TFN) for the agreed scale for rating the alternatives illustrated in Table 5 and Figure 5. The initial decision matrix X is constructed, where the m number of alternatives is assessed with respect to the n number of criteria. For generalisation, the number of decision-makers that are contributing to the assessment procedure of m alternatives is denoted as k regarding the n number of criteria. The fuzzy initial decision matrix is as follows: where, X (k) = x k l , x k m , x k u ; x l ≤ x m ≤ x u ; i = 1, 2, 3, . . . , m; j = 1, 2, 3, . . . , n; k = number of experts. Absolutely Important (AI) (7,9,9) Sustainability 2021, 13, x FOR PEER REVIEW 13 of 32

Modified Fuzzy MAIRCA Method
The main shortcoming of the traditional MAIRCA method is in responding to the preference of decision-makers by using crisp values and its inability to calculate the total aggregated performance of multiple decision-makers. This study proposes a modified fuzzy MAIRCA method for the alternative's assessment with the use of triangular fuzzy numbers and the aggregated performance method [92]. The selection of the most impacted province by overtourism is conducted based on the gap values. The ranking is accomplished to measure every region with highlighted objectives of the analysis. The following steps are required to apply the proposed fuzzy MAIRCA model: Step 1: In this step, the linguistic variables are translated into triangular fuzzy numbers (TFN) for the agreed scale for rating the alternatives illustrated in Table 5 and Figure 5. The initial decision matrix is constructed, where the number of alternatives is assessed with respect to the number of criteria. For generalisation, the number of decision-makers that are contributing to the assessment procedure of alternatives is denoted as regarding the number of criteria. The fuzzy initial decision matrix is as follows: Absolutely Important (AI) (7,9,9) where, ( ) = , , ; ≤ ≤ ; = 1,2,3, … , ; = 1,2,3, … , ; = number of experts. Step 2: The aggregated performance decision matrix provided by experts is calculated with the use of the geometrical mean and represented by Equation (13). Step 2: The aggregated performance decision matrix provided by experts is calculated with the use of the geometrical mean and represented by Equation (13).
Step 3: This step considers the normalisation, for each given criterion, of aggregated decision matrix X F as generated in Step 2. Assuming that this study does not hold a large number of incompatible criteria and the decision-makers do not take the type of the criteria (e.g., beneficial and non-beneficial) into account throughout the decision-making process. For the sake of ease in calculations, the normalisation process of the aggregated data (by geometric mean) holds a new technique unlike the one proposed by Pamučar [64], as given in Equation (14).
Next, the fuzzy normalised decision matrix is represented by Equation (15).
Step 4: Since the experts are indifferent towards the choice of the optimum alternative in the study, the present step ensures that each alternative is given an equal probability for selection as the ideal one. Hence, the preferences for each alternative can be defined by using Equation (16).
where m denotes the number of alternatives.
Step 5: Next, the elements of the fuzzy matrix of theoretical ponder T P X are computed by multiplying the preferences to the selection of alternatives and the fuzzy criteria weights taken from the application of fuzzy AHP method w i = w l i , w m i , w u i , where w l i = w m i = w u i . The fuzzy theoretical ponder matrix is computed as follows: Thus, the fuzzy theoretical ponder matrix looks like the following: where t P ij = 1 m × w j ; i = 1, 2, 3, . . . , m and j = 1, 2, 3, . . . , n.
Step 6: The fuzzy elements of the actual ponder matrix T r X is computed by multiplying the elements of the fuzzy normalised decision matrix to the elements of the fuzzy theoretical ponder matrix as represented in Equation (19): Step 7: To calculate the total gap matrix G, it is recommended to use the fuzzy Euclidean distance for each alternative for each criterion from the theoretical and actual ponder matrix. The reason for using this notion is to diminish the calculation steps and simply get the crisp gap value (de-fuzzification). The total gap matrix is represented by Equation (20).
where, the elements of total gap matrix G are represented by g ij .
Step 8: Final obtained value of the criteria functions is computed by summing the gap values for each alternative regarding each criterion as indicated in Equation (21). Then, the final values are ranked in ascending order.
As stated earlier, the best-ranked alternative is the one with the minimum gap value, whereas the worst-ranked alternative is the one with the maximum gap value.

Results
This study proposes a novel hybrid model, namely, fuzzy AHP-MAIRCA, adapted to assist the tourism policymakers to make better decisions in future and implement wideranging sustainable strategies. For good measure, the present paper aims at focusing on specific regions that suffered from overtourism before the COVID-19 outbreak and seeking to revive tourism activity after the crisis in Spain through a series of realistic applications.

Weights Calculated in F-AHP
Fuzzy AHP is applied in the present work to identify the relative fuzzy weights of the dimensions: environmental impact (C 1 ), social impact (C 2 ), and economic impact (C 3 ) as well as their corresponding criteria (C 11 , C 12 , C 13 , C 21 , C 22 , C 23 , C 24 , C 31 , C 32 , C 33 ). The overtourism complexity is solved by disintegrating it into a hierarchy entailing numerous sub-issues, which are scrutinised independently. Here, only the dimensions calculations are shown due to space constrictions. The same procedure is performed to compute the fuzzy weights of the 10 main criteria to obtain at the end of the local and global weights.
Step 1: The initial decision matrix D (1) is constructed with the support of the decisionmaker's assessment. The fuzzy pairwise comparison matrix of the dimensions is illustrated in Table 6. Step 2: Next, the initial decision matrix is aggregated by using Equation (8) as given in Table 7. Step 3: Using Equation (9), the consistency index (CI) of the fuzzy aggregated pairwise comparison matrix D = d ij is computed regarding the number of the criteria. In this case study, the value is CI = 0.0215. Then, the consistency rate is calculated by applying Equation (10) for the value of the random index. The value is CR = 0.0371. Similarly, all values of the consistency rate are less than 0.1 which means that the evaluation process is satisfactory and the collected data from experts is consistent. Table 8 shows the calculated hierarchical structure of local and global weights of the first decisionmaker. The global weights are obtained by multiplying the local weight of each dimension with the local weight of each criterion. Further, the results can be achieved via a similar computational way. Step 4: The fuzzy geometric mean of the triangular fuzzy numbers and the fuzzy weights for twelve decision-makers is computed according to Equation (11) and presented in Table 9.

Ranking Generated in F-MAIRCA
In the next step, the fuzzy obtained weights of the F-AHP model are integrated with the F-MAIRCA model to rank the alternatives according to the degree of overtourism impact. The assessment matrix provided by decision-makers is constructed regarding the scale given in Table 5. Only one of the decision-makers' assessment matrices is shown in Table 10 due to paper space limitations. Step 1: First, translate the linguistic variables obtained from experts into triangular fuzzy numbers using Table 5. Applying Equation (12), the initial decision matrix X 1 is developed with respect to ten criteria C j = (j = 1, 2, . . . , 10) and eight alternatives' vectors A i = (j = 1, 2, . . . , 8) as shown in Table 11.  (5,7,9) (1, 3, 5) (1, 1, 3) (5, 7, 9) (7,9,9) (5, 7, 9) (7,9,9) (1, 3, 5) (1, 1, 3) (7,9,9)  Step 2: The aggregated performance matrix is calculated using Equation (13) and presented in Table 12. Step 3: Next, the elements of the aggregated performance matrix are normalised applying Equation (14) as represented in the normalised decision matrix N X = n ij 8×10 in Table 13.  Step 4: In this step, the preference for each alternative is calculated using Equation (16). In the present work, the optimum value is 1 8 = 0.125. The matrix of preference is constructed for the preference value.
Step 5: Next, the elements of the fuzzy theoretical ponder matrix T P X = t P ij 8×10 are calculated by multiplying the obtained fuzzy weights from the fuzzy AHP method as shown in Equation (17), and are provided in Table 14. Step 6: The elements of the fuzzy actual ponder T r X = t r ij 8×10 are obtained using Equation (19) and are provided in Table 15. Step 7: The total gap matrix [G] 8×10 is computed by applying Equation (20) as given in Table 16. Step 8: Final criteria function values are calculated by using Equation (21) as provided in Table 17. As mentioned earlier, the municipality with the highest gap value is the most impacted one by overtourism and vice versa as per the proposed extended fuzzy MAIRCA methodology. The obtained results demonstrate that the Nerja municipality turned out to be the most impacted municipality by overtourism and Malaga capital is the least impacted one.

Discussion
In this section, the municipalities' rankings result obtained by applying the proposed model is validated. First, a comparative analysis of municipalities is furnished for the set of criteria. The ranking results are compared with various MCDM methods. Second, a sensitivity analysis is conducted by modifying the criteria weight coefficients to reveal the subsequent variations in municipalities' rankings.

Comparison Analysis
To validate the efficiency of the proposed model, a comparative analysis with some established techniques grounded in the tourism sector is performed. The comparative analysis will mainly focus on the MCDM approaches that have been applied in the tourism domain in general. A few of them are fuzzy VIKOR [94], fuzzy TOPSIS [95], fuzzy ARAS [96], and the fuzzy extension of the WASPAS method [97]. The obtained rankings are also compared with the fuzzy extension of the classical MAIRCA method developed by [64].
The ranking order of the municipalities A i (i = 1, 2, 3, . . . , 8) by different MCDM methods is presented in Table 18. The criteria weights obtained by fuzzy AHP methods in the present study are directly integrated with all other selected MCDM methodologies to obtain constant outcomes. The calculation performance steps are not presented in this sub-section since it is related to the comparative analysis of the ranking results. According to the ranking of municipalities by using different fuzzy MCDM methods, it is perceived that, each time, Malaga capital (A4) is identified as the less impacted municipality by overtourism phenomenon. Moreover, it is deduced that Nerja (A8) is identified as the second most impacted municipality by all other established fuzzy methods, except the FARAS method. Thus, the ranking results agree to middling with the original ranking order. Further, compared with the five methodologies, the proposed fuzzy hybrid model generates reasonable results. It is concluded that the tighter the municipalities' area is, the higher the overtourism impact is.
To check the reliability of the final decision results, it is preferable to compute Spearman's rank correlation coefficient to define the correlation of ranks generated by applying different approaches. According to Keshavarz Ghorabee et al. [98], the correlation is generally calculated to reveal the statistical significance of difference among the ranking generated throughout the application of various approaches. As per Table 19, it is observed that all the correlation coefficients are higher than 0.80, with an average value of 0.9131, which implies a strong correlation between the existing MCDM approaches and the proposed fuzzy hybrid model. Figure 6 also confirms the credibility of the proposed approach. by [64]. The ranking order of the municipalities ( = 1,2,3, … ,8) by different MCDM methods is presented in Table 18. The criteria weights obtained by fuzzy AHP methods in the present study are directly integrated with all other selected MCDM methodologies to obtain constant outcomes. The calculation performance steps are not presented in this sub-section since it is related to the comparative analysis of the ranking results. According to the ranking of municipalities by using different fuzzy MCDM methods, it is perceived that, each time, Malaga capital (A4) is identified as the less impacted municipality by overtourism phenomenon. Moreover, it is deduced that Nerja (A8) is identified as the second most impacted municipality by all other established fuzzy methods, except the FARAS method. Thus, the ranking results agree to middling with the original ranking order. Further, compared with the five methodologies, the proposed fuzzy hybrid model generates reasonable results. It is concluded that the tighter the municipalities' area is, the higher the overtourism impact is.
To check the reliability of the final decision results, it is preferable to compute Spearman's rank correlation coefficient to define the correlation of ranks generated by applying different approaches. According to Keshavarz Ghorabee et al. [98], the correlation is generally calculated to reveal the statistical significance of difference among the ranking generated throughout the application of various approaches. As per Table 19, it is observed that all the correlation coefficients are higher than 0.80, with an average value of 0.9131, which implies a strong correlation between the existing MCDM approaches and the proposed fuzzy hybrid model. Figure 6 also confirms the credibility of the proposed approach.

Sensitivity Analysis
In the MCDM methodology, the ranking results are directly linked to the criteria weight coefficients. Therefore, it is required to conduct sensitivity analysis by modifying the weight coefficients of the criteria to test the constancy and permanency of the new proposed approach. Generally, the variation in criteria weight coefficients leads to a slight change in the ranking order of the alternatives. Hence, it is essential to spot the sensitivity of the proposed model applied to rank the alternatives and determine at the same time the criteria that drive the important changes to alternatives rankings. As mentioned earlier, the weight coefficients in this study are performed based on decision-makers' and experts' relative importance, and the outcomes of possible deviance of these criteria weights should be evaluated. In various studies, sensitivity analysis is performed by adopting a slight modification in the original weights. For the sake of unity, the sensitivity analysis is executed by applying the following equation [99]: where, γ denotes the percentage of the variation of the original weights.
In this work, a total of 40 scenarios are assigned based on criteria C j (j = 1, 2, 3 . . . , 10) and their corresponding alternatives A i (i = 1, 2, . . . 8), and commonly divided into four stages, where each stage entails 10 different scenarios. For each scenario, the relative weights are decreased and increased with a variation of (10%, 30%, 50%, and 70%) for the total new sum as provided in Table 20. Additionally, Figure 7 shows the variation in municipalities rankings for various level of inconstancy for each criterion (a, b, c, d, e, f, g, h, i, j). After further observation, the ranking order remains steady unless rough variations are executed in the criteria weights. Ranking results

Municipalities
Proposed model 10% variation 30% variation 50% variation 70% variation According to the findings, it is noticed that Nerja municipality (A8) preserves its first position in 100% of the scenarios, while Estepona municipality (A6) maintains its second position in 80% of the scenarios. Furthermore, Nerja is considered as the most impacted municipality by overtourism and Estepona can be considered as the second most impacted municipality. The obtained results also imply an important correlation (90%) by applying the fuzzy adapted F-AHP and F-MAIRCA. This significant correlation indicates that the final raking results of Malaga region's municipalities using the proposed methodology are credible and valid.
Overtourism is a phenomenon that occurred before the pandemic. If the pre-pandemic regulations remain in the destinations, similar situations and processes are likely to emerge with renewed intensity [100]. The rhetorical question is whether one should return to the previous status as soon as possible or take advantage of this opportunity to modify, partially or totally, the economic and social conditions that led to this situation [101]. Spain is a typical example of current alterations in its tourism policy. These fluctuations are occurring due to several reasons as mentioned above. The findings revealed that the economic ground is worsening the tourism situation, followed by the social and environmental ground. There is urgent need to conduct resilient public policies that align with sustainability concerns. This study assessed the impact of overtourism from different backgrounds (economic, environmental, and social). Therefore, the sustainability dimensions must be examined interdependently and not independently. The findings of this study also demonstrate that overtourism cannot be examined by taking the tourism perspective alone. Instead, the tourism policymakers are required to consider management strategies in a wider scope to keep up with the fast-growing movement within the tourism sector. The proposed methodology allows local governments to take measures according to the needs of their destinations without going back to the mistakes of the past and seeking a balance among environment, society, and economy [102]. Sustainable interventions tend to be valued positively by tourists and, in addition, promote increased competitiveness of destinations and touristic companies, although their implementation may incur additional costs. These costs of adapting tourist destinations to the new times should not affect a process of differentiation by prices for the promotion of the destination.
Nerja is a clear example of a tourist destination in which results and interpretations about tourism can change in a short period (overtourism/undertourism). The attractiveness of this destination is the result of many voices of supply and demand, in which economic development must be combined with DMOs. Therefore, the solution of balance between this overtourism/undertourism is to make everyone understand the role of tourism and accept the tourism product as such, that spatial and seasonal distribution of tourism is carried out and continuous monitoring and development of the Smart Destination are carried out to allow a process of continuous improvement [103]. Additional measures could be aimed at: extending the length of stay of tourists, stocking up in closer or proximity markets, promoting the consumption of organic and seasonal products, moving towards a model with high added value, encouraging local benefit and not that of global foreign-owned platforms, abandoning carbon-intensive travel, reducing seasonality, preserving the social fabric of cities, and promoting public-public, private-private, and public-private collaboration [104]. This period can be seen as an opportunity to reconsider a transformation of the global tourism system, focusing it towards one that is more aligned with the United Nations Sustainable Development Goals for 2030 (SDG) [105]. As indicated by Pons et al. [106], tourism sustainability does not only imply providing a balanced solution to economic profitability but also to social well-being and ecological balance.

Conclusions
Tourism endures incessant metamorphosis and requires sustainable management for residents, tourists and visitors. The economic, social, and environmental negative effects of overtourism have to be corrected. Hence, there lies the importance of reflection on the perception of this problem, determining its impacts, and proposing answers. Currently, Spain is facing a global pandemic that has paralysed the tourism sector and has exposed its extreme volatility. Therefore, it is time to put the innovation and accumulated tourist knowledge at the service of society, especially in areas such as the Costa del Sol. In the post-COVID-19 scenario, sustainability will be the differentiating factor to be taken into account in decision making in tourist destinations. The local destinations should start working more forcefully and sustainably, betting on collaboration also between destinations, due to the competitiveness to obtain the attention of future tourists. The visitor will be significant and their choice must be defined based on the product, the excellence, and the experience offered by the destination and not based on the price. Along the course of the pandemic, until it reaches its end, many changes will surely appear. The tourist who comes will have new requirements, in line with the experiences in previous destinations, together with the possible trauma caused by the difficult environment created by COVID-19. The requirement of safety (health, economic, and social) linked to environmental and social sustainability will be the pillar of the new framework within the tourism sector once post-COVID-19 begins to emerge.
The fuzzy modified AHP method was applied in this study for calculating the weights of the overtourism dimensions and main criteria. On the other hand, the fuzzy modified MAIRCA method was used for spotting the most impacted municipality in Malaga province by the overtourism phenomenon in the pre-pandemic scenario. The obtained global weights via fuzzy AHP were utilised as inputs in the adapted fuzzy MAIRCA method for assessing the municipalities and designating the most and less impacted ones. The attained findings were acknowledged by the professionals and there was a strong agreement among them. Firstly, the fuzzy AHP proposes a consistency rate (CR) equation to verify the pairwise comparison matrix of the decision-makers. Further, it demonstrates the mutual relationship between the push factors of the overtourism phenomenon. This can help policymakers and managers to conduct local and global weights for each criterion to reveal the impacted municipality. The main influencing criteria are: "CO 2 emission/marine pollution (C 11 )", "inorganic waste (C 12 )", and "Economic dependence on tourism sector (C 33 )". Secondly, the fuzzy MAIRCA illustrates to us how to aggregate the experts' judgments and select the most impacted municipality by overtourism issue. The eightieth alternative (Nerja municipality) turned out to be the most impacted municipality in Malaga province. Thus, the tourism policymakers in the public sector and managers in the private sector would be capable to support the impacted municipalities by implementing sustainable tourism actions that go along with overtourism and undertourism in the post-COVID-19 scenario. Thus, the proposed framework may help as a reference for depicting and analysing the push factors as well as communities worldwide.
Andalusia has not been the only region hardest hit by the effects of the pandemic on the health of its population. On the contrary, it may be one of those that suffer the most from its economic consequences. The structural characteristics of the Andalusia economy, its lower density and business size, its sectorial specialisation, its low capacity to generate employment, and the trends observed in recent years of high relative entrepreneurship out of necessity, accompanied by some incipient experiences of entrepreneurship in High potential (Start-ups) very localised in very few geographical locations, make the consequences of the crisis particularly worrying in Andalusia. The impact has been high in consumer-oriented sectors such as tourism, hospitality, leisure, and retail with a high presence in the area. It has also been relatively important in micro-enterprises, self-employed and recently created companies, which have more relative weight in the region. It has especially hit women, youth, and people with less training and with temporary contracts. In this context, localities affected by overtourism can apply the eleven strategies proposed by the UNWTO [9] to deal with this phenomenon: promote the dispersion of visitors inside and outside of the city; promote the dispersion of visitors based on time; stimulate new itineraries and attractions; review and adapt regulation; improve visitor segmentation; ensure that local communities benefit from tourism; create urban experiences for residents and visitors; improve the infrastructure and facilities of the city; communicate and in-volve local stakeholders; communicate and attract visitors; and establish supervision and response measures.
The present study has few boundaries as well. Only 3 dimensions and 10 main criteria for overtourism phenomenon analysis have been acknowledged. Further dimensions and criteria (intuitional, technology, industrial) have not been acknowledged here. Further, it is an analysis of the current situation, with which the results cannot be extrapolated to other populations or the same ones analysed in other periods. The calculation steps of the proposed model are accomplished by considering the preference of the decisionmakers. Henceforth, it is important to achieve these evaluations prudently. Furthermore, the proposed framework findings are grounded on a sole area; therefore, numerous MCDM techniques like (grey cognitive maps, ANP, MABAC, and DEMATEL) can be practical to analyse the push factors of the overtourism phenomenon.
Supplementary Materials: The following is available online at https://www.mdpi.com/article/10 .3390/su13116394/s1. Table S1: the questionnaire example. Funding: The research in this paper was funded by the project "Overtourism in Spanish coastal destinations. Tourism degrowth strategies an approach from the social dimension" (RTI2018-094844-B-C33) financed by the Spanish Ministry of Science, Innovation and Universities (National Plan for I+D+i), the Spanish State Research Agency, and the European Regional Development Fund (ERDF).

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement: Not applicable.
Acknowledgments: The authors would like to thank the Spanish Ministry of Science Innovation and Universities (National Plan for I+D+i), the Spanish State Research Agency, and the European Regional Development Fund (ERDF) for making a collaborative program for carrying out the research work.

Conflicts of Interest:
The authors declare no conflict of interest. Table A1. Experts Identification.

Experts Expertise
DM 1 Former head of tourism and promotion of Malaga city council and expert in tourism delegation with 30 years of work experience DM 2 Dean of tourism faculty of the University of Malaga Spain with 25 years of work experience DM 3 Area advisor of tourism delegation of Estepona city council with 17 years of work experience DM 4 Coordinator of the tourism department of Fuengirola city council with 12 years of work experience DM 5 Head of the tourism service of the territorial delegation of tourism ministry (JUNTA DE ANDALUCÍA) in Malaga with 41 years of work experience DM 6 Expert supervisor of Marbella city tourism project with 25 years of work experience DM 7 Managing director of tourism and planning in Costa del Sol area with 15 years of work experience DM 8 Head of the tourism department of Mijas city council with 35 years of work experience DM 9 Councillor delegate for tourism of Nerja city council 27 years of work experience DM 10 Head of marketing section, promotion and tourism development of Malaga city council with 20 years of work experience DM 11 Responsible for the tourism delegation of Torremolinos city council with 20 years of work experience DM 12 Coordinator of the tourism area for foreigners of Benalmadena city council with 25 years of work experience