Research article 28 May 2021
Research article  28 May 2021
Wind inflow observation from load harmonics: initial steps towards a field validation
 Wind Energy Institute, Technische Universität München, 85748 Garching bei München, Germany
 Wind Energy Institute, Technische Universität München, 85748 Garching bei München, Germany
Correspondence: Carlo L. Bottasso (carlo.bottasso@tum.de)
Hide author detailsCorrespondence: Carlo L. Bottasso (carlo.bottasso@tum.de)
A previously published wind sensing method is applied to an experimental dataset obtained from a 3.5 MW turbine. The method is based on a loadwind model that correlates onceperrevolution blade load harmonics to rotorequivalent shears and wind directions. Loads measured during turbine operation are used to estimate online – through the loadwind model – the inflow at the rotor disk, thereby turning the whole turbine into a sort of generalized anemometer.
The experimental dataset consists of synchronous measurements of loads, from blademounted strain gages, and of the inflow, obtained from a nearby met mast. As the mast reaches only to hub height, a second independent method is used to extend the metmastmeasured shear above hub height to cover the entire rotor disk. Part of the dataset is first used to identify the loadwind model, and then the performance of the wind observer is characterized with the rest of the data.
Although the experimental setup falls short of providing a thorough validation of the method, it still allows for a realistic practical demonstration of some of its main features. Results indicate a good quality of the estimated linear shear both in terms of 1 and 10 min averages and of resolved time histories, with mean average errors around 0.04. A similarly accuracy is found in the estimation of the yaw misalignment, with mean errors typically below 3^{∘}.
This paper presents a first attempt at the field validation of a wind sensing method based on load harmonics. Wind sensing refers to the general concept of using the response of the turbine to estimate certain characteristics of the inflow, a task that can be accomplished in several different ways (Bottasso et al., 2010; Bottasso and Riboldi, 2014; Simley and Pao, 2016; Bottasso and Riboldi, 2015; Bertelè et al., 2017; Bottasso et al., 2018; Schreiber et al., 2020). Information on the inflow can support a variety of applications, including turbine and farmlevel control, lifetime assessment and fatigue consumption estimation, power and wind forecasting, and others (Schreiber et al., 2020).
A detailed knowledge of the wind inflow is today lacking in essentially all wind turbine installations, with the exception of experimental prototypes, certification tests, and turbines equipped with forwarding looking lidars (Scholbrock et al., 2013; Schlipf et al., 2014; Peña Diaz et al., 2014). In fact, standard production machines are typically equipped with nacellemounted anemometers and wind vanes. These sensors need to be carefully calibrated to eliminate effects caused by – among others – the rotor wake, blade passage, and flow distortions caused by the large bluff body represented by the nacelle. Even when these effects are properly accounted for, nacelle mounted sensors suffer from the unavoidable limitation of providing only pointwise measurements. As such, they are blind to flow features that are characterized by a variability in the wind field over the rotor swept area, namely horizontal and vertical shears, veer, or the presence of an impinging wake released by a turbine located upstream. Additionally, especially for today's very large turbines of ever increasing diameters, pointwise measurements might not fully reflect the actual conditions experienced at the rotor disk. At some wind plants, met masts are available and can in principle provide additional information on the wind characteristics at different heights above ground. However, here again this information is of limited use: only a small number of masts are typically available at production sites, and, clearly, these measurements are not colocated with the turbines. Lidars and radars (Lang and McKeogh, 2011; Scholbrock et al., 2013; Mikkelsen, 2014; Schlipf et al., 2014; Hirth et al., 2015; Valldecabres et al., 2018; Peña Diaz et al., 2014) are remote sensing devices that can be used to scan the flow, providing maps of wind characteristics in space and time. Such devices are, however, not yet routinely used on production machines because of cost, reliability, and availability issues.
Wind sensing was first proposed to address the need for a simple, lowcost way of measuring the inflow at the rotor disk during operation of a turbine, a capability that is today still lacking.
The wind sensing approach is based on two main observations.
The first observation is mainly an economic one and relates to the opportunity offered by sensor technology, in particular strain gages, accelerometers, and pressure sensors. While sensors have been and still are routinely used on prototypes and during certification and research tests, they are today becoming more commonly deployed also on production machines, for example for enabling loadreducing control, or for condition monitoring, digital twin applications, faulticingerosion detection, etc. Indeed, a growing number of use cases, improved technology, and reduced purchase and maintenance costs have led several original equipment manufacturers (OEMs) to equip their latest models with rotor sensors, while retrofitting solutions are becoming readily available on the market (Bachmann Electronic GmbH, 2021; fos4X GmbH, 2021; WindConsult GmbH, 2021). The question is then, can these same devices also be used for wind sensing? In principle, a positive answer to this question opens up very interesting opportunities: when a rotor is already sensorized, the estimation of the wind inflow could be a simple software upgrade, with no extra equipment needed, and therefore no extra purchase and maintenance costs. One could then, at no or very limited cost, turn each rotor in a farm into a sophisticated anemometer, this way providing a wealth of information on the actual flow conditions throughout the plant. But even if a rotor is not already equipped with sensors, a wind sensing technology based on existing proven sensors readily available on the market might still be very attractive and economically viable. Based on this first observation, the development of wind sensing has mainly revolved up to now on the idea of using information provided by blade load sensors (either standard bridgebased ones or optical ones based on fiber Bragg grating technology; Schubel et al., 2013), which is also the approach pursued in this paper. However, future work could try to exploit accelerometers, alone or fused together with load sensors.
The second observation at the heart of wind sensing is rather obviously that changes in the wind inflow will affect the rotor response. The basic scientific question is then, if one could measure the rotor response (through the rotor sensors discussed above), would it be possible to infer the wind inflow from such measurements? In other words, is there a windresponse relationship that can be inverted to estimate the inflow given a measurement of the response? A positive answer to this question was first given by Bottasso and Riboldi (2014), which showed that shear and misalignment (i.e., a relative angle between rotor axis and wind vector) do leave distinguishing effects on the onceperrevolution (1P) sine and cosine harmonic components of the blade response. The method, termed here harmonicbased approach, was then further elaborated and improved by Bottasso and Riboldi (2015), Cacciola et al. (2016a), and Bertelè et al. (2017, 2018, 2019). With time, a simpler method was developed (Bottasso et al., 2018; Schreiber et al., 2016, 2020) which uses blade load time histories to estimate the average wind over sectors of the rotor disk. This second method, termed here sectoreffective approach to distinguish it from the harmonicbased one, is capable of detecting shears and an impinging wake (Schreiber et al., 2020) but not wind directions.
The present paper uses an existing dataset, previously collected for other purposes, to attempt a first field validation of the harmonicbased wind inflow estimator. In short, this method is based on a correlation between the in and outofplane blade root bending 1P harmonics and four specific characteristics of the rotor inflow: namely, the horizontal and vertical shears and the lateral and vertical misalignment angles. Indeed, it can be shown (Bertelè et al., 2017) that each one of these inflow characteristics generates a specific response in the rotor at the 1P frequency. This is a very desirable feature because of the following reasons:

The 1P frequency is strongly dominated by these four “deterministic” characteristics of the wind and much less so by turbulent fluctuations (Bertelè et al., 2017). On the other hand, higher response frequencies are associated with smallerscale variations in the flow in space and time caused by turbulence (Bertelè et al., 2017).

The measurement of such low frequencies requires low sampling rates, which eases the requirements on the sensors.

There should be limited variability in such low frequencies among different installations of the same wind turbine type. Although this has not been demonstrated yet, it would imply that, once the method has been tuned on one machine, it should be applicable with minimum recalibration also on different turbines of the same type.

The lower frequencies of the response of a wind turbine should be reasonably well captured by existing simulation tools used for design and certification. This implies that the method can be tested in a simulation environment with the expectation of realistic results on its performance. This is clearly important for a number of reasons, not least the fact that in a simulation environment the actual inflow at the rotor disk is precisely known, something that is much harder to do in the field. Indeed, as shown later on, the incomplete knowledge of the actual inflow is one of the main limitations of the present experimental study, and of any field test.
If the turbine implements individual pitch control, the map correlating loads and wind characteristics is extended to include also blade pitch angle harmonics (Bertelè et al., 2018).
The loadharmonic method requires a training dataset consisting of measured rotor loads and corresponding measured wind characteristics. The dataset can be based on experimental measurements or be generated synthetically using a simulation model; these two approaches were respectively termed modelfree and modelbased in Bottasso and Riboldi (2014). Here we consider the former approach. Indeed, a simulation model with the necessary characteristics might not always be available, for example in cases when wind sensing is applied to a turbine without the support of the manufacturer. Even when a model is available, it might not have been fully validated, so a purely datadriven approach has a significant appeal.
It is envisioned that, in a practical application of the modelfree harmonicbased method, the training phase would be a oneoff activity performed at a test site equipped with a met mast or other wind measuring devices such as lidars (Mikkelsen, 2014; Peña Diaz et al., 2014). Indeed, hubtall met masts are routinely used during certification (IEC, 2017) and could be employed for the additional purpose of training the observer. After training, the method could be used on other installations of that same turbine type at normal production sites without necessitating met masts or other devices. This claim, however, still remains to be demonstrated in practice.
The principal goal of this paper is to present the application of the loadharmonic estimator to field test data collected at a test site on a 3.5 MW wind turbine and a nearby met mast (Schreiber et al., 2020; Bertelè and Bottasso, 2020). This experimental setup is a realistic representation of the scenario outlined above, in which a hubtall met mast is located in close proximity to a wind turbine for certification purposes. From this point of view, the present dataset provides opportunities not only for a first – partial – field demonstration of the method but also for addressing some important practical implementational aspects.
Specifically, the vertical shear requires special attention. In fact, a hubtall met mast with more than one anemometer can only measure the wind shear over the lower part of the rotor disk; on the other hand, the loadharmonic observer estimates a rotorequivalent shear (i.e., a shear over the entire rotor disk area). For large modern rotors, halfrotor or fullrotor shears are not necessarily equal (Murphy et al., 2020; Schreiber et al., 2020). Therefore, a way is needed to extend the measurement of the inflow above the met mast, possibly without resorting to extra windscanning equipment to reduce cost and complexity. This problem is solved here using the sectoreffective wind sensing method (Bottasso et al., 2018; Schreiber et al., 2016, 2020). This second approach uses blade loads to estimate the average local speed over sectors of the rotor disk; from these sectorequivalent wind speeds, one can then estimate shears, including a vertical shear defined over just the lower half of the rotor.
The sectoreffective speed and loadharmonic observers have distinct characteristics which make them somewhat complementary and applicable to different scenarios. In fact, the sectoreffective observer does not need to be trained with data before it can be used since it is derived from standard performance characteristics of the rotor (Schreiber et al., 2020). Although not indispensable, field data can optionally be employed to finetune the observer, as shown in Schreiber et al. (2020). The sectoreffective approach, however, can only reconstruct shears and not wind directions. The loadharmonic observer, on the other hand, can reconstruct both shears and directions but needs to be trained from data, which is a potential complication. Here, a novel threestep procedure is developed and demonstrated, in which the two observers are used in synergy combining some of their complementary features:

The lowerhalfrotor shear measured by the sectorequivalent speed method is tuned and validated with respect to the metmast reference.

The fullrotor shear is computed using the validated sectorequivalent speed method, extending the measurement of the inflow above the met mast.

This rotorequivalent shear is finally used for training the harmonicbased estimator.
Although the present setup allows for a first demonstration of this procedure, it also presents some limitations that hinder a real and complete validation of the method. First, the extension of the shear above the met mast is performed through the same rotor loads that are also used by the harmonicbased estimator. Clearly, a completely independent measurement of the inflow up to the tip of the rotor would be preferable for validation purposes. Second, the present met mast only includes a wind vane at hub height. This is a pointwise measurement, whereas the one provided by the observer – being obtained through the response of the rotor – is a rotoreffective quantity. Here again, it would be desirable to train and verify the method with an independently derived rotorequivalent quantity. Third, a met mast cannot really provide a true and absolute ground truth as it measures the flow away from the rotor disk (two and a half diameters away, in the present case). When the wind is not directly aligned with turbine and mast, the wind shear and direction may be slightly different on account of wind spatial variability because of orographic and vegetationinduced effects. These differences are indeed visible to some extent in the present dataset. Even when wind, mast, and turbine are aligned, the two measurements are not colocated and therefore not necessarily identical. Fourth, the met mast does not provide measurements for two of the four observed quantities, namely horizontal shear and upflow, for which, consequently, no comparison nor conclusion can be made. Clearly, a more precise characterization of the effective inflow experienced by the rotor disk would be desirable for validation purposes. A lidar scanning the inflow immediately in front of the disk plane – to ensure colocation of the measurements – might be a possible solution.
Although the present study clearly falls short of a true validation of the harmonicbased formulation of wind sensing, it still provides for a first field demonstration of this method, giving also a useful insight into some of its main characteristics.
The paper is organized as follows. Section 2 describes the overall methodology, including a brief review of the harmonicbased estimator in Sect. 2.2 and a description of the test site and the measurement of the inflow characteristics in Sect. 2.3. The analysis of the wind observer performance is presented in Sect. 3, while Sect. 4 concludes the paper.
2.1 Wind parametrization
The wind inflow is described by four parameters: the vertical linear shear κ_{v}, the horizontal linear shear κ_{h}, the vertical wind misalignment angle (or upflow) χ, and the horizontal (or yaw) misalignment angle ϕ. These quantities are illustrated in Fig. 1 and are expressed in a hubcentered nacelleattached reference frame, where x is parallel to the axis of rotation (and it is therefore inclined with respect to the ground because of uptilt), y is horizontal with respect to the ground and points left looking downstream, and z forms a righthanded triad. It should be noticed that the vertical shear is customarily defined with respect to the horizontal instead of the uptilt direction; additionally, its profile is typically either logarithmic or expressed as a power law instead of a linear function. As explained later, these choices are made here to exploit the rotational symmetry of the rotor (Bertelè et al., 2019); this is useful in the present context because it allows the lack of horizontal shear and upflow measurements in the available dataset to be overcome. Clearly, the four wind parameters, once estimated, can be readily transformed into a horizontal frame, if necessary. Furthermore, abandoning the rotational symmetry, the observer can be formulated in terms of a vertical nonlinear shear, as shown in Bertelè et al. (2017).
A linearly sheared wind speed W at the rotor disk is defined as
where V_{h} is the hubheight speed and R the rotor radius. By projecting the wind vector along the x, y, and z axes, respectively, the three nacelleattached velocity components u, v, and w are readily obtained as
where $\stackrel{\mathrm{\u0303}}{v}$ and $\stackrel{\mathrm{\u0303}}{w}$ are defined as
For notational simplicity, the four wind parameters are grouped together in the wind state vector $\mathit{\theta}=\mathit{\{}\stackrel{\mathrm{\u0303}}{v},{\mathit{\kappa}}_{\mathrm{v}},\stackrel{\mathrm{\u0303}}{w},{\mathit{\kappa}}_{\mathrm{h}}{\mathit{\}}}^{\mathrm{T}}$. Given θ, the misalignment angles can be readily computed by inverting Eq. (3) to get $\mathit{\chi}=\mathrm{arcsin}\stackrel{\mathrm{\u0303}}{w}$ and $\mathit{\varphi}=\mathrm{arcsin}\stackrel{\mathrm{\u0303}}{v}/\mathrm{cos}\mathit{\chi}$.
2.2 Wind observer formulation
The relationship between wind states and rotor loads is assumed in the linear form
where F and m_{0} are model coefficients. These coefficients depend on wind speed V and air density ρ on account of the different behavior, control, and deformation of the machine in different operating conditions. For example, under the push of the rotor thrust, the tower will bend backward, in turn slightly changing the rotor uptilt; if this effect is not accounted for, this deformationinduced uptilt will affect the observed wind upflow. The dependency on wind speed is taken into account by discretizing the wind speed range in nodal values and linearly interpolating the model based on the rotorequivalent wind speed (see Sect. 2.3.4), while density is accounted for as explained in Sect. 2.2.1. The load vector m is defined as
where m indicates the blade bending moment, subscripts (⋅)_{1s} and (⋅)_{1c}, respectively, indicate 1P sine and cosine harmonic amplitudes, and superscripts (⋅)^{OP} and (⋅)^{IP} indicate out and inplane load components, respectively. Harmonic components are obtained from measured blade loads using the Coleman transformation (Coleman and Feingold, 1958), followed by lowpass filtering.
The model coefficients F are not all independent because of the rotational symmetry of the rotor (Bertelè et al., 2019). In fact, neglecting the disturbance caused by the tower, the effects on loads caused by a horizontal shear are the same as the ones caused by an equal vertical shear with a phase delay of $\mathit{\pi}/\mathrm{2}$; the same holds true for the wind misalignment angles. This not only reduces the number of unknowns but also eases the identification of the model, especially when using longer time averages. In fact, while both vertical and horizontal shear undergo rapid changes due to spatial turbulence variability, it is easier to observe slower changes in vertical shear than in the horizontal one. In fact, vertical shear exhibits slow natural changes over a significant range, for example because of diurnal fluctuations. On the other hand, horizontal shear might exhibit slowscale variability because of orographic effects or in waked conditions, which – depending on the turbine – might or might not happen very frequently or be particularly pronounced. Similarly, whereas yaw misalignment changes significantly in normal operation because of the inability of the yaw system to immediately and exactly track rapid wind direction fluctuations, upflow changes little (except that for orographic winddirectiondependent effects). Therefore, by exploiting the rotational symmetry, a complete model can be identified simply from variable vertical shear and horizontal misalignment because the effects of the other two wind states are obtained by the symmetry of the coefficients.
Dropping the dependency on V and ρ for notational simplicity, the model coefficients T are identified by stacking side by side measured wind states θ into a matrix $\mathbf{\Theta}=[{\stackrel{\mathrm{\u203e}}{\mathit{\theta}}}_{\mathrm{1}},\phantom{\rule{0.125em}{0ex}}\mathrm{\dots},\phantom{\rule{0.125em}{0ex}}{\stackrel{\mathrm{\u203e}}{\mathit{\theta}}}_{N}]$, while the corresponding measured blade loads m are stacked into matrix $\mathbf{M}=[{\mathit{m}}_{\mathrm{1}},\phantom{\rule{0.125em}{0ex}}\mathrm{\dots},\phantom{\rule{0.125em}{0ex}}{\mathit{m}}_{N}]$, obtaining
The model coefficients are then computed by least squares as
Measured loads m_{M} are defined as
where m is given in Eq. (4) and r is the residual with covariance Q=E[rr^{T}]. Residuals are assumed to be zeromean and colored and are due to measurement noise and unmodeled physics (Jategaonkar, 2015). Given the model coefficients, a maximum likelihood (Strutz, 2016) estimate θ_{E} of the wind states can be computed online during turbine operation from the measured loads m_{M} from Eqs. (4) and (8) as follows:
2.2.1 Density correction
Aerodynamic moments can be written as
where $q=\mathrm{1}/\mathrm{2}\mathit{\rho}{V}^{\mathrm{2}}$ is the dynamic pressure, A=πR^{2} is the rotor disk area, and C is a nondimensional coefficient. A correction for density can be simply obtained as
where ρ_{ref} is a reference density and ρ_{i} the density corresponding to measurement ${m}_{{\mathrm{A}}_{i}}$.
However, blade load sensors measure not only aerodynamic loads but also the effects of inertia and gravity, which do not depend on air density. Inertial loads for a balanced rotor spinning at constant rotor speed do not generate rotating 1P harmonic components and hence do not appear in Eq. (4). On the other hand, gravitational terms generate 1P loads represented by the nonhomogeneous term m_{0} in that same equation. According to Bertelè et al. (2017), this term can be written as
The first term is a gravityinduced load due to the rotor deformation caused by aerodynamic loads; for example, if the blade bends under the push of thrust, the resulting deformation generates a nonnull moment arm for gravity with respect to the blade root where the load sensor is located, resulting in a 1P load. This term is proportional to dynamic pressure and can be corrected for density. The second term g accounts for inplane and outofplane gravityinduced loads, the latter being caused by blade precone, prebend, and rotor uptilt. This term does not depend on density, and hence it should be eliminated by the equations before a density correction can be applied. To this end, first the model coefficients of Eq. (4) were identified for a very low wind speed, just above cutin. Here the effects caused by qARC are negligible, and hence g≈m_{0}. Having first identified the gravity term g and then having eliminated it from model (Eq. 4), each measured load was finally corrected for density using Eq. (11).
2.3 Wind parametrization in the field
Before wind states can be estimated at run time from measured loads using Eq. (9), the model coefficients must be identified through the simultaneous measurements of wind states and associated loads using Eq. (7). This section presents a practical method to perform this task, based on the use of a standard IECcompliant (IEC, 2017) hubtall met mast. A similar procedure could be used to identify the observer for a specific wind turbine type. Having obtained the model coefficients, one should be able to use the same observer for other installations of that same wind turbine type. Although there is yet no direct demonstration of this assertion, it seems reasonable to assume that wind turbines of the same model will have a similar 1P response to shears and misalignment angles. Additionally, Bottasso and Riboldi (2015) showed that the method is robust to the typical changes occurring in some of the wind turbine parameters across different installations of the same wind turbine type, including changes in the stiffness of foundations, orographic effects, imbalance due to pitch misalignment, miscalibration of the load sensors, and changes in airfoil lift and drag due to soiling/erosion.
2.3.1 Test site
Data were measured at a test site between 19 October and 29 November 2017 for a campaign unrelated to the present study. Since the data were collected long before the beginning of this work, the data had to be used “as is” without the possibility of verifications, calibrations, or any other activity meant to improve the knowledge of the conditions in which the dataset was collected.
Figure 2 shows a panoramic view of the test site (Bromm et al., 2018), which is located in Germany a few kilometers inland from the Baltic Sea and is characterized by gentle hills, open fields, and forests. A 3.5 MW eno114 turbine is installed at the site. The machine (labeled WT1 in the figure) has a 92 m hub height and a rotor diameter of 114.9 m.
A met mast is situated at about 2.5 D (diameters) from the turbine. Wind direction was measured at a height above ground of 89.3 m with a Thies GmbH wind vane, while wind speed measurements were obtained with three cup anemometers produced from the same company and located at 89.3 m, 91.5 m, and at the lower tip of the rotor (about 34 m). All measurements obtained on the mast were shifted in time on account of the distance between turbine and met mast, the time delay being computed from the average wind speed.
A second turbine (labeled WT2) is also present on site, and its wake affects the met mast and WT1 for easterly and southeasterly winds. Similarly, the wake of WT1 affects the met mast for northern wind directions. All these conditions were discarded from the training dataset in addition to all other situations when WT1 was not in a normal power production state. A forest of 15–20 m tall trees is located 300 m east of WT1; as only wind directions $\mathrm{\Gamma}\in [\mathrm{180},\mathrm{340}]$^{∘} were considered in this work, this high roughness area was never in the inflow direction. On the other hand, the town of Brusow is located about 1 km to the west of the site, and its effects on the inflow are unknown. A test campaign conducted at the same site in the period July–November of the previous year revealed an almost equal distribution of unstable, neutral, and stable conditions, as measured by an eddy covariance station (Bromm et al., 2018).
Synchronized turbine and blade load data were sampled at 10 Hz on WT1. Blades 1 and 3 were equipped with strain gages, installed in close proximity to the blade roots and measuring both flapwise and edgewise bending components; unfortunately, however, the same load sensors were not installed on blade 2. To reconstruct the missing load components, the measurements of blades 1 and 3 were shifted by $\pm \mathrm{2}\mathit{\pi}/\mathrm{3}$, averaged together and then attributed to blade 2. This approximation assumes that neighboring blades experience the same loads when they are at the same azimuthal position, which is reasonable because loads and wind states are timeaveraged quantities linked by a steady loadwind model (cf. Eq. 4).
In general, sensors deployed in the field cannot be assumed to be always exactly calibrated, and they may suffer from a variety of issues that affect the quality of the measurements that they provide. To address this problem, it is useful to devise simple and practical ways to correct the measurements even when the root cause of the problem is unknown. Here, consistent mismatches between the longterm mean readings of the two blade load sensors were observed. To correct for this inconsistency, the signals were adjusted a posteriori by a factor s to enforce the same mean loads on the two blades. This was obtained by scaling the measurements as ${\stackrel{\mathrm{\u203e}}{m}}_{\mathrm{1}}(\mathrm{1}+s)={\stackrel{\mathrm{\u203e}}{m}}_{\mathrm{3}}(\mathrm{1}s)$, with s=0.0274. Clearly, this is different from a true calibration meant to ensure the correct reading of a known quantity. However, since the data had been collected prior to this study and no additional information was available, this is probably the only possible adjustment that can be applied. Additionally, the azimuth signal was corrected to account for sensor bias and dynamic effects, as explained in Schreiber et al. (2020). The yaw encoder signal was also corrected for an apparent inconsistency in its readings, as explained later in this section. The turbine onboard wind vane was found to correlate well with the signal provided by the mast after correcting for time delays due to their different locations. However, for coherence with the reference wind speed measurements, also the wind direction reference was taken as the one provided by the mast.
2.3.2 Wind shears
The met mast present at the test site reaches only up to hub height; this is also the typical case of IECcompliant met masts used for certification (IEC, 2017). The three anemometers at 34, 89, and 92 m can be used to estimate the vertical shear over the lower half of the rotor, which, however, in general differs from the shear computed over the whole rotor height.
To address this issue, the sectoreffective wind speed (SEWS) estimation method described in Schreiber et al. (2020) was employed to obtain a rotoreffective reference for the shears. In short, the method works as follows: the blades are used as local speed sensors that, scanning the rotor disk, provide average speeds over four rotor quadrants. By using the two lateral and the lower quadrants, the shear over the lower part of the rotor disk can be computed. This quantity is validated with respect to the shear measured by the met mast, assumed as a ground truth. Then, having verified a good correlation between the measured and estimated shears over the lower part of the rotor, the SEWSs for all four quadrants are used to calculate the wind shear over the whole rotor disk. A brief overview of the SEWS estimator is reported next, and the interested reader is referred to Schreiber et al. (2020) for further details.
The rotor cone coefficient is defined as
where β is the pitch angle, $\mathit{\lambda}=\mathrm{\Omega}R/V$ the tip speed ratio, Ω the rotor speed, m_{i} the outofplane bending load of the ith blade, and ψ_{i} its azimuthal position. The dependency of the coefficient on the azimuthal position of the blade is primarily dictated by the effects of gravity which for an uptilted rotor generate an outofplane bending moment that needs to be taken into account. Accuracy can be improved by considering the deformation of tower and rotor depending on operating conditions (Bottasso et al., 2018), an effect that was neglected here for simplicity. Coefficient C_{m} was computed from a complete aeroelastic model of the turbine, implemented with the code FAST (Jonkman and Jonkman, 2018). Inverting Eq. (13), a lookup table (LUT) is generated that returns the bladeeffective wind speed V_{i} given measured blade pitch angle, rotor speed, azimuthal blade position, bending moment, and density:
This way each individual blade is turned into a local wind speed sensor which scans the rotor disk. Since this local measurement is noisy, the rotor disk is divided into sectors of area A_{S}, and a sectorequivalent wind speed is computed as
Here the four sectors shown in Fig. 3 were used. This yields four measurements of the local speed at the rotor disk, namely above, below, and to the sides of the hub center. Bottasso et al. (2018) showed that, for a linear shear and a sector 90^{∘} wide, the SEWS corresponds to the inflow speed at a distance of approximately $\mathrm{2}/\mathrm{3}R$ from the hub center.
The rotoreffective horizontal linear shear can be computed inserting the SEWSs in Eq. (1) to get
The analysis of the vertical shears requires some care. In fact, the linear vertical shear estimated by the met mast and by the sectoreffective speeds is computed from measurements obtained at different heights above ground; as such, they are not directly comparable because shear has typically a nonlinear variability with height. To address this issue, a power law is first fitted to the measurements to accurately represent the shape of the wind speed gradient; once the power law parameters have been determined, linear shears are computed for mast and observer between the same two heights, resulting in comparable quantities. As previously mentioned, the vertical shear can be parameterized in various ways. In this work, a linear fit was chosen in order to match the linear definition of the horizontal shear because this avoids the need for horizontal shear measurements by using the rotor symmetry (Bertelè et al., 2019).
More precisely, the calculation of the linear shears is conducted as follows. The power law profile is defined as
where H is the height of the hub, V_{ref} the wind speed at that point, and α the power law exponent. Given n measurements V_{i} at z_{i}, the parameters of the power law are computed by the following best fit:
Notice that two measurements at two different heights are sufficient to estimate the power law since it depends on only the two free parameters V_{ref} and α. Having solved the fitting problem (Eq. 18), the linear shear κ_{v} between two generic heights z_{A} and z_{B} is computed as
Figure 4a shows the correlation between 10 min averages of the vertical shears obtained by the met mast and by the sectoreffective wind speeds on the lower half of the rotor. Only wind directions between 170 and 215^{∘} are considered, when the turbine and met mast are aligned. The power law for the met mast was obtained by using all three speed measurements, although the two at 89.3 and 91.5 m above ground are almost coincident. For the sectoreffective observer the power law was obtained by using the two estimates $({V}_{\mathrm{S},\mathrm{left}}+{V}_{\mathrm{S},\mathrm{right}})/\mathrm{2}$ at z=0 and V_{S,down} at $z=\mathrm{2}/\mathrm{3}R$ (although this latter value is strictly valid only for linear shears). For both cases, the power law coefficients were first computed using Eq. (18), and then the lowerhalfrotor linear shear was obtained from Eq. (19) using z_{A}=0 and ${z}_{B}=R$. The figure shows that there is a good correlation between the two lowerhalfrotor shears, resulting in a Pearson's coefficient of 0.91.
The figure also shows that the linear fit (dashed red line) has a different slope than the ideal match (solid black line). This could be due to a nonideal power law profile but also by a nonexact elimination of the effects of gravity, for example because of a different position of the load sensors in the model and reality or a slightly modified uptilt on account of tower deformation. Unfortunately, not enough information on the present experimental setup was available to determine the cause of this discrepancy with certainty. However, the results presented later in Sect. 3 were pragmatically corrected to account for this error: the slope deviation was evaluated from Fig. 4a, and the estimates were modified accordingly in order to generate results lying on the bisector.
For the same data points, Fig. 4b shows the correlation between the vertical shears obtained by the met mast and by the sectoreffective estimator over the complete rotor. Here again the power law for the met mast was obtained by using all three speed measurements. For the sectoreffective estimator the power law was obtained by using Eq. (18) with the three estimates V_{S,up} at $z=\mathrm{2}/\mathrm{3}R$, $({V}_{\mathrm{S},\mathrm{left}}+{V}_{\mathrm{S},\mathrm{right}})/\mathrm{2}$ at z=0, and V_{S,down} at $z=\mathrm{2}/\mathrm{3}R$, although here again the vertical coordinates are strictly valid only for a linear shear. For both cases, the fullrotor linear shear was computed from Eq. (19) using z_{A}=R and ${z}_{B}=R$. It should be noted that, since the height of the top anemometer reaches only up to hub height, for the met mast the calculation of the fullrotor shear implies a considerable extrapolation outside of the available measurements.
A comparison of Fig. 4a and b shows that, in the fullrotor case, there is a lower correlation between the met mast and the SEWS observer than in the lowerhalf rotor case. This indicates that the shear changes over the height of the rotor disk. In addition, as expected for a typical power law in which the profile gradient increases with height, the lowerhalfshear coefficient is typically higher than the fullrotor one.
Based on these results, it appears that the rotoreffective shear used for identifying the model of Sect. 2.2 would require a tall met mast or other wind measurement devices such as lidars or sodars capable of scanning the inflow reaching the top of the rotor. Here – as such a tall mast was not available – an alternative approach was adopted: the sectorequivalent wind speed method was used to virtually extend the metmast measurements to the required height. Based on the good correlation shown in Fig. 4a for the lowerhalfrotor shear, it was concluded that the two lateral and the lower sectorequivalent speeds are sufficiently accurate for the purpose of estimating shears. Since the top sector speed is based on exactly the same calculation procedure as the other ones, all four speeds were then used to estimate the fullrotor shear, which in turn was adopted as the reference for the identification of the model of Sect. 2.2.
Unfortunately a similar validation cannot be performed for the horizontal shear with the present met mast because of the lack of multiple lateral measurements. However, the horizontal shear is based on the same sectorequivalent wind speeds that estimate the vertical shear with good accuracy so that there is no reason to believe that Eq. (16) should not provide a similarly goodquality estimate. Additionally, the horizontal shear based on the two lateral sectoreffective wind speeds was shown in Schreiber et al. (2020) to track the movement of an impinging wake with remarkable accuracy.
2.3.3 Wind misalignment angles
The met mast is equipped with a single wind vane measuring the wind direction Γ at hub height. Unfortunately, this means that only a pointwise measurement is available instead of the rotorequivalent one that would be ideally necessary for the training of the loadharmonic method of Sect. 2.2. This is a limit of the current setup and of the present attempt at validating the approach. Nonetheless, a pragmatic choice was made here to use this signal as a proxy for the rotoreffective horizontal wind direction. The misalignment angle between turbine and wind was obtained by subtracting the absolute yaw angle of the nacelle from the metmastmeasured wind direction. The result was filtered with a 1 min moving average to remove the faster fluctuations. Figure 5a shows 10 min averages of the resulting metmast yaw misalignment angle Φ_{MM} plotted as a function of wind direction Γ. The clear trend visible in the plot is probably due to a miscalibration of the nacelle yaw encoder. Indeed, Bromm et al. (2018) also noticed a nonconstant offset when comparing the turbine supervisory control and data acquisition (SCADA) orientation with the one provided by a temporarily installed GPS system. This trend was removed using the first 10 d of data, excluding waked directions, obtaining Fig. 5b.
As the current setup does not provide for measurements of the upflow, the rotational symmetry of the rotor was used to compute the relevant model coefficients.
2.3.4 Wind speed and density
Since the loadwind model expressed by Eq. (4) depends on the operating conditions, a rotoreffective wind speed was computed with the torque balance equation (Ma et al., 1995; Van der Hooft and Engelen, 2004; Soltani et al., 2013; Schreiber et al., 2020) and used as the scheduling parameter of the wind observer. Figure 6 shows an excellent correlation for the 10 min averages of the computed rotoreffective wind speed and the metmast hubheight speed, with a Pearson coefficient of 0.988 and a root mean square (RMS) error of ϵ_{RMS}=0.418 m s^{−1}. Density was obtained from the ideal gas law based on temperature since no additional information was available, and it was used to rescale the load measurements.
3.1 Model identification
The observer coefficients were identified with Eq. (7) using the horizontal and vertical shears obtained from the sectoreffective wind speeds, and with the yaw misalignment angle computed from the metmast wind vane and the nacelle yaw encoder corrected according to Fig. 5. The upflow model coefficients were obtained from the rotational symmetry of the rotor behavior. Load measurements were corrected for density, the reference value being set to 1.238 kg m^{−3}.
The model coefficients were scheduled as functions of the rotoreffective wind speed computed from the torque balance equation. The wind speed nodes of the linear parameter varying (LPV) model (Eq. 4) were defined as $V=\mathit{\{}\mathrm{4},\phantom{\rule{0.125em}{0ex}}\mathrm{5},\phantom{\rule{0.125em}{0ex}}\mathrm{6.5},\phantom{\rule{0.125em}{0ex}}\mathrm{8},\phantom{\rule{0.125em}{0ex}}\mathrm{9},\phantom{\rule{0.125em}{0ex}}\mathrm{10},\phantom{\rule{0.125em}{0ex}}\mathrm{12},\phantom{\rule{0.125em}{0ex}}\mathrm{13.5}\mathit{\}}$ m s^{−1}. This means that model coefficients were computed at each of these wind speed nodes, while any speed within the range [4, 13.5] m s^{−1} – i.e., not necessarily at the nodes – was used for identification by linearly distributing its contributions to the two neighboring nodes. At run time, the coefficients were interpolated from the LPV model based on the current wind speed.
Table 1 shows the range covered by each parameter within the training dataset.
About 15 % of the available data was used for identification, leaving about 370 h of measurements for validation. In the following, the performance of the harmonic observer is evaluated solely based on the validation dataset, i.e., excluding all data points used for training.
3.2 Wind observer performance
Models were identified based on different time averages of the raw 10 Hz data. Here, the two cases of 1 and 10 min averages are presented because they correspond to the typical outputs of standard SCADA systems. In both cases, the raw data points were the same; this means that the 1 min model was identified on 10 times more loadstate pairs than in the case of the 10 min model.
An overview of the performance of the two models is given in Fig. 7 (for the 10 min case) and Fig. 8 (for the 1 min case). The figures report correlations between reference and observed parameters using the validation subset for wind speeds above 8 m s^{−1}. For each parameter, one per subplot, the reference state is shown on the x axis, whereas the observed one is on the y axis.
A comparison of the 10 and 1 min cases shows that results are essentially identical for the shears. For the misalignment angle, results are very slightly better using the longer time window notwithstanding the smaller number of loadstate pairs used for identification. Probably this is because longer time averaging alleviates the effects of outliers. Based on these results, the rest of the paper only considers the 10 min case.
Considering the shears, Fig. 7b and c show that the Pearson's correlation coefficient (R) is above 0.9, and the RMS error ϵ_{RMS} is of the order of 10^{−3}. The yaw misalignment angle is less accurate possibly because the reference is pointwise, whereas the estimate is rotoreffective. Indeed, investigations at the same site with a more complete setup including a lidar profiler reported significant veer at the inflow (Bromm et al., 2018). However, with a correlation coefficient of 0.85 and an ϵ_{RMS} of 1.9^{∘}, the matching is still good.
It is very interesting to observe that even a model trained only with 10 min averages is still able to provide for timeresolved estimates of the parameters. To illustrate this fact, Fig. 9 reports a 10 Hz time history of the vertical shears from the validation subset. The figure corresponds to about 2 d of operation during which the wind direction, Fig. 9c, was $\mathrm{\Gamma}\in [\mathrm{145},\mathrm{260}]$^{∘}. Turbine and met mast are roughly aligned for $\mathrm{\Gamma}\in [\mathrm{177.5},\mathrm{215}]$^{∘}; WT1 is in the wake of WT2 for approximatively $\mathrm{\Gamma}\in [\mathrm{120},\mathrm{170}]$^{∘}, the two directions being indicated in the plot with two dashed horizontal lines. Figure 9a shows the lowerhalfrotor shears measured at the met mast and by the sectorequivalent speeds. Although some discrepancies are present, the figure shows that the sectoreffective observer is capable of following the main changes in shear captured by the met mast. The main discrepancies can be found between 14:00 on 21 October and about 04:00 on 22 October when WT1 is in the wake of WT2 or in its close proximity. However, one should not forget that the two estimates correspond to two locations spaced 2.5 D apart and that the exact ground truth at the rotor disk – where the observers operate – is unknown. Figure 9b shows the rotorequivalent shear estimated by Eq. (9) based on rotor harmonics and its reference quantity obtained by the sectorequivalent speeds. The two vertical shears are in excellent agreement even with respect to relatively fast fluctuations.
To provide for a more complete statistical characterization of the observer performance, the 10 min data points were binned for the various relevant parameters. For each bin, the mean absolute error (MAE) between the estimated θ_{E} and reference θ_{R} wind parameters was computed as $\mathit{\u03f5}=\mathrm{1}/N\sum _{i}^{N}{\mathit{\theta}}_{{\mathrm{R}}_{i}}{\mathit{\theta}}_{{\mathrm{E}}_{i}}$.
Figure 10 shows the MAE ϵ for yaw misalignment (Fig. 10a) and vertical (Fig. 10b) and horizontal shear (Fig. 10d), plotted as functions of binned wind speed, for various binned turbulence intensity (TI) levels. The number of available hours of data is reported in Fig. 10c, to help determine the statistical significance of the results. Looking at the yaw angle results, it appears that the maximum error is about 3^{∘} and that accuracy tends to increase with higher wind speeds. Moreover, TI appears to have only a small effect on the results. As previously mentioned, this can be attributed to the fact that 1P harmonics are dominated by the four wind states and only modestly affected by turbulent fluctuations.
The error in the vertical shear includes the error between the met mast and the sectoreffective observer of Sect. 2.3.2. Even in this case the error is small, and effects of TI are present but relatively mild. The figure also reports the horizontal shear, whose error – although very small – might not be very indicative: since no reference value was available from the met mast for this quantity, only the error with respect to the sectoreffective observer of Sect. 2.3.2 could be quantified.
Figure 11 reports the results for varying binned air densities. The plots show that the density correction of Sect. 2.2.1 is not perfect probably because of an only approximate identification of the gravity term in Eq. (12).
Finally, Fig. 12 reports the results for varying wind directions. Looking at the vertical shear, the best results are obtained for wind directions between 170 and 210^{∘} when turbine and met mast are aligned, whereas the error increases significantly for other wind directions. When turbine and met mast are not aligned, the two can be subjected to slightly different inflows on account of orographic and vegetationinduced effects. This indicates once again that, as noted earlier on, the information provided by the reference met mast cannot be regarded as an absolute ground truth. The yaw misalignment angle seems to be less influenced by these local effects, which might induce stronger local changes in shear than in direction at this particular site.
This paper has presented the application of a previously published harmonicbased wind sensing method to an experimental dataset. The setup at the test site is not complete enough to provide for a true field validation of the method. However, it is representative of a practical scenario in which, by using a hubtall certification met mast, the method is trained for a given turbine model before being deployed on assets of that same type at other production sites. After having explained the methodology and described the test site, the paper has also formulated a new method to extend the shear measured by a hubtall mast to the tip of the rotor in order to compute a fullrotor shear.
Based on the results analyzed herein, and notwithstanding the limits of the present dataset, the following conclusions can be drawn:

There is a good correlation between met mast and estimated lowerhalfrotor shears, with Pearson's coefficients above 0.9 and RMS errors around $\mathrm{4}\times {\mathrm{10}}^{\mathrm{2}}$.

There is an excellent correlation between the fullrotor shear extended above the mast and the one estimated by harmonic loads, with Pearson's coefficients above 0.99 and RMS errors around $\mathrm{4}\times {\mathrm{10}}^{\mathrm{3}}$.

Training with 1 or 10 min averages produces shear estimates of a very similar quality, but there is a marginal improvement of the wind direction for the longer time window. This is probably due to the noisier nature of wind direction, which is measured here only at hub height.

Notwithstanding a training based on 10 min averages, the quality of the correlation between estimates and references does not only apply to 10 min quantities, but it also extends to timeresolved 10 Hz signals. In this sense, the observer seems capable of following relatively fast changes in shear. This might be useful for certain application scenarios, as for example the tracking of horizontal shears induced by wake interactions.

There is a nonnegligible effect of nonexact wind–mast–turbine alignment. In this sense, the actual quality of the correlation might be even better than what appears from the results shown here. This is in fact an intrinsic limit of field testing, for which an exact ground truth is in general difficult if not impossible to obtain. Realistic simulations and wind tunnel studies as the ones reported in Bertelè et al. (2017, 2018, 2019) – where the ground truth is known – may help in this sense.

Yaw misalignment is also estimated with reasonable quality, maximum errors being in general below 3^{∘}. However, the results here are less conclusive due to the fact that the metmast reference is a pointwise measurement that might not fully represent rotoreffective conditions.

There is only a modest effect of TI, which supports the hypothesis that 1P harmonics are mostly driven by “deterministic” wind characteristics and less affected by turbulent fluctuations.

Notwithstanding the complicated effect of gravity on harmonic load components, its presence can be eliminated with enough accuracy to allow for a reasonably precise density correction.
The main limits of the present dataset are as follows: independent reference measurements for horizontal shear and upflow were completely missing, yaw misalignment was measured only at a point instead of over the rotor disk, and the vertical shear had to be extended over the hub by the use of another estimation method. Although the utmost care was put into the reconstruction of the fullrotor vertical shear, this operation still had to rely on the same blade load measurements used by the harmonic estimator, which is clearly a weakness. Additionally, as the test campaign was performed prior to the present study, the dataset had to be used as is, without any possible verification, correction, or calibration of the sensors. Other less substantial limitations are also present, for example those caused by the missing load sensors on one of the blades.
A continuation of this work would greatly benefit from access to a more complete dataset. Multiple, independent rotoreffective measurements of the inflow in close proximity to the rotor disk would be necessary to establish an effective ground truth. This would enable a better characterization of the accuracy of this method and a way to study the effects induced by training with a standard hubtall mast. A remaining open point is the sensitivity of the method to phenomena like aging, soiling, and rotor imbalances. Indeed, any exogenous cause affecting the 1P response will be interpreted by the harmonic observer as a change in the wind states. Some reassuring results have already been reported by Bottasso and Riboldi (2015), although a more thorough experimental investigation is necessary. Finally, it remains to be shown that the method can indeed be trained on a turbine and then applied to another machine of that same model at another site; although this seems to be a very reasonable assumption, the evidence that this is indeed possible is lacking.
A  Rotor area 
C_{m}  Cone coefficient 
H  Height of the hub above ground 
m  Blade bending moment 
m  Vector of moment harmonics 
N  Number of available data points 
q  Dynamic pressure 
R  Rotor radius or Pearson's coefficient 
Q  Covariance matrix 
V  Wind speed 
V_{h}  Wind speed at hub height 
V_{PL}(z)  Power law wind speed profile 
V_{S}  Sectoreffective wind speed 
V_{TB}  Torquebalance rotoreffective wind speed 
$\stackrel{\mathrm{\u0303}}{v}$  Nondimensional horizontal cross flow at hub height 
$\stackrel{\mathrm{\u0303}}{w}$  Nondimensional vertical cross flow at hub height 
x, y, z  Hubcentered nacelleattached axes 
β  Pitch angle 
Γ  Wind direction 
ϵ  Mean absolute error 
θ  Wind state vector 
κ_{h}  Horizontal shear 
κ_{v}  Vertical shear 
λ  Tip speed ratio 
ρ  Air density 
ϕ  Yaw misalignment angle 
χ  Upflow angle 
ψ  Azimuth angle 
Ω  Rotor speed 
(⋅)^{T}  Transpose 
(⋅)^{IP}  Inplane component 
(⋅)^{OP}  Outofplane component 
(⋅)_{1c}  1P cosine amplitude 
(⋅)_{1s}  1P sine amplitude 
(⋅)_{E}  Estimated quantity 
(⋅)_{MM}  Metmast measurement 
(⋅)_{ref}  Reference quantity 
(⋅)_{RMS}  Root mean square 
1P  Once per revolution 
MAE  Mean absolute error 
Lidar  Light detection and ranging 
LPV  Linear parameter varying 
LUT  Lookup table 
OEM  Original equipment manufacturer 
RMS  Root mean square 
SCADA  Supervisory control and data acquisition 
SEWS  Sectoreffective wind speed 
Sodar  Sound detection and ranging 
TI  Turbulence intensity 
WT  Wind turbine 
The operational data and turbine model used in this research are property of Eno Energy Systems GmbH.
MB conducted the main research work and prepared a first draft of the manuscript. CLB developed the core idea of loadbased wind sensing, formulated the harmonic and sectoreffective methods, supervised the research, and contributed to the writing of the paper. JS assisted in the measurement postprocessing and analysis. All authors provided important input to this research work through discussions and feedback and by improving the manuscript.
The authors declare that they have no conflict of interest.
The authors express their gratitude to Stefan Bockholt and Alexander Gerds of eno energy systems GmbH, who granted access to the measurement data and turbine model, and to Marijn van Dooren, Anantha Sekar, and Martin Kühn of ForWind Oldenburg, who shared insight on the data.
This work has been supported by the CompactWind II project (FKZ: 0325492G), which receives funding from the German Federal Ministry for Economic Affairs and Energy (BMWi).
This work was supported by the German Research Foundation (DFG) and the Technical University of Munich (TUM) in the framework of the openaccess publishing program.
This paper was edited by Katherine Dykes and reviewed by two anonymous referees.
Bachmann Electronic GmbH: Hardware und installation, available at: https://www.bachmann.info/produkte/conditionmonitoringsystem/hardwareundinstallation/, last access: 12 January 2021. a
Bertelè, M. and Bottasso, C. L.: Initial results from the field testing of the rotorasasensor concept, J. Phys.: Conf. Ser., 1452, 012074, https://doi.org/10.1088/17426596/1452/1/012074, 2020. a
Bertelè, M., Bottasso, C. L., Cacciola, S., Daher Adegas, F., and Delport, S.: Wind inflow observation from load harmonics, Wind Energ. Sci., 2, 615–640, https://doi.org/10.5194/wes26152017, 2017. a, b, c, d, e, f, g, h
Bertelè, M., Bottasso, C. L., and Cacciola, S.: Simultaneous estimation of wind shears and misalignments from rotor loads: formulation for IPCcontrolled wind turbines, J. Phys.: Conf. Ser., 1037, 032007, https://doi.org/10.1088/17426596/1037/3/032007, 2018. a, b, c
Bertelè, M., Bottasso, C. L., and Cacciola, S.: Brief communication: Wind inflow observation from load harmonics – wind tunnel validation of the rotationally symmetric formulation, Wind Energ. Sci., 4, 89–97, https://doi.org/10.5194/wes4892019, 2019. a, b, c, d, e
Bottasso, C. L. and Riboldi, C. E. D.: Estimation of wind misalignment and vertical shear from blade loads, Renew. Energ., 62, 293–302, https://doi.org/10.1016/j.renene.2013.07.021, 2014. a, b, c
Bottasso, C. L. and Riboldi, C. E. D.: Validation of a wind misalignment observer using field test data, Renew. Energ., 74, 298–306, https://doi.org/10.1016/j.renene.2014.07.048, 2015. a, b, c, d
Bottasso, C. L., Croce, A., and Riboldi, C. E. D.: Spatial estimation of wind states from the aeroelastic response of a wind turbine, in: The Science of Making Torque from Wind (TORQUE 2010), 28–30 June 2010, Heraklion, Crete, Greece, 2010. a
Bottasso, C. L., Cacciola, S., and Schreiber, J.: Local wind speed estimation, with application to wake impingement detection, Renew. Energ., 116, 155–168, https://doi.org/10.1016/j.renene.2017.09.044, 2018. a, b, c, d, e
Bromm, M., Rott, A., Beck, H., Vollmer, L., Steinfeld, G. and Kühn, M.: Field investigation on the influence of yaw misalignment on the propagation of wind turbine wakes, Wind Energy, 21, 1011–1028, https://doi.org/10.1002/we.2210, 2018. a, b, c, d
Cacciola, S., Bertelè, M., and Bottasso, C. L.: Simultaneous observation of wind shears and misalignments from rotor loads, J. Phys.: Conf. Ser., 753, 05200218, https://doi.org/10.1088/17426596/753/5/052002, 2016. a
Coleman, R. P. and Feingold, A. M.: Theory of selfexcited mechanical oscillations of helicopter rotors with hinged blades, Technical Report, NACA TN 1351, available at: https://digital.library.unt.edu/ark:/67531/metadc56181/m2/1/high_res_d/19930084623.pdf (last access: 19 May 2021), 1958. a
fos4X GmbH: retrox, available at: https://fos4x.com/en/solutions/retrox/, last access: 12 January 2021. a
Hirth, B. D., Schroeder, J. L., Gunter, W. S., and Guynes, J. G.: Coupling Doppler radarderived wind maps with operational turbine data to document wind farm complex flows, Wind Energy, 18, 529–540, https://doi.org/10.1002/we.1701, 2015. a
IEC – International Electrotechnical Commission: DV IEC 61400121, Technical Report, IEC, Geneva, Switzerland, 2017. a, b, c
Jategaonkar, R. V.: Flight Vehicle System Identification – A TimeDomain Methodology, American Institute of Aeronautics and Astronautics Inc., Reston, Virginia, USA, https://doi.org/10.2514/4.102790, 2015. a
Jonkman, J. and Jonkman, B.: FAST v7, available at: https://nwtc.nrel.gov/FAST7 (last access: 26 June 2020), 2018. a
Lang, S. and McKeogh, E.: Lidar and sodar measurements of wind speed and direction in upland terrain for wind energy purposes, Remote Sens., 3, 1871–1901, https://doi.org/10.3390/rs3091871, 2011. a
Ma, X., Poulsen, N., and Bindner, H.: Estimation of Wind Speed in Connection to a Wind Turbine, Technical Report, Informatics and Mathematical Modelling, Technical University of Denmark, Copenhagen, Denmark, 1995. a
Mikkelsen, T. K.: Lidarbased research and innovation at DTU Wind Energy – A review, J. Phys.: Conf. Ser., 524, 012007, https://doi.org/10.1088/17426596/524/1/012007, 2014. a, b
Murphy, P., Lundquist, J. K., and Fleming, P.: How wind speed shear and directional veer affect the power production of a megawattscale operational wind turbine, Wind Energ. Sci., 5, 1169–1190, https://doi.org/10.5194/wes511692020, 2020. a
Peña Diaz, A., Sjöholm, M., Mikkelsen, T. K., and Hasager, C. B.: Inflow characterization using measurements from the SpinnerLidar: the ScanFlow experiment, J. Phys.: Conf. Ser., 1037, 052027, https://doi.org/10.1088/17426596/1037/5/052027, 2018. a, b, c
Schlipf, D., Fleming, P., Haizmann, F., Scholbrock, A., Hofsäß, M., Wright, A., and Cheng, P. W.: Field testing of feedforward collective pitch control on the CART2 using a nacellebased lidar scanner, J. Phys.: Conf. Ser., 555, 012090, https://doi.org/10.1088/17426596/555/1/012090, 2014. a, b
Scholbrock, A., Fleming, P., Fingersh, L., Wright, A., Schlipf, D., Haizmann, F., and Belen, F.: Field testing lidarbased feedforward controls on the NREL controls advanced research turbine, American Institute of Aeronautics and Astronautics, Reston, Virginia, USA, https://doi.org/10.2514/6.2013818, 2013. a, b
Schreiber, J., Cacciola, S., Campagnolo, F., Petrović, V., Mourembles, D., and Bottasso, C. L.: Wind shear estimation and wake detection by rotor loads – First wind tunnel verification, J. Phys.: Conf. Ser., 753, 032027, https://doi.org/10.1088/17426596/753/3/032027, 2016. a, b
Schreiber, J., Bottasso, C. L., and Bertelè, M.: Field testing of a local wind inflow estimator and wake detector, Wind Energ. Sci., 5, 867–884, https://doi.org/10.5194/wes58672020, 2020. a, b, c, d, e, f, g, h, i, j, k, l, m, n
Schubel, P. J., Crossley, R. J., Boateng, E. K. G., and Hutchinson, J. R.: Review of structural health and cure monitoring techniques for large wind turbine blades, Renew. Energ., 51, 113–123, https://doi.org/10.1016/j.renene.2012.08.072, 2013. a
Simley, E. and Pao, L. Y.: Evaluation of a wind speed estimator for effective hubheight and shear components, Wind Energy, 19, 167–184, https://doi.org/10.1002/we.1817, 2016. a
Soltani, M., Knudsen, T., Svenstrup, M., Wisniewski, R., Brath, P., Ortega, R., and Johnson, K.: Estimation of rotor effective wind speed: a comparison, IEEE Trans. Control Syst. Technol., 21, 1155e1167, https://doi.org/10.1109/tcst.2013.2260751, 2013. a
Strutz, T.: Data Fitting and Uncertainty (A practical introduction to weighted least squares and beyond), Springer Vieweg, Wiesbaden, Germany, ISBN 9783658114558, 2016. a
Valldecabres, L., Nygaard, N., VeraTudela, L., von Bremen, L., and Kühn, M.: On the use of dualdoppler radar measurements for very shortterm wind power forecasts, Remote Sens., 10, 1701, https://doi.org/10.3390/rs10111701, 2018. a
Van der Hooft, E. L. and Engelen, T.: Estimated wind speed feed forward control for wind turbine operation optimisation, in: European Wind Energy Conference & Exhibition (EWEC 2004), 22–25 November 2004, London, UK, 2004. a
WindConsult GmbH: Loads, available at: https://www.windconsult.de/lasten/, last access: 12 January 2021. a