Assessment of Coastal Locations Safety Using a Fuzzy Analytical Hierarchy Process-Based Model

: Worldwide tourist beaches have been an ideal destination for people searching for a recreational place to visit; however, several conditions could generate risks for these visitors. At present, the main efforts to assess these risks and prevent people from danger are essentially focused on monitoring tidal height in the zone, ignoring other risk sources. Therefore, this work generates an overall evaluation that considers the tidal height, bathymetry, temperature, solar radiation, and wind speed, establishing a relationship between parameters and safe beach conditions using a fuzzy logic approach. In addition, this paper presents the design and implementation of a computational model, based on a fuzzy analytical hierarchy process (FAHP), to evaluate coastal and climatological parameters involved in tourist safety, which can be continuously monitored. Tourist beaches in South Baja California, M é xico, with diverse environmental and coastal conditions, were assessed with the proposed model, providing several safety scenarios that contrast the results between sites and demonstrate the capacity of the model to evaluate them. The evaluation, obtained as a result of the computational model, presents information about the safety conditions in the assessed zone, considering the possible risks for all the parameters analyzed, which could be presented to beach visitors to prevent dangerous situations and avoid accidents.


Introduction
In several countries, tourism has become an essential economic income, and beaches are considered an important part of the tourist sector.Particularly in México, this industry grew by 12.7% in 2017 with the arrival of 3.4 million international tourists [1]; 73% of the economic income associated with tourism proceeded from beach areas in this country [2].Notwithstanding the attractive beach and weather conditions, occasional natural events, such as red tide, sea, or environmental conditions, could represent a risk for tourism activities.At present, most of the efforts to assess security in coastal zones are limited to monitoring the height of the waves, and few metrics define if the conditions are safe for tourist activities.One of the most popular ways of assessing a beach location is the flag color-coding.Unfortunately, in a large number of tourist destinations, safety signs do not exist; moreover, some tourists ignore them, or do not know how to interpret them when they do exist, and this leads to many accidents being reported in these zones.Consequently, there is a need for a continuous assessment of safety conditions in tourist beach zones that considers the physical characteristics of the zone and its weather conditions.In this sense, some questions can be raised to determine if it is possible to improve the assessment of beach conditions using different oceanological parameters instead of just tidal height.Likewise, can different levels of hazardous conditions be assessed to alert visitors?To this end, we hypothesize that a model based on monitoring parameters with different importance levels and preferences can help determine good or bad beach conditions to provide real-time information about better and safer places for visitors.
A study of the existing risks in coastal locations in Australia and their signage is presented in [3]; in this paper, results of beach user interviews are presented, demonstrating that almost half the visitors did not even see the signals.Hammerton [4], by using interviews, determined that 82% of the participants knew someone who had drowned and that a significant number of the total drownings took place in the wave-dominated coastline of the analyzed cases in Ghana.Another study presented by Silva-Cavalcanti [5] analyzed the locations of warning signals and the locations of drowning incidents, finding that the safety conditions did not always correspond with the warning signs.Different criteria have been defined to evaluate the safety level of areas for tourist activities.Nevertheless, there is a need to create tools for risk detection and the evaluation of suitable conditions for tourist activities in coastal zones, based on climate and physical parameters, and present them in an easily understandable way.Some evaluations of coastal areas have been carried out from different perspectives.Considering beach quality, Leatherman [6] presented a ranking of beaches, taking into account 50 physical, biological, and human parameters evaluated quantitatively, adapting the evaluation methodology from [7] designed for river scenery.Looking also for beach quality, Micallef [8] described a quantitative assessment in the region of Andalusia, Spain, that considers five beach parameters considered relevant, but only one of them forecasts safety.The evaluation of this one parameter contemplates the existence of infrastructure features and neglects the environmental or sea conditions.In the results of the works just mentioned, the evaluator perspective and preferences had importance, since the appraisal of some parameters comes from subjective judgments.Campuzano [9] made characterizations of the Bahia Blanca in Argentina, employing data analysis and numerical modeling to improve the understanding of the dominant hydrodynamic processes in the area, providing useful information for the management of the area.
Moreover, taking into account coastal landscape assessment, Ergin, [10] defined 26 parameters to evaluate coastal scenery through surveys of tourists and landscape perception experts in Turkey and the UK.For each parameter, five ranges of values or attributes were established.The parameters are integrated using a fuzzy logic approach (FLA), providing a general assessment that overcomes subjectivity and uncertainty.The mentioned methodology was applied to assess the landscape in coastal zones from the Western Black Sea, and the results were published in Ergin [11].
In México, the analysis for coastal safety assessment is limited to the observation of the tide, calculation of the amount of stool, and the Blue Flag certification [12], but there are few incipient models to integrate this study or provide results with important variables related to the safety of tourism in coastal zones.A first approach was proposed in [13], where a safety index for coastal areas was proposed by analyzing tide, bathymetry, temperature, and wind speed; the parameters are integrated by employing an analytical hierarchy process approach for importance weight assessment.However, the parameters themselves do not define a safety condition as a whole, which is needed to determine when site conditions are dangerous or safe, using a site and priority classification.
On the other hand, multi-criteria decision making (MCDM) models are utilized in problems with multi-criteria dimensions, where each criterion can be quantitative or qualitative, and could present conflict with other input.Therefore, the output must be a compromise [14].For these properties, MCDM models have been applied with success in several disciplines, such as medicine, business, government, and environmental management.Particularly, in this last discipline, MCDM has solved complex problems; for example, in [15], Deveci used a CODAS based model to identify the best rehabilitation solution for abandoned mines, and in [16], the same author made a modification to the WASPAS model, applying Einstein norms in a fuzzy environment over triangular fuzzy numbers to assign priorities to climate change mitigation strategies in urban mobility planning.Furthermore, the most popular variety of MCDM model, especially in environmental science, is the AHP methodology [17], which is flexible in considering a group of factors to combine the opinion of many experts, reducing complex problems into simple pairwise comparisons between elements [18].
This paper presents the design of a computational model that assesses safety in coastal zones for the practice of tourist activities considering physical, oceanographic, and climatological variables, where the importance of each parameter is pondered and with the aim of providing information to visitors in an easy and understandable format that other evaluations do not consider, improving the actual safety evaluations and addressing the flaws.The relative importance of each evaluated parameter is assigned with a fuzzy analytical hierarchy process (FAHP) methodology because of its capacity to consider a group of factors and the opinion of many experts simultaneously, as well as prioritize the elements that have more impact on the overall evaluation.Then, all the parameters are integrated with a fuzzy logic approach, obtaining a Beach Safety Index (BSI) as a result.This paper considers the region of Los Cabos as a case study; this destination is located in the south of the Baja California peninsula, México, and it is one of the most important tourist destinations in the country, where tourism is the main source of income for the region.The climate conditions and favorable tide for surfing attract thousands of tourists each year, but not all beaches are adequate for bathing.The data provided from the existing infrastructure for the monitoring of physical and climate conditions, in addition, to the information obtained from previous research works [19], are entered as inputs to the model.Temperature, solar radiation, wind speed, bathymetry, and tide height were the parameters selected for analysis and safety evaluation.
The rest of the paper has been structured as follows: Section 2 explains the considerations for coastal safety conditions in tourist locations, as well as the model building process.Section 3 describes the implementation of the proposed model and the data used for assessing specific beaches, and Section 4 presents the discussions and suggests further work.

Literature Review
Several works have provided methodologies for beach assessment.In Table 1, a summary of these works found in the literature is presented as it relates to the problem of the safety assessment of the beach zones and the MCDM models, as well as their principal contributions.

Matthews et al. (2014) [3]
Analysis of interviews with beach visitors; almost half of them did not see safety signals.

Hammerton et al. (2013) [4]
Analysis of interviews with local people; 82% of them knew someone who had drowned.A significant number of drownings occurred on the wave-dominated coast.
Silva-Cavalcanti et al. (2018) [5] Study of incidents location.Safety conditions do not correspond to the warning signals.

Micallef et al. (2011) [8]
Quantitative assessment of beach quality.Safety is evaluated considering the infrastructure and response capacity for an emergency.

Campuzano et al. (2014) [9]
Characterization of Bahia Blanca, Argentina, employing numerical modeling and data analysis to improve understanding of the dominant hydrodynamic process.

Sampling Site
The regions of Los Cabos and La Paz, located in the state of South Baja California, México, are annually visited by thousands of international and domestic tourists.The most popular destinations in this region are Cabo San Lúcas, San Jose del Cabo, La Paz, and Los Barriles.These sites are also the most populated in the state (besides Los Barriles, with only 1674 residents): La Paz, with 250,141, Cabo San Lúcas, with 202,694, and San José del Cabo, with 136,285 residents, according to the local government statistics [20]; for these reasons, they have been considered as sampling sites for this study.Facilities in these locations provide information about weather and ocean measurements, and the closeness of different tourist sites is ideal for our proposal.The beaches in South Baja California are attractive to tourists because of their enjoyable scenery, hot and dry climate, and turquoise waters.However, several differences exist between the beaches of the region, from low deep bays, with low tide in La Paz, to open sea beaches, where swimming or bathing activities are extremely dangerous, such as in San Jose City (Figure 1).Different activities can be practiced depending on the characteristics of a particular zone; some places are ideal for swimming and sunbathing, while others are better for the practice of aquatic sports.Climate and oceanographic conditions at these sites are continuously monitored, providing data that can be used to evaluate the safety of the area.The sampling sites were chosen to take into account the variety of conditions between sampling sites and the importance of the destinations for the tourism industry.

Factor Analysis
Beach safety assessment depends on several factors, where the behaviors of the involved parameters each have a different impact because the risk involved in the behavior of some parameters is greater than others.In this work, three principal factors are considered to assign the relative importance between parameters, but in contrast with the methods in [13], factors taken into consideration are the principal risks in beach zones.The first contemplated factor is the drowning risk.The next factor is the environmental risk, where environmental conditions are not desirable, such as high wind speed or high solar radiation.The third analyzed factor is red tide risk, the presence of algal blooms that cause skin irritation.In this sense, factors are denominated as follows: can be practiced depending on the characteristics of a particular zone; some places are ideal for swimming and sunbathing, while others are better for the practice of aquatic sports.Climate and oceanographic conditions at these sites are continuously monitored, providing data that can be used to evaluate the safety of the area.The sampling sites were chosen to take into account the variety of conditions between sampling sites and the importance of the destinations for the tourism industry.In this work, five parameters have been identified as common to the studied factors and the most relevant to compute a safety evaluation.Their importance lies in the fact that they can be automated and monitored for data acquisition purposes.These parameters are tide, temperature, wind speed, solar radiation, and bathymetry.The tide and bathymetry are related to the drowning risk [21][22][23], the temperature and wind speed are associated with the red tide risk, and wind speed and solar radiation are mainly associated with the environmental risk.Although in the literature, a high number of parameters are assessed, they correspond to static observations, such as sand type, relief, and beach material, among others [6,10].This work has the aim of monitoring parameters that present dynamic behaviors and can change the safety conditions on the beach; this set will be continuously monitored.
Different categories and optimal ranges of each parameter have been defined in Table 2.These classifications are based on the international scales for sea conditions defined by the World Meteorological Organization (WMO) [24], Encyclopedia Britannica [25], and information collected from research works [26][27][28], and the health recommendations come from and World Health Organization (WHO) [29].These categories include different status types presented in beach locations dedicated to tourism.

Model Architecture
In general, computational models for beach assessments have been developed only for aesthetic purposes [8,11].For this reason, this work proposes a safety indicator through a fuzzy logic approach which, due to this technique, is well adapted to environmental problems in which uncertainty is used for modeling parameters with different ranges of validity.
The relative importance between parameters and factors is assigned by employing the FAHP methodology, aiming to eliminate subjectivism in the weighting procedure.The model architecture is shown in Figure 2 and is formed by four principal stages.The first stage obtains the data inputs for each parameter.The next stage is the fuzzy assessment, where each input is evaluated in its optimal scale, retrieving a normalized output.In the third stage, each fuzzy output is multiplied by an importance weight obtained through the FAHP methodology.The last stage integrates all the outputs, retrieving as output a normalized index that represents the safety label in the evaluated coastal zone.
x FOR PEER REVIEW 6 of 20

Model Architecture
In general, computational models for beach assessments have been developed only for aesthetic purposes [8,11].For this reason, this work proposes a safety indicator through a fuzzy logic approach which, due to this technique, is well adapted to environmental problems in which uncertainty is used for modeling parameters with different ranges of validity.
The relative importance between parameters and factors is assigned by employing the FAHP methodology, aiming to eliminate subjectivism in the weighting procedure.The model architecture is shown in Figure 2 and is formed by four principal stages.The first stage obtains the data inputs for each parameter.The next stage is the fuzzy assessment, where each input is evaluated in its optimal scale, retrieving a normalized output.In the third stage, each fuzzy output is multiplied by an importance weight obtained through the FAHP methodology.The last stage integrates all the outputs, retrieving as output a normalized index that represents the safety label in the evaluated coastal zone.

Fuzzy Logic Approach
Usually, the modeling of environmental parameters is complex, and the development of a mathematical equation can be a hard task.Fuzzy logic has overcome this issue, classifying parameters according to a certain grade against Boolean evaluations, where only one true or false range can be defined as a value that varies from a bad to an excellent condition.In this case, all parameters are normalized in more understandable ranges, which can be used as inputs on the same scale, being easier to process.The purpose of the fuzzy assessment is to evaluate each parameter according to its scale and

Fuzzy Logic Approach
Usually, the modeling of environmental parameters is complex, and the development of a mathematical equation can be a hard task.Fuzzy logic has overcome this issue, classifying parameters according to a certain grade against Boolean evaluations, where only one true or false range can be defined as a value that varies from a bad to an excellent condition.In this case, all parameters are normalized in more understandable ranges, which can be used as inputs on the same scale, being easier to process.The purpose of the fuzzy assessment is to evaluate each parameter according to its scale and optimal level, transforming the measured values into a [0, 1] range.This normalizes all parameter ranges and scales [30].Fuzzy expressions were defined for the normalization of each parameter, setting the limits based on the permissible ranges, as detailed in Table 2.A letter was assigned to refer to each of the parameters: d refers to the tide, t to temperature, s to wind speed, r to solar radiation, and b to bathymetry; the expressions are assigned as follows:

Importance Comparison Rules
To assign the relative importance weights with the FAHP methodology, verbal comparisons of importance are made between factors and independently between parameters, considering each factor, and then these comparisons are transformed into fuzzy triangular functions using Table 3, as defined by Wang [31].An important consideration for the analysis is that if a first element is considered more important at any level of the scale than a second element, the comparison of the second element against the first is the reciprocal function.

Fuzzy Analytical Hierarchy Process Modelling
The safety evaluation of beaches can be defined according to the importance of each oceanographic parameter.For this purpose, an analysis of parameters with a high impact on the safety assessment is proposed.The fuzzy analytic hierarchy process methodology (FAHP) was proposed by Laarhoven and Pedryez [32], extending Saaty's AHP theory [33], establishing the importance of comparisons using fuzzy functions, represented by a set of numbers denoted as (l, m, u) and denominated triangular fuzzy numbers.The membership function of a fuzzy triangular number M in R, µ M (x) : R → [0, 1] is described by the equation: where l ≤ m ≤ u, l and u represent the upper and lower value of the triangular function, respectively, and m represents the value with the maximum membership.The range of M is delimited by the set elements, consequently {x ∈ R | l < x < u}.If l = m = u, the set is not a fuzzy number by convention (Table 3).Several approaches have been proposed to obtain the importance weights using the FAHP methodology.This work employs the extent analysis method developed by Chang [34], and detailed in [35], computing the mathematical operations as described by Deng [36].The first step is to order the importance comparisons, in the form of triangular fuzzy numbers, in a positive reciprocal pairwise matrix (A = M j g i ), with n × n dimensions, as follows: The value of the synthetic extent for the ith object is calculated using Equation (8); M 1 g i , M 2 g i , . . .M m g i are defined as values of the extent analysis of the ith object of m objectives.
The possibility degree of M 1 ≥ M 2 represents the membership value for the intersection of the fuzzy functions M 1 and M 2 .This definition is shown in Figure 3 and described mathematically in the next expression: 2, 13, x FOR PEER REVIEW 8 of 20 The possibility degree of  1 ≥  2 represents the membership value for the intersection of the fuzzy functions 1 and  2 .This definition is shown in Figure 3 and described mathematically in the next expression: Since M1 and M2 are convex fuzzy numbers: The possibility degree for a convex fuzzy number to be greater than k convex fuzzy numbers can be assumed as: where k =1, 2, …, n;  ≠ .The  ′ (  ) values for all the objectives are appended in a vector  ′ .The final weight vector W is obtained by normalizing the values of  ′ : Since M 1 and M 2 are convex fuzzy numbers: Sustainability 2022, 14, 5972 9 of 21 The possibility degree for a convex fuzzy number to be greater than k convex fuzzy numbers can be assumed as: where k = 1, 2, . . ., n; k = i.The d (A i ) values for all the objectives are appended in a vector W .The final weight vector W is obtained by normalizing the values of W :

Weights Calculation
Defining importance criteria is not an easy task, and it must be done based on the judgment of experts on the subject.In this work, the opinions of 12 experts on beach locations were collected through questionaries.Considerations about relationships of importance among factors and parameters for the assessment of safety in coastal tourist zones were provided by 4 academics from ITES (Instituto Técnológico de Estudios Superiores), Los Cabos [37], dedicated to beach conservation from tourist and civil protection profiles, with expertise between 10 and 15 years; 4 academics from CICIMAR (Centro Interdisciplinario de Ciencias Marinas) [38] focused on marine biology and oceanography, with experience between 8 and 13 years; and 4 lifesavers that belong to the Mexican Fund for Tourism Promotion (Fondo Nacional de Fomento al Turismo) [39] and the Federal Marine Zone (Zona Federal Marítima) [40], with 5 years of experience.Although the 12 experts' work is focused on beach zones, they have different profiles, and their judgments are established from different perspectives of the problem, which are contained in the overall statements of safety, avoiding bias or tendency to only one opinion or perspective; furthermore, several studies have been developed with smaller groups of experts [41].The questionnaire was applied in four stages: in the first, the risk danger is determined by employing verbal pairwise comparisons for Saaty's scale [33].In the following three stages, the relative impact of each parameter against the others, considering one risk at the time, is also defined with verbal pairwise comparisons.The results for the questionnaires are presented as pairwise fuzzy numbers shown in Tables 4-7, transforming the linguistic comparisons into fuzzy functions according to the criteria defined by [31].In this case, the comparisons between factors are established in Table 4.The relationships of importance between parameters are presented in Tables 5-7, defining each of the factors as decision criteria, respectively.According to the drowning risk factor (F1), bathymetry has been considered as the most important parameter because beaches that present high inclinations are non-optimal for any visitor, and some swimming requirements are needed to have a minimal risk; the second-highest importance is assigned to the tide, as shown in Table 5.In the case of environmental risk factor (F2), solar radiation is assigned with the greatest weight of relative importance, as can be seen in Table 6.The second more important parameter is the wind speed.Finally, in the comparison taking into account red tide risk (F3), wind speed and temperature are the most important parameters due to their influence on the appearance of red tide, as seen in Table 7.The general weight priorities are calculated by combining the importance weights of Tables 4-7 according to the sum of the weight multiplication, where the resulting importance weights show the unity.The result is presented in Table 8.

Beach Safety Index
The final step is to create a method of assessing safety according to an understandable score.In this sense, the beach safety index (BSI) provides a final value between the [0, 1] range.The BSI provides a score of "1" when all the parameter conditions are in their optimal range; contrastingly, if one or more parameters present non-optimal values, the score decreases, reaching "0" when the safety conditions of the zone are completely unsuitable for the practice of tourist activities.The BSI index involves the complete parameter set evaluations, integrating partial scores into a global result.The FAHP weights amplify or attenuate the negative impact of each parameter according to its importance level, generating a score concerning the safety condition in that location.According to this, the integration of the evaluated parameters in the BSI index is described by: where k refers to a parameter, µ i (k) is the fuzzy function to normalize the measured values for the parameter k, defined in Section 3.1, and W k i is the importance weight calculated using the FAHP methodology for the parameter k.By replacing each parameter in Equation ( 14), the following expression results: where d refers to the tide, t is the temperature, s is wind speed, r is radiation, and b is bathymetry.Replacing the calculated weights in Equation ( 15), the safety index can be rewritten as follows:

Data Collection
To show the safety index (SI) performance, a database of collected measurements was created using measurements obtained from meteorological stations (Vantage Vue Weather Station) at sampling sites located on Médano Beach at Cabo San Lúcas, Costa Azul Beach at San Jose del Cabo, Malecon Beach at La Paz, and Barriles Beach at Los Barriles.Likewise, measurements obtained from Meteoblue [42] and Windy [43] websites complemented the required data information.The proposed database contains 336 register samples for each place, with measurements at every hour from the 7th to the 20th of October 2019.Each register includes values for temperature, wind speed, solar radiation, tide, and bathymetry.

Parameter Analysis
In Figure 4, database information is represented as signals, where parameter dynamics can be observed.According to the permissible limits defined in Table 1, the wind speed was below the optimal range at almost all the measurement times and for all the locations.This is a condition that increases the probability of red tide presence, representing one of the assessed risks for the practice of tourist activities.The higher mean value of wind speed during the sampling period was registered on Malecon Beach (Figure 4a).Temperature presented a fluctuating behavior, which can be explained due to sunset and sunrise conditions (Figure 4b).Only a few samples per day exceed the permissible limits.Since the sampling sites are in the same region, the temperature signals were measured with similarity.Solar radiation (Figure 4c) showed similar behavior in the four places almost all the time, usually reaching values out of the optimal range at midday.The tide is presented in Figure 4d, where higher waves were measured at the Costa Azul and Médano sites.In some cases, the high waves were out of the optimal range.For Malecon and Barriles beaches, the tide was low, which means a good condition for tourist activities.On some days, the obtained values from these sites increased the tide safety conditions; however, they were never cataloged as dangerous.In the sampling sites, the maximum depth allowed for tourists to swim was defined as 2 m.

Fuzzy Assessment
The results of assessing the measured values for each parameter, using their respective fuzzy expressions, Equations ( 1)-( 5), are shown in Figure 5.The fuzzy assessment for wind speed indicated that the conditions of this parameter, in a great part of the sampled period, were bad because it was almost always out of the optimal range (Figure 5a). Figure 5b shows a highly variable behavior that corresponds with the measured signal of temperature changes during the sampled period; since high temperatures are out of the optimal range, these moments were considered as bad conditions of temperature for safety.Similarly, in the solar radiation fuzzy assessment (Figure 5c), the daily variations can be observed and the high radiation values are evaluated as a bad condition for safety.In the case of the tide, the fuzzy assessment scored

Fuzzy Assessment
The results of assessing the measured values for each parameter, using their respective fuzzy expressions, Equations ( 1)-( 5), are shown in Figure 5.The fuzzy assessment for wind speed indicated that the conditions of this parameter, in a great part of the sampled period, were bad because it was almost always out of the optimal range (Figure 5a). Figure 5b shows a highly variable behavior that corresponds with the measured signal of temperature changes during the sampled period; since high temperatures are out of the optimal range, these moments were considered as bad conditions of temperature for safety.Similarly, in the solar radiation fuzzy assessment (Figure 5c), the daily variations can be observed and the high radiation values are evaluated as a bad condition for safety.In the case of the tide, the fuzzy assessment scored 1 in a great portion of the sampled period for Malecon and Barriles, while in the other two sites, the evaluation varied between values from 0 to 0.73 (Figure 5c).Fuzzy evaluations are integrated into the BSI index, and the final assessments are computed and explained in the following section.
Sustainability 2022, 13, x FOR PEER REVIEW 13 of 20 1 in a great portion of the sampled period for Malecon and Barriles, while in the other two sites, the evaluation varied between values from 0 to 0.73 (Figure 5c).Fuzzy evaluations are integrated into the BSI index, and the final assessments are computed and explained in the following section.

Beach Safety Index (BSI)
The BSI integrates the fuzzy evaluations of all the parameters with the importance weights obtained from the FAHP analysis using Equation (16).Since the weights and the parameter inputs are normalized, the obtained output is a number between the [0, 1] range, 0 being the worst case of safety conditions, and 1 the optimal safety conditions.If all parameters are in the optimal range, then the score will be 1; otherwise, it will decrease when one or more parameters fall out of the optimal range.If the BSI score is a value between 0.85 and 1, the safety conditions are considered excellent; between 0.7 and 0.84,

Beach Safety Index (BSI)
The BSI integrates the fuzzy evaluations of all the parameters with the importance weights obtained from the FAHP analysis using Equation (16).Since the weights and the parameter inputs are normalized, the obtained output is a number between the [0, 1] range, 0 being the worst case of safety conditions, and 1 the optimal safety conditions.If all parameters are in the optimal range, then the score will be 1; otherwise, it will decrease when one or more parameters fall out of the optimal range.If the BSI score is a value between 0.85 and 1, the safety conditions are considered excellent; between 0.7 and 0.84, good; from 0.5 to 0.69, regular; from 0.25 to 0.49, bad; and values from 0 to 0.24 are considered very dangerous safety conditions.The negative effect of undesirable conditions in temperature, wind speed, and radiation can be noticed in the BSI score; however, those changes have little influence on the safety index response.In contrast, the behavior of the index is affected in great measure by tide and bathymetry variations, because the weights obtained in the FAHP process for these parameters are significantly higher than those calculated for the other parameters.The safety conditions during the testing dates at Malecon and Barriles beaches were from regular to excellent conditions almost all the time good; the minimum value was assessed on October 11th and the better conditions were registered on October 14 and 15; Malecon Beach safety conditions were better than the other sites because the tide height was lower, as can be seen in Figure 4d).Costa Azul and Médano beaches show similar behavior in the BSI evaluation, registering conditions as regular and bad in the sampled period.The minimum values for both sites were on October 10; on this day, the measured values were out of the optimal ranges, yielding a low score in the fuzzy assessment (Figure 5).The beach safety index variations for the sampled sites are displayed in Figure 6; the mean, maximum, and minimum values for each site are included.
Sustainability 2022, 13, x FOR PEER REVIEW 14 of 20 good; from 0.5 to 0.69, regular; from 0.25 to 0.49, bad; and values from 0 to 0.24 are considered very dangerous safety conditions.The negative effect of undesirable conditions in temperature, wind speed, and radiation can be noticed in the BSI score; however, those changes have little influence on the safety index response.In contrast, the behavior of the index is affected in great measure by tide and bathymetry variations, because the weights obtained in the FAHP process for these parameters are significantly higher than those calculated for the other parameters.The safety conditions during the testing dates at Malecon and Barriles beaches were from regular to excellent conditions almost all the time good; the minimum value was assessed on October 11th and the better conditions were registered on October 14 and 15; Malecon Beach safety conditions were better than the other sites because the tide height was lower, as can be seen in Figure 4d).Costa Azul and Médano beaches show similar behavior in the BSI evaluation, registering conditions as regular and bad in the sampled period.The minimum values for both sites were on October 10; on this day, the measured values were out of the optimal ranges, yielding a low score in the fuzzy assessment (Figure 5).The beach safety index variations for the sampled sites are displayed in Figure 6; the mean, maximum, and minimum values for each site are included.To validate the proposed approach, an index was generated using the methodology proposed by [6], taking into account only the four parameters considered in this study.This index is denominated as the statistical index (SI).Simultaneously, the proposed index is also compared against the approach presented by Dávila-Lamas [13], which also considers four parameters.For the three indexes, the fuzzy inputs are the same.The evolution of the BSI for Médano, Costa Azul, and Malecon beaches are shown in Figure 7, respectively, and are compared with the results of the SI and AHP models.The differences in the behavior of both indexes correspond to the differences in weights assigned for each model, while in the SI, all the parameters have the same weight.The proposed fuzzy AHP methodology calculates different weights for each parameter according to its relative importance, which is different from that calculated with the AHP approach because of the methodology differences and the changes in factor definitions and considerations.Moreover, the solar radiation is not considered in any of the comparison models.The statistical index retrieves a lower score for several samples, due to the fact that most of the parameters are out of their optimal ranges, and the weights assigned for all the parameters are equal.The BSI FAHP score is higher in these cases because the weight for those parameters is lower than the weights assigned for bathymetry and tide.Because of the above, the statistical model presented lower minimums and mean values, while the To validate the proposed approach, an index was generated using the methodology proposed by [6], taking into account only the four parameters considered in this study.This index is denominated as the statistical index (SI).Simultaneously, the proposed index is also compared against the approach presented by Dávila-Lamas [13], which also considers four parameters.For the three indexes, the fuzzy inputs are the same.The evolution of the BSI for Médano, Costa Azul, and Malecon beaches are shown in Figure 7, respectively, and are compared with the results of the SI and AHP models.The differences in the behavior of both indexes correspond to the differences in weights assigned for each model, while in the SI, all the parameters have the same weight.The proposed fuzzy AHP methodology calculates different weights for each parameter according to its relative importance, which is different from that calculated with the AHP approach because of the methodology differences and the changes in factor definitions and considerations.Moreover, the solar radiation is not considered in any of the comparison models.The statistical index retrieves a lower score for several samples, due to the fact that most of the parameters are out of their optimal ranges, and the weights assigned for all the parameters are equal.The BSI FAHP score is higher in these cases because the weight for those parameters is lower than the weights assigned for bathymetry and tide.Because of the above, the statistical model presented lower minimums and mean values, while the maximum values are almost equal to those obtained with BSI FAHP.In the case of the AHP index, the values obtained are generally greater than the ones calculated with the SI and BSI FAHP.Nevertheless, dangerous situations related to temperature, wind speed, and radiation are mostly ignored, since these parameters were assigned with a low importance weight.The BSI FAHP proposed in this work assigns weights of relative importance, in contrast with the SI, and considers the most common risks that can affect a beach visitor as criteria for defining the priority for each parameter, avoiding the risks associated when the parameters with lower importance are underestimated, contrary to the AHP index, and taking solar radiation into account.
Sustainability 2022, 13, x FOR PEER REVIEW 15 of 20 maximum values are almost equal to those obtained with BSI FAHP.In the case of the AHP index, the values obtained are generally greater than the ones calculated with the SI and BSI FAHP.Nevertheless, dangerous situations related to temperature, wind speed, and radiation are mostly ignored, since these parameters were assigned with a low importance weight.The BSI FAHP proposed in this work assigns weights of relative importance, in contrast with the SI, and considers the most common risks that can affect a beach visitor as criteria for defining the priority for each parameter, avoiding the risks associated when the parameters with lower importance are underestimated, contrary to the AHP index, and taking solar radiation into account.The standard deviation and variance were higher in the statistical index (σ SI and σ 2 SI ), in comparison with the values obtained for the BSI FAHP model (σ FAHP and σ 2 FAHP ), which are slightly greater than the ones obtained using the AHP approach (σ AHP and σ 2 AHP ). Figure 8 shows the correlation between the SI and BSI FAHP indexes, while Figure 9 presents the correlation between BSI FAHP and BSI AHP; these plots are obtained by graphing the values of each comparison model against the values obtained with the BSI FAHP.In Table 9, some numerical examples of the more significant differences in the values obtained with both indexes are detailed, along with their input parameters.ustainability 2022, 13, x FOR PEER REVIEW 16 20 The standard deviation and variance were higher in the statistical index (  and  2  ), in comparison with the values obtained for the BSI FAHP model (  and  2  ) which are slightly greater than the ones obtained using the AHP approach (  and  2  ). Figure 8 shows the correlation between the SI and BSI FAHP indexes, while Figure 9 presents the correlation between BSI FAHP and BSI AHP; these plots are obtained by graphing the values of each comparison model against the values obtained with the BS FAHP.In Table 9, some numerical examples of the more significant differences in the values obtained with both indexes are detailed, along with their input parameters.Coastal studies are mainly focused on assessing tourist considerations about beach scenery [10,11] or beach quality [6], and while safety is one of the factors evaluated by some works [8], it is evaluated by subjective judgments, or by considering the presence or absence of personnel and infrastructure for attending emergencies.These evaluations neglect the natural events that could represent a risk for beach visitors.In contrast, the approach proposed in this work evaluates the safety conditions based on the analysis, study, and integration of physical, meteorological, and oceanographic parameters that can be monitored continuously.While in the work proposed by Micallef [8], safety is defined considering the capacity of response in case of emergency, in the approach presented in this

Conclusions
The computational model for the calculation of the BSI index proposed in this paper evaluates the safety conditions in beach zones, considering three principal risks for the tourist as factors, and five climatological and physical variables for the definition and calculation of each risk.The results of its implementation in four sites in Baja California Sur, México, are presented and discussed, successfully detecting the differences between the particular safety conditions of each of the sampling sites.Furthermore, continuous monitoring enables the model to reveal variations in the safety conditions during the sampling period.The proposal given in this work has been developed to provide information about excellent or bad conditions for tourist activities on beaches, taking into consideration the oceanographic and meteorological conditions at the site.The principal limitations of the case study detailed in this work are the variety in safety conditions, limited by the characteristics of the sampling sites and sampling period.Although the model generated in this research was applied to assess the safety of Mexican tourist beaches, it can be used in any tourist coastal zone in the world, simply by replacing input data with those of the locations of interest.This will allow for the testing of the model in a wider variety of conditions.For future work, more parameters have been proposed to be studied and analyzed, such as harmful wildlife, or ocean aspects; however, to consider more parameters in the FAHP methodology, the whole analysis must to be performed from the start.On the other hand, the integration of other techniques and computational models for better analysis is also proposed, such as the geospatial analysis, which could help improve the accuracy and granularity of the results, geographically differentiating the results in a beach zone.This model could be implemented to display to beach visitors the resultant evaluation in the form of a color-coded score, helping visitors to have a better understanding of the safety conditions in the location, and as a consequence, to avoid dangerous moments and areas, preventing accidents and providing a safe visit to tourist beaches.Safety and risk management has proven to be an efficient strategy of local administrations to improve the tourism offerings of their destinations, and consequently, to increase their regional economies.

Figure 1 .
Figure 1.Sampling sites monitored: Médano Beach in Cabo San Lúcas, Costa Azul Beach in San Jose del Cabo, Malecon Beach in La Paz, and Barriles Beach in Los Barriles, all located at Southern Baja California, México.

Figure 2 .
Figure 2. Flowchart of the proposed model used for computing the BSI index.

Figure 2 .
Figure 2. Flowchart of the proposed model used for computing the BSI index.

Figure 3 .
Figure 3. Definition of possibility degree functions M 1 and M 2 .

Figure 4 .
Figure 4. Parameter measurements were used for assessing safety conditions in México: Costa Azul Beach in San Jose del Cabo (SJC), Médano Beach in Cabo San Lúcas (CSL), and Malecon Beach in La Paz (LP).a) Wind speed b) Temperature c) Solar radiation d) Tide.

Figure 4 .
Figure 4. Parameter measurements were used for assessing safety conditions in México: Costa Azul Beach in San Jose del Cabo (SJC), Médano Beach in Cabo San Lúcas (CSL), and Malecon Beach in La Paz (LP).(a) Wind speed (b) Temperature (c) Solar radiation (d) Tide.

Figure 6 .
Figure 6.Safety assessment using the BSI index for the Costa Azul, Médano, Malecon, and Barriles beaches.

Figure 6 .
Figure 6.Safety assessment using the BSI index for the Costa Azul, Médano, Malecon, and Barriles beaches.

Figure 7 .
Figure 7. Safety assessment using the BSI FAHP index compared against the SI and the BSI AHP.(a) Costa Azul Beach at San José del Cabo.(b) Médano Beach at Cabo San Lúcas.(c) Malecon Beach at La Paz.(d) Barriles Beach at Los Barriles.

Figure 7 .
Figure 7. Safety assessment using the BSI FAHP index compared against the SI and the BSI AHP.(a) Costa Azul Beach at San José del Cabo.(b) Médano Beach at Cabo San Lúcas.(c) Malecon Beach at La Paz.(d) Barriles Beach at Los Barriles.

Figure 8 .
Figure 8. Correlation between the BSI FAHP and the statistical index.(a) Costa Azul Beach at San José del Cabo.(b) Médano Beach at Cabo San Lúcas.(c) Malecon Beach at La Paz.(d) Barriles Beach at Los Barriles.

Figure 8 .
Figure 8. Correlation between the BSI FAHP and the statistical index.(a) Costa Azul Beach at San José del Cabo.(b) Médano Beach at Cabo San Lúcas.(c) Malecon Beach at La Paz.(d) Barriles Beach at Los Barriles.

Table 2 .
Permissible ranges for safety parameters.

Table 3 .
Linguistic scale conversion to triangular fuzzy functions.

Table 4 .
Comparison levels between factors.

Table 5 .
Parameter comparison according to F1: drowning risk.

Table 6 .
Parameter comparison according to F2: environmental risk.

Table 7 .
Parameter comparison according to F3: red tide risk.

Table 8 .
Parameter comparison considering all factors and weights.

Table 9 .
Numerical comparison between evaluations for the BSI and SM indexes.