EXTENDED ABSTRACT Introduction Bora is a strong, temporally and spatially transient downslope wind. Its mean wind velocities rarely exceed 30 m s-1, while one of Bora main characteristics is the gustiness when the wind velocity may reach up to five times the mean wind velocity, up to 69 m s-1 (e.g., Bajić and Peroš, 2005). Bora is more common in winter when it can last for several days, (e.g., Poje, 1992). It is more frequent at the eastern Adriatic and weakens seaward from the shore in a way that it is rarely stormy over the western Adriatic (e.g., Grisogono and Belušić, 2009). Depending on the triggering system, Bora can be cyclonic, bringing clouds with a high possibility for precipitation, anticyclonic, accompanied with clear weather, or frontal (e.g., Jurčec and Visković, 1994). Previous meteorological and geophysical studies laid groundwork regarding Bora large- and mesoscale motions, but further research is still required to fully understand and resolve Bora turbulence in a form usable for engineers. Bora is characteristic for the eastern Adriatic Coast but also for other dynamically similar places around the world thus making the Bora research globally applicable, (e.g., Jurčec, 1981). For example, Bora affects lives of people and may create significant problems to traffic and engineering structures (e.g., Kozmar et al., 2014, 2015). Due to Bora, ferry traffic among the coastal ports and the islands usually stops, major roads and highways remain closed to traffic for several days. Bora can cause vehicles to roll over or slide away from the road, blow off roofs and building façades. It can damage cable-stayed bridges as well. For example, the Dubrovnik cable-stayed bridge was considerably damaged in several severe Bora episodes in 2005 and 2006. Wind turbines cannot operate during strong Bora events because they experience structural loading considerably different than the wind turbines in other areas without Bora. Bora blows down the trees, causing larger areas or even entire counties to lose electricity due to collapsing poles. There are not many records of building collapsing due to the severe Bora wind because buildings in coastal Croatia are generally low-rise and made of concrete that makes them not very prone to wind effects. In the summer, wildfires are significantly enhanced when Bora blows. For these reasons, further research is still needed to fully elucidate Bora characteristics in a form useful to engineers studying wind loads on structures. Therefore, an analysis of vertical wind profiles, time series of turbulence intensity, Reynolds shear stress, turbulence length scales, power spectral density and other relevant parameters, was carried out and the results were compared to international standards. Furthermore, velocity power spectra and thermal stratification of the low-level atmosphere during Bora was analysed. Methodology The measured data considered in this study were obtained from April 2010 until June 2011 at the meteorological tower on the Pometeno brdo (43°36'28,9''N and 16°28'37,4''E) in the hinterland of Split, Croatia on the lee side of the central Dinaric Alps. This rather long measurement period allows for making quite robust and reliable conclusions. Northern, eastern and vertical wind velocity components as well as the ultrasonic temperature were measured using ultrasonic anemometers at the sampling frequency of 5 Hz. That frequency was at the time the finest sampling rate ever applied to study Bora turbulence characteristics. Anemometers were mounted on a 60 m high meteorological tower at 10, 22 and 40 m heights placed along the main Bora direction. The vegetation surrounding the measurement site is mostly shrubs and bushes lower than 3 m and this area is known for frequent and strong Bora occurrences (e.g., Makjanić, 1978, Jurčec, 1981). Prior to the analysis of the measured data, it is necessary to check the data quality and remove the outliers (e.g., Pandžić, 2002). Afterwards, the suspicious data that was removed (rather small amount less than approximately 0,004 % of the total dataset), this missing portion of the data was linearly interpolated, and Bora episodes were accordingly selected. A data set is considered as one Bora episode if the wind direction is from the northeast with negative longitudinal and lateral wind velocity components and if these conditions lasted continuously for a minimum of three hours. During 18 months, the period in which the measurements were carried out, there were in total 296 Bora episodes that lasted from a few hours up to 127 hours. In total, 119 Bora episodes lasting longer than ten hours were selected and analysed. As expected, the wintertime Bora episodes generally last longer and have larger mean and maximum wind velocities than the summertime episodes, in agreement with Enger and Grisogono (1998). Initial calculations and analyses were based on one summer Bora episode that lasted a total of 62 h between 24 and 27 July 2010 (Magjarević et al., 2011). During this period, the wind was relatively strong and long-lasting with mean hourly values of wind velocities above 15m s-1. To check whether the trends obtained for this summer episode are valid for other 118 analysed Bora episodes, the same analyses were performed for each individual Bora occurrence. The studied parameters include thermal stratification, the mean Bora wind velocity, turbulence intensity, Reynolds shear stress, turbulence length scales and velocity power spectra. To investigate the applicability of the power- and logarithmic laws on Bora mean wind velocity profiles, a suitable coordinate system was adopted with the x-axis aligned along the mean wind direction. Therefore, the original data for the two horizontal velocity components were adjusted to this new coordinate system. To test the data sensitivity to the choice of the time averaging period, the time-averaged wind velocity for the summertime Bora episode was calculated using moving average with four different time averaging periods. Periods of 5 min, 8 min, 17 min and 20 min were selected because it appears that, for this particular episode, the quasi-periodic Bora pulsations occur at the time period of about 8 min and the time period of 17 ± 3 min represents a suitable turbulence averaging scale - the scale that separates turbulence from the mean flows at large scales (Magjarević et al., 2011). The analysis for the remaining 118 Bora episodes was accordingly performed using moving average with the time averaging period of 17 min. Thermal stratification is commonly described using the dimensionless stability parameter ζ (e.g., Stull, 1988). Positive ζ implies statically stable, negative ζ implies statically unstable and ζ equal to zero implies statically neutral thermal stratification. In the past, the power law proposed by Hellman (1916) proved to represent well the mean wind velocity profile throughout the entire atmospheric boundary layer, while the logarithmic law (e.g., Thuillier and Lappe, 1964), is considered valid within the surface layer. While it is hypothesized that a typical Bora surface layer is almost always near a neutral stratification, it is reasonable to ignore the thermal stability correction term in the logarithmic law and use the simplified expression in further analysis. The applicability of the power law on the Bora wind velocity profile is analysed by calculating the power-law exponent, α, and subsequently fitting the power law to the mean wind velocity profile that comprises of measured data at three studied heights. In particular, the time-averaged mean wind velocities in the x-direction at three levels are normalized using the time-averaged mean wind velocity in the x-direction at 40 m height. To test the logarithmic law suitability for the mean Bora wind velocity profile, friction velocity, u*, and aerodynamic surface roughness length, z0, were calculated in two different ways, i.e., by data adjustment to the logarithmic law, and by directly applying the logarithmic law to a layer between 10 m and 40 m (Lepri et al., 2014, 2015). Turbulence intensity is defined as the root-mean-square of the fluctuating velocity components normalized with the time-averaged wind velocity (e.g., Simiu and Scanlan, 1996). Reynolds shear stress was calculated by using the fluctuating velocity correlations. The integral length scales of turbulence in the x-direction related to u’, v’, w’ fluctuations were calculated using autocorrelation functions and assuming the validity of the Taylor frozen turbulence hypothesis (e.g., ESDU 74030, 1976). To remove the long-wave meandering from the time history prior to calculating Bora turbulence parameters, the moving average was removed from the velocity time history because this long-wave meandering is not considered to be a part of Bora turbulence. The power spectral density of longitudinal, lateral and vertical Bora velocity fluctuations on all three height levels were calculated using the fast Fourier transform (FFT) (e.g., Stull, 1988). The results obtained for turbulence intensity, Reynolds shear stress, turbulence length scales and power spectrum densities were compared with values recommended in major international wind engineering standards and codes. Results and concluding remarks Due to strong mechanical mixing, Bora is near-neutral thermally stratified. The differences in day to night trends are negligible for all parameters considered in this study (Lepri et al., 2017). Trends in the power-law exponent α, friction velocity, u*, and aerodynamic surface roughness length, z0, were studied dependent on the mean Bora wind velocity. For the first time, it was established that the observed profiles of the mean wind velocity along the dominant Bora direction agree well with the power law and logarithmic law approximations. The logarithmic law and the power law fit the measured data better when performing the analysis using the median rather than using the arithmetic mean wind velocity due to a drifting of the mean aerodynamic surface roughness length during the recording time. This indicates that a more robust median may be a more suitable parameter to determine the aerodynamic surface roughness length than the arithmetic mean value. For the friction velocity, the arithmetic mean proved to be independent of the time-averaging period length, while for the power-law exponent neither the arithmetic mean nor the median are influenced by the time-averaging period. Another interesting feature is a decrease in α and z0 and an increase in u* with increasing Bora mean wind velocity, and vice versa. This indicates an urban-like velocity profile for smaller wind velocities and a rural-like velocity profile for larger wind velocities, which is due to a stronger increase in absolute velocity at each of the heights observed as compared to the respective velocity gradient (difference in mean wind velocity among two different heights), (Lepri et al., 2014). Variations in velocity profiles at the same site during different wind periods are important to be noted because in the engineering community it has been commonly adopted that the aerodynamic characteristics at a particular site remain the same during various wind regimes. For the first time, Reynolds shear stress during Bora is assessed directly from suitable data sets. The observed values for turbulence intensity and Reynolds shear stress remain within the same range during the time record if the mean wind velocities exceed 5 m s-1. When the mean wind velocities are lower than 5 m s-1, the spread of these respective values increases considerably. This means that the wind velocity of ~5 m s-1 may be considered as the critical velocity for Bora. Turbulence intensity, unlike Reynolds shear stress, agrees well with the values recommended in international standards for the respective terrain type. Turbulence length scales increase with increasing mean wind velocity and vice versa. With increasing height from the ground, turbulence intensity and absolute Reynolds shear stress decrease, while turbulence length scales increase, all in agreement with characteristics of the typical atmospheric boundary layer. Longitudinal turbulence length scales generally agree well with the values recommended in international standards for the respective terrain type, while turbulence length scales related to the lateral and vertical velocity fluctuations are much larger than the standard values. The choice of the time-averaging period proved not to influence the observed trends, while the median and arithmetic mean values are nearly the same for all three parameters. Spectral distribution of Bora turbulence kinetic energy differs from ESDU 85020 (1985) in a way that it is distributed across a wider range of frequencies compared to ESDU 85020 (1985). Bora spectrum is flatter, without a distinct peak, and there is a strong energy content at high frequencies. To understand the real importance of the Bora wind characterization from the engineering point of view, further work is still required with respect to Bora turbulence for other terrain types, seasons and elevations from the ground, as well as the extreme value distributions in comparison with cyclonic events.