This study is a follow up on a previous one carried out within the frame of the French project SMARTEOLE, during which, a ground-based scanning LiDAR measurement campaign was conducted in the onshore wind farm of Sole du Moulin Vieux. That previous study focused on the wakes of two wind turbines that experienced different degrees of interaction depending on the incoming wind direction, through the processing of LiDAR measurements. The measurement duration (7 months) ensured the statistical convergence of the ensemble-averaged flow fields obtained after holding a categorisation process based on the wind speed at hub height, wind direction, and atmospheric stability, where only near-neutral stability conditions were considered. The present study focuses on integrating the operational data of the wind turbines through SCADA processing to complement the LiDAR wake field observations and to be used as an input for analytical wake models. First, the correlation between the atmospheric stability, deduced from MERRA-2 dataset, and the free-stream turbulence intensity, measured by the wind turbines’ anemometers, is studied for different wind speed ranges. It is observed that the turbulence intensity tends towards a consistent value as the atmospheric stability approaches near-neutral stability conditions, giving confidence into the applied strategy of data categorisation based on MERRA-2 outputs. The influence of the degree of wake interaction on the wake added turbulence, the velocity and power deficits between both turbines is assessed. Clear trends between the wake added turbulence and both the velocity and power deficits are detected. Consequently, two fitting laws are proposed. Then, different analytical wake models and wake superposition methods are fed with the operational data deduced from the processed SCADA data, and are used for predicting the evolution of the velocity deficit within the wake. Some statistical metrics are used for error quantification of the different engineering wake models compared to the scanning LiDAR measurements, used as reference, and Blondel and Cathelain produces the closest results to the field measurements.