Journal cover Journal topic
Atmospheric Chemistry and Physics An interactive open-access journal of the European Geosciences Union
Journal topic
Atmos. Chem. Phys., 19, 2713-2723, 2019
https://doi.org/10.5194/acp-19-2713-2019
Atmos. Chem. Phys., 19, 2713-2723, 2019
https://doi.org/10.5194/acp-19-2713-2019

Research article 01 Mar 2019

Research article | 01 Mar 2019

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

Characterization of flow recirculation zones
Robert Menke1, Nikola Vasiljević1, Jakob Mann1, and Julie K. Lundquist2,3 Robert Menke et al.
• 1Technical University of Denmark – DTU Wind Energy, Frederiksborgvej 399, Building 118, 4000 Roskilde, Denmark
• 2Department of Atmospheric and Oceanic Sciences, University of Colorado, Boulder, CO, USA
• 3National Renewable Energy Laboratory, Golden, CO, USA
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 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 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.

1 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 . 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 , such as flow acceleration and channeling effects, the formation of lee waves and flow recirculation (Stull2012), which are not captured well by frequently applied computer models (e.g., the Wind Atlas Analysis and Application Program; ) 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 , 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 (Stull2012). Recirculation occurs behind abrupt drops in terrain elevation such as cliffs but also in the lee of gentle slopes of a certain steepness . In addition, small terrain features and the presence of forest canopy 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 . 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 , 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 . 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 and large-eddy simulations (LESs) as presented in and . 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 . 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 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 , 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.

2 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.

Figure 1(a) Elevation map of the Perdigão site. Measurement transects are indicated by the solid blue (transect 1), dashed red (transect 2) and dotted black (transect 3) lines. Black disks: WindScanner positions; numbers: WindScanner station numbers as in 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.

## 2.1 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-ridge-top 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 .

Table 1Positions 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.

## 2.2 WindScanner measurements

We deployed six long-range WindScanners, scanning Doppler lidars , 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 upwards 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 $<-\mathrm{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 consistently through the measurement campaign .

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).

Figure 2Elevation profiles of measurement transects. The transect colors correspond to Fig. 1. The vertical scale is exaggerated by a factor of 2. The origin (range of 0) refers to the positions of WindScanners 105, 102 and 106 for transects 1, 2 and 3, respectively. The coordinates of the transects are available as CSV files in the Supplement.

Figure 3Distribution of the 1613 available 10 min periods with wind directions perpendicular to the ridges (214–254 and 34–74). The observational period lasted for 68 days starting on 8 April and ending on 14 June 2017. The intensive observation period of the campaign started on day 24 (1 May 2017). The availability data can be found in the Supplement.

Table 2Parameters 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 hSW and hNE. ASW and ANE is the half width at half height of the faces inside the valley.

3 Methods

## 3.1 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 in :

$\begin{array}{}\text{(1)}& {\mathit{Ri}}_{\mathrm{G}}=\frac{\frac{g}{\mathrm{\Theta }}\frac{\partial \stackrel{\mathrm{‾}}{\mathrm{\Theta }}}{\partial z}}{\left[{\left(\frac{\partial \stackrel{\mathrm{‾}}{u}}{\partial z}\right)}^{\mathrm{2}}+{\left(\frac{\partial \stackrel{\mathrm{‾}}{v}}{\partial z}\right)}^{\mathrm{2}}\right]},\end{array}$

where Θ is the mean temperature measured at 100 m, $\partial \stackrel{\mathrm{‾}}{\mathrm{\Theta }}$ 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 $\stackrel{\mathrm{‾}}{\mathrm{\Theta }}\approx T+\left(g/{C}_{\mathrm{p}}\right)\cdot z$, where gCp = 0.0098 K m−1 (Stull2012). 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 . 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

$\begin{array}{}\text{(2)}& N=\sqrt{\frac{g}{\mathrm{\Theta }}\frac{\partial \stackrel{\mathrm{‾}}{\mathrm{\Theta }}}{\partial \stackrel{\mathrm{‾}}{z}}},\end{array}$

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.

Figure 4Detection of recirculation zones at the second transect (lateral distance of 0 is the position of WindScanner 102), view from the southeast. (a) Schematic of the overlay of the two RHI scans at the transect. Blue and red areas are corresponding to the RHI scans and the grey area shows the terrain. (b) Lidar measurements along the transect and detected recirculation zone for one 10 min period with RiG=0.016. Positive LOS velocities indicate flow towards the northeast (from left to right in the plot). The terrain is plotted in grey and the contour line defining the recirculation zone is plotted in black inside the velocity field. Time is in UTC.

## 3.2 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 angles 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.

Figure 5Flow 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”.

4 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 assess 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.

Figure 6Recirculation occurrence as a function of the Richardson number. The results are presented in bins with a width of 0.2. Data from all three transects are presented together, segregated by the wind direction.

## 4.1 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 features suppress the formation of recirculation along transect 3. Other terrain parameters and land use/land cover characteristics only vary insignificantly among the transects.

Figure 7Occurrence of recirculation depending on the mean wind speed measured by the 100 m sonic on the upstream mast. The histograms are separated by transect and wind direction.

### 4.1.1 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 , 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 $\mathit{Ri}=-\mathrm{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.

### 4.1.2 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.

Figure 8Zero contour line (black line) defining the upper border of recirculation zones of the mean velocity field of all recirculation periods per transect and flow direction. The flow speed is positive in the direction of the mean flow field and the transect origins (lateral distance of 0) refer to the positions of WindScanners 105, 102 and 106 for transects 1, 2 and 3, respectively.

### 4.1.3 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.

Figure 9Observational periods with recirculation (colored markers) and without recirculation (grey cross marks) under stratified conditions as a function of mean lee-side slope (hA) and Nhu. The round markers refer to recirculation at the southwest ridge and the square markers to northeast ridge. The colored diamond markers show the mean value of Nhu per transect during periods with recirculation and the grey markers the mean for periods without recirculation. The solid black line shows where $NA/u=\mathit{\pi }$.

### 4.1.4 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 (Baines1998). The occurrence of these behaviors is sensitive to the downstream slope hA of the obstacle and Nhu. Steeper downstream slopes support the formation of flow separation behind an obstacle and stronger stratification and higher values (separation occurs when $NA/u<\mathit{\pi }$) of Nhu suppress separation, respectively. The relation of these mechanisms is well summarized in Fig. 5.8 in , and we present our findings in the same non-dimensional framework in Fig. 9. The observations made in Perdigão, recirculation and non-recirculation periods, are falling almost entirely in an area where $NA/u<\mathit{\pi }$ in which separation is expected. However, looking at the mean value for recirculation and non-recirculation periods, it becomes noticeable that the mean of non-recirculation periods of Nhu 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 Nhu.

Figure 10Maximum reverse flow speed divided by the inflow wind speed over the Richardson number. The color scaling shows the wind speed measured at the upstream mast 100 m above ground.

Table 3Comparison 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={\mathit{\sigma }}_{U}{\stackrel{\mathrm{‾}}{U}}^{-\mathrm{1}}$, where $\stackrel{\mathrm{‾}}{U}$ is the mean wind speed and σU the standard deviation of U. Turbulent kinetic energy is calculated as $e=\frac{\mathrm{1}}{\mathrm{2}}\left[\stackrel{\mathrm{‾}}{{{u}^{\prime }}^{\mathrm{2}}}+\stackrel{\mathrm{‾}}{{{v}^{\prime }}^{\mathrm{2}}}+\stackrel{\mathrm{‾}}{{{w}^{\prime }}^{\mathrm{2}}}\right]$, 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.

Table 4Changes in the wind characteristics at the downwind ridge.

## 4.2 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 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.

## 4.3 Impact on potential turbine sites

Ridges and hills are advantageous sites for wind turbines, due to the terrain-induced speed-up . 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.

### 4.3.1 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).

### 4.3.2 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 dividing with the measurements taken at the upstream ridge:

$\begin{array}{}\text{(3)}& \mathrm{\Delta }\left(U\right)=\frac{\left({U}_{\mathrm{downwind}}-{U}_{\mathrm{upwind}}\right)}{{U}_{\mathrm{upwind}}}\cdot \mathrm{100}.\end{array}$

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.

5 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.

Data availability
Data availability.

All datasets measured during the Perdigão 2017 campaign are publicly available through dedicated web portals of the University of Porto (https://windsp.fe.up.pt/, ) and NCAR/EOL (https://data.eol.ucar.edu/master_list/?project=PERDIGAO, ). The lidar dataset measured by DTU is additionally available under the following DOI: https://doi.org/10.11583/DTU.7228544.v1 ().

Supplement
Supplement.

Author contributions
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
Competing interests.

The authors declare that they have no conflict of interest.

Special issue statement
Special issue statement.

Acknowledgements
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

References

Aitken, M. L., Banta, R. M., Pichugina, Y. L., and Lundquist, J. K.: Quantifying Wind Turbine Wake Characteristics from Scanning Remote Sensor Data, J. Atmos. Ocean. Tech., 31, 765–787, https://doi.org/10.1175/JTECH-D-13-00104.1, 2014. a

Baines, P.: Topographic Effects in Stratified Flows, Cambridge Monographs on Mechanics and Applied Mathematics, Cambridge University Press, 1998. a, b

Berg, J., Troldborg, N., Menke, R., Patton, E. G., Sullivan, P. P., Mann, J., and Sørensen, N. N.: Flow in complex terrain – a Large Eddy Simulation comparison study, Journal of Physics: Conference Series, 1037, p072015, 2018. a

Burns, S. P., Sun, J., Lenschow, D. H., Oncley, S. P., Stephens, B. B., Yi, C., Anderson, D. E., Hu, J., and Monson, R. K.: Atmospheric Stability Effects on Wind Fields and Scalar Mixing Within and Just Above a Subalpine Forest in Sloping Terrain, Bound.-Lay. Meteorol., 138, 231–262, https://doi.org/10.1007/s10546-010-9560-6, 2011. a

Chow, F. K. and Street, R. L.: Evaluation of Turbulence Closure Models for Large-Eddy Simulation over Complex Terrain: Flow over Askervein Hill, J. Appl. Meteorol. Clim., 48, 1050–1065, https://doi.org/10.1175/2008JAMC1862.1, 2009. a

Earth Observing Laboratory: Perdigão data archive, available at: http://data.eol.ucar.edu/master_list/?project=PERDIGAO, last access: 26 February 2019. a

Fernando, H., Mann, J., Palma, J., Lundquist, J., Barthelmie, R., BeloPereira, M., Brown, W., Chow, F., Gerz, T., Hocut, C., Klein, P., Leo, L., Matos, J., Oncley, S., Pryor, S., Bariteau, L., Bell, T., Bodini, N., Carney, M., Courtney, M., Creegan, E., Dimitrova, R., Gomes, S., Hagen, M., Hyde, J., Kigle, S., Krishnamurthy, R., Lopes, J., Mazzaro, L., Neher, J., Menke, R., Murphy, P., Oswald, L., Otarola-Bustos, S., Pattantyus, A., Rodrigues, C. V., Schady, A., Sirin, N., Spuler, S., Svensson, E., Tomaszewski, J., Turner, D., van Veen, L., Vasiljević, N., Vassallo, D., Voss, S., Wildmann, N., and Wang, Y.: The Perdigão: Peering into Microscale Details of Mountain Winds, B. Am. Meteorol. Soc., online first, https://doi.org/10.1175/BAMS-D-17-0227.1, 2018. a

Finnigan, J. and Belcher, S.: Flow over a hill covered with a plant canopy, Q. J. Roy. Meteor. Soc., 130, 1–29, 2004. a

Hunt, J. C. R. and Snyder, W. H.: Experiments on stably and neutrally stratified flow over a model three-dimensional hill, J. Fluid Mech., 96, 671–704, https://doi.org/10.1017/S0022112080002303, 1980. a

Iungo, G. V., Wu, Y.-T., and Porté-Agel, F.: Field Measurements of Wind Turbine Wakes with Lidars, J. Atmos. Ocean. Tech., 30, 274–287, https://doi.org/10.1175/JTECH-D-12-00051.1, 2013. a

Kaimal, J. C. and Finnigan, J. J.: Atmospheric boundary layer flows: their structure and measurement, Oxford university press, 1994. a

Käsler, Y., Rahm, S., Simmet, R., and Kühn, M.: Wake Measurements of a Multi-MW Wind Turbine with Coherent Long-Range Pulsed Doppler Wind Lidar, J. Atmos. Ocean. Tech., 27, 1529–1532, https://doi.org/10.1175/2010JTECHA1483.1, 2010. a

Kutter, E., Yi, C., Hendrey, G., Liu, H., Eaton, T., and Ni-Meister, W.: Recirculation over complex terrain: Recirculation Over Complex Terrain, J. Geophys. Res.-Atmos., 122, 6637–6651, https://doi.org/10.1002/2016JD026409, 2017. a, b, c

Lange, J., Mann, J., Angelou, N., Berg, J., Sjöholm, M., and Mikkelsen, T.: Variations of the Wake Height over the Bolund Escarpment Measured by a Scanning Lidar, Bound.-Lay. Meteorol., 159, 147–159, https://doi.org/10.1007/s10546-015-0107-8, 2016. a

Lange, J., Mann, J., Berg, J., Parvu, D., Kilpatrick, R., Costache, A., Chowdhury, J., Siddiqui, K., and Hangan, H.: For wind turbines in complex terrain, the devil is in the detail, Environ. Res. Lett., 12, 094020, https://doi.org/10.1088/1748-9326/aa81db, 2017. a

Mann, J., Angelou, N., Arnqvist, J., Callies, D., Cantero, E., Arroyo, R. C., Courtney, M., Cuxart, J., Dellwik, E., Gottschall, J., Ivanell, S., Kÿhn, P., Lea, G., Matos, J. C., Palma, J. M. L. M., Pauscher, L., Peña, A., Rodrigo, J. S., Söderberg, S., Vasiljevic, N., and Rodrigues, C. V.: Complex terrain experiments in the New European Wind Atlas, Philos. T. Roy. Soc. A, 375, 2016101, https://doi.org/10.1098/rsta.2016.0101, 2017. a, b, c

Menke, R., Vasiljevic, N., Hansen, K., Hahmann, A. N., and Mann, J.: Does the wind turbine wake follow the topography? – A multi-lidar study in complex terrain, Wind Energ. Sci., 3, 681–691, https://doi.org/10.5194/wes-3-681-2018, 2018a. a

Menke, R., Mann, J., and Vasiljevic, N.: Perdigão-2017: multi-lidar flow mapping over the complex terrain site, https://doi.org/10.11583/DTU.7228544.v1, 2018b. a

Porto University: WindsP data archive, available at: https://windsp.fe.up.pt/, last access: 26 February 2019.  a

Rotach, M. W. and Zardi, D.: On the boundary-layer structure over highly complex terrain: Key findings from MAP, Q. J. Roy. Meteor. Soc., 133, 937–948, https://doi.org/10.1002/qj.71, 2007. a

Sathe, A., Mann, J., Barlas, T., Bierbooms, W., and van Bussel, G.: Influence of atmospheric stability on wind turbine loads: Atmospheric stability and loads, Wind Energy, 16, 1013–1032, https://doi.org/10.1002/we.1528, 2013. a

Silva Lopes, A., Palma, J. M. L. M., and Castro, F. A.: Simulation of the Askervein flow. Part 2: Large-eddy simulations, Bound.-Lay. Meteorol., 125, 85–108, https://doi.org/10.1007/s10546-007-9195-4, 2007. a

Sogachev, A., Rannik, Ü., and Vesala, T.: Flux footprints over complex terrain covered by heterogeneous forest, Agr. Forest Meteorol., 127, 143–158, https://doi.org/10.1016/j.agrformet.2004.07.010, 2004. a

Stull, R. B.: An Introduction to Boundary Layer Meteorology, Springer Science & Business Media, 2012. a, b, c, d

Taylor, P. A. and Teunissen, H. W.: Askervein'82: Report on the September/October 1982 Experiment to Study Boundary Layer Flow over Askervein, South Uist, Meteorological Services Research Branch, Atmospheric Environment Service, 1983. a

Troen, I. and Petersen, E. L.: European wind atlas, Published for the Commission of the European Communities, Directorate-General for Science, Research, and Development, Brussels, Belgium by Risø National Laboratory, Roskilde, Denmark, 1989. a

Vasiljevic, N., Lea, G., Courtney, M., Cariou, J.-P., Mann, J., and Mikkelsen, T.: Long-range WindScanner system, Remote Sensing, 8, 896, https://doi.org/10.3390/rs8110896, 2016. a

Vasiljević, N., L. M. Palma, J. M., Angelou, N., Carlos Matos, J., Menke, R., Lea, G., Mann, J., Courtney, M., Frölen Ribeiro, L., and M. G. C. Gomes, V. M.: Perdigão 2015: methodology for atmospheric multi-Doppler lidar experiment Atmos. Meas. Tech., 10, 3463–3483, https://doi.org/10.5194/amt-10-3463-2017, 2017. a, b

WASP: Wind Atlas Analysis and Application Program, available at: http://www.wasp.dk, last access: January 31 2019. a

Wood, N.: The onset of separation in neutral, turbulent flow over hills, Bound.-Lay. Meteorol., 76, 137–164, https://doi.org/10.1007/BF00710894, 1995. a

Xu, X. and Yi, C.: The influence of geometry on recirculation and CO2 transport over forested hills, Meteorol. Atmos. Phys., 119, 187–196, https://doi.org/10.1007/s00703-012-0224-6, 2013. a