Characterization of flow recirculation zones at the Perdigão site using multi-lidar measurements

. Because ﬂow recirculation can generate signiﬁcant amounts of turbulence, it can impact the success of wind energy projects. This study uses unique Doppler lidar observations to quantify occurrences of ﬂow recirculation on lee sides of ridges. An extensive dataset of observations of ﬂow over complex terrain is available from the Perdigão 2017 ﬁeld campaign over a period of 3 months. The campaign site was selected because of the unique terrain feature of two nearly parallel ridges with a valley-to-ridge-top height difference of about 200 m and a ridge-to-ridge distance of 1.4 km. Six scanning Doppler lidars probed the ﬂow ﬁeld in sev-eral vertical planes orthogonal to the ridges using range– height indicator scans. With this lidar setup, we achieved vertical scans of the recirculation zone at three positions along two parallel ridges. We construct a method to identify ﬂow recirculation zones in the scans, as well as deﬁne characteristics of these zones. According to our data analysis, ﬂow recirculation, with reverse ﬂow wind speeds greater than 0.5 m s − 1 , occurs over 50% of the time when the wind direction is perpendicular to the direction of the ridges. Atmospheric conditions, such as atmospheric stability and wind speed, affect the occurrence of ﬂow Flow recirculation occurs more frequently during periods with wind speeds above 8 m s − 1 . Recirculation within the valley affects the mean wind and turbulence ﬁelds at turbine heights on the downwind ridge in magnitudes signiﬁcant for wind resource assessment.

Abstract. Because flow recirculation can generate significant amounts of turbulence, it can impact the success of wind energy projects. This study uses unique Doppler lidar observations to quantify occurrences of flow recirculation on lee sides of ridges. An extensive dataset of observations of flow over complex terrain is available from the Perdigão 2017 field campaign over a period of 3 months. The campaign site was selected because of the unique terrain feature of two nearly parallel ridges with a valley-to-ridge-top height difference of about 200 m and a ridge-to-ridge distance of 1.4 km.
Six scanning Doppler lidars probed the flow field in several vertical planes orthogonal to the ridges using rangeheight indicator scans. With this lidar setup, we achieved vertical scans of the recirculation zone at three positions along two parallel ridges. We construct a method to identify flow recirculation zones in the scans, as well as define characteristics of these zones. According to our data analysis, flow recirculation, with reverse flow wind speeds greater than 0.5 m s −1 , occurs over 50 % of the time when the wind direction is perpendicular to the direction of the ridges. Atmospheric conditions, such as atmospheric stability and wind speed, affect the occurrence of flow recirculation. Flow recirculation occurs more frequently during periods with wind speeds above 8 m s −1 . Recirculation within the valley affects the mean wind and turbulence fields at turbine heights on the downwind ridge in magnitudes significant for wind resource assessment.

Introduction
Traditional wind turbine siting relies on wind measurements from a single mast or a small number of masts deployed at the site of interest. These measurements are extrapolated over the entire site using a linearized flow model to provide a general assessment of wind resources (Troen and Petersen, 1989). This method proved suitable for flat terrain sites, gently sloping terrains and offshore where the flow is close to being homogeneous. In complex terrain, however, where wind fields are strongly affected by the topography and roughness changes, the uncertainties of wind resource estimations are increased because models struggle to predict flow under these conditions, despite continuous advancement of flow models. A vast number of flow phenomena occur at sites with complex geometry (Rotach and Zardi, 2007), such as flow acceleration and channeling effects, the formation of lee waves and flow recirculation (Stull, 2012), which are not captured well by frequently applied computer models (e.g., the Wind Atlas Analysis and Application Program; WASP, 2019) in the wind farm planning process. These effects have high spatial and temporal variability and are thus difficult to investigate with sparse measurements. However, these phenomena can induce high loads on wind turbines which will decrease turbine lifetime (Sathe et al., 2013), and therefore such effects must be characterized and understood in order to be able to predict them.
In this study, we investigate flow recirculation occurring in the lee of ridges. Flow recirculation can be identified by flow moving reversely to the prevailing flow direction. In this region, both turbulence and wind shear are increased (Stull, 2012). Recirculation occurs behind abrupt drops in terrain elevation such as cliffs but also in the lee of gentle slopes of a certain steepness (Wood, 1995;Xu and Yi, 2013;Kutter et al., 2017). In addition, small terrain features (Lange et al., 2017) and the presence of forest canopy (Sogachev et al., 2004;Finnigan and Belcher, 2004) can modify the occurrence of recirculation drastically. Furthermore, the stratification of the atmosphere influences the occurrence of recirculation; for instance, recirculation was prevalent during neutral and unstable atmospheric conditions during the observational study of Kutter et al. (2017). The length of flow separation zones is mainly dependent on the hill shape and the surface roughness. Increasing downwind slopes and larger roughness are causing longer separation zones. In Kaimal and Finnigan (1994), a usual extent of the separation zones for naturally shaped hills of 2-3 times the hill height is stated.
Many studies that are describing recirculation or in general flow over complex terrain rely on experimental data measured in the field. Probably the best known field experiment took place at Askervein Hill (Taylor and Teunissen, 1983). The measured data are heavily used to validate flow simulations such as linear models but also more sophisticated code using Reynolds-averaged Navier-Stokes (RANS) methods (Silva Lopes et al., 2007) and large-eddy simulations (LESs) as presented in Street (2009) andBerg et al. (2018). Another more recent field experiment is the Bolund Hill campaign, which took place in 2007-2008. A blind test of different microscale flow models for the Bolund case revealed major uncertainties in predicting the influence of the Bolund escarpment on the flow field. The experiment also demonstrated the strength of lidars, in their case continuous-wave Doppler lidars, to measure the influence of the Bolund escarpment on the flow field (Lange et al., 2016). Here, we investigate measurements from the Perdigão 2017 campaign where, in addition to sonic anemometers on 100 m high masts, pulsed Doppler lidars observed the flow over the two parallel ridges. To surpass the capabilities of point measurements from networks of masts, Doppler lidars have, not only in the Bolund campaign, proved capable of measuring flow phenomena such as wind turbine and wind farm wakes (Käsler et al., 2010;Iungo et al., 2013;Aitken et al., 2014;Menke et al., 2018a) and the effect of roughness changes on wind fields in different heights and nearshore flows . In addition to the advancements in the measurement technology, the orography of the Perdigão site can be approximated by two two-dimensional ridges as the Askervein's orography with a three-dimensional hill of about half the height of the Perdigão ridges.
The present work uses multi-lidar measurements to observe and classify recirculation that occurs in the lee of ridges at Perdigão. The flow field is measured by six scanning lidars at three different transects across the ridges. The scans cover the entire recirculation zone, enabling us to characterize the dimensions of recirculation zones and to investigate their spatial variability along the ridges. Moreover, we relate these recirculation zones to resulting changes in turbulence and wind shear measured by sonic anemometers on measurement masts at potential wind turbine sites downwind of recirculation zones. The measurement campaign, Perdigão 2017 Fernando et al., 2018), is presented in Sect. 2 of this paper. In Sect. 3, we describe how the lidar measurements are used to detect recirculation zones. The occurrences of recirculation at the Perdigão site are analyzed with respect to mean wind speed and atmospheric stability in Sect. 4. In addition, we assess the impact of recirculation at potential turbine sites in this section. Conclusions are given in the last section.

Field campaign overview and measurement data
The experiment is part of a series of measurement campaigns within the New European Wind Atlas with the objective to provide an extensive dataset of high-quality measurements of flow over complex terrain . Starting in January 2017, a network of 50 measurement masts was deployed at the site, and over the next few months, 19 scanning Doppler lidars were added. For this study, we use the data from two 100 m masts equipped with 3-D ultrasonic anemometers (Gill WindMaster Pro) and temperature sensors (NCAR SHT75) located on the tops of the ridges as well as six scanning lidars ( Fig. 1 and, for exact positions, Table 1). Here, we focus on the period of 68 days from 8 April to 14 June during which long-range WindScanners and masts operated in parallel.

Measurement site
The measurement site, located in central Portugal near the Spanish border, was selected because of the unique terrain feature of two nearly parallel ridges with a valley-to-ridgetop elevation difference of about 200 m and a ridge-to-ridge distance of 1.4 km (Fig. 1). The ridges run from the northwest to the southeast and are hereafter referred to as SW and NE ridges, respectively. The terrain roughness is characterized by a patchy vegetation of pine and eucalyptus trees with heights of up to 15 m. A detailed site description can be found in Vasiljević et al. (2017).

WindScanner measurements
We deployed six long-range WindScanners, scanning Doppler lidars (Vasiljevic et al., 2016), on top of the ridges at the Perdigão site. All scanners performed range-height indicator (RHI) scans along three transects which are almost perpendicular to the ridges. The transects are orientated along the SW-NE direction ( Fig. 1; see Table 2 for the exact directions). The scans with a maximal range of 3 km covered the area within the valley. Outside the valley, the scans covered areas to a maximum range of 3 km and up to 1600 m above ground level. Elevation profiles along the transects are shown in Fig. 2. We measured continuously during the up-  Table 1; black + symbol: wind turbine position; white X symbols: position of measurement masts used in this investigation. PT-TM06/ETRS89 coordinate system, height above sea level. (b) Distribution of 30 min averaged wind measurements taken by the 100 m anemometer on the mast located on the northeast ridge over the period from 1 April to 15 June 2017. Measurements are in m s −1 . Table 1. Positions of instruments used in this study. Easting and northing are in the PT-TM06/ETRS89 coordinate system and elevation is above sea level. Elevation of the masts (first two rows) refers to the mast base and elevation of the WindScanners (remaining rows) refers to the scanner head position. LRWS indicates long-range WindScanners. wards movement of the scanners with an averaging elevation angle of 0.75 • over a range of 25 • . Range gates were placed every 15 m, starting at a range of 100 m, extending to 3000 m. With range gates, we refer to the time intervals used for determining the wind speed from the back-scattered light. These time intervals correspond to spatial intervals along the line of sight (LOS) for which the LOS wind speed is evaluated. They translate into a weighting function along the LOS, which in this case has a full-width half maximum of approximately 30 m. The RHI scans were carried out twice per hour for 10 min by all lidars. During one 10 min period, 23 scans were performed; a single scan took about 24 s to complete. We averaged the data over each 10 min period, after filtering the data by the carrier-to-noise ratio (CNR < −24 dB) and removing hard target influences. To ensure that the pointing accuracy of a minimum of 0.05 • is maintained during the entire measurement period, the WindScanner's leveling was checked and the known hard targets were mapped consis-tently through the measurement campaign (Vasiljević et al., 2017). For the analysis, we selected periods within two 40 • wind direction sectors centered at a line perpendicular to the ridge (214-254 and 34-74 • ). The direction is measured by the 100 m sonic anemometer on the mast located on the southwest ridge. In total, 1619 10 min periods of lidar measurements are available for the analysis in these sectors (Fig. 3).

Atmospheric stability
To define atmospheric stability, we calculated the gradient Richardson number at the mast located on the northeast ridge ( Fig. 1). The SW ridge mast could not be used due to missing temperature measurements for a major part of the observational periods. The gradient Richardson number is defined as Table 2. Parameters describing the measurement transects. Peak height gives the highest points at the southwest (SW) and northeast (NE) ridges, respectively. The average terrain elevation is measured along the transects in between the peaks. The maximum height of the ridges relative to the terrain height at the transect before and after the ridges is given by h SW and h NE . A SW and A NE is the half width at half height of the faces inside the valley.   in Stull (2012): where is the mean temperature measured at 100 m, ∂ the potential temperature difference between 10 and 100 m, g = 9.81 m s −2 is the gravitational acceleration, and u and v are the two horizontal components of the mean wind vector. The potential temperature is approximated by ≈ T +(g/C p )·z, where g/C p = 0.0098 K m −1 (Stull, 2012). In order to reduce the effect of terrain-induced flow on the stability estimates, we calculate the difference between the wind speed measured by the 100 m sonic and at ground level, 0 m, where the wind speed is assumed to be zero. A similar approach has been used by Burns et al. (2011). All calculations are based on 30 min averages. In total, we could assess the atmospheric stability for 82 % of the 10 min periods due to missing temperature data from the meteorological mast. Moreover, the Brunt-Väisälä frequency, which is defined as is estimated from the mast measurements for stably stratified conditions. We use again the potential temperature difference between the 10 and 100 m levels at the mast on the northeast ridge.

Detection of recirculation zones
We develop a method to identify recirculation zones in measurements taken by lidars along the measurement transects (Fig. 5). This method is based on measurements by two lidars located on the same transect. By using two lidars, the area covered by the scans within the valley can be maximized. The RHI scans provide LOS wind speed measurements in a polar coordinate system with the lidar at the origin. Consequently, two lidars which are measuring in the same vertical plane will provide data in polar coordinate systems with different origins. Therefore, the LOS wind speeds measured by the lidars are interpolated to a Cartesian coordinate system with the abscissa orientated along the measurement transect and a grid size of 10 m in both the vertical and horizontal directions. In overlap regions, only measurements of the lidar further downwind are used. Additional regions that are sampled by the second lidar upwind are filled in a second step. This procedure prevents discontinuities between the two scans from interfering with estimates within the recirculation zone (Fig. 4a). We do not attempt to combine LOS velocities measured by the lidars to get the true horizontal wind along the measurement transect, since the elevation an-  gles are too small to measure the vertical component precisely. Moreover, the assumption of zero vertical wind speed, as often applied in flat terrain, is not valid in complex terrain.
Having the measurements combined, the zero contour line of the velocity field behind the upwind ridge is determined using linear interpolation. The area bounded by this contour and the terrain beneath describes the dimensions of the recirculation zone (Fig. 4). Recirculation is defined to be presented for reverse flow speeds greater than 0.5 m s −1 . Also, the length of these zones is required to exceed five grid cells or 50 m in the horizontal dimension. To increase the robustness of this method, we require that at least one valid measurement is available below the contour line. A potential disadvantage is that recirculation appearing close to the ground cannot be captured since measurement range gates of the lidars that reach into the ground or vegetation must be discarded due to the hard target returns.

Results and discussion
The method to detect the recirculation zones is applied to all available periods so that we can analyze the occurrences of recirculation under different atmospheric conditions in combination with the mast measurements. Specifically, we as-sess the impact of flow direction, mean flow speed and atmospheric stability. The measurement period from April to mid-June was generally very hot with maximum temperatures above 40 • C. Days were normally cloud free, while on some mornings fog formed inside the valley. Therefore, stably stratified atmospheric conditions at night and unstable conditions during the day are prevailing.

Occurrence of recirculation
On average, recirculation occurs frequently, in 52 % of all 10 min time periods examined. This ratio changes for the individual transects. Along transect 1 (the furthest northwest, with the lowest point in the valley and the shortest distance between ridges), recirculation occurs most frequently 69 % of the time, while transect 2 observes recirculation 51 % of the time. Transect 3 (the furthest southeast, with the highest low point in the valley and the longest peak-to-peak distance) only observes recirculation in 32 % of the available periods. These variations of recirculation occurrence may be related to the transects' elevation profiles within the valley. The average elevation along transects 1 and 2 is the same, whereas the average elevation along transect 3 is 15 m higher (Table 2). Additionally, the elevation profile of transect 3 shows a hill at the center of the valley (Fig. 2). We assume that these fea- Figure 5. Flow diagram of recirculation zone detection process. Recirculation is detected for periods that are flagged with "1"; no recirculation is detected for periods flagged with "0". tures suppress the formation of recirculation along transect 3. Other terrain parameters and land use/land cover characteristics only vary insignificantly among the transects.

Dependence on atmospheric stability
Recirculation is more likely to occur during unstable or neutral atmospheric conditions (Ri ≤ 0) than for stable conditions (Ri > 0) for both SW and NE winds (Fig. 6), as expected for flow over steep geometries such as the Perdigão ridges. These results resonate with the findings of Kutter et al. (2017), who found, in a study of recirculation at a single forested hill, that recirculation is less likely to occur under stable conditions. In fact, they found only one period with recirculation during stable conditions, whereas we have several such cases of recirculation under stable conditions.
Only the bin centered at Ri = −0.2 shows an opposite dominance of recirculation versus no recirculation for the two flow directions. Recirculation occurs in this bin mainly for northeasterly flow (Fig. 6b), whereas the opposite is true for southwesterly flow (Fig. 6a). For stability calculations, only measurements from the mast on the northeast ridge are used since the availability of temperature measurements at the other mast is very low. Thus, turbulent mixing is increased at the mast location for southwesterly flow which could affect the Richardson number and thereby these results.

Dependence on mean wind direction
Recirculation occurs more often for southwesterly wind directions (Fig. 7). A possible explanation for this difference could be that the northeast face of the southwest ridge (the downwind face in southwesterly flow) has steep escarpments close to the ridge top and the average slope is slightly steeper compared to the southwest face of the northeast ridge (the downwind face in northeasterly flow). Both the higher steepness and the escarpments make flow separation more likely. At transect 1, recirculation occurs predominantly, especially for southwesterly flow. The ratio of periods with recirculation occurrence to no occurrence is balanced at transect 2, and at transect 3 recirculation can only be observed 32 % of the time. For all transects, the chance of recirculation is higher for southwesterly wind directions.

Dependence on mean wind speed
More differences between the transects emerge when we group the data by wind speed. Transects 1 and 2 show a very high chance of recirculation to occur for wind speeds above 8 m s −1 (Fig. 7). For wind speeds above 8 m s −1 , fewer cases of recirculation occur for northeasterly wind directions, whereas recirculation is present in almost all cases for southwesterly wind directions. Moreover, the wind speed bin of 6 m s −1 shows interesting results in the inter-transect comparisons: at transect 1, recirculation dominates this bin, whereas a contrary trend can be observed at the other transects.

Characteristics under stratified conditions
The findings above show that recirculation occurs as expected, prevailingly under unstable atmospheric conditions, but occurrences under stably stratified conditions are not absent. In general, three types of behaviors can be expected when flow encounters an obstacle (e.g., a hill or a ridge): (i) complete attachment of the flow, where separation downstream is suppressed by strong stratification; (ii) downstream separation of flow direction directly behind the obstacle; and (iii) post-wave separation (Baines, 1998). The occurrence of these behaviors is sensitive to the downstream slope h/A of the obstacle and N h/u. Steeper downstream slopes support the formation of flow separation behind an obstacle and stronger stratification and higher values (separation occurs when N A/u < π ) of N h/u suppress separation, respectively. The relation of these mechanisms is well summarized in Fig. 5.8 in Baines (1998), and we present our    findings in the same non-dimensional framework in Fig. 9. The observations made in Perdigão, recirculation and nonrecirculation periods, are falling almost entirely in an area where NA/u < π in which separation is expected. However, looking at the mean value for recirculation and nonrecirculation periods, it becomes noticeable that the mean of non-recirculation periods of N h/u is higher consistently at all transects. The higher mean value of non-recirculation periods aligns with the expectations described above to find attached flow or post-wave separation (which is not necessarily detected by our detection method) for higher values of N h/u.

Recirculation zone characteristics
The observed recirculation zones can be described by their shape in terms of extent and height and the magnitude of reverse flow speed. On average, the observed recirculation Table 3. Comparison of wind characteristics during the period of 16:30-16:40 UTC on 5 May 2017 measured by the sonic anemometers at 100 m height on the upwind SW and downwind NE ridges. Turbulence intensity is defined as I = σ U U −1 , where U is the mean wind speed and σ U the standard deviation of U . Turbulent kinetic energy is calculated as e = 1 2 u 2 + v 2 + w 2 , where u , v and w are the fluctuating parts of the wind vector components over a 10 min average as measured by the sonic anemometers. areas extend to 697 m in the horizontal for both wind directions, which is almost exactly half the average ridge-to-ridge distance at the transects. At the southeastern transect (transect 3, the longest transect with the least frequent recirculation), the average extent is 22 % smaller than the mean of all transects (Fig. 4). The average extent at the northwestern transect (transect 1, the shortest transect) is 9 % higher than the average extent. The average depth of observed recirculation zones is 157 m above the valley bottom height and changes insignificantly among the transects, whereas the maximum depth shows larger variation among the transects. We found that recirculation zones at transect 1 reach up to 300 m above the valley bottom, and 215 and 189 m at transects 2 and 3, respectively. Maximum reverse flow speeds of above 4 m s −1 are observed. The median reverse flow speeds at the transects are 2.0, 1.5 and 1.4 m s −1 for transects 1, 2 and 3, respectively. For all transects and both flow directions combined, a trend for higher reverse flow speeds for Richardson numbers close to zero is found (Fig. 10). For all transects and both flow directions combined, the relation of the reverse flow speeds to the Richardson number is presented in Fig. 10. The reverse flow speed measurements are divided by the inflow wind speed measured at the meteorological masts. For unstable atmospheric conditions (Ri < 0), no clear trend for the relation of reverse wind speed and inflow wind speed can be observed. Ratios from less than 0.1 to more than 0.5 can be observed. During conditions with Ri > 0, ratios are generally lower and appear to decrease with increasing Ri. For high wind speeds (greater than 12 m s −1 ), ratios of less than 0.3 are observed and they occur, as expected, for small or slightly positive Richardson numbers.

Impact on potential turbine sites
Ridges and hills are advantageous sites for wind turbines, due to the terrain-induced speed-up (Hunt and Snyder, 1980  However, these stronger winds can come at a cost, due to higher turbulence levels and gustiness which increase loads and decrease the lifetime of turbines. In this section, we assess differences in wind characteristics at the downwind ridge for the entire dataset as well as a specific case study of the period shown in Fig. 4.

Case study of recirculation under neutral stratification
On 5 May, from 16:30 to 16:40 UTC, winds were 10.4 m s −1 from the southwesterly direction (254 • ). A strong and distinct recirculation zone occurred behind the southwest ridge within the valley at the second transect (Fig. 4). The zone has a total length of 807 m, spreading out beyond the valley center. It reaches over the ridge-peak height to 192 m above the valley bottom. Reverse flows inside the recirculation zone have a magnitude of up to 1.8 m s −1 or 17 % of the inflow. For this case with a strong recirculation zone, we observe a reduction in mean wind speed and increased turbulence intensity at the northeast ridge. Measurements at the meteorological masts show the impact at the downwind ridge (Table 3). At the 100 m level, the wind speed is reduced 50 % as compared to the upwind 100 m mast. Similarly, the turbulence intensity increases from 12 % at the upwind ridge to 31 % at the downwind ridge. The wind veer observed from 60 to 100 m equals 14 • at the downwind ridge versus 2 • at the upwind ridge. Moreover, the lidar scan shows that the wake of the upwind ridge reaches to approximately 250 m above the northeast peak height (Fig. 4).

General impact of recirculation at downwind locations
The detailed analysis of the period in the previous section shows significant changes of the wind characteristics at the downwind ridge. To generalize this observation, we analyze all periods available from transect 2 (on which the measurement masts are located), or 455 (260) periods to consider for northeasterly (southwesterly) flow. The measurements of 10 min mean wind speed U (we drop the bar for simplicity) at the downwind ridge are normalized for each period by di-viding with the measurements taken at the upstream ridge: At both ridges, data from the 100 m sonic anemometers are used. For both wind directions, a decrease in mean wind speed and an increase in turbulence intensity when recirculation is present are found at the downwind ridge (Table 4). For southwesterly wind directions which tended to have larger and more vigorous recirculation zones, these changes are more pronounced. Notice that the mean of I for recirculation periods shows a increase of I at the downwind ridge compared to upwind ridge of 7.8 % (1.9 %) for southwesterly (northeasterly) wind directions. During periods without recirculation, I only increases by 3.1 % (0.8 %) for southwesterly (northeasterly) wind directions.

Conclusions
Recirculation zones at two parallel ridges have been analyzed using line-of-sight measurements from six long-range WindScanners and two 100 m measurement masts equipped with sonic anemometers. The data were collected during the Perdigão 2017 measurement campaign in spring and summer 2017 in central Portugal. A method is developed to detect recirculation zones from RHI scans performed at three transects perpendicular to the ridges. Atmospheric stability is characterized using measurements of sonic anemometers and temperature sensors on a 100 m mast on the ridges. For flow perpendicular to the ridges, recirculation occurs frequently, in 52 % of the scans, 55 % for southwesterly wind and 49 % for northeasterly wind. The occurrence of recirculation depends on mean wind speed and atmospheric stability. Recirculation is more likely for wind speeds above 8 m s −1 (measured 100 m above the ridges) and is less likely during stable atmospheric conditions. An intercomparison of the recirculation occurrence per transect revealed significant differences. The northwestern transect shows a 69 % probability of recirculation to occur, compared to a 32 % probability at the southeasternmost transect, perhaps due to the topographic variations between the transects.
Finally, recirculation affects the wind characteristics which are of importance for wind energy generation at the downwind ridge. During recirculation periods, the mean wind speed is lower at the downwind ridge than when recirculation is not present. Further, turbulence intensity is increased at the downwind ridge. These increases in turbulence intensity were not symmetric: a increase of the mean turbulence intensity of 7.8 % (3.1 %) is found at the downwind ridge for southwesterly (northeasterly) flow compared to the upwind ridge. These differences can result from differences in the orography: slopes at the southwest ridge are steeper compared to the northeast ridge and the southwest ridge is higher than the northeast ridge. These findings are of importance for wind energy projects in complex terrain. Developers may reduce loads on turbines, and thereby increase turbine lifetimes, by considering the spatial extent of recirculation zones and positioning the turbines accordingly. Because assessment with multiple scanning lidars, as presented in this study, is yet not feasible for commercial projects, efforts should be made to account for recirculation in the flow models used for wind resource assessment.
Author contributions. RM performed the analysis and wrote the main body of the manuscript. JKL, JM and NV gave scientific advice. RM, JM and NV designed the lidar part of the field campaign. All authors contributed critical feedback to this paper.
Competing interests. The authors declare that they have no conflict of interest.
Special issue statement. This article is part of the special issue "Flow in complex terrain: the Perdigão campaigns (ACP/WES/AMT inter-journal SI)". It is not associated with a conference.
Acknowledgements. We acknowledge the work of everyone involved in the planning and execution of the campaign; in particular, we would like to thank Stephan Voß, Julian Hieronimus (ForWind, University of Oldenburg), Per Hansen and Preben Aagaard (DTU Wind Energy) for their help with the installation of the WindScanners, and Steven Oncley and Kurt Knudson (National Center for Atmospheric Research) for their help with the installation of the meteorological masts. We are also grateful for the contribution of three WindScanners to the campaign by ForWind. Moreover, only the intensive negotiations of José Carlos Matos (INEGI) with local landowners about specific locations for our WindScanners made this research possible. We are grateful to the municipality of Vila Velha de Ródão, landowners who authorized installation of scientific equipment in their properties, the residents of Vale do Cobrão, Foz do Cobrão, Alvaiade and Chão das Servas and local businesses who kindly contributed to the success of the campaign. The space for the operational center was generously provided by Centro Sócio-Cultural e Recreativo de Alvaiade in Vila Velha de Rodão. We thank the Danish Energy Agency for funding through the New European Wind Atlas project (EUDP 14-II). Julie K. Lundquist's effort is supported by the US National Science Foundation under grant AGS-1565498.
Edited by: Jose Laginha Palma Reviewed by: three anonymous referees