Monitoring, cataloguing, and weather scenarios of thunderstorm outflows in the northern Mediterranean

High sampling rate (10 Hz) anemometric measurements of the “Wind, Ports, and Sea” monitoring network in the northern Tyrrhenian Sea have been analysed to extract the thunderstorm-related signals and catalogue them into three families according to the different time-scale of each event, subdivided among 10 min, 1, and 10 h events. Their characteristics in terms of direction of motion and seasonality/daily occurrence have been analysed: the results showed that most of the selected events come from the sea and occur from 12:00 to 00:00 UTC during the winter season. In terms of peak wind speed, the strongest events all belonged to the 10 min family, but no systematic correlation was found between event duration and peaks. Three events, each one representative of the corresponding class of duration, have been analysed from the meteorological point of view, in order to investigate their physical nature. According to this analysis, which was mainly based on satellite images, meteorological fields obtained from GFS analyses related to convection in the atmosphere, and lightning activity, the thunderstorm-related nature of the 10 min and 1 h events was confirmed. The 10 h event turned out to be a synoptic event, related to extra-tropical cyclone activity.

Abstract. High sampling rate (10 Hz) anemometric measurements of the "Wind, Ports, and Sea" monitoring network in the northern Tyrrhenian Sea have been analysed to extract the thunderstorm-related signals and catalogue them into three families according to the different time-scale of each event, subdivided among 10 min, 1, and 10 h events. Their characteristics in terms of direction of motion and seasonality/daily occurrence have been analysed: the results showed that most of the selected events come from the sea and occur from 12:00 to 00:00 UTC during the winter season. In terms of peak wind speed, the strongest events all belonged to the 10 min family, but no systematic correlation was found between event duration and peaks.
Three events, each one representative of the corresponding class of duration, have been analysed from the meteorological point of view, in order to investigate their physical nature. According to this analysis, which was mainly based on satellite images, meteorological fields obtained from GFS analyses related to convection in the atmosphere, and lightning activity, the thunderstorm-related nature of the 10 min and 1 h events was confirmed. The 10 h event turned out to be a synoptic event, related to extra-tropical cyclone activity.

Introduction
A primary aim of politics, science, and technology is to pursue the safety and cost-efficiency of built-up environments exposed to natural hazards by forecasting catastrophic events, assessing their precursors, and evaluating and managing the risks they cause. Wind is the most destructive natural phenomenon, over 70 % of the damage and deaths caused by nature are due to the wind (Tamura and Cao, 2012;Ulbrich et al., 2013), so the evaluation of its actions is crucial to guarantee the safety and limit the cost of structures, infrastructure and territory management, thus representing a cornerstone in the field of civil engineering and atmospheric sciences and a societal need.
A climatological condition in which wind phenomena of different natures coexist, for example, extra-tropical and tropical cyclones, monsoons, tornadoes, downslope winds, and thunderstorms, is referred to as a mixed wind climate (Gomes and Vickery, 1978). Extra-tropical cyclones are the most typical events that strike the mid-latitude areas. Tropical cyclones, monsoons, tornadoes, and downslope winds are key features of the wind climate of specific zones. Thunderstorms occur almost everywhere. The European wind climate and that of many countries at the mid-latitudes is dominated by extra-tropical cyclones (Deroche et al., 2014) and thunderstorms .
The polar front theory (Bjerknes and Solberg, 1922) explains and describes the genesis and life cycle of extratropical cyclones. These are synoptic phenomena that develop in a few days on a few thousand kilometres (Fig. 1a).
The surface velocity field in extra-tropical cyclones is characterized by a mean wind profile in equilibrium with an atmospheric boundary layer whose depth is of the order of magnitude of 1-3 km. Here, within time intervals between 10 min and 1 h, turbulent fluctuations are stationary and Gaussian. Davenport (1961) identified the most intense wind events with the extra-tropical cyclones and introduced a model, based on this hypothesis, to determine the wind loading  of structures. After over half century from the Davenport's model definition, it is still a foundation of wind engineering (Simiu and Scanlan, 1996;Holmes, 2015).
The modern study of thunderstorms started when Byers and Braham (1949) proved that these events are mesoscale phenomena that develop over a few kilometres (Fig. 1b). They consist of convective cells that evolve in about 30 min through three stages in which an updraft of warm air is followed by a downdraft of cold air. Fujita (1985) showed that the transient downdraft that impinges on the ground produces radial outflows (Fig. 2) that, even though in most cases do not produce really damaging winds, can be as high as about 75 m s −1 (Fujita, 1990). The whole of these air movements is called downburst and is referred to as a macroburst or a microburst depending on whether the downdraft diameter is greater or smaller than 4 km, respectively. Radial outflows exhibit non-stationary and non-Gaussian wind speed properties, and a vertical "nose profile" that increases up to about 50-100 m height, then decreases above. These studies gave rise to an extraordinary fervour of research in atmospheric science (Goff, 1976;Wakimoto, 1982;Hjelmfelt, 1988).
However, despite this huge amount of research, this matter is still dominated by large uncertainties; furthermore, there is not yet a shared model of thunderstorm outflows and their actions on structures like that developed by Davenport (1961) for extra-tropical cyclones. This happens because the complexity of the thunderstorm downbursts makes it difficult to establish physically realistic and simple engineering schemes, their short duration and small size means few data are available, and a large gap exists between wind engineering and atmospheric science. It follows that the wind loading of structures is still evaluated by the Davenport's model without any concern for the real nature and the properties of the meteorological event that causes the loading. This is nonsense because extra-tropical cyclones and thunderstorm outflows are different phenomena that need separate assessments (Solari, 2014).
The research on thunderstorm outflows carried out by the Wind Engineering and Structural Dynamics (WinDyn) Research Group (www.windyn.org, last access: 29 August 2018) at the University of Genoa takes cue from two European Projects, i.e. "Wind and Ports" (WP) (2009)(2010)(2011)(2012)   and "Wind, Ports and Sea" (WPS) (2013-2015) , financed by the European Cross-border Cooperation Program "Italy-France Maritime 2007-2013". They handled the problem of the safe management and risk assessment of North Tyrrhenian seaport areas with respect to strong wind conditions through a joint cooperation between the Windyn group -the unique scientific partner in these projects -and the port authorities of Genoa, Savona -Vado Ligure, La Spezia, Livorno (Italy) and Bastia -L'Île-Rousse (France). In this framework, a wide in situ wind monitoring network has been created that is generating an unprecedented amount of high quality wind measurements. In addition to their institutional role of supporting the activities of port authorities, they represent an unlimited source of information for carrying out scientific research in several different fields.
The analysis of these data shows the presence of recordings due to wind phenomena of a different nature, namely extra-tropical cyclones, thunderstorms outflows, and intermediate events (Kasperski, 2002;Zhang et al., 2018). Thus, in order to focus on the study of intense thunderstorm outflows, a semi-automatic procedure was implemented to recognize and extract these phenomena (De Gaetano et al., 2014). This approach is consistent with previous procedures developed and calibrated in order to process a huge amount of data, based on a few synthetic elements, derived from the sole anemometric recordings (Riera and Nanni, 1989;Twisdale and Vickery, 1992;Choi and Tanurdjaja, 2002;Kasperski, 2002;Duranona et al., 2006;Lombardo et al., 2009), without carrying out systematic and prohibitive meteorological surveys of the weather scenarios out of which they took place. According to this criterion, an extensive set of records labelled as thunderstorms has been gathered and subjected to probabilistic signal analyses aiming at evaluating their main properties relevant to the wind loading of structures (Solari et al., 2015a;Zhang et al., 2018). These properties have formed the base on which two novel methods have been proposed to determine the structural response to thunderstorm outflows (Solari et al., 2015bSolari, 2016).
Despite its inherent advantages and merits, this approach has a shortcoming in that it misses the knowledge of the weather scenarios that occur during events classified as thunderstorms, without recognizing their actual meteorological nature. In order to take the first step towards filling this gap, the thunderstorm downburst that occurred on 1 October 2012 over Livorno, Italy, was selected as a reference test case (Burlando et al., 2017a). Detailed analyses were carried out of the wind speed and direction records detected by the WP and WPS network. In parallel, the atmospheric conditions concurrent with this event were studied in great detail by gathering all the meteorological data available in this area, which included model analyses, standard in situ measurements (stations and radio-soundings), remote sensing techniques (radar and satellite), proxy data (lightning), and direct observations (from the European Severe Weather Database, Dotzek et al., 2009). This information led to reconstructing the weather scenario, classifying this event as a wet downburst, determining its space-time evolution, and embedding in this framework signal analyses that aim to extract the key parameters for determining the wind loading of structures.
From this point of view, the above study may become a reference model to carry out comprehensive analyses of the major thunderstorm events recorded by the WP and WPS monitoring network and by any network in general. However, this is only relevant when repeating such investigations for several events and elaborating the results in a probabilistic framework aiming to construct suitable models of the wind field of thunderstorm outflows truly related to different classes of meteorological events. Such a probabilistic model should be based on the recognition of clear meteorological precursors of these phenomena and their main parameters, and on the statistical analysis of their occurrence in terms of touch-down position, size, duration, intensity, motion pattern, and background flow (Mohr et al., 2017). These data are in turn the key information for implementing advanced thunderstorm models addressed to hazard and risk analyses for a broad spectrum of applications.
Unfortunately, the framework depicted above may represent a utopian prospect due to the burden of collecting so many data from different sources and performing their joint analysis for several events. Finding a reasonable balance between expeditionary evaluations based on the sole wind records detected and studies that encapsulate the above information within detailed meteorological surveys is a very difficult and ambitious aim. This paper represents a first step and a pilot attempt in this direction. Section 2 describes the WP and WPS wind monitoring network and dataset. Section 3 illustrates the separation and classification procedure applied in order to gather a rich sub-dataset of records labelled as thunderstorm outflows. Section 4 provides a new procedure to extract and catalogue thunderstorm outflows with reference to their duration and intensity, these parameters being the key features for evaluating the wind loading of structures. Section 5 furnishes additional elements on the direction, seasonality, and hour of daily occurrence of these events. Section 6 introduces a synthetic meteorological survey, coherent with but easier than the general methodology described by Burlando et al. (2017a), of the weather scenarios corresponding to three events preliminarily labelled as "thunderstorms" and characterized by different lifetime scales. Section 7 summarizes the main conclusions and highlights some prospects for future research.

Monitoring network and dataset
The wind monitoring and measurement through direct-and remote-sensing techniques is essential for a broad range of scientific disciplines including atmospheric physics, meteorology, climatology, wind, and civil and environmental engineering. Besides, the analysis of wind data is often carried out in the framework of a broad range of applications con-cerning both the interpretation of wind phenomena that occurred in the past, i.e. in climatological analyses, as well as to forecast future scenarios, e.g. in climate change surveys (Ortego et al., 2014;Dawkins et al., 2016) or in the study of wind hazards, damage and risks to which buildings, infrastructure, and human activities in contact with Earth's atmosphere are exposed.
The WP project gave rise to a wind monitoring network made up of 23 ultrasonic anemometers (yellow circles in Fig. 4) distributed in the ports of Genoa (2), La Spezia (5), Livorno (5), Savona -Vado Ligure (6), and Bastia (5). WPS, which represents the continuation and development of WP, enhanced and enlarged this network by means of five ultrasonic anemometers installed in the ports of Savona (1), La Spezia (1), Livorno (1), and L'Île-Rousse (2) (orange triangles in Fig. 4), three weather stations,each one including an additional ultra-sonic anemometer, a thermometer, a barometer, and a hygrometer, in the ports of Genoa (1), Savona (1), and Livorno (1) (blue diamonds in Fig. 4), and three lidar (light detection and ranging) wind profilers, again in the ports of Genoa (1), Savona (1), and Livorno (1) (red squares in Fig. 4), which detect the vertical wind profile at 12 heights from 40 to 250 m above the ground level (a.g.l.). Other sensors installed autonomously by individual port authorities are now in the stage of becoming integral parts of the WP and WPS network.
The ultra-sonic anemometric stations consist of bi-axial or three-axial sensors that detect the wind speed and direction with a precision of 0.01 m s −1 and 1 • , respectively. Their sampling rate is 10 Hz, with the exception of one sensor installed in Savona, which samples at 1 Hz, and the sensors installed in the ports of Bastia and L'Île-Rousse, whose sampling rate are 2 Hz. Their position was chosen in order to cover homogeneously the port areas involved in these projects and to register undisturbed wind velocity histories. The instruments are mounted on high-rise towers or on antenna masts at the top of buildings, paying attention to avoid local effects that may contaminate the quality of measurements. Table 1 illustrates the most important features of the 31 anemometers that the current WP and WPS network is made up of, with h being their height a.g.l. that varies from 10 to 84 m.
A set of local servers placed in each port authority's headquarters receive the data recorded by anemometers and lidars in its own port area and elaborates the basic statistics on 10 min of averaging period, namely the mean and peak wind velocities and the mean wind direction. Each server automatically sends this information to the central server at the University of Genoa, which stores the data in a database through a four-step automatic procedure: 1. raw data are systematically checked and validated; 2. 10 min statistics, including turbulence intensity and gust factor, are evaluated from raw data and stored in the database; 3. 1 min statistics, including turbulence intensity and gust factor, are evaluated and stored in the database; and 4. an automatic check report is produced and sent to port authorities daily.
These results are the main outcome transferred to port authorities for their institutional activity. At the same time, they represent the starting point for a broad range of researches carried out by the WinDyn Research Group. In this regard, starting from the original database, a new one has been created that collects further statistical parameters of the anemometric measurements. In particular, one record for each subsequent T = 10 min period is stored that gathers such parameters into three groups: 1. 1 s peak wind velocityV , mean wind velocity V m10 , mean wind direction α m10 , gust factor G 10 =V /V m10 , turbulence intensity I 10 , skewness γ 10 , and kurtosis κ 10 in the interval T ; 2. mean wind velocity V m60 , gust factor G 60 =V /V m60 , turbulence intensity I 60 , skewness γ 60 , and kurtosis κ 60 in the 1 h time interval centred around T ; 3. maximum mean wind velocity averaged over 1 min V m1 and gust factor G 1 =V /V m1 in T .
This dataset represents the starting point for the separation and classification procedure described in the next section. Similar analyses are currently in progress with regard to lidar measurements (Burlando et al., 2017b).

Separation and classification procedure
A thorough examination of the huge amount of data collected during the WP and WPS Projects reveals that intense wind events can be separated and classified into three families characterized by different statistical properties (De Gaetano et al., 2014).
1. Stationary Gaussian events with relatively large mean wind velocities and small gust factors; these usually correspond to synoptic neutral atmospheric conditions and are here referred to as extra-tropical cyclones or depressions ( Fig. 5).
2. Non-stationary non-Gaussian events with relatively small mean wind velocities, large and quite isolated peaks, and high gust factors; these are here referred to as thunderstorms outflows (Fig. 6).
3. Stationary non-Gaussian events with relatively small mean wind velocities, large and repeated peaks, and moderately high gust factors; these are here referred to as intermediate events ( Fig. 7) or gust fronts (Kasperski, 2002). While waiting to carry out a systematic meteorological survey and interpretation of these events, it seems to be reasonable to advance the hypothesis that they are associated to strongly unstable atmospheric conditions, or downslope winds, or recirculating vortices as well.
The separation of intense wind phenomena into homogeneous families is a key topic to interpret the events of engineering interest and to deal with them by models coherent   with their physical reality. This is possible by inspecting each event and the weather scenario in which it occurs by merging the analysis of the anemometric recordings with that of the meteorological data detected in the same area and concurrent with the event considered (Burlando et al., 2017a). This operation is clearly not possible when many instruments, many years of measurements, and many different events have to be examined as in this case. To achieve the separation and classification of different intense wind events as easily and efficiently as possible, De Gaetano et al. (2014) developed a semi-automated procedure applied to the records with peak wind velocity (V ) greater than 15 m s −1 . This choice is coherent with thunderstorm analyses carried out by other authors (Choi, 2000(Choi, , 2004Duranona et al., 2006) and with the tradition of evaluating the parameters of synoptic events by collecting all records that satisfy the requirement of neutral atmospheric conditions (Solari and Piccardo, 2001;Solari and Tubino, 2002), including several phenomena of limited engineering interest. The alternative approach of restricting analyses to thunderstorm outflows with higher peak values (Geerts, 2001;Lombardo et al., 2014) improves the information related to the phenomena of major engineering interest, but reduces the statistical representativeness of the results.
This semi-automated procedure involves a suitable mix of systematic quantitative controls and qualitative expert judgments. The quantitative controls are based on the comparison between the detected values of the gust factors, G 60 , G 10 , G 1 , and their reference values, G 0 60 , G 0 10 , G 0 1 , evaluated by means of numerical simulations (ESDU, 1993;Burlando et al., 2007Burlando et al., , 2013 assuming that intense wind speeds occur in neutral atmospheric conditions during synoptic extra-tropical cyclones with stationary Gaussian properties. Conceptually, a wind event is labelled as an extra-tropical cyclone if the ratios G 60 /G 0 60 , G 10 /G 0 10 , G 1 /G 0 1 are small. Instead, it is labelled as a thunderstorm outflow or an intermediate event when the ratio G 10 /G 0 10 is large. In this way, De Gaetano et al. (2014) automatically identified over 99.5 % of the records related to intense wind events as synoptic extra-tropical cyclones. The remaining ones were submitted to qualitative ex-pert judgments. These analyses were based on the wind speed and direction raw data detected by nine anemometers in the period 2011-2012.
Figures 5, 6, and 7 show three typical 1 h records registered by the anemometer 3 of the port of La Spezia that correspond to an extra-tropical cyclone, a thunderstorm outflow and an intermediate event, respectively. Panels (a) and (b) show the time-series of the wind speed and direction raw data, respectively, and they also provide their mean values over 1 h, V m60 and α m60 (solid lines), and 10 min periods, V m10 and α m10 (dotted lines). In panel (a), the 1 s peak wind speed (red circles) is also shown, which is obviously smaller than the instantaneous peak. Panel (c) shows the ratios G 60 /G 0 60 , G 10 /G 0 10 , and G 1 /G 0 1 over subsequent 10 min periods.
The extra-tropical cyclone record shown in Fig. 5 is characterized by a relatively high mean wind speed (V m10 = 12.04 m s −1 , V m60 = 12.21 m s −1 ) and gust peak (V = 20.46 m s −1 ). The gust factor (G 10 = 1.70, G 60 = 1.68) is rather high but typical of synoptic neutral atmospheric conditions. Likewise, the wind speed and also the wind direction involve stationary features. The ratios G 60 /G 0 60 = 0.80, G 10 /G 0 10 = 0.85, and G 1 /G 0 1 = 0.85 are well below 1. The thunderstorm outflow record shown in Fig. 6 is characterized by a relatively low mean wind speed (V m10 = 6.60 m s −1 , V m60 = 2.99 m s −1 ), a relatively high gust peak (V = 19.61 m s −1 ) and a very high gust factor (G 10 = 2.97, G 60 = 6.56). The ratio G 60 /G 0 60 = 3.69 exhibits a sudden increase in correspondence of the gust peak whereas the ratios G 10 /G 0 10 = 1.78 is well above 1 and G 1 /G 0 1 = 0.91 is larger than 0.80. The wind direction changes of almost 180 • as usually occurs when a downburst passes over the anemometer (Orwig and Schroeder, 2007).
The intermediate event record shown in Fig. 7 is characterized by a relatively low mean wind speed (V m10 = 9.71 m s −1 , V m60 = 11.85 m s −1 ) and a rather high gust peak (V = 23.75 m s −1 ); the gust factor (G 10 = 2.45, G 60 = 2.00) is much greater than the typical values in neutral atmospheric conditions, this being confirmed by the ratios G 60 /G 0 60 = 1.22, G 10 /G 0 10 = 1.58, and G 1 /G 0 1 = 1.17. Both the wind speed and the direction exhibit quite regular trends without apparent transient features.
It is worth noting that the records shown in Figs. 5-7 have clear trends that do not imply doubtful decisions. Unfortunately, this ideal condition does not always occur and the extraction of the thunderstorm outflows carried out by De Gaetano et al. (2014) may give rise to controversial choices. In contrast, though not availing of records detected with high sample rates, Duranona (2015) made a fine selection of the convective events that take place in Uruguay by inspecting 10 min mean and peak wind speeds over a period 10 h. This study inspired a re-calibration of the procedure applied by De Gaetano et al. (2014) based on 10 min, 1, and 10 h records. Joined with the possibility of processing a more extensive set of measurements, this approach leads to a novel extraction and cataloguing procedure described in the next section.

Thunderstorm extraction and cataloguing
As mentioned in the previous section, the above procedure is firstly extended to the data recorded by 14 ultrasonic anemometers in the period 2011-2015, including the nine anemometers previously inspected in the period 2011-2012. Secondly, following the method suggested by Duranona (2015), the previous analyses have been improved by performing the qualitative selection of thunderstorm outflows based not only on 10 min and 1 h records as in De Gaetano (2014), but also inspecting the 10 h records centred around the peak wind speed. Table 2 shows a general framework of the anemometers, periods and data herein examined. NTE = 198 and NTR = 277 are the number of events and the number of records labelled as thunderstorms, respectively; NTR is always greater than NTE as the same thunderstorm event may be detected simultaneously by more than one anemometer.
It is worth noting that the percentage and the total amount of valid data are quite different at each experimental site. This depends first on the successive installation of sensors, then on the different periods in which measurements were not carried out due to accidents or malfunctions of instruments, including some cases in which they have not been restored yet; there are also periods in which measurements have been judged not reliable enough to be examined and have been disregarded (Cook, 2014a, b). It is also worth noting that the databases of several anemometers have not been analysed yet, that the monitoring network continuously produces new data and that new anemometers are progressively added to make the network richer and richer.
The analysis developed here involves not only the assemblage of a more comprehensive and controlled thunderstorm outflow dataset but, more importantly, a major advance in understanding the time-scale of transient events. Accordingly, it favours the classification and cataloguing of such events based on the duration of the transient peak and on the wind speed itself. This information is crucial to investigate the loading and response of structures.
In this regard, transient records have been separated into three families depending on whether the presence of a rampup and the transient peak are clearly detectable in 10 min (Fig. 8), 1 h (Fig. 9), or 10 h (Fig. 10) records; for sake of simplicity they are referred to as 10 min, 1 h, and 10 h events. It was found that 50.9 % of the extracted transient records are detectable on 10 min periods, 38.3 % of them can be recognized on 1 h periods whereas only 10.8 % are pointed out by inspecting 10 h records. This aspect reflects on the duration of the ramp-up and has a key engineering role. As demonstrated by Kwon and Kareem (2009), the structural response increases after reducing the length of the impulse related to   the passage of the gust front. In Figs. 8, 9, and 10, panels (a, b), (c, d) and (e, f) refer to 10 min, 1, and 10 h records, respectively, centred around the gust peak; panels (a, c, e) and (b, d, f) correspond to wind speed and direction, respectively. In all these diagrams the time variation of the wind direction reflects the time variation of the wind speed, exhibiting a change that may be perceived at the same time-scale in which the ramp-up and the transient peak are perceived. However, not all the gathered records have the same property. The diversity between this approach and typical meteorological surveys is apparent. For instance, according to the Federal Meteorological Handbook No. 1 (NOAA, 2005), "the beginning of a thunderstorm is to be reported as the earliest time: (1) thunder is heard; (2) lightning is observed at the station when the local noise level is sufficient to prevent hearing thunder; or (3) lightning is detected by an automated sensor". Conversely, "the ending of a thunderstorm shall be reported as 15 minutes after the last occurrence of any of the above criteria". The gathering of high-quality wind speed records makes it possible to introduce a diverse duration criterion that represents a key advance for hazard analyses, wind loading modelling, structural response and sensitivity to thunderstorm events. In addition to the above classification criterion, transient records are separated into four groups as a function of the peak wind speed (in m s −1 ), namely 15 ≤V < 20, 20 ≤V < 25, 25 ≤V < 30, and 30 ≤V < 35. This classification is aimed at recognizing the existence of any correlation between the duration of the gust front and its wind speed; shortduration events with high wind speed are the most relevant hazard from a structural viewpoint. Figure 11 shows the example of a thunderstorm record for each wind speed class. Table 3 shows the results of intersecting the two classification criteria based on the duration and the wind speed pointing out, at least for the available data, no systematic correlation between the duration of the most intense part of the record and the peak wind speed. However, it is worth noting that the four events whose peak wind speed exceeds 30 m s −1 are characterized by a limited duration. In contrast, the 21 events whose peak wind speed exceeds 25 m s −1 are quite uniformly distributed over different durations.

Thunderstorm direction, seasonality, and hour of daily occurrence
As shown by Figs. 5-10, the definition of the thunderstorm outflow direction is quite controversial. A thunderstorm cell may be classified as stationary or non-stationary depending on whether it has a translational motion. In the case of a stationary event the direction of the thunderstorm outflow, of radial nature, strictly depends on the position of the axis of the downdraft, usually assumed as vertical, with regard to the position of the sensor. In the case of a non-stationary event, the velocity and direction detected by the sensor is the vector composition of the velocity and direction of the thunderstorm cell, dealt with as stationary, and its translational components (Holmes and Oliver, 2000). The situation becomes more complex in the frequent case in which a thunderstorm cell is embedded in a background larger-scale boundary layer flow field, usually of synoptic type, or ever more in the case in which multiple downdraft are generated by single or multiple thunderstorm cells. In principle, all these situations may be treated by vector compositions of the velocity and direction of component flows; in reality, no proof exists that this approach is physically and mathematically suitable. Lombardo et al. (2009) identified the thunderstorm outflow direction with its average value in a 5 s period centred around the peak wind speed. Solari et al. (2014) defined it as the average value in a 30 s period centred around the gust peak. In this paper, the period in which the wind direction is averaged is increased to the 1 min centred around the peak. Figure 12 shows the distribution of the thunderstorm records reported in Table 2 with reference to the day of occurrence shown on the polar radius (1 relates to 1 January) and to the wind direction (0 • refers to the north). The background of every diagram is the map of the port, roughly illustrating the relationship between the directions of intense outflows and geographic conditions. Most of the events (68 %) occur in the months between September and January. Besides, most of the events (78 %) are characterized by wind directions coming from the sea. In the port of Savona only one year of reliable data has been gathered, so the number of thunderstorm records is small and no definite trend is identified. However, it is possible that the spread of these results may be due not only to the scarcity of data but also to the fact that the port of Savona contains two different port areas: Savona itself and that of Vado Ligure, though these areas are rather close, they have different features. Also, the presence of intense downslope winds in the Vado Ligure area (Burlando et al., 2017c) may contribute to the increase of the spread. Figure 13 shows the distribution of thunderstorm records (NTR) corresponding to the three families defined above, namely 10 min, 1, and 10 h events, reported in Table 2 with respect to the day time. Overall, for the area under study, 20 % of the thunderstorms occur between 00:00 and 06:00 UTC, 24 % occur in the morning between 06:00-12:00 UTC, 30 % occur in the afternoon from 12:00 to 18:00 UTC and 26 % occur in the evening from 18:00 until 00:00 UTC. Therefore, at least in the monitored area, it seems that thunderstorm events are slightly more likely to occur in the warmer hours of the day, possibly because of the role of solar heating in triggering thermals from the earth's surface. The percentage of 1 h events grows in the afternoon with respect to the 10 min events as well.
However, it is worth noting that the occurrence of thunderstorms during night hours is an open topic that needs further study. These phenomena are most likely to form when the temperature of the air decreases with height pretty rapidly, i.e. when it is hot at the ground and cold aloft. Thunderstorms that form at night occur in the absence of heating at the ground by the sun, so that they are due to different forcing mechanisms. One of the possible explanations is that during night, sea surface is warmer than land due to the higher thermal capacity of water. Therefore, the land-to-sea breezes bring cold air from land to sea and the advected air becomes statically unstable thus resulting in convection. That convection, if strong enough, should be one of the main contributors to the development of thunderstorms over sea surface.

Short meteorological survey and weather scenarios
As previously noted, the procedure depicted by Burlando et al. (2017a) may represent a reference model to carry out comprehensive analyses of the major thunderstorm events.
However, its burden is so high that it makes this procedure realistically unusable for systematic analyses of historical series of such events. Hence, the need arises, or at least the objective, to develop a faster approach that integrates the data provided by an anemometric network such as the WP and WPS ones, with few essential meteorological information that may qualify, albeit preliminarily, the transient intense wind events detected by the network.
In this spirit, making treasure of the experience matured during the detailed analysis of several downbursts, in particular the one that stroke Livorno on 1 October 2012 (Burlando et al., 2017a), this section describes simplified meteorological surveys and preliminary reconstructions of the weather scenarios that occurred during the three events presented in Sect. 4, Figs. 8-10, and already described from the signal analysis viewpoint. In the following three subsections, the meteorological conditions that brought about all these intense wind events, whose characteristic lifetime-scales are 10 min, 1, and 10 h, are briefly reported one by one. The analysis is performed on two spatial scales: the meteorological conditions at the synoptic-scale are firstly inspected in order to evaluate the pattern of cyclones, anticyclones, and fronts that determined instability and cloudiness in the atmosphere; the phenomena possibly occurring at the meso-scale that developed over the areas of interest are then investigated in or-der to understand the specific convective structures present during the events under consideration. The whole analysis is based primarily on the following data: the Global Forecast System (GFS) analyses, obtained from the National Center for Environmental Prediction (NCEP) through the National Centers for Environmental Information database; the cloud top height, obtained from the cloud analysis performed by Eumetsat (EUMETSAT, 2013;Derrien et al., 2013) based on infrared measurements collected by SEVIRI (Spinning Enhanced Visible & Infrared Imager) on board Meteosat Second Generation (MSG) satellites; the lightning strikes, obtained from the Blitzortung database.
GFS analyses are available worldwide on a 0.5 • by 0.5 • geographical grid with 6 h time step; these spatial and temporal resolutions are suitable to evaluate the movement of the large-scale structures that determine the evolution of the weather conditions. The satellite measurements of Meteosat 10, which is the one used here, are available every 15 min as a full disk imagery centred at 0 • of longitude and latitude, with a spatial resolution over Europe of a few kilometres. The spatial and temporal resolution of such data allow for distinguishing the convective structures at the meso-scale and their evolution in time with sufficient accuracy also in the case of thunderstorms, whose typical scales are of the order of 10 km and 1 h. Finally, the lightning strikes are used to confirm whether the convective activity can be associated to cumulonimbus clouds, which are typical of thunderstorms.
The meteorological analysis based on such data cannot provide precise information concerning the very small-scale structure of the thunderstorm considered, like the shape of its gust front or the area covered by heavy precipitation. This information can be obtained using higher-resolution measurements like radar imagery or finer monitoring networks (Burlando et al., 2017a), and possibly high-resolution numerical simulations (Lompar et al., 2017). However, as shown in the following sections, the meteorological analysis described above may be sufficient to state the nature, i.e. convective or synoptic, of the transient intense wind events investigated, which is the main focus of the present paper.

Weather analysis of the 10 m event on 25 October 2011 in La Spezia
On 25 October 2011, the exceptionally large Anticyclone Ulla (names given by the Institute of Meteorology of the Freie Universität Berlin, Germany) located over Eastern Europe, with pressure maxima greater than 1030 hPa, determined the blocking of the zonal shift of Cyclone Meeno, which remained approximately stationary to the west of Ireland, extending its low-pressure minimum of 985 hPa northward to Iceland, as shown in Fig. 14a. The cold front of Cyclone Meeno, which extended meridionally to northern Africa passing over the Alps, was moving slightly westward over the Mediterranean during the day. The warm sector uplift ahead of Meeno's cold front, due to the south-easterly winds forced over southern Italy and the Adriatic Sea, by the extension of the influence of Ulla to the Balkans Region, determined a wide area of cloudiness all over northern and central Italy, shown in Fig. 14b. The distribution of cloud top heights in this figure shows that the prevailing southerly flow of warm and humid air from the Mediterranean favoured the development of strong instability, associated with deep convection, which in turn locally determined very intense thunderstorm events. The advection and uplift of moist air from the southern quadrants in the lee of the Alps is confirmed by values of relative humidity higher than 100 % over a large part of the Padan plain, eastern Liguria and southward along a convergence line that corresponds approximately to the frontal zone beneath. According to the GFS analysis, Fig. 15a shows that saturated air conditions, i.e. RH values (shaded contours) equal to 100 %, were occurring exactly above La Spezia, at 18:00 UTC. Moreover, the values of mean storm motion in the 0-6000 m a.g.l. in the same area were between 10 and 15 m s −1 from south or southwest. According to satellite measurements, Fig. 15b shows that a deep convective cloud with the top at about 12 000 m above the sea level (a.s.l.) occurred at around 15:45 UTC over La Spezia. However, during the morning of 25 October, a series of thunderstorms developed slightly to the northwest of La Spezia because of the orographic forcing determined by the south-eastern flow, as measured by the anemometers in the port of La Spezia (see Fig. 8e, f). In that area, many meteorological stations of the Liguria region (e.g. Monterosso, Serò di Zignago, Levanto San Gottardo, Brugnato) reported precipitation rates larger than 50 and even 100 mm h −1 from 09:00 to 15:00 UTC. The thunderstorm over La Spezia at 15:40 UTC was the last of this series, just before the wind shifted to the north when the cold front overcame the Western Alps. The sudden change in wind speed and direction reported in Fig. 8 represents the transition between thus higher-to-lower wind speed regimes. It is worth noting that, as reported in Zhang et al. (2018), the transition between these kinds of wind regime is usually slower, so that these events are often classified as 1 h events. In the present case, the dynamics of the transition is particularly fast so that the event is classified as a 10 min one.
As far as the cloud-ground strikes measured by the Blitzortung network are concerned, this convective cell showed only a moderate lightning activity, which confirms the convective nature of this phenomenon. The lightning occurrence on 25 October 2011, shown in Fig. 16, in the interval from 15:25 to 15:55 UTC (half an hour centred with respect to the maximum wind speed recorded during this event) was 164 strikes, with a slightly increasing frequency during the first 25 min. The meteorological situation over Europe was dominated by the presence of the anticyclone Netti, with its high-pressure maximum of 1025 hPa situated over Ukraine and Russia, indicating a blocking condition (Rex, 1950). The cyclone Quirin, which began on 2 October to the north of the Pyrenees, where a strong thermal contrast between a colder maritime Atlantic air mass to the north and the much warmer and moist tropical air to the south occurred, had slightly moved its low-pressure minimum north-eastward over northern France on 3 October and Belgium on 4 October. According to GFS analyses, at 06:00 UTC the tropopause height showed an abrupt discontinuity to the north of the Alps with a roughly meridional gradient ranging from a maximum 14 214 m to a minimum 9119 m, denoting the existence of a frontal zone beneath, consisting of a warm core of tropical air southward (red contour reported in Fig. 17a) and a colder core of Atlantic air northward (blue contour in Fig. 17a). Cyclone Quirin, which was relatively small but very active because of the strong thermal gradients between Atlantic and tropical air masses, led to intense local thunderstorms and heavy precipitation in southern France and northern Italy during the night of 3 and 4 October. The distribution of cloud top heights shows clearly the presence of the cyclone Quirin, with its surface low-pressure minimum located over Belgium at 06:00 UTC on 4 October (Fig. 17b) and the occluded front extending southward to the Adriatic Sea.
In the lee side of the western Alps, the presence of a cold front extending from the Gulf of Genoa to the aforementioned occluded front in the northern Adriatic Sea is revealed  Figure 18a shows that RH values (shaded contours) close to saturation that occurred along a narrow band in the northern Tyrrhenian Sea, according to the GFS analysis at 06:00 UTC. Vectors in Fig. 18a, which represent the mean storm motion in the 0-6000 m a.g.l., show that a storm developing along the frontal boundary would eventually move from the sea towards the coast of Tuscany at approximately 5 m s −1 or less. This is the case of the thunderstorm shown in Fig. 9 that developed on 4 October early in the morning off the coast of Livorno. The sequence of three satellite images reported in Fig. 18b-d resembles the development of a single-cell thunderstorm caused by a cumulonimbus cloud that started developing its cumulus stage around 04:30 UTC (Fig. 18b) over Livorno, then reached the mature stage with a cloud top height of about 12 000 m a.s.l. at 05:15 UTC (Fig. 18c), and finally gradually dissipated while moving slowly farther inland (Fig. 18d). Because of the low storm advection from sea to land, the anemometric signals recorded in the port of Livorno have caught the whole three stages of the thunderstorm evolution on a time-scale, i.e. 1 h, which is approximately coincident to the typical order of magnitude of a single-cell thunderstorm life-cycle.
The intense convective activity that occurred in the surroundings of Livorno because of this thunderstorm is confirmed by the great number of lightning strikes registered by the Blitzortung network, as reported in Fig. 19. The lightning occurrence from 05:00 to 05:30 UTC was almost 1500 strikes, quite regularly distributed during the 30 min (see the frequency histogram in the bottom right corner). Note that over Livorno the whitish symbols, which correspond to times closer to 05:30 UTC, are shifted slightly eastward with respect to the reddish ones, according to the slow thunderstorm advection from sea to land.

Weather analysis of the 10-hour event on 21
November 2013 in Genoa The low-pressure system known as Quentin was born on 18 November 2013 in the Baffin Bay, between Canada and Greenland, and moved zonally to the north of England. During the 20 November, it moved south-eastward under the influence of anticyclone Susanne I, located at mid-latitudes over the Atlantic Ocean, and approached Belgium on 21 November at about 00:00 UTC. After this, it moved farther to the south and determined a low-pressure minimum of 995 hPa in the lee side of the Alps that was over the Ligurian Sea at 12:00 UTC, as shown in Fig. 20a. The cold core aloft of Quentin, i.e. tropopause heights (shaded contours) below 8000 m, over France determined strong atmospheric instability and caused high precipitation rates in Belgium, Holland, and France while advecting meridionally during the day. However, in the lee of the Alps, the south-eastward motion of the cold front of Quentin did not induce very deep convection over Liguria, as demonstrated by the relatively low values of the cloud top shown in Fig. 20b, which are between 8000 and 10 000 m a.s.l. The position of the cold front at the surface in the lee of the Alps at 12:00 UTC is indicated by the high-RH values that extends as an arc-shaped band from the eastern Liguria to the west of Sardinia Island, shown in Fig. 21a. On 21 November in the morning, the front had just passed over the Alps and the secondary pressure minimum aloft determined the wind rotation from north to southwest over the Ligurian Sea. The mean storm motion (vectors in Fig. 21a), which in this case corresponds approximately to the mean flow in the lower half of the troposphere as the directional wind shear from 0 to 6000 m a.s.l. was rather low, was about 20-25 m s −1 from southwest, indeed. The strong forcing aloft was likely the main reason for the relatively sudden increase of wind speed recorded by anemometer 2 of the port of Genoa (Fig. 10) rather than some deep convective phenomenon that did not seem to occur according to the top height of clouds obtained from satellite data (Fig. 21b). This is also confirmed by the Blitzortung network, which did not record any strike in this area in between 06:00 and 14:00 UTC.

Conclusions and perspectives
The wind monitoring network realized for the European Projects "Wind and Ports" and "Wind, Ports and Sea" is an inexhaustible source of measurements that highlight the speed and frequency of transient events, likely of convective nature, often disregarded, especially in the past, from classical wind engineering. It is now recognized that they are crucial with respect to hazard assessments aimed at the safety of construction, infrastructure, and territory, but they are often still undistinguished from extra-tropical cyclones in most cataloguing procedures (Stucki et al., 2014).
In the first part of this paper (Sects. 2-5) a description of the main properties of the anemometric network and the database generated by it is provided. It then illustrates the procedure used to separate the records associated with different wind phenomena based on information, stationary and Gaussian features, typical of signal analysis but lacking meteorological contents. Thanks to this method, events labelled as thunderstorm outflows are selected and classified into three families according to the time-scale (10 min, 1 h, 10 h) over which the transient part of the wind speed develops. In addition, analyses related to speed, direction, seasonality, and hour of daily occurrence are also presented.
The second part of this paper takes cue from a detailed meteorological survey of the wet downburst that occurred on 1 October 2012 over Livorno (Burlando et al., 2017a). As the burden of this approach prevents its realistic application within systematic analyses of historical series of similar phenomena, an expeditionary procedure is codified and pro- posed herein to integrate the anemometric records with few essential meteorological features that at least qualify the convective or synoptic nature of different detected events. This procedure has been applied to three sample events referred to as 10 min, 1, and 10 h records.
The analysis of the two shorter events has confirmed their convective nature. However, these thunderstorms have different triggering mechanisms. The 10 min event was determined by the mechanical lift of maritime air exerted by the orography; its transition between higher-to-lower wind speed regimes corresponds to the passage of the cold front. The 1 h event was most probably brought about by the mechanical lift due to the cold front's southeastward movement. These remarks may provide preliminary motivations to the different temporal scales of fast transient events and stimulate research towards the comprehension of this delicate issue.
Conversely, the longer event, i.e. the 10 h one, turned out to be a synoptic phenomenon, endowed with a rapid evolu- tion, initially misclassified as a potential thunderstorm, and it should be likely catalogued as an extra-tropical cyclonerelated windstorm instead (Roberts et al., 2014). Even if this result cannot be generalized for the whole group of 10 h intense wind events, it raises the question whether in some particular cases these phenomena can really have a convective genesis. This question remains open and will deserve further and more systematic investigations in the future.
The meteorological analysis of the three events considered here may be helpful for clarifying the relation between the shape of anemometric signals, discussed by Zhang et al. (2018), and the underlying meteorological phenomena. However, it is worth noting that each particular shape could be determined in principle by more than one phenomenon, especially if different locations are considered. Therefore, from this perspective this analysis should become more systematic and should be repeated for different databases of recordings taken at different latitudes and in geographical contexts different from the coastal area considered here. In order to make the meteorological analysis simpler, it would be beneficial in the future to adopt more advanced techniques of data sharing between raw data stored by different institutions. A possible strategy is to create a so-called interplanetary file system between sparse databases, which is a protocol designed to create a content-addressable, peer-to-peer method of storing and sharing hypermedia in a distributed file system, connecting different computing devices with the same system of files. Alternatively, an easier solution is represented by the creation of a program to automatically access databases that contain meteorological data and download the desired data automatically for the investigated events.
Data availability. The "Wind and Ports" database is not publicly accessible. Part of this database concerning selected thunderstorm events will be available soon through a portal linked to the website of THUNDERR project, which is currently under construction. 2016. It is also funded by "Compagnia di San Paolo" for the project "Wind monitoring, simulation, and forecasting for the smart management and safety of port, urban and territorial systems" (grant no. 2015.0333, ID ROL: 9820), by the Italian Ministry of Instruction and Scientific Research (PRIN 2015), with regard to the Project "Identification and diagnostic of complex structural systems" (grant no. 2015TTJN95), and by the 111 project "Innovation on mitigating wind-induced disaster of infrastructures sensitive to wind" supported by the Ministry of Education, China, at the Beijing Jiaotong University. The data exploited for this research have been recorded by the monitoring network realized in the course of the European Projects "Winds and Ports" and "Wind, Ports and Sea", funded by the European Territorial Cooperation Objective, cross-border program Italy-France Maritime, 2007-2013. Satellite images are based on level 1 data recorded by SEVIRI instrument on-board Meteosat Second Generation satellites, operated by EUMETSAT.
Edited by: Maria-Carmen Llasat Reviewed by: two anonymous referees