Next Article in Journal
DEM Generation from GF-7 Satellite Stereo Imagery Assisted by Space-Borne LiDAR and Its Application to Active Tectonics
Previous Article in Journal
Center-Ness and Repulsion: Constraints to Improve Remote Sensing Object Detection via RepPoints
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Unified Formulation for the Computation of the Six-Degrees-of-Freedom-Motion-Induced Errors in Floating Doppler Wind LiDARs

by
Andreu Salcedo-Bosch
1,
Joan Farré-Guarné
1,
Marcos Paulo Araújo da Silva
1 and
Francesc Rocadenbosch
1,2,*
1
CommSensLab, Department of Signal Theory and Communications, Universitat Politècnica de Catalunya (UPC), E-08034 Barcelona, Spain
2
Institut d’Estudis Espacials de Catalunya (IEEC), Universitat Politècnica de Catalunya, E-08034 Barcelona, Spain
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(6), 1478; https://doi.org/10.3390/rs15061478
Submission received: 3 February 2023 / Revised: 3 March 2023 / Accepted: 4 March 2023 / Published: 7 March 2023
(This article belongs to the Section Atmospheric Remote Sensing)

Abstract

:
This work presents an analytical formulation to assess the six-degrees-of-freedom-motion-induced error in floating Doppler wind LiDARs (FDWLs). The error products derive from the horizontal wind speed bias and apparent turbulence intensity. Departing from a geometrical formulation of the FDWL attitude and of the LiDAR retrieval algorithm, the contributions of the rotational and translational motion to the FDWL-measured total error are computed. Central to this process is the interpretation of the velocity–azimuth display retrieval algorithm in terms of a first-order Fourier series. The obtained 6 DoF formulation is validated numerically by means of a floating LiDAR motion simulator and experimentally in nearshore and open-sea scenarios in the framework of the Pont del Petroli and IJmuiden campaigns, respectively. Both measurement campaigns involved a fixed and a floating Z e p h I R TM 300 LiDAR. The proposed formulation proved capable of estimating the motion-induced FDWL horizontal wind speed bias and returned similar percentiles when comparing the FDWL with the fixed LiDAR. The estimations of the turbulence intensity increment statistically matched the FDWL measurements under all motional and wind scenarios when clustering the data as a function of the buoy’s mean tilt amplitude, mean translational-velocity amplitude, and mean horizontal wind speed.

1. Introduction

The more homogeneous and stronger wind fields found on open-sea environments make offshore wind energy a dependable source of electricity generation [1,2]. In the last few decades, the wind energy industry has made important investments in offshore wind energy [3]. Specifically, a cumulative offshore wind capacity of 35 GW was deployed around the world in 2020, with Europe and China standing out in the market [4]. However, it is still one of the most expensive sources of renewable energy [5], and cost reduction is the main goal of the industry [6,7]. One of the main considerations is to obtain reliable data to assess the viability of a future wind farm based on the location. Usually, offshore wind farms are deployed in shallow waters, but there is a tendency to go further offshore and to deeper waters [8], where a higher wind capacity can be obtained. At those locations, the traditional technology used for atmospheric assessment, namely meteorological masts planted on the seabed, is no longer an option. This motivates the use of alternative remote sensing sensors such as satellites [9,10], radars [11], sodars [12,13], and combined methods [14,15]. Nevertheless, due to the high requirements of the industry in terms of the measurement accuracy, temporal and spatial resolution, and data availability, Doppler Wind LiDARs (DWLs) have become the most trusted devices for wind remote sensing [7].
In general, two types of DWLs are used for wind remote sensing: pulsed and focusable continuous-wave (CW) DWLs [16]. Pulsed DWLs use the time of flight to build up a vertically resolved wind profile, whereas focusable CW DWLs achieve the same by focusing the emitted beam on a particular height by changing its optics. Among pulsed DWLs, two main detection techniques are used: direct detection and heterodyne detection. Direct detection utilizes an optical frequency analyzer, typically an interferometer, and the receiver directly measures the signal photoelectrons [17], whereas in heterodyne detection, the backscattered signal is mixed with the frequency-offset version of the beam of a local oscillator laser, and the receiver measures the resultant beating frequency [18].
When placed over floating platforms or buoys, DWLs are able to measure the wind profile in a cost-effective way [7,19,20]. Moreover, in contrast to metmasts, they can easily be redeployed and thus cover large areas [21]. On the other hand, FDWLs suffer 6 degrees of freedom (DoF) motion, induced by the waves [22,23,24], which increases the variance on the reconstructed wind vector by the LiDAR [25,26]. However, in wind energy standard averaging periods, typically 10 or 30 min, the motion-induced error on the retrieved mean wind vector can be neglected, as it is averaged out [25,27,28,29].
The motion-induced variance in FDWL measurements is encountered at higher wind statistical moments, such as the turbulence intensity ( T I ) [28]. The T I is defined as the ratio between the standard deviation of the horizontal wind speed ( H W S ) and the mean H W S . In the presence of the wave motion, FDWLs measure an apparent added turbulence due to the motion-induced variance, which corrupts the T I measurements [25,26]. T I is one of the parameters of main importance for wind farm design and operation, as wrong estimates of T I could lead to turbine overdesign or malfunctioning [30] and thus higher costs.
In the state of the art, the study of the motion-induced error on FDWLs is twofold: error estimation and error compensation. This paper addresses the former topic, specifically, error estimation in focusable CW FDWLs, which is carried out using a thorough analytical formulation of the system mechanics and the wind vector retrieval algorithm. The error estimation in pulsed FDWLs has been addressed elsewhere [31,32]. The latter topic, i.e., error compensation, using either mechanical [27] or signal processing techniques [25,28,33,34] is outside the scope of the present work.
So far, multiple studies have addressed the study of the motion-induced error on focusable CW FDWLs on both the scan-time and 10 min levels:
On the scan-time level, Tiana-Alsina et al. [35] proposed a basic CW FDWL motion simulator, in which the roll and pitch motions were simulated by means of Euler’s rotation matrices. This simulator provided insights on the H W S error in relation to the wind direction ( W D ) and the buoy tilt but was unable to simulate the effects of the translational motion on the H W S error. Bischoff et al. [36] presented a simulation environment to estimate the uncertainties in the measurements by CW FDWLs using a simplified buoy model. However, the error estimates could not match the measurements for the high-resolution (1 Hz) data. Kelberlau et al. [28] provided a thorough geometrical description of the FDWL dynamics taking into account the 6 DoF motion of the FDWL buoy and provided a method for compensating the apparent T I measured by the SEAWATCH CW FDWL buoy by Fugro TM . This geometrical description permitted a preliminary study of the 6-DoF-motion-induced error in a LiDAR scan but was limited by the oversimplification of assuming a constant value for the initial scan phase. Departing from this geometrical description, Salcedo et al. numerically simulated the 6 DoF CW DWL measurement error [37] and provided a first 6 DoF motion correction algorithm using an unscented Kalman filter [33]. The simulator [37] modeled each of the six DoF as a sinusoidal signal with a given amplitude, frequency, and phase and enabled full understanding of the motion-induced error in a LiDAR scan through the principle of error superposition. Most of the CW FDWL motion-related error studies in the literature resorted to numerical simulations due to the inherent complexity of the velocity–azimuth display (VAD) algorithm, because it involves a least squares fit as a nonlinear operation. These simulations are computationally expensive and difficult to implement [35,36,37,38].
On the 10 min level, Kelberlau et al. [38] recently presented a study on the motion-induced error in the measurement of the mean H W S using a motion simulator. The study theoretically proved the reduced sensitivity of the mean H W S error to the LiDAR motion and its relation to the buoy motion frequency. However, the study was limited to the mean H W S and provided error values up to one order of magnitude smaller than the typically measured values. The study by Bischoff et al. [36] was also limited to studying the bias in the H W S mean, and difficulties in validating the estimated values with the experimental data were encountered. Gutiérrez-Antuñano et al. [25] presented a method to estimate the T I measurement error based on a 2 DoF (roll and pitch) motion simulator as a way to compensate for this error. Overall, it showed good estimation of the motion-induced error, but the performance of the method largely depended on the motion scenario because only two DoF were taken into account. In addition, the simulation was computationally demanding.
The above considerations motivate the study of a unified 6 DoF (rotational and translational) method to quantify the motion-induced error in focusable CW FDWLs: In this work, we provide, for the first time, a complete analytical formulation and calculus of the 6-DoF-motion-induced error on the H W S retrieved by the focusable CW FDWL in a scan. The formulation relies on the VAD LiDAR wind retrieval algorithm and motion parameterization using the characteristic amplitude, frequency, and phase of each of the six DoF (i.e., three rotational and three translational components). Methods to statistically quantify the 10 min H W S bias and the T I increment in response to the floating LiDAR motion are also derived. The framework is the focusable continuous-wave Z e p h I R TM 300 LiDAR, using experimental data gathered during the “Pont del Petroli” and IJmuiden measurement campaigns.
This paper is structured as follows: Section 2 presents the “Pont del Petroli” and IJmuiden measurement campaigns and describes the instrumentation used; Section 3 describes the analytical motion model of the FDWL and presents the motion estimation method; Section 4 validates the analytical model by means of the simulation [37] and the experimental measurements from the “Pont del Petroli” and IJmuiden campaigns; Section 5 gives some concluding remarks.

2. Materials

2.1. Pont del Petroli Campaign

The Pont-del-Petroli (PdP) campaign took place in June 2013 at the “Pont del Petroli” pier (Badalona, Barcelona, Spain) [39]. It aimed to assess the wind measurement capabilities of a proof-of-concept FDWL buoy in the context of the NEPTUNE project by comparing its performance against a fixed LiDAR [27]. Both the fixed and the FDWL were ZephIR 300 focusable continuous-wave DWLs. The fixed LiDAR was situated at the PdP pier in a standstill configuration. The FDWL was deployed offshore 50 m away from the fixed LiDAR. Both ZephIR 300s were identical instruments and were calibrated onshore one beside the other for a period of 3 h. They were configured to measure the wind at a height of 100 m above sea level (a.s.l.), and the calibration ensured identical measurements at 1 s and 10 min intervals. The experiment surroundings were dominated by an open-sea topology in the south and east directions and an urban topology in the west and north directions. The coastline followed the southwestern to northeastern direction.
ZephIR 300 DWL. The ZephIR 300 is a continuous-wave focusable Doppler LiDAR specially conceived for offshore operation [40]. It is able to measure the wind at user-defined heights from 10 m up to 200 m in steps of 1 m [41]. In order to retrieve the wind vector for a particular measurement height, the ZephIR 300 measures the Doppler radial speed along 50 lines of sight (LoSs) every second over a scanning cone of a 30 deg aperture from the zenith. A wedge prism rotating with a uniform circular motion of one revolution per second is used to steer the laser beam in the directions of the 50 LoSs. In uniform wind, the projection of the wind vector along the 50 LoSs, as a function of the azimuth scanning angle (i.e., the so called VAD function), takes the shape of a cosine. The LiDAR uses the VAD algorithm to retrieve the wind vector from the LoSs by fitting a cosine function to the LoS measurements [16]. The ZephIR 300 is a homodyne detection LiDAR, which means that it cannot discern the sign of the Doppler shift. As a result, the VAD function takes the form of a rectified cosine wave, creating an ambiguity of ±180 deg for the W D measurement [42]. The ZephIR 300 is equipped with a wind vane to disambiguate the W D .
FDWL buoy. The FDWL buoy is a proof-of-concept experimental buoy specially conceived to host a ZephIR 300 LiDAR in an offshore operation [43,44]. It weighs 3 tons and has a width of 3.77 m. It consists of a four-floater structure designed to conform with the wind energy measurement requirements as well as to monitor multiple sea-related parameters. Apart from the ZephIR 300, it hosts a MicroStrain 3DM-GX3-45 inertial measurement unit (IMU) that was used to monitor the 6 DoF motion of the lidar.
3DM-GX3-45 IMU. The IMU consists of a GPS unit, an accelerometer, a magnetometer, and a gyro [45], measuring the LiDAR’s tilt and velocity at a 10 Hz sampling frequency. The IMU was located on the bottom of the FDWL buoy structure. We defined, at the location of this FDWL buoy bottom, an inertial right-handed clockwise north–east–down (NED) frame of reference. The gyro separately measured the angular accelerations around the north, east, and down axes, and the accelerometer measured the translational accelerations along these axes. In order to avoid incurring an accumulated error when integrating the accelerations, the 3DM-GX3-45 uses an extended Kalman filter, which permits the monitoring of the buoy’s attitude by considering the GPS and magnetometer measurements. The 3DM-GX3-45 is able to measure the buoy’s velocity and tilt with an accuracy of ±0.1 m/s and ±0.35 deg, respectively.

2.2. IJmuiden Campaign

The IJmuiden campaign took place between 1 April and 1 June 2015 at the IJmuiden metmast test facilities (North Sea, 85 km offshore of the Netherlands coast). The goal of the campaign was to validate the precommercial FDWL buoy EOLOS TM against the reference IJmuiden metmast, located 300 m apart. The EOLOS TM FDWL was similar to the FDWL used in the PdP campaign but prepared for harsher scenarios and higher endurance. The FDWL buoy hosted a ZephIR TM 300 DWL and a 3DM-GX3-45 IMU as well. The FDWL was configured to sense the wind profile at four different heights (25, 38, 56, and 83 m) coincidental with the measurement heights of the metmast anemometers. An identical ZephIR TM 300 DWL was installed on the metmast platform (21 m) sensing the wind above the metmast top from 90 m up to 315 m every 25 m.

3. Methods

3.1. Basic Definitions

The instantaneous wind vector u ( u = U + u , where U denotes the mean wind and u the turbulent component) is defined as the three-dimensional vector ( u 1 , u 2 , u 3 ) computed here from the H W S , W D , and vertical wind speed ( V W S ) as
u = [ H W S · cos ( W D ) , H W S · sin ( W D ) , V W S ] .
In wind energy, the 10 min mean wind vector is the usual standard, and hereafter, H W S , W D , and V W S refer to the 10 min mean values unless otherwise indicated.
The T I is indicative of the H W S variations with respect to the mean H W S , and it is computed as
T I = σ H W S H W S ,
where σ H W S is the 10 min H W S standard deviation.
Instantaneous error. To study the motion-induced FDWL measurement error in the retrieved H W S at the i-th lidar scan, we defined the “instantaneous” error (in practice, the 1 s error considering a scan rate of 1 scan/s) as
ϵ H W S , i = H W S i ^ H W S i ,
where H W S ^ denotes the H W S “estimated” by the VAD algorithm, and H W S i denotes the “true” H W S at the i-th scan.
10 min error. The wind energy industry is interested in quantifying both the H W S measurement bias and the T I increment due to the wave-induced motion on FDWLs at a 10 min level.
The bias during a sequence of N measurement scans (i.e., the sample statistical interval; usually N = 600 for 10 min measurements when measuring at a single height) is defined as the difference between the mean H W S measured by the FDWL and the mean H W S measured by the reference fixed LiDAR. Formally,
b i a s m e a s = 1 N i = 1 N H W S ^ i 1 N i = 1 N H W S i = 1 N i = 1 N ϵ H W S , i .
Statistically, it can be expressed as
b i a s e s t = E [ H W S ^ ] H W S = E [ ϵ H W S , i ] ,
where E [ · ] denotes the expectation operator, H W S ^ is the 10 min horizontal wind speed estimated by the FDWL, and H W S is the “true” 10 min horizontal wind speed.
The T I increment, Δ T I m e a s , is defined as the T I measured by the FDWL minus the T I measured by the reference fixed LiDAR. Formally,
Δ T I m e a s = σ H W S F D W L H W S ¯ F D W L | b i n σ H W S f i x e d H W S ¯ f i x e d | b i n ,
where the overbar denotes the “mean” over the 10 min H W S samples in the calculus bin.
On the other hand, the T I increment can be estimated as
Δ T I e s t = σ ϵ H W S , i H W S ,
where σ ϵ H W S , i is the standard deviation of the instantaneous error (Equation (3)).

3.2. Reconciling the Estimated and the Measured TI

Next, we related Equation (7) to Equation (6) from a statistical and probabilistic theoretical perspective. Toward this goal and in order to aid the mathematical notation, we introduce the short-hand notation F-R-W as follows: we changed the mnemonics “FDWL” into “F” and “fixed” into “R” (representing “reference” fixed LiDAR). We also introduce the subscript “W” to denote the wind that a nonmoving DWL would measure at the same place as the moving FDWL. In other words, subscript “W” denotes the “true” wind measurement but includes the 1 s average and spatial smoothing inherent to the DWL scanning procedure. To simplify the notation, the “true” wind HWS, H W S W , is simply denoted as H W S .

3.2.1. On the Estimated TI

According to probability theory, the variance in the difference of two random variables (i.e., ϵ F , W = H W S F H W S W , which reproduces Equation (3)) is equal to the sum of each of their variances minus twice their covariance [46],
σ ϵ F , W 2 = σ F 2 + σ W 2 2 ρ σ F σ W ,
where ρ is the correlation coefficient between the “F” and “W” variables, and c o v ( F , W ) = ρ F , W σ F σ W .
At this point, we introduce the assumption that the FDWL-measured and the “true” H W S (subscripts “F” and “W”, respectively) are linearly correlated so that ρ F , W 1 . By experiment, Tiana-Alsina et al. [35] reported correlation coefficients of ρ F , R = 0.9 , ρ F , R = 0.86 , and ρ F , R = 0.66 between the co-located moving and reference DWLs in a controlled-motion experiment for tilt motion amplitudes of 10, 16, and 25 deg, respectively. In past experimental campaigns involving different FDWL buoy topologies [23,25,28,38], the most frequent buoy tilts were reported to range from 0 up to 10 deg. These low tilts are also supportive of high F-to-W correlation coefficients, ρ F , W 1 , at the FDWL location. Inserting ρ F , W = 1 in Equation (8) above and rearranging the terms, we obtain
σ ϵ F , W 2 ( σ F σ W ) 2 .
By taking the square root of Equation (9) above and substituting it into Equation (7), the estimated turbulence intensity can be written as
Δ T I e s t σ F σ W H W S .

3.2.2. On the Measured TI

As in the preceding subsection, one can assume that the 10 min FDWL- and the reference-LiDAR-measured H W S are virtually identical so that H W S F D W L H W S f i x e d H W S [20,29]. Introducing this approximation into Equation (6) and using the shorthand notation F-R-W introduced above, it can be rewritten as
Δ T I m e a s = σ F σ R H W S .
If the FDWL and the fixed LiDAR carried out exactly co-located measurements, the estimated T I (Equation (10)) and the measured T I (Equation (11)) would be identical (i.e., σ R = σ W ). In practice, the fixed and the FDWL were located 50 m apart in the PdP campaign and 300 m apart in the IJmuiden campaign (see Section 2); therefore, small differences up to 1% [33] were expected to arise between the estimated and the measured T I increment when different wind flows were measured by the reference and floating LiDARs.

3.3. FDWL Geometrical Model

We defined the right-handed Cartesian XYZ “moving body” coordinate system of the LiDAR buoy (hereafter, the “moving” coordinate system) and the right-handed Cartesian north–east–down (NED) “fixed” inertial frame of reference (hereafter, the “fixed” system, Figure 1). The IMU (see Section 2) measured the rotation and translation of the “moving” coordinate system with respect to the “fixed” system. The vectors n ^ , e ^ , and d ^ were the unitary vectors aligned with the N, E, and D axes of the fixed coordinate system, respectively. x ^ , y ^ , and z ^ were the unitary vectors aligned with the X, Y, and Z axes, respectively, of the moving coordinate system. h ^ was a unit vector in the outbound direction of the LiDAR scanning cone axis and is defined as the opposite of z ^ . The half-angle cone aperture was θ 0 = 30 deg.
The unitary vector r ^ defined the LiDAR beam pointing direction in each LoS measurement during the LiDAR scan. In what follows, r ^ is given in spherical coordinates relative to the moving coordinate system. Specifically, r ^ is given by the zenith angle relative to h ^ , which is denoted θ 0 (i.e., the cone aperture), and by the azimuth angle relative to x ^ , which is denoted ϕ (Figure 1). As mentioned, θ 0 is the constant value of 30 deg corresponding to the scanning-cone aperture. ϕ is a time-variable angle representing the LiDAR scanning phase (positive towards y ^ ).
During the LiDAR scan, the vector r ^ rotated around h ^ with a uniform circular motion at a rate of 1 revolution per second, i.e., ϕ ( t ) = 2 π × t [rad]. In what follows, for simplicity, we use the variable ϕ instead of the usual time variable t. The variable ϕ ranged from 0 to 2 π [rad] during a LiDAR scan. Therefore, r ^ can be formulated in the moving coordinate system as follows:
r ^ ( ϕ ) = sin ( θ 0 ) · cos ( ϕ ϕ 0 ) x ^ + sin ( θ 0 ) · sin ( ϕ ϕ 0 ) y ^ cos ( θ 0 ) z ^ ,
where ϕ 0 is the LiDAR initial scan phase or, hereafter, the “initial phase”.

3.4. The VAD Algorithm as a First-Order Fourier Series

Focusable CW DWLs use the VAD algorithm to retrieve the wind vector from the measured along-LoS-projected wind velocity during the LiDAR scan, which is computed as the so-called VAD function,
f ( ϕ ) = u · r ^ ( ϕ ) .
In a motionless uniform nonturbulent wind scenario, the wind vector can be decomposed as u = H W S ρ ^ + w h ^ , where ρ ^ = c o s ( W D ) x ^ + s i n ( W D ) y ^ , with W D as positive counterclockwise. Then, the ideal VAD function f ( ϕ ) takes the shape of the model cosine function [16,47],
f m o d e l ( ϕ ) = A · c o s ( ϕ B ) + C ,
where A is the amplitude, B is the phase, and C is the offset term. From Equation (14) above, the H W S , W D , and V W S components of the wind vector are obtained as [28]
H W S = A s i n ( θ 0 ) , W D = B , V W S = C c o s ( θ 0 ) .
In Equation (15) above, the H W S , W D , and V W S refer to the values retrieved at the i-th scan.
The VAD algorithm uses the least squares (LSQ) algorithm to fit the model cosine function f m o d e l (Equation (14)) to the measured function f ( ϕ ) in order to obtain the sought-after parameters A, B, and C:
m i n | A , B , C | | f ( ϕ ) [ A · c o s ( ϕ B ) + C ] | | 2 .
The LSQ fit is usually carried out by means of numerical methods. Yet, they do not allow the explicit analytical formulation of solving coefficients A, B, and C with a view to understanding the propagation of the motion-induced error. A way out of this problem is to treat the error-norm minimization of Equation (16) in the Hilbert space of the square integrable functions L 2 [ 0 , 2 π ] , in which 1 , cos ( m ϕ ) , sin ( n ϕ ) , m , n = 1 form an orthogonal basis. Equivalently, the Fourier series is the LSQ projection of a target function to an orthogonal space defined by the sum of the sine and cosine functions [48]. Therefore, the solution of Equation (16) is the first-order Fourier series of f ( ϕ ) in the basis function 1 , cos ( ϕ ) , sin ( ϕ ) . The first-order Fourier series can be rewritten in the form of Equation (14), in which the A, B, and C parameters can be expressed as a function of the Fourier coefficients as [49]
A = a 1 2 + b 1 2 , B = a r c t a n b 1 a 1 , C = a 0 2 ,
where a 0 , a 1 , and b 1 are the first-order Fourier coefficients of f ( ϕ ) , obtained as [50]
a 0 = 2 P 0 P f ( ϕ ) d ϕ , a 1 = 2 P 0 P f ( ϕ ) cos ( ϕ ) d ϕ , b 1 = 2 P 0 P f ( ϕ ) sin ( ϕ ) d ϕ ,
where P = 2 π .
From Equations (15) and (17) above, it emerges that the Fourier coefficients a 1 and b 1 propagate errors to the H W S and W D (i.e., to the wind horizontal components), while the coefficient a 0 propagates errors only to the retrieved V W S .
As mentioned, homodyne-detection LiDARs, as in the case of the ZephIR 300, are only able to measure the magnitude of the along-LoS radial velocities but not the sign. This is equivalent to saying that f ( ϕ ) takes the shape of a rectified cosine, i.e., the absolute value of f ( ϕ ) [28]. When applying the VAD algorithm to | f ( ϕ ) | , parameters A and C retain the same value as for the heterodyne case (VAD applied to f ( ϕ ) ), but B has a ± π rad ambiguity. The latter can be disambiguated by the wind vane co-located with the LiDAR buoy.

3.5. Estimation Error Methodology

In the absence of buoy motion, the X, Y, and Z axes of the moving coordinate system would coincide with the N, E, and D axes of the fixed system [23]. In the presence of waves, the FDWL experiences angular and translational motion around the three NED axes. In this section, we formulate analytical expressions for the motion influence on the VAD-retrieved wind vector.
In order to simplify the obtained analytical expressions, the following assumptions are considered: (i) the motion in each of the six DoF is considered to be a zero-mean simple harmonic motion; (ii) the rotational and translational motions are considered to be independent error sources; and (iii) under motion, the V W S contribution to the H W S error is considered to be null as compared to the H W S contribution [51].

3.5.1. Rotational Motion Model

The FDWL buoy undergoes roll, pitch, and yaw motions around the three fixed coordinate axes, i.e., N, E, and D, respectively. As observed experimentally [23], roll and pitch motions can be modeled using a first-order small-angle approximation as simple harmonic motions. The yaw motion, which shows larger and more constant values with time, can be modeled by the constant tilt.
We define roll and pitch tilts (in units of radians) as sinusoidal signals of the form
Ω x ( ϕ ) = A x · sin ( f x ϕ α x ) ,
where A x , f x , and α x denote the amplitude, frequency, and motional phase, respectively. Subindex x = r o l l , p i t c h denotes the r o l l and p i t c h motion, respectively. Note that, for convenience, Equation (19) is expressed as a function of the LiDAR scan-phase variable ϕ ( ϕ = 2 π t ) instead of time variable t.
As mentioned, the yaw tilt is considered a constant value Ω y a w ( ϕ ) = A y a w . This is conducted without limiting the general formulation of the problem. Thus, the reader could extend the formulation of the yaw angle to the harmonic case by defining the yaw motion as a sinusoidal signal (Equation (19)) plus an offset term equal to its mean value during the interval under study, and by following the formulas in Section 3.5.1 and Appendix A.1.
First, in order to compute the rotated LiDAR pointing direction, r ^ r o t , in the fixed NED coordinate system given the roll, pitch, and yaw rotations, we used the Euler rotation matrix. Euler’s rotation theorem states that any rigid-body rotation can be defined by three rotation angles [52]. There are many different conventions regarding the rotation angles and their order of application. Here, we considered the NED convention, i.e., three elemental rotations were carried out sequentially: first, a rotation around the N axis (roll motion); second, a rotation around the E axis (pitch motion); and finally, a rotation around the D axis (yaw motion). The rotations were positive counterclockwise. Therefore, the Euler rotation matrix can be written as [53]
R = R D · R E · R N ,
where R N , R E , and R D are the elemental rotation matrices around the N, E, and D axes, respectively. They can be formulated as follows:
R D = cos ( Ω y ) sin ( Ω y ) 0 sin ( Ω y ) cos ( Ω y ) 0 0 0 1 , R E = cos ( Ω p ( ϕ ) ) 0 sin ( Ω p ( ϕ ) ) 0 1 0 sin ( Ω p ( ϕ ) ) 0 cos ( Ω p ( ϕ ) ) , R N = 1 0 0 0 cos ( Ω r ( ϕ ) ) sin ( Ω r ( ϕ ) ) 0 sin ( Ω r ( ϕ ) ) cos ( Ω r ( ϕ ) ) ,
where subindexes y, p, and r stand for the “yaw”, “pitch”, and “roll”, respectively.
The small-angle approximation for roll and pitch angles [23,28] translates into a first-order Taylor approximation for the sine and cosine functions (i.e., s i n ( x ) x and c o s ( x ) 1 ). These first-order approximations were incorporated into Equation (A1), Appendix A.1, in order to simplify the rotation matrix R .
The second step was to compute the rotated version of the LiDAR pointing vector, r ^ r o t ( ϕ ) , in the NED fixed coordinate system:
r ^ r o t ( ϕ ) = R · r ^ ( ϕ ) .
Third, the FDWL-measured wind under rotational motion was computed by projecting the wind vector over the rotated LiDAR pointing vector for each LoS. This was expressed by means of the dot product as follows:
f r o t ( ϕ ) = u · r ^ r o t ( ϕ ) .
Eventually, the motion-corrupted FDWL wind vector was derived by applying the VAD algorithm to the function f r o t ( ϕ ) above. With this aim, the first-order Fourier coefficients of f r o t ( ϕ ) (i.e., a 0 r o t , a 1 r o t , and b 1 r o t ) were computed through Equation (18) by substituting f r o t ( ϕ ) in place of f ( ϕ ) . The H W S retrieval error was computed as the difference between the “true” H W S (i.e., without motion influence) and the FDWL-retrieved H W S under rotational motion:
ϵ H W S , i r o t = H W S , i 1 sin ( θ 0 ) ( a 1 r o t ) 2 + ( b 1 r o t ) 2 ,
where a 1 r o t and b 1 r o t are the Fourier coefficients of f r o t ( ϕ ) . See Appendix A.3 for the mathematical results in expanded form.

3.5.2. Translational-Motion Model

Waves also induce translational motion to the FDWL in the N, E, and D directions, i.e., surge, sway, and heave motions, respectively. Similar to the derivation of the rotational motion components in Equation (19), sinusoidal variation was assumed for each translational-motion component. They were formulated as follows:
t x ( ϕ ) = A x · sin ( f x ϕ α x ) ,
where A x , f x , and α x , are the amplitude, frequency, and phase of the translational motion x, where x denotes the surge, sway, or heave components. We also defined the translational-velocity vector t ( ϕ ) in the NED fixed coordinate system as the three-component vector,
t ( ϕ ) = [ t s u ( ϕ ) , t s w ( ϕ ) , t h e ( ϕ ) ] ,
where the subscripts s u , s w , and h e referred to the surge, sway, and heave, respectively.
First, to compute the FDWL-measured wind under translational motion, we computed the apparent wind vector. The apparent wind vector measured by the FDWL was the difference between the wind vector and the translational-velocity vector:
u t r a n s ( ϕ ) = u t ( ϕ ) .
Second, the translational-velocity contribution to the LiDAR-measured LoS, f t r a n s ( ϕ ) , was obtained by projecting the apparent wind vector, u t r a n s , on the LiDAR pointing direction r ^ ( ϕ ) as follows:
f t r a n s ( ϕ ) = u t r a n s ( ϕ ) · r ^ ( ϕ ) .
The LiDAR-retrieved wind vector was the result of applying the VAD algorithm over the motion-corrupted LoS, f t r a n s ( ϕ ) . The Fourier coefficients a 0 t r a n s , a 1 t r a n s , and b 1 t r a n s were obtained through Equation (18) by substituting f ( ϕ ) by f t r a n s ( ϕ ) . The Fourier coefficients are given in expanded form in Appendix A.4 of Appendix A.
Finally, the H W S measurement error due to the FDWL’s translational motion ϵ H W S , i t r a n s became
ϵ H W S , i t r a n s = H W S i 1 sin ( θ 0 ) ( a 1 t r a n s ) 2 + ( b 1 t r a n s ) 2 .

3.5.3. Total Error Model

In order to estimate the total error, we departed from the assumption that the rotational and translational motion were independent error sources (see Section 3.5). Therefore, the total error in the retrieved H W S in a LiDAR scan is the superposition of errors:
ϵ H W S , i ( u , ϕ 0 ) = ϵ H W S , i r o t ( u , ϕ 0 ) + ϵ H W S , i t r a n s ( u , ϕ 0 ) .
In Equation (30), the total error, ϵ H W S , i ( u , ϕ 0 ) , is parameterized by the instantaneous wind vector u and initial phase ϕ 0 .
Figure 2a illustrates the estimated H W S total error in a LiDAR scan parameterized as a function of the wind direction ( W D ranging from 0 to 360 deg) and the initial phase ϕ 0 in a roll-only motion scenario ( H W S = 10 m/s, V W S = 0 m/s). As can be observed, the H W S error was highly dependent on both input parameters (WD and initial phase), ranging from −1.5 up to 1.5 m/s (the H W S errors for other input H W S s are scalable through simple direct proportionality). If the W D is known, the H W S error in a LiDAR scan is dependent only on the initial phase. This is represented in Figure 2b, where the H W S error in the LiDAR scan in Figure 2a is represented as a function of the initial phase ϕ 0 for the constraint settings H W S = 10 m/s, W D = 275 deg, and V W S = 0 m/s.
However, the initial phase is the manufacturer’s undisclosed parameter for the ZephIR 300 LiDAR. As the instantaneous value of the initial phase at each LiDAR scan cannot be known nor, consequently, the total error, we assumed the initial phase was a random variable with a uniform distribution between 0 and 2 π rad [25]. As a result, H W S error moments rather than instantaneous values were derived.
The 10 min H W S bias (Equation (5)) was computed as the expectation (first raw moment) of the H W S total error function (Equation (30)), ϵ H W S , i ( u , ϕ 0 ) , with respect to the initial phase, ϕ 0 , constrained to the mean wind vector U as follows:
b i a s e s t = E ϕ 0 [ ϵ H W S , i ( U , ϕ 0 ) ] .
In practice, the b i a s e s t is simply the mean value of ϵ H W S , i ( U , ϕ 0 ) computed as b i a s e s t = 1 N i = 1 N ϵ H W S , i ( U , ϕ 0 , i ) , where the discrete random variable ϕ 0 is defined by a finite list { ϕ 0 , 1 , , ϕ 0 , N } of equally likely outcomes (vertical line in Figure 2a).
Analogously, the 10 min TI increment of Equation (7) is estimated as follows:
Δ T I e s t = S T D ϕ 0 [ ϵ H W S , i ( U , ϕ 0 ) ] H W S ^ ,
where S T D stands for the standard deviation, which is the square root of the variance of the total error ( σ ϵ H W S , i ), and H W S ^ is the 10 min mean H W S measured by the FDWL.

3.6. Bias and TI-Increment Estimation Procedure

Figure 3 summarizes the procedure described in Section 3.5 to compute the estimated 10 min bias and T I increment. The 10 min mean wind vector ( U ), as well as the characteristic amplitude, frequency, and phase of each of the 6 DoF, were the inputs. In practice, the yaw frequency and phase were set to zero. The Fourier coefficients a 1 r o t and b 1 r o t were obtained through Equations (A3) and (A4), respectively, which estimated the H W S measurement error due to the rotational motion ( ϵ H W S , i r o t , Equation (24)). Similarly, the Fourier coefficients a 1 t r a n s and b 1 t r a n s were obtained through Equations (A5) and (A6), respectively, which estimated the H W S measurement error due to the translational motion ( ϵ H W S , i t r a n s , Equation (29)). The following step was to sum the rotational- and translational-error expressions to yield the total error as a function of the wind vector and initial scan phase (Equation (30)). Then, it was constrained to the 10 min mean wind vector, U , by substituting the u in Equation (30) with U . Finally, the bias and T I increment were estimated by means of Equations (31) and (32), respectively.
Alternatively, for a specific 1 s error model validation (refer to Section 4.1), the total error (Equation (30)) was obtained as a function of the instantaneous wind vector and the LiDAR initial scan phase.

3.7. Sinusoidal Characterization of the Measured Motion Time Series

In order to estimate the FDWL motion-induced error with the presented methodology based on Equations (19) and (25), there is a requirement that the motion time series in each of the DoF be modeled as a sinusoidal signal with a characteristic amplitude, frequency, and phase (A, f, and α , respectively). The estimation procedure for these parameters is described next through a case example.
Figure 4a shows a record of the IMU-measured roll time series during the PdP campaign. As can be noticed, the signal was composed of multiple frequency components. To analyze the signal spectrum, the power spectral density (PSD) was computed in batches of 10 min segments using the Blackman–Tukey method [54]. In Figure 4b, it emerges that the signal had a broad spectrum with relevant frequency components from approximately 0.25 to 0.5 Hz. The characteristic frequency f was computed as the frequency corresponding to the peak PSD maximum (≃0.3 Hz).
In order to estimate the characteristic amplitude A, we made use of the fact that the amplitude of a sinusoidal signal is related to its power, P, through the relationship P = A 2 / 2 [54]. Moreover, the mean power P ¯ of a stochastic signal s ( t ) (i.e., the IMU-measured motion) with duration T is computed as [54]
P ¯ = 0 T | s ( t ) | 2 d t .
The mean power, P ¯ , takes into account all the frequency components of the signal spectrum [55]. Therefore, the characteristic amplitude was estimated from the mean signal power as A = 2 P ¯ . Finally, the characteristic phase α was obtained from the first-order Fourier decomposition of s ( t ) in a similar fashion to the calculation of phase term B in Equation (17), given the first-order coefficients a 1 and b 1 in Equation (18) (see Section 3.4).
Figure 4a shows the IMU-measured roll time series against the “fitted” sinusoidal signal with the characteristic amplitude, frequency, and phase of A = 1.3 deg, f = 0.3 Hz, and α = 1.1 deg, respectively. It can be observed that the fitted sinusoidal signal reproduced the measured time series under a first-harmonic approximation with reasonable accuracy.

3.8. A Note on Appendix A and the Supplementary Materials Math Formulations

Provided in Appendix A is a formulation compendium to compute the 6-DoF-motion-induced error (Section 3.5.3). A MATLAB TM R2020a code to compute the total error (Section 3.5.3) is also included as part of the Supplementary Materials of this manuscript.
The mathematical derivations of this Appendix A were checked both algebraically and numerically. Algebraically, the manually derived math expressions were validated using the MATLAB TM R2020a symbolic toolbox. Numerically, these expressions were coded in MATLAB and compared with the outputs of the FDWL simulator [35] (Figure 5, Figure 6 and Figure 7) for quality assurance.

4. Results and Discussion

4.1. Error Model Validation

In order to validate the presented analytical-error formulation (Section 3.5 and Appendix A; in what follows, “the unified-error formulation”), the estimated total error in the FDWL-retrieved H W S (Equation (30)) was compared with the LiDAR motion simulator as a reference [37]. The study was carried out by inputting the same reference wind vector ( H W S = 10 m/s and V W S = 0 m/s) into both the unified formulation and the simulator under three different motion scenarios: (i) rotational motion only, (ii) translational motion only, and (iii) 6 DoF motion. In all cases, f x = 0.3 Hz and α x = 0 deg, x = r o l l , p i t c h were used. In each scenario, the performance of the analytically estimated error was evaluated as a function of the W D and initial scan phase ϕ 0 .
In roll-only scenario #1 (Figure 5), the analytically estimated error (Figure 5a, Equation (24)) was virtually identical to the simulator-estimated error (Figure 5b) over the whole span of the W D s and the initial-phase values. From Figure 5c, all the error differences between panels (a) and (b) were lower than 0.3 m/s, which is evidence of the satisfactory accuracy attained using the first-order approximation in the formulation. The R M S E between the error estimates of Figure 5a,b (equivalently, the square root of the mean of the squares of the error values shown in Figure 5c) was 0.04 m/s.
When addressing the translational-only scenario #2 (Figure 6), null differences (0 m/s) were found in Figure 6c between the analytically estimated and the simulator-estimated error. The exact 0 m/s error difference was explained by the fact that no approximations were used when deriving the analytical expressions for the translational-motion component.
Opposite to other floating remote-sensing devices such as high-frequency surface-wave radars, which measure the horizontal propagation of waves and thus are not sensitive to heave-induced Doppler effect [56,57], the ZephIR TM 300 FDWL is sensitive to all three translational motion components due to its inherent conical scanning pattern [16].
Finally, the successful performance of the unified formulation presented for the analytically estimated error was re-encountered in Figure 7, where both the translational and rotational motion were combined into the 6 DoF motion scenario #3. Despite the greater complexity of this scenario, the error differences (Figure 7c) remained below 0.7 m/s. The R M S E was 0.22 m/s.

4.2. Experimental Results

The unified-error formulation (Section 3.5 and Appendix A) was also validated by comparing the analytically estimated H W S mean bias and T I increment (Equations (31) and (32), respectively) against their experimental values (Equations (4) and (6), respectively). The experimental values were computed from the 10 min errors between the FDWL and the reference fixed LiDAR used in both the nearshore (“Pont del Petroli”) and open-sea (IJmuiden) scenarios. For the IJmuiden campaign, the measurements of the FDWL at 83 m and the measurements of the metmast DWL at 92 m were considered. Towards this purpose, 848 data records (6 to 30 June 2013) from PdP campaign and 3893 data records (1 April to 7 May 2015) from the IJmuiden campaign were used (see Section 2). The data were filtered for quality assurance following the outlier-rejection criteria [33]: the measured H W S values outside the 1–80 m/s range, the rain-flagged data, the spatial variation (SV) values higher than 0.2, and the backscatter coefficients lower than 0.02 were removed. The SV is an indicator of the goodness of the VAD fitting of the measured LoS in a LiDAR scan [58]. The backscatter coefficient is indicative of the power of the received LiDAR echo.

4.2.1. The Performance according to the Estimation of the H W S Bias

In the PdP campaign, using statistical analysis, the measured 10 min H W S bias (Equation (4)) showed virtually nil values in accordance with the previous results in nearshore locations in the state of the art [27,28]. The median of the measured 10 min H W S bias was −0.02%, with 25th and 75th percentiles of −0.89% and 0.76%, respectively. On the other hand, the analytically estimated 10 min H W S bias (Equation (5)) showed values of the same order of magnitude as the measured ones, thus validating the method. Thus, the median of the estimated bias was 0.06%, and the 25th and 75th percentiles were −0.66%, and 0.64%, respectively. The median bias was in accordance with the results obtained in other state-of-the-art studies [38].
In the Ijmuiden campaign, the measured 10 min H W S showed higher values than the observed ones in the “Pont del Petroli” campaign: the median of the measured 10 min H W S bias was −0.53%, with 25th and 75th percentiles of −2.30% and 1.08%, respectively (and notice the asymmetrical distribution of these percentiles around 0%). These values were in accordance with the previous open-sea measurement campaigns [32,59]. This is due to the fact that in an open-sea scenario there is stronger wave motion, which induces a higher error in the FDWL measurements. The abovementioned asymmetry on the measured biases accounts for the wind shear between the fixed DWL and FDWL measurement heights (92 and 83 m, respectively). The analytical model estimated biases of the same order of magnitude as the measured ones with a median of 0.05%, and 25th and 75th percentiles of −1.48% and 1.47%, respectively, which further validated the analytical model.

4.2.2. Performance according to the Estimation of the T I Increment (I): Case Example

As part of the validation study, in the next case example, the analytically estimated 10 min T I increment, Δ T I e s t (Equation (32)), was compared with the measured one, Δ T I m e a s (Equation (7)). The different motion conditions studied are described as a function of the mean tilt amplitude and mean translational-velocity amplitude of the LiDAR buoy.
The mean tilt amplitude was computed as the 10 min mean of the quadratic sum of roll and pitch tilts (simple harmonic motion, Equation (19)) [28]. It can be formulated as follows:
A t i l t ¯ = k = 1 N Ω r o l l ( k ) 2 + Ω p i t c h ( k ) 2 N ,
where Ω x , x = r o l l p i t c h is the IMU-measured roll/pitch tilt, N = 6000 is the number of samples in a 10 min interval at a 10 Hz IMU-sampling frequency, and k is a reminder of the discrete time t k .
The mean translational-velocity amplitude was computed as the 10 min mean norm of the FDWL translational-velocity vector,
A v t r a n s ¯ = k = 1 N | | t ( k ) | | N ,
where t ( k ) is the translational-velocity vector defined by Equation (26).
Figure 8a compares the estimated T I -increment time series ( Δ T I e s t ) with the measured one ( Δ T I m e a s ) in the context of the PdP campaign. The motion-induced error on the FDWL-measured T I manifested as T I increments between 0.5 and 5%. The highest T I increments occurred during the daytime (from 10 a.m. to 4 p.m., approximately), corresponding to time periods with high HWSs. These periods were also related to the high wave motion [27,60]. This was corroborated by Figure 8b, where the mean tilt amplitude and mean translational-velocity amplitude of the floating LiDAR are depicted as time series. It can be observed that the high TI-increment values in Figure 8a were linked to high-motion scenarios, i.e., to the time periods with high tilt and velocity amplitudes in Figure 8b.
As shown in Figure 8a, the estimated T I increment exhibited similar values to the measured ones over the period under study. The 1 h averaged time series further showed the goodness of the TI-increment estimates. However, a small underestimation of the peak values occurred. Thus, while the measured T I increment peaked at 5% close to noon, the estimated one only reached 3.5%. Larger differences were observed in one-to-one comparisons of the 10 min estimates with the 10 min measurements. A suitable explanation for this is that in the unified formulation presented in Section 3.5.3, the wind flow during 10 min periods was assumed to be uniform, and only the motion influence was considered. As the LiDARs were located 50 m apart, wind scenarios with high spatial variability (e.g., turbulent winds coming from the urban area, Section 2) usually led to different instantaneous measurements between the floating and the reference LiDAR. Therefore, the correlation coefficients ρ F , R < 1 (Section 3.2.1) may limit the validity of the estimates, Δ T I e s t . This effect may also occur in low-motion scenarios [33].

4.3. The Performance of the Estimation of the T I Increment (II): Statistical Analysis

The performance of the analytical-error formulation when estimating the 10 min T I increment was also assessed with reference to the FDWL-measured T I increment under nearshore and open-sea motional scenarios (the PdP and IJmuiden campaigns, respectively). The statistical sample consisted of 848 data records (from 6 to 30 June 2013) for the PdP campaign and 3893 data records (from 1 April to 7 May of 2015) for the IJmuiden campaign. The estimated and the measured T I increments ( Δ T I e s t and Δ T I m e a s , respectively) for the whole campaign were clustered as a function of the wind intensity, using the HWS as the binning variable, and as a function of the type of motion (i.e., rotational or translational), using the mean tilt amplitude and the mean translational-velocity amplitude as the respective binning variables.
The estimation accuracy was quantitatively assessed according to the root mean square error ( R M S E ), which is formulated as follows:
R M S E = p = 1 N [ Δ T I e s t ( p ) Δ T I m e a s ( p ) ] 2 N ,
where N is the number of samples in the bin under study, Δ T I e s t ( p ) is the p-th T I -increment estimate in the bin, and Δ T I m e a s ( p ) is the corresponding p-th measurement.
Figure 9 depicts the statistical results for the PdP (left panels) and IJmuiden (right panels) campaigns. In all panels, it can be observed that the T I increment, Δ T I , was always positive, evidencing the motion-induced additive turbulence [26,28]. Moreover, it increased with increasing wave motion, i.e., with increasing mean tilt amplitude (panels (a), (b)) and mean translational-velocity amplitude (panels (c), (d)). Δ T I ranged from median values of Δ T I 1 % in low-motion scenarios to Δ T I 2 % in high-motion ones. The green bars illustrating the statistically significant bins (i.e., bins containing greater than or equal to 5% of the population of the PdP campaign, i.e., ≃50 samples) showed mean tilt amplitudes between 1 and 3.5 deg for the PdP campaign and between 1 and 5.5 deg for the IJmuiden campaign (notice the different X-axis ranges used). The observed mean translational-velocity amplitudes ranged between 0.1 and 0.2 m/s for the PdP campaign and between 0.1 and 0.65 m/s for the IJmuiden campaign (notice the different X-axis ranges used). The H W S ranged between 2 and 8 m/s in the PdP campaign and between 2 and 22 m/s in the IJmuiden campaign (notice the different X-axis ranges used). The higher motion and wind magnitudes evidenced the harsher scenario experienced by the FDWL during the IJmuiden campaign.
On one hand, Figure 9a,b plot the T I increment as a function of the mean tilt amplitude in 0.5 deg bins, which was representative of the rotational-motion component. The median of the estimated Δ T I values virtually matched that of the measurements, thus validating the overall analytical-error formulation. In the PdP campaign (Figure 9a), the error bars showed a larger dispersion for the measured values ( Δ T I m e a s ) than for the estimated ones ( Δ T I e s t ). This is attributable to the urban topology surrounding the experimental area, which creates a high spatial variation in the wind field. In contrast, this was not found in IJmuiden. Thus, Figure 9b,d,f shows estimated error bars virtually coincidental with the measured ones due to the more homogeneous wind fields found over the ocean.
On the other hand, Figure 9c,d represent the T I increment as a function of the FDWL mean translational velocity in 0.025 m/s bins for the PdP campaign and 0.05 m/s bins for the IJmuiden campaign, respectively. Similar to Figure 9a,b, the median of the estimated Δ T I values almost ideally matched the measurements in all the significant bins (dark green), Figure 9c,d. The largest difference was Δ T I m e a s Δ T I e s t = 0.26 % for the 0.175 m/s bin. Similar comments apply to those for Figure 9a,b regarding the measurement error bars.
Finally, Figure 9e,f depict the T I increment as a function of the mean H W S in 1 m/s bins for the PdP campaign and 2 m/s for the IJmuiden campaign, respectively. Similar concurrent results were observed, with the estimated Δ T I median values matching the measurements for all H W S bins except for the 2 m/s one. This was due to the fact that the ZephIR TM 300 H W S measurements below 3 m/s were tagged as “unreliable”, according to the manufacturer’s specs for the LiDAR [40]. Similar discussion comments apply for the error bars shown in Figure 9e,f.
The R M S E between the estimated and the measured Δ T I (Equation (36)) is depicted in red for each of the tilt, velocity, and H W S bins. Similar values were obtained in all the statistically significant bins, ranging from R M S E = 0.004 (0.4%) up to R M S E = 0.013 (1.3%), which was comparable to (therefore, consistent with) the size of the measurement error bars obtained for Δ T I m e a s . As can be observed, the R M S E was slightly correlated with the tilt, so that the larger the buoy tilt, the larger the R M S E in rough approximation, which in turn accounted for the first-order approximation used in Equation (A1) (refer to Appendix A.1).

5. Conclusions

A unified analytical formulation for the computation of the 6-DoF-motion-induced error in focusable CW FDWLs was presented (Figure 3). The total error in the retrieved H W S was computed as the superposition of the rotational- and translational-motion errors. The formulation proved to be capable of estimating the H W S bias (Section 4.2.1) and T I increment (Figure 8 and Figure 9).
The analytical model departed from a thorough formulation of the FDWL buoy geometry and dynamics in order to derive the rotational- and translational-motion influence on the LiDAR-measured LoS measurements in a CW DWL scan. The well-known VAD algorithm was computed as a first-order Fourier series, which allowed derivation of the sought-after analytical expressions describing the FDWL-measured wind-vector error as a function of the buoy attitude and the “true” wind vector. The first-order approximation was retained when computing the buoy rotational matrix. The assumption of a uniformly distributed random initial scan phase permitted the estimation of the H W S measurement bias and R M S E over 10 min time intervals.
The wave-induced motion in each of the six DoF was modeled as a simple harmonic motion with a characteristic amplitude, frequency, and phase recomputed every 10 min for each DoF. A method to estimate these variables based on spectral analysis of the LiDAR-buoy motional time series was presented in Section 3.7.
The proposed formulation was numerically validated by comparing the H W S error figures obtained as a function of the W D and the initial scan phase with those output by the FDWL motion simulator [37]. A virtually perfect match between both methods was found except for the second-order approximation errors.
The methodology was further validated using experimental data from the “Pont del Petroli” (Badalona, Barcelona) and IJmuiden (North sea, 85 km offshore the Netherlands) campaigns, in which an FDWL and a reference fixed LiDAR were used. The estimated 10 min bias using the analytical formulation presented here yielded similar results to the measured one. Thus, for the PdP campaign, the median values obtained for the estimated and the experimental H W S bias were −0.02% and 0.06%, respectively, and the 25th and 75th percentiles were −0.66% (estimated) and −0.89% (measured) and 0.76% (measured) and 0.64% (estimated), respectively. For the IJmuiden campaign, median values of −0.53% (measured) and 0.05% (estimated) were found, and the 25th and 75th percentiles were −2.30% (measured) and −1.48% (estimated) and 1.08% (measured) and 1.47% (estimated), respectively. The inherent asymmetry between the positive and negative biases was due to the wind shear between the fixed DWL and FDWL measurement heights. The slight differences between the measured and estimated values were inconsequential and were due to the harsher sea conditions at the IJmuiden site. With regard to the estimation of the T I increment, the analytical estimates matched the measured values in statistical terms under all motion and wind scenarios for both the PdP and IJmuiden campaigns (Figure 9). Thus, the root mean square differences between the medians of Δ T I e s t and Δ T I m e a s were 0.16% (PdP) and 0.11% (IJmuiden) for the rotational motion (Figure 9a,b), 0.13% and 0.16% for the translational motion (Figure 9c,d), and 0.14% and 0.21% for the H W S (Figure 9e,f), giving mean values of about 0.015 (1.5%). The experimental validation was limited by the fact that the compared LiDARs were a minimum of 50 m apart and subject to wind flows that were not always uniform. This caused moderate dispersion in the turbulence values measured by the two LiDARs, which deteriorated in one-to-one comparison of the 10 min T I .
All in all, the unified formulation presented here proved to be a straightforward and accurate tool for quantitatively assessing the 6-DoF-motion-induced error in focusable CW FDWLs in terms of both the H W S -bias and T I -increment estimation, in addition to the real one. This information is increasingly being recognized as critical to the acceptance of FDWLs and ship-borne LiDARs in the Wind Energy industry as compared to fixed LiDARs and reference metmasts [61,62]. To the best of our knowledge, the proposed unified analytical formulation is the first in the state of the art able to estimate the T I increment at a 10-min basis from solely FDWL measurements under virtually all motional and wind scenarios. In addition, the novel analytical expressions obtained are less computationally demanding than numerical simulators and unveil a key understanding of the involved FDWL attitude and LiDAR-retrieval-algorithm error sources. These findings have important implications for the accurate measurement of wind speed and turbulence in FDWLs, particularly when they are being deployed as standalone measurement instruments. M a t l a b TM R2020a codes are also provided to the reader as part of the Supplementary Materials of this manuscript. Further steps could include the study of multimodal motion as well as the wind flow variability in a LiDAR scan. Finally, experimental validation would benefit from tests at different LiDAR sounding heights.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/rs15061478/s1.

Author Contributions

Derivation of analytical expressions, J.F.-G. and A.S.-B.; software development, A.S.-B. and M.P.A.d.S.; analyses and figures, A.S.-B. and J.F.-G.; writing—original draft preparation, A.S.-B., J.F.-G. and M.P.A.d.S.; review and editing, F.R.; funding acquisition, F.R.; conceptualization support, F.R. This work was developed as part of A.S.-B.’s doctoral thesis, supervised by F.R. and J.F.-G.’s master thesis co-advised by F.R. and A.S.-B. All authors have read and agreed to the published version of the manuscript.

Funding

This research project was part of the project PID2021-126436OB-C21 funded by the Ministerio de Ciencia e Investigación (MCIN)/Agencia Estatal de Investigación (AEI)/10.13039/501100011033 y FEDER “Una manera de hacer Europa”. The work of A. Salcedo-Bosch was supported by grant 2020 FISDU 00455 funded by Generalitat de Catalunya—AGAUR. The work of M.P Araújo da Silva was supported under Grant PRE2018-086054 funded by MCIN/AEI/10.13039/501100011033 and FSE “El FSE invierte en tu futuro”. The European Commission collaborated under projects H2020 ATMO-ACCESS (GA-101008004) and H2020 ACTRIS-IMP (GA-871115). The European Institute of Innovation and Technology (EIT), KIC InnoEnergy project NEPTUNE (call FP7), supported the measurement campaigns.

Data Availability Statement

Data are available from the authors.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CWcontinuous wave
DoFdegrees of freedom
DWLDoppler wind LiDAR
FDWLfloating Doppler wind LiDAR
HWShorizontal wind speed
IMUinertial measurement unit
LDlinear dichroism
LoSline of sight
LSQleast squares
NEDnorth–east–down
metmastmeteorological mast
MDPIMultidisciplinary Digital Publishing Institute
PdPPont del Petroli
PSDpower spectral density
RMSEroot mean square error
SVspatial variation
TIturbulence intensity
VADvelocity–azimuth display
VWSvertical wind speed
WDwind direction

Appendix A. Formulation Compendium

This Appendix provides exhaustive details for the formulation of the 6-DoF motion-induced error. Specifically, Appendix A.1 addresses the first-order approximation of the Euler rotation matrix, R , in connection with Equations (20) and (21). Appendix A.2 provides the calculation of the LoS-projected wind vector under rotational motion to obtain the VAD function f r o t ( ϕ ) . The first-order Fourier coefficients under rotational and translational motion are formulated in Appendix A.3 and Appendix A.4, respectively, by using the auxiliary integrals given in Appendix A.5. The Supplementary Materials provide the MATLAB TM R2020a code to compute the total error (Section 3.5.3).
Refer to Section 4.1 and Figure 5, Figure 6 and Figure 7 for quality assurance.

Appendix A.1. First-Order Approximation of the Rotation Matrix

The Euler rotation matrix, R , is the result of elemental rotations about the N, E, and D axes, R N , R E , and R D , respectively, representing the roll, pitch, and yaw rotations. Inserting Equations (21) into (20) and retaining the first-order terms yields:
R cos ( Ω y ) sin ( Ω y ) sin ( Ω y ) · Ω r ( ϕ ) + cos ( Ω y ) · Ω p ( ϕ ) sin ( Ω y ) cos ( Ω y ) cos ( Ω y ) · Ω r ( ϕ ) + sin ( Ω y ) · Ω p ( ϕ ) Ω p ( ϕ ) Ω r ( ϕ ) 1 ,
where Ω x , x = r , p , y is the tilt angle already defined in Equation (19). r , p , y is shorthand notation for “roll”, “pitch”, and “yaw”.

Appendix A.2. Wind-Vector Projection over the Rotated LiDAR Pointing Vector r ^ r o t (ϕ)

In rotational motion, the FDWL-measured wind as a function of the scan phase, f r o t ( ϕ ) , is the projection of the wind vector u on the LiDAR pointing vector r ^ r o t ( ϕ ) . Inserting Equations (12) and (A1) into Equation (22) gives r ^ r o t ( ϕ ) , and substituting this vector and definitional Equation (1) for the wind vector into Equation (23) gives the VAD function for the measured wind. The result can be expressed as follows:
f r o t ( ϕ ) = u · r ^ r o t ( ϕ ) = H W S · [ cos ( W D ) · ( sin ( θ 0 ) · [ cos ( ϕ ϕ 0 ) · r 11 + sin ( ϕ ϕ 0 ) · r 12 ] cos ( θ 0 ) · r 13 ) + sin ( W D ) · ( sin ( θ 0 ) · [ cos ( ϕ ϕ 0 ) · r 21 + sin ( ϕ ϕ 0 ) · r 22 ] cos ( θ 0 ) · r 23 ) ] + V W S · ( sin ( θ 0 ) · [ cos ( ϕ ϕ 0 ) · r 31 + sin ( ϕ ϕ 0 ) · r 32 ] cos ( θ 0 ) · r 33 ) ,
where coefficients r i j , i , j = 1 3 are the entries in the i-th row and j-th column of the rotation matrix R (Equation (A1)).

Appendix A.3. Fourier Coefficients for the Rotational Motion Model

In the rotational motion of the FDWL (Section 3.5.1), the first-order Fourier coefficients a 1 r o t and b 1 r o t of Equation (24) that give the H W S rotational error are obtained in analytical form by inserting Equation (A2) above into definitional Equation (18). They take the following form:
a 1 r o t = 1 π 0 2 π f r o t ( ϕ ) cos ( ϕ ) d ϕ = 1 π H W S · [ cos ( W D ) · ( sin ( θ 0 ) · [ cos ( y ) I 11 sin ( y ) I 12 ] cos ( θ 0 ) · [ A r sin ( y ) I 1 , r + A p cos ( y ) I 1 , p ] ) + sin ( W D ) · ( sin ( θ 0 ) · [ sin ( y ) I 11 + cos ( y ) I 12 ] cos ( θ 0 ) · [ A r cos ( y ) I 1 , r + A p sin ( y ) I 1 , p ] ) ] ,
and
b 1 r o t = 1 π 0 2 π f r o t ( ϕ ) sin ( ϕ ) d ϕ = 1 π H W S · [ cos ( W D ) · ( sin ( θ 0 ) · [ cos ( y ) I 21 sin ( y ) I 22 ] cos ( θ 0 ) · [ A r sin ( y ) I 2 , r + A p cos ( y ) I 2 , p ] ) + sin ( W D ) · ( sin ( θ 0 ) · [ sin ( y ) I 21 + cos ( y ) I 22 ] cos ( θ 0 ) · [ A r cos ( y ) I 2 , r + A p sin ( y ) I 2 , p ] ) ] ,
where coefficients I m n and I m , x with m , n = 1 , 2 and x = r , p are the auxiliary integrals of Table A1. The mnemonic subscript m = 1 ( m = 2 ) refers to the Fourier coefficient a 1 ( b 1 ) , and as mentioned, subscripts r and p denote “roll” and “pitch”, respectively.

Appendix A.4. Fourier Coefficients for the Translational-Motion Model

In the translational motion, the FDWL-measured wind as a function of the scan phase, f t r a n s ( ϕ ) (Equation (28)), is the projection of the apparent wind vector (Equation (27)) on the LiDAR pointing vector (Equation (12)). The first step towards computing the Fourier coefficients is, therefore, computation of the apparent wind by substituting translational-velocity vector Equation (26) and definitional Equation (1) for the wind vector into Equation (27). The second step is substituting Equation (27) for the measured wind and Equation (12) for the LiDAR pointing vector into Equation (28), which gives the VAD function for the measured wind, f t r a n s ( ϕ ) .
Finally, the first-order Fourier coefficients a 1 t r a n s and b 1 t r a n s in Equation (29), giving the H W S translational error, are derived in analytical form by inserting Equation (28) into definitional Equation (18). This gives
a 1 t r a n s = 1 π 0 2 π f t r a n s ( ϕ ) cos ( ϕ ) d ϕ = 1 π · ( sin ( θ 0 ) · [ H W S · ( cos ( W D ) I 11 + sin ( W D ) I 12 ) A s u · I 11 , s u A s w · I 12 , s w ] + A h · cos ( θ 0 ) · I 1 , h ) ,
and
b 1 t r a n s = 1 π 0 2 π f t r a n s ( ϕ ) sin ( ϕ ) d ϕ = 1 π · ( sin ( θ 0 ) · [ H W S · ( cos ( W D ) I 21 + sin ( W D ) I 22 ) A s u · I 21 , s u A s w · I 22 , s w ] + A h · cos ( θ 0 ) · I 2 , h ) ,
where coefficients I m n , I m , x , and I m n , x with m , n = 1 , 2 and x = s u , s w , h are the auxiliary integrals of Table A1. As mentioned, the mnemonic subscript m = 1 ( m = 2 ) refers to the Fourier coefficient a 1 ( b 1 ) . The subscripts s u , s w , and h denote the “surge”, “sway”, and “heave”, respectively.

Appendix A.5. Auxiliary Integrals

Table A1. Auxiliary integrals for computation of the rotational- and translational-motion Fourier coefficients given by Equations (A3)–(A6), respectively. Note that not all possible subscript combinations are simultaneously used.
Table A1. Auxiliary integrals for computation of the rotational- and translational-motion Fourier coefficients given by Equations (A3)–(A6), respectively. Note that not all possible subscript combinations are simultaneously used.
Coef. g ( ϕ ) 0 2 π g ( ϕ )   d ϕ
I 11 cos ( ϕ ) · cos ( ϕ ϕ 0 ) π · cos ( ϕ 0 )
I 12 cos ( ϕ ) · sin ( ϕ ϕ 0 ) π · sin ( ϕ 0 )
I 21 sin ( ϕ ) · cos ( ϕ ϕ 0 ) π · sin ( ϕ 0 )
I 22 sin ( ϕ ) · sin ( ϕ ϕ 0 ) π · cos ( ϕ 0 )
I 1 , x cos ( ϕ ) · sin ( f x ϕ α x ) f x f x 2 1 · [ cos ( α x 2 π f x ) cos ( α x ) ]
if f x = 1 π sin ( α x )
I 2 , x sin ( ϕ ) · sin ( f x ϕ α x ) 1 f x 2 1 · [ sin ( α x 2 π f x ) sin ( α x ) ]
if f x = 1 π cos ( α x )
I 11 , x cos ( ϕ ) · cos ( ϕ ϕ 0 ) · sin ( f x ϕ α x ) 1 f x ( f x 2 4 ) [ f x sin ( ϕ 0 ) A x ( f x 2 2 ) cos ( ϕ 0 ) B x ]
if f x = 2 π 2 sin ( α x ϕ 0 )
I 12 , x cos ( ϕ ) · sin ( ϕ ϕ 0 ) · sin ( f x ϕ α x ) 1 f x ( f x 2 4 ) [ f x cos ( ϕ 0 ) A x + ( f x 2 2 ) sin ( ϕ 0 ) B x ]
if f x = 2 π 2 cos ( α x ϕ 0 )
I 21 , x sin ( ϕ ) · cos ( ϕ ϕ 0 ) · sin ( f x ϕ α x ) 1 f x ( f x 2 4 ) [ f x cos ( ϕ 0 ) A x + 2 sin ( ϕ 0 ) B x ]
if f x = 2 π 2 cos ( α x ϕ 0 )
I 22 , x sin ( ϕ ) · sin ( ϕ ϕ 0 ) · sin ( f x ϕ α x ) 1 f x ( f x 2 4 ) [ f x sin ( ϕ 0 ) A x + 2 cos ( ϕ 0 ) B x ]
if f x = 2 π 2 sin ( α x ϕ 0 )
where:being:
x [ r = r o l l , p = p i t c h A x = sin ( α x ) ( 1 cos ( 2 π f x ) )
s u = s u r g e , s w = s w a y , h = h e a v e ] + cos ( α x ) sin ( 2 π f x )
B x = sin ( α x ) sin ( 2 π f x )
cos ( α x ) ( 1 cos ( 2 π f x ) )

References

  1. Global Wind Energy Council. Global Wind Report 2018; Technical Report; Global Wind Energy Council: Brussels, Belgium, 2019. [Google Scholar]
  2. Barthelmie, R.; Courtney, M.; Højstrup, J.; Larsen, S. Meteorological aspects of offshore wind energy: Observations from the Vindeby wind farm. J. Wind. Eng. Ind. Aerodyn. 1996, 62, 191–211. [Google Scholar] [CrossRef]
  3. Ramírez, L.; Fraile, D.; Brindley, G. Offshore Wind in Europe Key Trends and Statistics 2019; Technical Report; WindEurope: Brussels, Belgium, 2020. [Google Scholar]
  4. Global Wind Energy Council. Global Wind Report 2021; Technical Report; Global Wind Energy Council: Brussels, Belgium, 2021. [Google Scholar]
  5. Kost, K.; Mayer, J.N.; Thomsen, J.; Hartmann, N.; Senkpiel, C.; Philipps, S.; Nold, S.; Lude, S.; Saad, N.; Schlegl, T. Levelized Cost of Electricity Renewable Energy Technologies; Technical Report; Fraunhofer Institut for Solar Energy Systems (ISE): Freiburg, Germany, 2013. [Google Scholar]
  6. Timmons, D. 25—Optimal Renewable Energy Systems: Minimizing the Cost of Intermittent Sources and Energy Storage. In A Comprehensive Guide to Solar Energy Systems; Letcher, T.M., Fthenakis, V.M., Eds.; Academic Press: Cambridge, MA, USA, 2018; pp. 485–504. [Google Scholar] [CrossRef]
  7. Carbon Trust. Carbon Trust Offshore Wind Accelerator Roadmap for the Commercial Acceptance of Floating LiDAR Technology; Technical Report; Carbon Trust: London, UK, 2018. [Google Scholar]
  8. Roland Berger. Offshore Wind toward 2020: On the Pathway to Cost Competitiveness; Technical Report; Roland Berger: Munich, Germany, 2013. [Google Scholar]
  9. Barthelmie, R.; Pryor, S. Can Satellite Sampling of Offshore Wind Speeds Realistically Represent Wind Speed Distributions? J. Appl. Meteorol. 2003, 42, 83–94. [Google Scholar] [CrossRef]
  10. Chang, R.; Zhu, R.; Badger, M.; Hasager, C.B.; Zhou, R.; Ye, D.; Zhang, X. Applicability of Synthetic Aperture Radar wind retrievals on offshore wind resources assessment in Hangzhou Bay, China. Energies 2014, 7, 3339–3354. [Google Scholar] [CrossRef] [Green Version]
  11. Hirth, B.D.; Schroeder, J.L.; Gunter, W.S.; Guynes, J.G. Measuring a utility-scale turbine wake using the TTUKa mobile research radars. J. Atmos. Ocean. Technol. 2012, 29, 765–771. [Google Scholar] [CrossRef] [Green Version]
  12. Vogt, S.; Thomas, P. SODAR—A useful remote sounder to measure wind and turbulence. J. Wind. Eng. Ind. Aerodyn. 1995, 54, 163–172. [Google Scholar] [CrossRef]
  13. Lang, S.; McKeogh, E. LIDAR and SODAR Measurements of Wind Speed and Direction in Upland Terrain for Wind Energy Purposes. Remote Sens. 2011, 3, 1871–1901. [Google Scholar] [CrossRef] [Green Version]
  14. International Energy Association. State of the art of Remote Wind Speed Sensing Techniques Using Sodar, Lidar and Satellites; Technical Report; International Energy Association: Paris, France, 2007. [Google Scholar]
  15. Sempreviva, A.M.; Barthelmie, R.J.; Pryor, S.C. Review of Methodologies for Offshore Wind Resource Assessment in European Seas. Surv. Geophys. 2008, 29, 471–497. [Google Scholar] [CrossRef]
  16. Slinger, C.; Harris, M. Introduction to Continuous-Wave Doppler Lidar; Summer School in Remote Sensing for Wind Energy: Boulder, CO, USA, 2012. [Google Scholar]
  17. Lolli, S.; Delaval, A.; Loth, C.; Garnier, A.; Flamant, P.H. 0.355-micrometer direct detection wind lidar under testing during a field campaign in consideration of ESA’s ADM-Aeolus mission. Atmos. Meas. Tech. 2013, 6, 3349–3358. [Google Scholar] [CrossRef] [Green Version]
  18. Reitebuch, O.; Werner, C.; Leike, I.; Delville, P.; Flamant, P.H.; Cress, A.; Engelbart, D. Experimental Validation of Wind Profiling Performed by the Airborne 10-μm Heterodyne Doppler Lidar WIND. J. Atmos. Ocean. Technol. 2001, 18, 1331–1344. [Google Scholar] [CrossRef]
  19. Schuon, F.; González, D.; Rocadenbosch, F.; Bischoff, O.; Jané, R. KIC InnoEnergy Project Neptune: Development of a Floating LiDAR Buoy for Wind, Wave and Current Measurements. In Proceedings of the DEWEK 2012 German Wind Energy Conference, Bremen, Germany, 7–8 November 2012. [Google Scholar]
  20. Silva, M.P.A.d.; Rocadenbosch, F.; Farré-Guarné, J.; Salcedo-Bosch, A.; González-Marco, D.; Peña, A. Assessing Obukhov Length and Friction Velocity from Floating Lidar Observations: A Data Screening and Sensitivity Computation Approach. Remote Sens. 2022, 14, 1394. [Google Scholar] [CrossRef]
  21. Pichugina, Y.; Banta, R.; Brewer, W.; Sandberg, S.; Hardesty, R. Doppler Lidar–Based Wind-Profile Measurement System for Offshore Wind-Energy and Other Marine Boundary Layer Applications. J. Appl. Meteorol. Climatol. 2012, 51, 327–349. [Google Scholar] [CrossRef]
  22. Gottschall, J.; Wolken-Möhlmann, G.; Lange, B. About offshore resource assessment with floating lidars with special respect to turbulence and extreme events. J. Phys. Conf. Ser. 2014, 555, 12–43. [Google Scholar] [CrossRef]
  23. Salcedo-Bosch, A.; Rocadenbosch, F.; Gutiérrez-Antuñano, M.A.; Tiana-Alsina, J. Estimation of Wave Period from Pitch and Roll of a Lidar Buoy. Sensors 2021, 21, 1310. [Google Scholar] [CrossRef] [PubMed]
  24. Salcedo-Bosch, A.; Gutierrez-Antunano, M.A.; Tiana-Alsina, J.; Rocadenbosch, F. Motional Behavior Estimation Using Simple Spectral Estimation: Application to the Off-Shore Wind Lidar. In Proceedings of the 2020 IEEE International Geoscience and Remote Sensing Symposium (IGARSS-2020), Waikoloa, HI, USA, 26 September–2 October 2020. [Google Scholar]
  25. Gutiérrez-Antuñano, M.; Tiana-Alsina, J.; Salcedo, A.; Rocadenbosch, F. Estimation of the Motion-Induced Horizontal-Wind-Speed Standard Deviation in an Offshore Doppler Lidar. Remote Sens. 2018, 10, 2037. [Google Scholar] [CrossRef] [Green Version]
  26. Salcedo-Bosch, A.; Gutierrez-Antunano, M.A.; Tiana-Alsina, J.; Rocadenbosch, F. Floating Doppler Wind Lidar Measurement of Wind Turbulence: A Cluster Analysis. In Proceedings of the 2020 IEEE International Geoscience and Remote Sensing Symposium (IGARSS-2020), Waikoloa, HI, USA, 26 September–2 October 2020. [Google Scholar]
  27. Gutiérrez-Antuñano, M.A.; Tiana-Alsina, J.; Rocadenbosch, F. Performance evaluation of a floating lidar buoy in nearshore conditions. Wind. Energy 2017, 20, 1711–1726. [Google Scholar] [CrossRef] [Green Version]
  28. Kelberlau, F.; Neshaug, V.; Lønseth, L.; Bracchi, T.; Mann, J. Taking the Motion out of Floating Lidar: Turbulence Intensity Estimates with a Continuous-Wave Wind Lidar. Remote Sens. 2020, 12, 898. [Google Scholar] [CrossRef] [Green Version]
  29. Gottschall, J.; Gribben, B.; Stein, D.; Würth, I. Floating lidar as an advanced offshore wind speed measurement technique: Current technology status and gap analysis in regard to full maturity. Wiley Interdiscip. Rev. Energy Environ. 2017, 6, e250. [Google Scholar] [CrossRef]
  30. Mücke, T.; Kleinhans, D.; Peinke, J. Atmospheric turbulence and its influence on the alternating loads on wind turbines. Wind Energy 2011, 14, 301–316. [Google Scholar] [CrossRef]
  31. Désert, T.; Knapp, G.; Aubrun, S. Quantification and Correction of Wave-Induced Turbulence Intensity Bias for a Floating LIDAR System. Remote Sens. 2021, 13, 2973. [Google Scholar] [CrossRef]
  32. Gottschall, J.; Wolken-Möhlmann, G.; Viergutz, T.; Lange, B. Results and conclusions of a floating-lidar offshore test. Energy Procedia 2014, 53, 156–161. [Google Scholar] [CrossRef]
  33. Salcedo-Bosch, A.; Rocadenbosch, F.; Sospedra, J. A Robust Adaptive Unscented Kalman Filter for Floating Doppler Wind-LiDAR Motion Correction. Remote Sens. 2021, 13, 4167. [Google Scholar] [CrossRef]
  34. Salcedo-Bosch, A.; Rocadenbosch, F.; Sospedra, J. Enhanced Dual Filter for Floating Wind Lidar Motion Correction: The Impact of Wind and Initial Scan Phase Models. Remote Sens. 2022, 14, 4704. [Google Scholar] [CrossRef]
  35. Tiana-Alsina, J.; Gutiérrez, M.A.; Würth, I.; Puigdefàbregas, J.; Rocadenbosch, F. Motion compensation study for a floating doppler wind lidar. In Proceedings of the Geoscience and Remote Sensing Symposium Proceedings, Milan, Italy, 26–31 July 2015. [Google Scholar]
  36. Bischoff, O.; Yu, W.; Gottschall, J.; Cheng, P.W. Validating a simulation environment for floating lidar systems. J. Phys. Conf. Ser. 2018, 1037, 052036. [Google Scholar] [CrossRef] [Green Version]
  37. Salcedo-Bosch, A.; Farré-Guarné, J.; Sala-Álvarez, J.; Villares-Piera, J.; Rocadenbosch, F.; Tanamachi, R.L. Floating Doppler wind lidar simulator for horizontal wind speed measurement error assessment. In Proceedings of the 2021 IEEE International Geoscience and Remote Sensing Symposium (IGARSS-2021), Brussels, Belgium, 11–16 July 2021. [Google Scholar]
  38. Kelberlau, F.; Mann, J. Quantification of Motion-Induced Measurement Error on Floating Lidar Systems. Atmos. Meas. Tech. Discuss. 2022, 2022, 1–31. [Google Scholar] [CrossRef]
  39. Sospedra, J.; Cateura, J.; Puigdefàbregas, J. Novel multipurpose buoy for offshore wind profile measurements EOLOS platform faces validation at ijmuiden offshore metmast. Sea Technol. 2015, 56, 25–28. [Google Scholar]
  40. Scientific, C. ZephIR 300; Technical Report; Campbell Scientific: Edmonton, AB, Canada, 2016. [Google Scholar]
  41. Pitter, M.; Burin des Roziers, E.; Medley, J.; Mangat, M.; Slinger, C.; Harris, M. Performance Stability of Zephir in High Motion Enviroments: Floating and Turbine Mounted; Technical Report; ZephIR: Ledbury, UK, 2014. [Google Scholar]
  42. Knoop, S.; Bosveld, F.C.; de Haij, M.J.; Apituley, A. A 2-year intercomparison of continuous-wave focusing wind lidar and tall mast wind measurements at Cabauw. Atmos. Meas. Tech. 2021, 14, 2219–2235. [Google Scholar] [CrossRef]
  43. Gutiérrez Antuñano, M. Doppler Wind LIDAR Systems Data Processing and Applications: An Overview towards Developing the New Generation of Wind Remote-Sensing Sensors for Off-Shore Wind Farms. Ph.D. Thesis, Universitat Politècnica de Catalunya, Barcelona, Spain, 2019. [Google Scholar]
  44. KIC InnoEnnergy. Neptune Project Leaflet; Technical Report; KIC InnoEnergy: Eindhoven, The Netherlands, 2015. [Google Scholar]
  45. MicroStrain. 3DM-GX3® -45 Theory of Operation; Technical Report; 459 Hurricane Lane: Williston, VT, USA, 2012. [Google Scholar]
  46. Barlow, R.J. Statistics: A Guide to the Use of Statistical Methods in the Physical Sciences; Manchester Physics Series; Wiley: Chichester, UK, 1989. [Google Scholar]
  47. Pẽnha, A.; Hasager, C.; Merete, B.; Barthelmie, R.; Ferhat, B.; Jean-Pierre, C.; Stefan, E.; Frandsen, S.; Michael, H.; Ioanna, K.; et al. Remote Sensing for Wind Energy; DTU Wind Energy: Roskilde, Denmark, 2015. [Google Scholar]
  48. Tolstov, J.P. Fourier Series; Prentice-Hall Inc.: Englewood Cliffs, NJ, USA, 1962. [Google Scholar]
  49. Stearns, S.D.; Hush, D.R. Digital Signal Processing with Examples in MATLAB, 2nd ed.; Taylor & Francis Inc.: Bosa Roca, FL, USA, 2011. [Google Scholar]
  50. Stade, E. Fourier Analysis, 1st ed.; Wiley-Interscience: Hoboken, NJ, USA, 2015. [Google Scholar]
  51. Kogaki, T.; Sakurai, K.; Shimada, S.; Kawabata, H.; Otake, Y.; Kondo, K.; Fujita, E. Field Measurements of Wind Characteristics Using LiDAR on a Wind Farm with Downwind Turbines Installed in a Complex Terrain Region. Energies 2020, 13, 5135. [Google Scholar] [CrossRef]
  52. Palais, B.; Palais, R.; Rodi, S. A Disorienting Look at Euler’s Theorem on the Axis of a Rotation. Am. Math. Mon. 2009, 116, 892–909. [Google Scholar] [CrossRef] [Green Version]
  53. Roithmayr, C.M.; Hodges, D.H. Dynamics: Theory and Application of Kane’s Method; American Society of Mechanical Engineers Digital Collection; NASA-Langley Research Center: Hampton, VA, USA, 2016.
  54. Proakis, J.; Manolakis, D. Digital Signal Processing, 4th ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2006. [Google Scholar]
  55. Plancherel, M.; Leffler, M. Contribution À ĽÉtude de la reprÉsentation D’une fonction arbitraire par des intÉgrales dÉfinies. Rend. Del Circ. Mat. Palermo (1884–1940) 1910, 30, 289–335. [Google Scholar] [CrossRef]
  56. Yao, G.; Xie, J.; Huang, W. HF Radar Ocean Surface Cross Section for the Case of Floating Platform Incorporating a Six-DOF Oscillation Motion Model. IEEE J. Ocean. Eng. 2021, 46, 156–171. [Google Scholar] [CrossRef]
  57. Wan, B.; Wu, X.; Yue, X.; Zhang, L.; Wang, L. Calibration of Phased-Array High-Frequency Radar on an Anchored Floating Platform. Remote Sens. 2022, 14, 2174. [Google Scholar] [CrossRef]
  58. Wagner, R.; Mikkelsen, T.; Courtney, M. Investigation of Turbulence Measurements with a Continuous Wave, Conically Scanning LiDAR; Technical Report; DTU: Roskilde, Denmark, 2009. [Google Scholar]
  59. Gottschall, J.; Courtney, M.S.; Wagner, R.; Jørgensen, H.E.; Antoniou, I. Lidar profilers in the context of wind energy–a verification procedure for traceable measurements. Wind Energy 2012, 15, 147–159. [Google Scholar] [CrossRef]
  60. Jeffreys, H.; Taylor, G.I. On the formation of water waves by wind. Proc. R. Soc. Lond. Ser. A Contain. Pap. Math. Phys. Character 1925, 107, 189–206. [Google Scholar] [CrossRef] [Green Version]
  61. IEA Wind TCP. Available online: https://iea-wind.org/ (accessed on 3 March 2023).
  62. Task 52—Large-Scale Deployment of Wind Lidar|IEA Wind TCP. Available online: https://iea-wind.org/task52/ (accessed on 3 March 2023).
Figure 1. The motion geometry of the FDWL buoy.
Figure 1. The motion geometry of the FDWL buoy.
Remotesensing 15 01478 g001
Figure 2. Example of the estimation of the H W S total error in a LiDAR scan ( H W S = 10 m/s, V W S = 0 m/s, and a W D ranging from 0 to 360 deg, roll-only motion). (a) Color plot of the H W S error as a function of the wind direction, W D , and initial phase, ϕ 0 . (b) H W S error for W D = 275 deg (vertical black line in panel (a)) as a function of the initial phase. The blue horizontal line indicates the mean H W S error.
Figure 2. Example of the estimation of the H W S total error in a LiDAR scan ( H W S = 10 m/s, V W S = 0 m/s, and a W D ranging from 0 to 360 deg, roll-only motion). (a) Color plot of the H W S error as a function of the wind direction, W D , and initial phase, ϕ 0 . (b) H W S error for W D = 275 deg (vertical black line in panel (a)) as a function of the initial phase. The blue horizontal line indicates the mean H W S error.
Remotesensing 15 01478 g002
Figure 3. Block diagram depicting the Section 3.5 estimation procedure for the bias and T I increment.
Figure 3. Block diagram depicting the Section 3.5 estimation procedure for the bias and T I increment.
Remotesensing 15 01478 g003
Figure 4. The FDWL motion characterization as a sinusoidal signal: roll motion example (PdP, 24 June 2013, 00:00 LT). (a) The 10 min roll time series (blue trace) and fitted sinusoidal signal (black trace). (b) PSD of the panel (a) time series.
Figure 4. The FDWL motion characterization as a sinusoidal signal: roll motion example (PdP, 24 June 2013, 00:00 LT). (a) The 10 min roll time series (blue trace) and fitted sinusoidal signal (black trace). (b) PSD of the panel (a) time series.
Remotesensing 15 01478 g004
Figure 5. Validation of the analytical-error formulation in the roll-only scenario #1: HWS error (parameters: A r o l l = 10 deg, f r o l l = 0.3 Hz and α r o l l = 0 deg. Input wind: H W S = 10 m/s, V W S = 0 m/s). (a) The analytically estimated H W S error as a function of the wind direction (WD) and the initial scan phase ( ϕ 0 ). (b) The simulator-estimated H W S error. (c) The difference of (a) minus (b) (in absolute value).
Figure 5. Validation of the analytical-error formulation in the roll-only scenario #1: HWS error (parameters: A r o l l = 10 deg, f r o l l = 0.3 Hz and α r o l l = 0 deg. Input wind: H W S = 10 m/s, V W S = 0 m/s). (a) The analytically estimated H W S error as a function of the wind direction (WD) and the initial scan phase ( ϕ 0 ). (b) The simulator-estimated H W S error. (c) The difference of (a) minus (b) (in absolute value).
Remotesensing 15 01478 g005
Figure 6. Validation of the analytical-error formulation in the translational-only scenario #2: HWS error (parameters: A s u r g e = 10 deg). (a) The analytically estimated H W S error as a function of the wind direction (WD) and the initial scan phase ( ϕ 0 ). (b) The simulator-estimated H W S error. (c) The difference of (a) minus (b) (in absolute value).
Figure 6. Validation of the analytical-error formulation in the translational-only scenario #2: HWS error (parameters: A s u r g e = 10 deg). (a) The analytically estimated H W S error as a function of the wind direction (WD) and the initial scan phase ( ϕ 0 ). (b) The simulator-estimated H W S error. (c) The difference of (a) minus (b) (in absolute value).
Remotesensing 15 01478 g006
Figure 7. Validation of the analytical-error formulation in the 6 DoF scenario #3: HWS error (parameters: A r o l l = A p i t c h = 10 deg, and A s u r g e = A s w a y = A h e a v e = 2 m/s). (a) The analytically estimated H W S error as a function of the wind direction (WD) and the initial scan phase ( ϕ 0 ). (b) The simulator-estimated H W S error. (c) The difference of (a) minus (b) (in absolute value).
Figure 7. Validation of the analytical-error formulation in the 6 DoF scenario #3: HWS error (parameters: A r o l l = A p i t c h = 10 deg, and A s u r g e = A s w a y = A h e a v e = 2 m/s). (a) The analytically estimated H W S error as a function of the wind direction (WD) and the initial scan phase ( ϕ 0 ). (b) The simulator-estimated H W S error. (c) The difference of (a) minus (b) (in absolute value).
Remotesensing 15 01478 g007
Figure 8. Validation of the analytical-error formulation in a real case scenario (PdP campaign, 25 June 2013, 00:00 UTC to 27 June 2013, 00:00 UTC): TI increment. (a) The analytically estimated vs. the measured TI-increment time series. The black dots represent the 10 min analytically estimated T I increment, Δ T I e s t (Equation (7)). The blue dots represent the 10 min measured T I increment, Δ T I m e a s (Equation (6)). (b) The mean tilt amplitude (Equation (34)) and the mean translational-velocity amplitude (Equation (35)) time series. In both panels, the solid trace represents the 1 h averaged time series.
Figure 8. Validation of the analytical-error formulation in a real case scenario (PdP campaign, 25 June 2013, 00:00 UTC to 27 June 2013, 00:00 UTC): TI increment. (a) The analytically estimated vs. the measured TI-increment time series. The black dots represent the 10 min analytically estimated T I increment, Δ T I e s t (Equation (7)). The blue dots represent the 10 min measured T I increment, Δ T I m e a s (Equation (6)). (b) The mean tilt amplitude (Equation (34)) and the mean translational-velocity amplitude (Equation (35)) time series. In both panels, the solid trace represents the 1 h averaged time series.
Remotesensing 15 01478 g008
Figure 9. Validation of the analytical-error formulation: Global statistics for the PdP campaign (6 to 30 June 2013) and IJmuiden campaign (1 April to 1 June 2015). The left panels show the results obtained for the PdP campaign, and the right panels show the results obtained for the IJmuiden campaign, with the comparison between the measured and the estimated T I increment as a function of different motion conditions clustered by (panels (a,b)) mean tilt amplitude (rotational motion), (panels (c,d)) mean translational-velocity amplitude, and (panels (e,f)) mean HWS. The error bars indicate the 25th and 75th percentiles. The bar graphs represent the number of occurrences in the data for each category of bin. The solid red lines represent the R M S E between Δ T I e s t and Δ T I m e a s for each category of bin. The bins containing less than 5% of the PdP population are gray shaded.
Figure 9. Validation of the analytical-error formulation: Global statistics for the PdP campaign (6 to 30 June 2013) and IJmuiden campaign (1 April to 1 June 2015). The left panels show the results obtained for the PdP campaign, and the right panels show the results obtained for the IJmuiden campaign, with the comparison between the measured and the estimated T I increment as a function of different motion conditions clustered by (panels (a,b)) mean tilt amplitude (rotational motion), (panels (c,d)) mean translational-velocity amplitude, and (panels (e,f)) mean HWS. The error bars indicate the 25th and 75th percentiles. The bar graphs represent the number of occurrences in the data for each category of bin. The solid red lines represent the R M S E between Δ T I e s t and Δ T I m e a s for each category of bin. The bins containing less than 5% of the PdP population are gray shaded.
Remotesensing 15 01478 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Salcedo-Bosch, A.; Farré-Guarné, J.; Araújo da Silva, M.P.; Rocadenbosch, F. A Unified Formulation for the Computation of the Six-Degrees-of-Freedom-Motion-Induced Errors in Floating Doppler Wind LiDARs. Remote Sens. 2023, 15, 1478. https://doi.org/10.3390/rs15061478

AMA Style

Salcedo-Bosch A, Farré-Guarné J, Araújo da Silva MP, Rocadenbosch F. A Unified Formulation for the Computation of the Six-Degrees-of-Freedom-Motion-Induced Errors in Floating Doppler Wind LiDARs. Remote Sensing. 2023; 15(6):1478. https://doi.org/10.3390/rs15061478

Chicago/Turabian Style

Salcedo-Bosch, Andreu, Joan Farré-Guarné, Marcos Paulo Araújo da Silva, and Francesc Rocadenbosch. 2023. "A Unified Formulation for the Computation of the Six-Degrees-of-Freedom-Motion-Induced Errors in Floating Doppler Wind LiDARs" Remote Sensing 15, no. 6: 1478. https://doi.org/10.3390/rs15061478

APA Style

Salcedo-Bosch, A., Farré-Guarné, J., Araújo da Silva, M. P., & Rocadenbosch, F. (2023). A Unified Formulation for the Computation of the Six-Degrees-of-Freedom-Motion-Induced Errors in Floating Doppler Wind LiDARs. Remote Sensing, 15(6), 1478. https://doi.org/10.3390/rs15061478

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop