Rip current hazard assessment on a sandy beach in Liguria, NW Mediterranean

Rip currents are one of the most significant environmental hazards for beachgoers and are of interest to coastal scientists. Several studies have been conducted to understand rip current dynamics, and several approaches for rip hazard assessment have been proposed. In general, the purpose is to provide knowledge and tools to support authorities and lifeguards in rip current risk prevention. This study proposes the application of an expeditious methodology to evaluate rip current hazard and risk, based on probability theory. The tested area was located along the Alassio beach, a renowned tourist destination located on the western Ligurian coast (NW Italy). A coastal video-monitoring system was used for rip currents individuation, whereas wave data were collected thanks to an oceanographic buoy managed by Regione Liguria. In detail, a yearly analysis was performed to identify the correspondence between rip currents and wave parameters data. The results showed that rip currents occur, in the study area, under moderate wave conditions (0.5≤Hs≤1.34\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.5 \le H_s \le 1.34$$\end{document} m; 4.7≤Tm≤7.0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4.7\le T_m \le 7.0$$\end{document} s; 150∘N≤θm≤227∘\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$150^{\circ }\,\hbox {N} \le \theta _m \le 227^{\circ }$$\end{document} N). Based on this analysis, an easy application of the probability theory was applied to evaluate the level of hazard. Moreover, considering the official tourist data, we also perform an expeditious rip currents risk evaluation. The results showed that the hazard level is considered high at annual time scale and moderate during the tourist season; the risk is related to seasonal presences. The study can propose a tool to support authorities and lifeguards in water safety planning and management.


Introduction
In Europe, more than 600 million people per year go to the beach during the summer period, when most of the 27.000/year drownings in coastal waters occur (Funari et al. 2016). In Italy, the death rate is attested to be between 50 and 100 cases. In 2014, 13 fatal drownings occurred along the Ligurian coast: three of those were local beachgoers and ten 1 3 were tourists (Funari et al. 2016). Italy has not official database on accidents and rescue. Additional information on accident dynamics is not described. The beachgoers risk exposure theme has a key role in coastal management practice . Even though official data are rare, it is obvious many people are potentially exposed to natural hazards along the Ligurian beaches-NW Italy. Among other natural hazards in this environment, rip currents are the greatest natural hazard for beachgoers, especially along wave-dominated coasts (Austin et al. 2012). These currents "are narrow and concentrated seaward-directed flows that extend from close to the shoreline and through the surf zone and varying distances beyond" (Castelle et al. 2016). Rip current behaviour and the related risk to bathers are, to a certain extent, predictable and directly related to the nearshore dynamics (MacMahan et al. 2006). For example, intermediate sandy beaches, i.e. beaches with transverse bar morphology subjected to medium-energy wave conditions, are considered particularly conducive to rip current activity (Brander 1999;Wright and Short 1984). Rip currents can also occur near rocky shoals and coastal structures, as well as headland-embayed beaches (Scott et al. 2014). Beachgoers are exposed to drowning risk related to rip currents. This is due to a range of factors such as: unfamiliarity with the environment and associated hazards (Ballantyne et al. 2005;Clifford et al. 2018;Moran and Ferner 2017), poor swimming abilities (Drozdzewski et al. 2015;Williamson et al. 2012) and inattention when on vacation (Clifford et al. 2018;Wilks and Pendergast 2010). Every year, hundreds of people drown and tens of thousands more are rescued from rip currents globally [e.g. da F. Klein et al. (2003), Hartmann (2006), Gensini and Ashley (2010), Brewster (2010), Brighton et al. (2013), Scott et al. (2011), Kumar and Prasad (2014), Arozarena et al. (2015) and Barlas and Beji (2016)]. Rip currents are an hazard responsible for swimmers' catch and drag out seaward (Shepard et al. 1941;Short and Hogan 1994). They can quickly move swimmers spatially to regions of deeper water and greater hazard and have globally documented as a significant natural hazard (Shepard et al. 1941;Short and Hogan 1994;Short 1999). Following Leatherman (2017), rip current hazard was primarily evaluated utilizing United States Lifesaving Association (USLA) rescue statistics. These are available online and show the rescue number from rip currents for each beach USLA lifeguards is present on (United States Lifesaving Association 2013). According to the USLA, 80 per cent of all lifeguard rescues are associated with rip currents. In the Netherlands, a rip current prediction model system is being used on the basis of the bathymetry measurement along with application of the XBeach hydrodynamic model (Roelvink et al. 2009), to predict current speeds and directions in the context of Egmond beach (Van Ormondt et al. 2012). Other work Austin et al. (2012) developed a rip currents operational tool using field measurements, a twodimensional horizontal (2DH) non-stationary model and ARGUS video sensing. Certainly, Italian beaches show minor risk factors compared to the Ocean coastlines. However, particular conditions such as the presence of structures and/or particular morphologies can increase rip current hazard (Brander and MacMahan 2011). In Tuscany, part of the "Perla" project (Cooperation Program border Italy-France "Maritime" 2007-2013 founded by EU) has started an awareness process addressed to beachgoers based on the characterization of the coasts and based on beach morphological features. The University of Florence has developed a methodology to define the sectors with different safety conditions for bathing including warning signs where rip current hazard are more probable (Pranzini 2015). Rip currents events affecting Mediterranean coastlines have been described by Benassai et al. (2017), via a model chain (WaveWatchIII + XBeach), sediment grain size analysis and unmanned aerial vehicle (UAV), highlighting rip events occurrence under wave conditions with offshore wave height between 0.3 and 1.4 m and peak periods from 2.7 to 5.8 s.

3
The aim of this study is to propose an expeditious method to evaluate the rip currents hazard and risk. In detail, we proposed a variant to the method proposed previously , based on the probability theory. The method is tested on the Alassio beach, located on the Ligurian Western coast (Ligurian Sea, NW Italy), which is one of the most popular seaside tourist location of the continent. The stretch of coast is managed by 20 beach resorts, and beach safety is guaranteed by resorts via seasonal lifeguards. The use of remote sensing video was considered with the purpose of recording rip current flows on the Alassio beach from April 2012 to March 2013 (Murray et al. 2013). These data, together with Capo Mele buoy wave dataset, defined the most common sea conditions to trigger rip currents in Alassio beach. Moreover, two rip events (among those observed) were simulated using 2DH hydrodynamic model (XBeach). XBeach is a well-known tool for rip currents modelling (Austin et al. 2012;Van Ormondt et al. 2012), and we used it to obtain a further qualitative confirmation regarding the rip currents genesis in the study area. The results show how the proposed approach for rip currents hazard and risk assessment is easily adoptable and customizable for single study cases, and how the probabilistic approach could be easily integrated in a marine forecasts system. Moreover, the proposed method follows the basic principles for hazard and risk assessment. This research is also a useful support for local authorities and coastal managers that deal with beach safety and risk prevention.

Study area
Alassio town is one of the most popular seaside tourism destinations in Italy and Europe, and it is located on the Ligurian Western coast (NW Italy) (Fig. 1). The main shore ( Fig. 1) can be described as a sandy beach on a 3-km long, NE-SW oriented coastline; the backshore is spatially limited by a promenade and a pedestrian wharf is situated in the central sector of the coastline (Fig. 1). The nearshore area, from the shoreline to 10 m water depth, has a slope of approximately 1.5 • . The zone from 10 to 20 m of water depth is of 4.5 • (Bowman et al. 2007). Alassio beach is classified as dissipative ( > 20) following (Vacchi et al. 2014) using surf scaling index supplied by (Jackson et al. 2005). Rivers feeding the coastline belong to small catchments characterized by a torrential regime and mainly located on the western Liguria flyschoid Units composed of limestone, marly limestone, marls and marly siltstone (Cortesogno and Vanossi 1984;Ferrari et al. 2014).
The investigated beach is exposed to the typical Western Liguria wave climate with storms from the South (Vietti et al. 2010). The more frequent storms come from SW, but Alassio beach is partially protected by capo Mele promontory that generates strong wave diffraction. Instead, marine weather conditions from the south-west are less frequent but more incisive in the beach processes showing a significant cross-shore sediment transport ).

Field survey and sampling
Dry beach profiles were surveyed using a Differential Global Position System (DGPS) in November 2011, while submerged beach morphology was acquired using single-beam echo sounder, mounted on a small boat and positioned by a Trimble real-time cinematic (RTK) GPS system. Topographic and bathymetric data have been collected along regularly spaced cross-shore transects from the emerged beach to 20 m water depth. Singlebeam data were corrected considering wave conditions as well as astronomical and baric tide. Bathymetries from 20 to 80 m in depth were supplied by the Istituto Idrografico della Marina Italiana. The morphology of the promontories was added to the dataset using the LIDAR data, acquired in 2008 (Italian Ministry of the Environmental Land and Sea, 2008-2011). These data have been subsequently merged together and interpolated onto a regular grid through GIS software, producing a Digital Terrain Model (DTM). Furthermore, wave data have been collected from Capo Mele Buoy, which is located 6.12 km far from the study area in correspondence to bathymetric 80 m (

Video monitoring
Images were collected by the web camera located in the northern side of the area and placed on the top of a building 25 m a.m.s.l. The web camera provides a view of 350 m of Alassio beach northern side, with 720 × 567 pixels resolution (Fig. 1). The system collected 120 images per hour that were successively rectified and elaborated using the Timex method. Through this study, a hydro-morphodynamic behaviour of the beach was analysed through a video monitoring system (Aarninkhof et al. 2005;Holland et al. 1997;Jiménez et al. 2007;Turner et al. 2004). Image datasets were analysed using Beachkeeper plus software (Brignone et al. 2012;Schiaffino et al. 2015), an image management tool providing Time-Exposure, Variance, Day-Timex, as suggested by the extensive coastal video monitoring literature (Aarninkhof et al. 2005;Alexander and Holman 2004;Davidson et al. 2006;Holman et al. 2003). The aim of these observations is to identify rip current flows occurring between April 2012 and March 2013 along the Alassio beach northside. We analysed 3600 Timex images, elaborated on 481800 snapshots for the year. In details, 12 Timex/day were elaborated during the period with major daylight (2160 images), while 8 Timex/day was used during the period with minor daylight (1440 images) to recognize hydrodynamic conditions with potential patterns. Flow detection from images was based on the physical consideration that the rip current produces a visible colour contrast, that is easy to detect in the surfzone when light and picture quality are adequate (Pugliano et al. 2019).

Rip currents modelling using XBeach
XBeach model was used to obtain a further qualitative confirmation for rip currents occurrence in the study area. XBeach concurrently solves the time-dependent short wave action balance, roller energy equations, nonlinear shallow water equations of mass and momentum, sediment transport formulations and bed update on the scale of wave groups (Roelvink et al. 2009). A computational grid was generated merging bathymetric and topographic data described in Sect. 3.1. The ultimate computational grid was generated by optimizing the spatial resolution across the study area, which is 765 m for cross-shore and 1230 m for alongshore directions. As shown in Fig. 1, an irregular mesh was used obtaining a higher resolution on the surf zone and a lower resolution for areas far from the coastline. The grid resolution was 5 m on the beach and 25 m off the coast. Offshore wave boundary conditions were implemented using the JONSWAP spectrum (Hasselmann et al. 1980) considering the Capo Mele buoy wave data for seaward flow simulation. XBeach has a number of free parameters used to calibrate the model. For this application, parameters were modified as shown in Table 1, whereas all other parameters have been set to their recommended default value. The breaking index ( ) has been calculated using a linear relationship with the product of local wave number and water depth (Ruessink et al. 2003); D50 and D90 values have been obtained from sedimentological analysis; bed friction has been computed using the White-Colebrook formulation, based on the relation between the D90 of the top bed layer and the geometrical roughness described previously (Nikuradse 1950) https ://xbeac h.readt hedoc s.io.

Rip currents hazard assessment
Hazard can be defined as "the probability of occurrence of a potentially damag-ing phenomenon within a specific period of time and within a given area" (Varnes 1984). To evaluate the rip hazard level, we propose a method based on probability theory . We considered three different hazard descriptors: (1) Significative wave height ( H s ), (2) Wave period ( T m ) and (3) Wave direction ( • m ). The hazard descriptors were combined using basic principles of probability theory (Spiegel 2008). Therefore, based on the Ferrari approach (Ferrari et al. 2019), the following equations can be written: gives the probability that all the events will occur simultaneously, while Eq. 2 expresses the probability of a single event occurrence. The probability of occurrence of the opposite event X (i.e. probability of non-occurrence) satisfies the following equation: Consequently, for each combination of the considered hazard descriptors, the rip hazard (RH) was calculated by using the following formula: Each hazard descriptor was quantified via three different probability values: 0.01 (low), 0.5 (moderate) and 0.99 (high). These values were assigned according to occurrence of the wave conditions (i.e. H s , T m , • m ) suitable for rip currents generation (days/ year or days/tourist season) ( Table 2). The choice to consider two different time scales (e.g. yearly and tourist season) allowed a complete assessment of the probability of rip currents occurrence (i.e. rip hazard). The yearly time scale allowed one to consider a wide range of wave boundary conditions and tourist season time scale is the most relevant in terms of beach safety management. The RH-index furnishes information about the rip current hazard degree on a given beach. Finally, three classes of hazard were determined based on the following conditions (Table 3):  probability value and the other one has low probability value, up to a condition where one hazard descriptor has low probability value, one has moderate probability value, and one has high probability value. • High hazard (RH > 0.995): at least two hazard descriptors have high probability value.

Rip currents risk assessment
Risk can be defined as a measure of the probability and severity of an adverse effect to health, property or the environment (Fell et al. 2008). Thus, if we try to obtain a rip currents risk analysis, we must identify and quantify the elements (i.e. bathers) at risk. Unfortunately, it is not possible to know the exact number of beachgoers that are on the beach at a given moment. Both for the study area and for the Italian coasts in general, control systems for beach tourists fluxes are not available at the moment. Thus, in order to obtain a quantitative evaluation of the beachgoers presence (i.e. elements at risk), we adopted the dataset of the tourist presences in 2012 for the Alassio municipality. A variable, identified as risk conditioning variable ( X 0 ), was defined based on tourist presences during the year (Fig. 5). Usually, tourist activities are divided into three seasons: Low, Medium and High, related to statistical tourism attendance (Fig. 5). Risk conditioning variable ( X 0 ) was introduced to quantify three different probability values for each season, as reported in Table 4. According to Ferrari et al. (2019), RH-index and risk conditioning variable ( X 0 ) were combined to obtain an evaluation of rip currents risk (RR) on the Alassio beach (Eq. 5).
Three RR classes (low risk, moderate risk and high risk) were obtained from the intersection of RH and X 0 . In this way, it was possible to obtain a RR-matrix (Fig. 9), a common approach to conducting a subjective risk assessment (Markowski and Mannan 2008) (Fig. 6).
(5) RR = P(X 0 )RH Table 4 Risk conditioning variable X 0 related to low, medium and high Season The numerical coefficients describe the probability to find people on the beach, and it is necessary to combining X0 and RH-index (through Eq. 5)

Web camera dataset and climate wave
Web camera images analysis was conducted for twelve months on the investigated beach. The correspondence between climate wave and recorded rip currents genesis was evaluated collecting wave parameters from Capo Mele buoy yearly data, described in Table 5. The  table reports mean values

XBeach simulations
As described in the previous section, to better analyse the rip currents features recorded by web camera, two cases (July 2012 and December 2012) were simulated using the 2DH coastal model XBeach. The output from XBeach model provides the reliable results in both considered cases as they correspond to the rip currents recorded by web camera. The results describe incident wave height and flow direction along the beach (Fig. 8). Figure 8 refers to the time steps elaborated by the XBeach model at 9 AM and at 11 AM on 22 July 2012, at 9 AM and 11 AM on 7 December 2012. Figure 8 illustrates the presence of flow rotation close to the coastline during the previously mentioned time steps, highlighting how wave height is lower in the rip neck than in the feeder zone. Velocities display mean values between 0.3 and 0.8 ms −1 in both simulations.

Fig. 8
Focus on hydrodynamic XBeach results for the investigated area: July (a1, a2) and December (b1, b2) 2012. Rip currents behaviour is indicated by black arrows and its magnitude is proportional to the arrows dimension; blue continuous line represents coastline

Rip currents hazard and risk
We applied the proposed approach for rip current hazard assessment on the Alassio beach, along the western Ligurian coast. Wave parameters and rip currents correspondence, detected by camera dataset and described in the previous section, were used to evaluate when rip currents could develop during the year (Figs. 2, 3, 4). A yearly investigation has reported that H s values, among 0.5 and 1.34 m, occur at 203 days, while T m among 4.7 and 7.0 s and • m among 151 • N and 227 • N occur for 90 and 121 days. Using the same range described above, tourist seasonal wave data analysis revealed 60 days of H s , 62 days for T m and 41 days for • m . Based on our approach, both H s and T m (hazard descriptors) show a moderate probability of occurrence during the tourist season. Term • m has a low probability of occurrence. Thus, Alassio beach displays a moderate level of rip hazard. Differently, considering the annual dataset of wave parameters, T m and • m show a moderate probability of occurrence whereas H s shows a high probability of occurrence. Therefore, the level of rip hazard at the annual scale is high. Combining RH with tourist presences ( X 0 ), we also obtained an evaluation of the rip currents risk (RR) on the Alassio beach. RR is low in lowtourist season, moderate in middle-tourist season and high in high-tourist season (Fig. 9).

Rip current monitoring
The first part of this study was focused on rip currents behaviour in relation to coastal hydrodynamics conditions. Camera-recorded images showed that rip currents are a common phenomenon occurring on the north side of Alassio beach. Thanks to the coastal video-monitoring, 28 cases of rip currents genesis were identified on time laps between the 1 April 2012 and the 31 March 2013. In detail, rip currents flows were identified along filmed area by means of Timex images. Offshore wave analysis conducted using the Capo Mele Buoy database shows a clear correspondence between wave parameters and rip current genesis. In accordance with the literature (Brander 1999;Scott et al. 2011), they occur under moderate-energy conditions, nevertheless, in the Mediterranean environment, the Fig. 9 Matrix of risk (RR-matrix) applied in the Alassio beach. RH-index is different by considering tourist season dataset (a) or yearly dataset (b). Anyhow, risk levels are the same tidal level factor is not considered decisive on rip current circulation assessment. Moreover, following Murray et al. (2013), Castelle et al. (2016), investigated rip currents can be classified as "flash rip": they occur when the beach was generally uniform alongshore in a shore-attached terrace state with dominant plunging breakers. This, is a near-normal wave incidence, and the rip currents are episodic and unpredictable, therefore, an increasing risk of drowning for bathers.
The web camera analysis used here is a suitable tool for the study of rip currents even though cameras need to be placed in an area where albedo effects and, more in general, weather conditions do not affect the recording. The use of snapshots is often, is insufficient to identify rip currents: the employment of timex methods is more useful given that it can discern rip currents flow. Data collected were also used for a rip currents evaluation by means of XBeach modelling. The model results from July and December 2012 show the genesis of rip currents with velocity values around 0.5 ms −1 located where they were observed by video monitoring in the first place. Rip currents velocity values are comparable to those observed by Wright and Short (1984) and Bowman et al. (1988), respectively, in the first, the beach can be classified as intermediate transverse bar under medium wave condition. In the second, Herzliyya littoral (Israel) is transverse bar beach with intermediate/transitional state and wave conditions show values between 0.3 and 1.1 m for Hs and between 3 and 5 s for Ts. The highest velocities are detected in the rip neck, while the lowest velocities are in the rip head when the flow expands and decreases in intensity (McKenzie 1958;Sonu 1972;Short and Hogan 1994;Masselink and Hegge 1995;Aagaard et al. 1997;Short 1999). The XBeach model application produced good results for the rips development study giving a qualitative evaluation of the hydrodynamic processes.

Rip hazard and risk assessment
According to Bowman et al. (1988), Benassai et al. (2017) and Dusek et al. (2011), the results show that rip currents occur under moderate conditions, when a significant wave height ( H s ) is between 0.50 and 1.34 m with mean direction ranging from 150 • and 227 • N. Moderate wave conditions increase swimmers' risk because their self-confidence causes them to underestimate the hazard level. Rip currents are one of the main hazards for bathers, Brighton et al. (2013), Castelle et al. (2016) and Ferrari et al. (2019), and thus, the development of a rip currents hazard assessment and prediction method is a well-known challenge for coastal scientists (Kumar and Prasad 2014;Austin et al. 2012;Brander and MacMahan 2011;Engle 2002). However, many proposed approaches are based on rescue and/or accidents data and are not suitable for a quick hazard assessment. Moreover, reliable accident statistics are often limited if not absent. For example, along the Italian coasts, which are among the most popular tourist beaches in the world, no reliable statistics dataset are available. We can find data on the fatal accidents, but they are not indicative because they often include all kinds of accidents on the beach (e.g. heat stroke, heart attack). Data on drownings usually do not supply additional information on accidents dynamics (e.g. sea conditions) and place (e.g. sea, river, lake or swimming pool). In Funari et al. (2016), an accidents statistic on drownings in the sea waters is provided. However, data are reported at the regional scale, and they are not useful as a support for lifeguard activity. Thus reason, the rip hazard assessment method proposed here is not based on accident statistics but on wave boundary conditions. Thanks to a rip currents video monitoring conducted on an annual time scale, we found that wave boundary conditions are responsible for rip currents development. In this way, considering H s , T m , and • m , we could define the RH-index for the considered beach. An interesting point is the different hazard levels observed considering annual time scale (high) and tourist season time scale (moderate). Indeed, even if the tourist season is generally the reference value concerning beach safety, it is important to highlight how Ligurian beaches (and Mediterranean beaches in general) are also usually frequented in low and medium season by bathers and by those who practice water sports. This aspect is important in terms of beach safety and management. The use of wave parameters allows a quick connection between the rip hazard assessment and marine forecasts. Indeed, the proposed approach offers a quick evaluation of RHindex based on forecasted wave conditions. The RH-index can also be used to evaluate rip currents risk (RR). Therefore, as described in the methods section, a risk evaluation is possible using the tourist presence as an index of possible elements at risk. The matrix of risk obtained shows how risk increases with increasing the tourist presence. This point appears banal but is the essence of risk definition (Fell et al. 2008); this point supports the validity of this proposed approach.

Beach management implications
The large and rapidly expanding tourism activities in the coastal environment impose an urgent need of both beach safety and sea bathing risk zoning ). Following (Brander 1999;Pranzini 2015), local managers should provide information to beachgoers about the rips hazard to reduce hazard level for Alassio beach water users. Safety signs in correspondence to beach access points should be used in patrolled and unpatrolled beaches. Signalling is the cheapest action to prevent risk on beaches (Teixeira 2014; Mucerino et al. 2019). For a higher efficiency, the symbols used should outline different situations according to ISO 201712 international directives and considering local regulations. Rip currents often occur in moderate sea conditions, precisely where and when beachgoers feel more confident. Many visitors from inland areas and international tourists lack a certain sea experience and knowledge. In addition, under certain circumstances, they may not be able to read beach warning signs or understand the hazard flag system (Leatherman 2017). These features must be considered because rip currents often occur in tourist season when the beaches are most crowded. Moreover, in Italy, beach safety skills reached a high level, and lifeguards need new tools to obtain a further improvement in drowning risk prevention (Funari and Giustini 2011). The proposed approach could provide suitable application to evaluate rip hazard and beachgoers exposure to rip risk. Indeed, the proposed approach can be implemented into a weather forecasting system in order to provide rip currents hazard predictions. This paper also highlighted the necessity to create a national database of incidents caused by rip currents along the Italian coasts as already conducted in the UK (Scott et al. 2014) to allow a more reliable classification of rip currents risk.

Conclusions
In this paper, an expeditious rip hazard and risk assessment is proposed, based on probability theory, that considers wave boundary conditions and field data observation. In order to evaluate rip currents hazard, web camera data analysis was performed using timex elaboration on yearly observations and compared with offshore wave conditions. The results present a clear correspondence between rip currents development and wave parameters range.

3
The XBeach model was computed to obtain an indicative evaluation of the velocities magnitude of seaward flow in two cases during summer and winter.
Based on this analysis, it was possible to apply the probabilistic calculation, described in this work to identify the level of hazard. Furthermore, the risk for swimmers was determined considering the tourist presence. On Alassio beach, the hazard level (RH) is considered high on an annual time scale, but it is moderate during tourist season. In any case, the risk (RR) for bathers seems to be proportional to tourist presence. Unfortunately, due to the absence of reliable data, it was not possible to obtain a more detailed rip current risk evaluation. This method is suitable because the beach safety tool-a relevant component of beach management-can be employed by beach managers as a support to lifeguards. Official websites, apps and signs can be useful in providing information on rip current risks to local population and summer tourists. In Alassio town, this aspect becomes more important due to the high seaside tourist vocation of the area, and since many summer tourists are not aware of the danger they are actually exposed to. Thus, the proposed approach can be considered a valid tool for hazard and risk evaluation of the rip currents in the Mediterranean environment with a high potential of replication and usability in other contexts.