A wind flow sensing system for determining turbulence quantities of wind flow at a set of different altitudes above a terrain from a set of measurements of radial velocities at each of the altitudes along the surface of a cone and along a center line of the cone, comprising: an input interface configured to receive the set of measurements of radial velocities at line-of-sight points above the terrain for each of the altitudes for a set of time steps; a processor configured to: estimate velocity fields for each of the altitudes based on data assimilation of the velocity fields above the terrain to fit the measurements of radial velocities, wherein the velocity fields are estimated for the set of time steps; and determine turbulence quantities based on unbiased horizontal velocities for each time step and the average of the unbiased horizontal velocity; and an output interface configured to render the turbulence quantities; characterized in that the processor is further configured to: estimate the unbiased horizontal velocities at each of the altitudes and for each time step as a horizontal projection of the corresponding radial velocities, and correct the unbiased horizontal velocities with corresponding horizontal derivatives of vertical velocities of the estimated velocity fields determined for the corresponding time step; determine, at each of the altitudes, an average of the unbiased horizontal velocities for a period of time including the set of time steps; wherein at each of the altitudes, turbulence quantities are determined based on the unbiased horizontal velocities for each time step and the average of the unbiased horizontal velocity; and the output interface is configured to render the turbulence quantities at each of the altitudes,wherein the data assimilation is performed by simulating computational fluid dynamics, CFD, of the wind flow or by analytical fluid mechanics approximation using potential flow approximation. an input interface configured to receive the set of measurements of radial velocities at line-of-sight points above the terrain for each of the altitudes for a set of time steps; a processor configured to: estimate velocity fields for each of the altitudes based on data assimilation of the velocity fields above the terrain to fit the measurements of radial velocities, wherein the velocity fields are estimated for the set of time steps; and determine turbulence quantities based on unbiased horizontal velocities for each time step and the average of the unbiased horizontal velocity; and an output interface configured to render the turbulence quantities; characterized in that the processor is further configured to: estimate the unbiased horizontal velocities at each of the altitudes and for each time step as a horizontal projection of the corresponding radial velocities, and correct the unbiased horizontal velocities with corresponding horizontal derivatives of vertical velocities of the estimated velocity fields determined for the corresponding time step; determine, at each of the altitudes, an average of the unbiased horizontal velocities for a period of time including the set of time steps; wherein at each of the altitudes, turbulence quantities are determined based on the unbiased horizontal velocities for each time step and the average of the unbiased horizontal velocity; and estimate velocity fields for each of the altitudes based on data assimilation of the velocity fields above the terrain to fit the measurements of radial velocities, wherein the velocity fields are estimated for the set of time steps; and determine turbulence quantities based on unbiased horizontal velocities for each time step and the average of the unbiased horizontal velocity; and an output interface configured to render the turbulence quantities; characterized in that the processor is further configured to: estimate the unbiased horizontal velocities at each of the altitudes and for each time step as a horizontal projection of the corresponding radial velocities, and correct the unbiased horizontal velocities with corresponding horizontal derivatives of vertical velocities of the estimated velocity fields determined for the corresponding time step; determine, at each of the altitudes, an average of the unbiased horizontal velocities for a period of time including the set of time steps; wherein at each of the altitudes, turbulence quantities are determined based on the unbiased horizontal velocities for each time step and the average of the unbiased horizontal velocity; and the output interface is configured to render the turbulence quantities at each of the altitudes,wherein the data assimilation is performed by simulating computational fluid dynamics, CFD, of the wind flow or by analytical fluid mechanics approximation using potential flow approximation.