Digital twins of the mooring line tension for floating offshore wind turbines to improve monitoring, lifespan, and safety

The number of installed floating offshore wind turbines (FOWTs) has doubled since 2017, quadrupling the total installed capacity, and is expected to increase significantly over the next decade. Consequently, there is a growing consideration towards the main challenges for FOWT projects: monitoring the system’s integrity, extending the lifespan of the components, and maintaining FOWTs safely at scale. Effectively and efficiently addressing these challenges would unlock the wide-scale deployment of FOWTs. In this work, we focus on one of the most critical components of the FOWTs, the Mooring Lines (MoLs), which are responsible for fixing the structure to the seabed. The primary mechanical failure mechanisms in MoLs are extreme load and fatigue, both of which are functions of the axial tension. An effective solution to detect long-term drifts in the mechanical response of the MoLs is to develop a Digital Twin (DT) able to accurately predict the behaviour of the healthy system to compare with the actual one. Moreover, we will develop another DT able to accurately predict the near future axial tension as an effective tool to improve the lifespan of the MoLs and the safety of FOWT maintenance operations. In fact, by changing the FOWT operational settings, according to the DT prediction, operators can increase the lifespan of the MoLs by reducing the stress and, additionally, in the case where FOWT operational maintenance is in progress, the prediction from the DT can serve as early safety warning to operators. Authors will leverage operational data collected from the world’s first commercial floating-wind farm [the Hywind Pilot Park (https://www.equinor.com/en/what-we-do/floating-wind/hywind-scotland.html.)] in 2018, to investigate the effectiveness of DTs for the prediction of the MoL axial tension for the two scenarios depicted above. The DTs will be developed using state-of-the-art data-driven methods, and results based on real operational data will support our proposal.


Introduction
Floating wind is one of the fastest-growing sectors within the Offshore Renewable Energy Industry and internationally recognised as one of the most promising renewable energy sources to satisfy a significant proportion of global energy demands (Carbon Trust 2015). The ability to economically deploy floating offshore wind turbines (FOWTs) in deepwater areas, that were previously unfeasible for development using fixed-bottom turbines, is one of the fundamental driving forces behind the success of floating wind (Carbon Trust 2018). In fact, deepwater areas are often characterised by higher average wind speeds and consequently an higher average capacity factors that could improve the economic viability of offshore wind energy (Hannon et al. 2019). However, floating wind is still an emerging market, and only a limited number of pilots have been deployed, so there is still a significant amount work required to address the unique aspects to floating wind: monitoring the system's integrity, extending the lifespan of the components, and maintaining FOWTs safely at scale (Carbon Trust 2015. Nevertheless, due to the success of Pilot FOWTs around the world, the industry is now focused on addressing these remaining challenges before deploying FOWTs at scale in deepwater (Carbon Trust 2020). In particular, the use, monitoring, and maintenance of the station-keeping devices (i.e., the Mooring Lines-MoLs) devoted to anchor the FOWT structure in place pose some of the most prominent challenges to overcome (Carbon Trust 2020).
Experience from the Oil&Gas industry indicates that the demanding environmental conditions such as the corrosive salt water and forceful waves, combined with the isolation of the deployment sites, are particularly damaging to the MoLs and may pose issues for checking and maintaining the integrity of FOWTs (Butterfield et al. 2007). In addition, unlike floating production storage and offloading (FPSO) vessels, where there are typically between 12 and 24 MoLs, economic drivers within the renewables sector tend to produce designs with no redundancy (Fugro 2020). For this reason, each MoL is critical to the FOWT structure, and failure is catastrophic (Ma et al. 2020). Nevertheless, the extreme conditions in which MoLs operate call for regular inspections to prevent disastrous failures and downtime. Inspections are usually performed in two ways: (1) close visual inspection by divers and (2) through remotely operated vehicles (ROVs) (Brown et al. 2005). When considering (1), diver inspections are not favoured due to the risk related to operating in deep water and human bias (Angulo et al. 2017). Similarly, with (2), ROVs are employed with limited success, since they need to be re-calibrated between successive measurements, and such may present errors in transmitting their positioning, affecting the localisation of damage (ABSG Consulting 2015). Moreover, both options remain uneconomical (ABSG Consulting 2015). For these reasons, within the Oil and Gas industry, alternative strategies for monitoring MoLs integrity are actually employed, like using load cell to detect rapid failure by extreme load (Lu 2016) or using GPS-based devices (e.g., LifeLine JIP (Marin LifeLine 2021)) to detect an unusual excursion caused by the loss of a single mooring (Marin LifeLine 2021).
Regarding FOWT failures and maintenance activities, since this industry is in its infancy, there are very little publicly available data describing the expected failure rates of FOWTs or the required maintenance schedules (Santos et al. 2016). According to (Röckmann et al. 2017) for a fixedbottom offshore wind farm with 200 turbines, it may be necessary to perform up to 3000 maintenance visits per year. Since the mechanics of FOWTs are similar to the fixedbottom turbines, it is safe to consider Floating-Wind farms will require (at least) a similar number of maintenance vis-its. With this in mind, the ability to schedule maintenance by forecasting the health status of FOWTs becomes fundamental to deploy FOWTs maintenance at scale (Ren et al. 2021). Similar considerations should also be made towards the safety of the maintenance operations. In fact, since a failure could be catastrophic (Ma et al. 2020), developing an early warning system able to predict a failure early enough to secure a safe working environment is essential to ensure the safety of maintenance operations.
To forecast MoL failures (short term) or health conditions (long term), it is required to develop computational models of the mooring systems. Contrarily to visual inspections, computational models can ensure an economic, reliable, and effective monitoring tools for the integrity of the MoLs (Angulo et al. 2017).
In the literature, it is possible to find a large body of work addressing MoL failures (Ma et al. 2013;Davidson and Ringwood 2017;Brown et al. 2005;Al-Solihat and Nahon 2018;Aranha and Pinto 2001;Bhinder et al. 2015;Borg et al. 2014;Hsu et al. 2014Hsu et al. , 2015Jin and Kim 2018;Qiao et al. 2020;Maroju et al. 2013;de Pina et al. 2013;Prislin and Maroju 2017;Jaiswal and Ruskin 2019;Li et al. 2018). The primary mechanical failure mechanisms in mooring systems are extreme load and fatigue, both of which are functions of the axial tension (Davidson and Ringwood 2017). By exploiting this relationship, it is possible to use the MoL axial tension as an indicator for the health status of the Mooring System. In addition, (Brown et al. 2005) acknowledge MoL failures are often focused at the top of the chain, as these are the points of highest tension and the most stressed locations. The works available in the literature can be grouped into two main families: Physical Models (PMs) and Data-Driven Models (DDMs). PMs and DDMs can be used to detect long-term drifts in the mechanical response of the MoLs by developing a Digital Twin (DT) able to accurately predict the behaviour of the healthy system to compare with the actual one. A DT is a specific type of model that embodies a precise digital copy of a physical system (Oneto et al. 2018). Consequently, DTs are an effective method of forecasting the future behaviour of a system after the DT has learned the behaviour from historical examples (Coraddu et al. 2021). Moreover, as an effective tool to improve the lifespan of the MoLs and safety of FOWT maintenance operations, it is possible to develop DTs to predict the near future axial tension.
Concerning the PMs, they are primarily based on the mechanistic knowledge of the MoLs behaviour. Borg et al. (2014) provided a good overview of the mooring dynamics PMs, ranging from simple, quick but approximate methods (linear force-displacement or force-displacement-velocity models), to the high fidelity, but computationally expensive methods (finite-element approach), comparing the main characteristics, as shown in Table 1. Focusing on failure modes, the linear model is not suitable, since it is valid only within a limited range of platform displacements (therefore not suitable for extreme loads), and by definition, it does not capture the non-linear dynamics of the tension and such, can under/overestimate the tension cycles. Moreover, it only models the total tension of all the mooring lines; therefore, no information about the tension in the single lines can be inferred (Cevasco et al. 2018). The quasi-static model is capable of modelling the single line tension and the non-linearity of the tension in the MoLs. It is considered suitable to estimate the global response motion of the platform, but may not be as accurate as a fully dynamic model (such as Multibody or Finite Elements) at predicting single line tension (Karimirad 2013). In particular, while at static level, a quasi-static model can be in very good agreement with multibody or finite-element models. However, when considering a stochastic load condition (e.g., random waves from a wave spectrum), the dynamic tension oscillation may be overestimated, since the inertia of the lines and, more importantly, the hydrodynamic viscous damping loads acting on the mooring lines are not captured (Hall et al. 2014). The multibody mooring model (also called the "lumped mass approach") is able to capture the effects of the line inertial loads and the hydrodynamic viscous forces acting on the mooring lines, and it is also suitable to estimate the extreme and oscillatory tension loads at a preliminary design stage, since it is characterised by a sufficient accuracy but still within an acceptable computational cost (Cevasco et al. 2018). The finite-element method is suggested for very accurate analyses, but this is usually reserved for a limited set of load cases/simulations, typically at an advanced design stage, due to its substantially higher computational costs (Cevasco et al. 2018;Hall et al. 2014).
Instead, DDMs allow us to develop computationally aware real-time monitoring systems for MoLs by learning the input-output behaviour of a system from historical examples without any a-priori knowledge about the MoLs. DDMs require a single intensive learning phase (i.e., model construction) and benefit from a computationally inexpensive forward phase (i.e., model used as a predictor) (Coraddu et al. 2020), making them well suited to develop DTs. Maroju et al. (2013);de Pina et al. (2013); Prislin and Maroju (2017); Jaiswal and Ruskin (2019) have all proposed using DDMs to monitor other marine energy systems (Floating platforms, FPSO vessels, etc.), and noted the potential for predictive models to reduce the cost of operational monitoring. In particular, de Pina et al. (2013) propose a DDM to approximate the MoL tension in place of a computationally expensive PM (specifically a Finite-Element Analysis) by use of a nonlinear autoregressive exogenous model. This model could accurately simulate the tension response over a period of multiple hours; however, the algorithm was effective only in forecasting the tension in the short term (0.1s). Fatigue damage of MoLs has been investigated by (Li et al. 2018) using an Artificial Neural Network (ANN), where they observe the significant influence of the environmental conditions (waves, winds, and currents) on the tension bands. Additionally, (Hsu et al. 2015(Hsu et al. , 2017 demonstrate the potential benefits of using DDMs for FOWT MoLs when scaling health monitoring tools site-wide.
To the best of the authors' knowledge, no one has yet developed DTs of the FOWT MoLs to detect long-term drifts in the mechanical response or to accurately predict the near future axial tension. Moreover, the previous literature on monitoring FOWTs has been focused only on synthetic data and scenarios. Developing strategies to address these problems and testing the solutions with real data could provide the necessary information to address the outlined challenges sufficiently. Specifically, the DT prediction could be used to change the FOWT operational settings to increase the lifespan of the MoLs by reducing the stress and, additionally, in the case where FOWT operational maintenance is in progress, serve as an early safety warning to operators.
For these reasons, to detect long-term drifts in the mechanical response of the MoLs, we will develop a DT able to accurately predict the behaviour of the healthy system to compare with the actual one. Moreover, as an effective tool to improve the lifespan of the MoLs and safety of FOWT maintenance operations, we will develop a secondary DT to accurately predict the near future axial tension. To these aims, we will leverage operational data collected from the world's first commercial floating-wind farm (the Hywind Pilot Park) to investigate the effectiveness of DTs for the prediction of the MoL axial tension for the two scenarios depicted above. The DTs will be developed using state-of-the-art data-driven methods (Shalev-Shwartz and Ben-David 2014; Oneto 2020) and results based on real operational data will support our proposal.
The rest of the paper is organised as follows. Section 2 summarises the problem we aim to address and the related available dataset. Section 3 describes the proposed modelling approach, based on DDMs, to address the problem summarised in Section 2 using the previously described dataset. Section 4 reports the experimental results showing the effectiveness of our proposal. Finally, Section 5 concludes the paper.

Problem and dataset description
In this section, we will summarise the problems under investigation and the available data that can be exploited for building models able to address them.
In this work, we aim to face two problems: (i) detect longterm drifts in the mechanical response of the MoLs and (ii) predict the near future MoL axial tension.
For what concerns the first problem, as described in Sect. 1, we will develop a DT able to infer the expected behaviour of the MoLs in healthy conditions to compare with the actual one. Drifts in the differences between the actual and the predicted (healthy) behaviour are an indicator of decay in the condition of the MoL (Coraddu et al. 2019). This DT will take as inputs the current status of the factors influencing the MoL behaviour (i.e., the motions of the turbine and the environmental conditions) and as output the MoL axial tension.
When it comes to the second problem, as described in Sect. 1, we will develop another DT able to predict the near future axial tension utilising past and present factors (i.e., not only the motions of the turbine and the environmental conditions, but also the MoL tension), in combination with selected forecasted future factors (i.e., the environmental conditions) influencing the MoL behaviour.
Note the fundamental difference between the two DTs. In the first case, we just use instantaneous information (i.e., turbine's motions and the environmental conditions) to predict the axial tension as we want to compare the current response of the MoL with the (predicted) one that it would have in perfectly healthy conditions. Whereas, in the second case, we also use the "context" information (i.e., past, present, and forecasted data), since we want to develop an high fidelity system which is able to predict (both in healthy and decayed states) the responses of the MoL for safety concerns. In fact, also using past and forecasted data in the first scenario would compromise our ability to develop a system able not to "follow" the decaying behaviour of the MoLs (i.e., learn the expected healthy behaviour) (Coraddu et al. 2019).
To develop our DTs, we will exploit the publicly available Hywind dataset (Catapult ORE 2020). The dataset describes the weather, motion, and mooring line axial tension of an SWT-6.0-154 turbine floating in depths of between 95 and 120m, 25 km off the coast of Scotland. The dataset is composed of 11 intervals of the turbine motion and response, in 30 min windows over the course of 2018. Table 2 describes the features of the data set. Table 3 reports the data collection periods together and the key operational information (i.e., the significant wave height H s and the peak waves period T p ).
A schematic of the turbine and the relevant sensor placement is seen in Fig. 1.

Proposed approach
In the proposed context, predicting the axial tension of the FOWT MoLs, a general modelisation framework can be defined, characterised by an input space X ⊆ R d , an output space Y ⊆ R b , and an unknown relation μ : X → Y to be learned (Shalev-Shwartz and Ben-David 2014;Hamilton 1994). For what concerns this work, X is composed by the the motions of the turbine and the environmental conditions (i.e., D e , D n , θ N , φ N , θ T , φ T , ψ T , v, and ψ dir in Table 2), while the output space Y refers to the axial tensions of the  FOWT MoL bridles (i.e., T ML 1 B 1 , T ML 1 B 2 , T ML 2 B 1 , T ML 2 B 2 , T ML 3 B 1 , and T ML 3 B 2 in Table 2). In this context, the authors define the model h : X → Y as an artificial simplification of μ. Now, the model h represents a DT of the FOWT MoL bridles. In our work, we aim to develop two kinds of DT (see Sect. 2): • the first one to infer the expected behaviour of the MoLs in healthy conditions to compare with the actual one; • the second one able to predict the near future axial tension utilising past and present factors.
The two DTs differ in the input space X and the output space (see Fig. 2).
For what concerns the first DT, the input space is composed by instantaneous information at time t (i.e., D e , D n , Table 2), while the output space is the axial tensions of the FOWT MoL bridles at time t (i.e., T ML 1 B 1 , T ML 1 B 2 , T ML 2 B 1 , T ML 2 B 2 , T ML 3 B 1 , and T ML 3 B 2 in Table 2). For what concerns the second DT, the input space is composed by past information during the time Table 2) and the forecasted information (t, t + + ] (i.e., D e , D n , and v, in Table 2) of the environmental conditions provided by a weather service, while the output space is the axial tensions of the FOWT MoL bridles at time t + + . − represents how much history of the different available data we want to exploit to make predictions. + , instead, represents the horizon of prediction. − is an hyperparameter for which an optimal value exists: too much history (too large − ) will make us face with the curse of dimensionality, while too little history (too small − ) will limit our ability to make accurate predictions (Shalev-Shwartz and Ben-David 2014; Oneto 2020; Hamilton 1994). + , instead, depends on the application: usually, the larger the better but the more we try to forecast deep into the future, the less accurate the prediction will be (Shalev-Shwartz and Ben-David 2014; Oneto 2020; Hamilton 1994). The model h, as described in Sect. 1 can be obtained with different kinds of techniques, for example, requiring some physical knowledge of the problem, as in PMs, or the acquisition of large amounts of data, as in DDMs. In this paper, we will use a state-of-the-art DDM for the reasons described in Sect. 1. Between the DDMs, it is possible to identify two families of approaches (Shalev-Shwartz and Ben-David 2014;Goodfellow et al. 2016). The first one, comprising traditional Machine Learning methods, needs an initial phase where the features must be defined a-priori from the data via feature engineering or implicit or explicit feature mapping (Shalev-Shwartz and Ben-David 2014; Zheng and Casari 2018;Shawe-Taylor and Cristianini 2004). The second family, which includes deep learning methods, automatically learns both the features and the models from the data Goodfellow et al. (2016). For small cardinality datasets and outside particular applications (e.g., computer vision and natural language processing), Deep Learning does not perform well, since they require huge amount of data to be reliable and to outperform traditional Machine Learn-ing models (Fernández-Delgado et al. 2014;Wainberg et al. 2016).
Machine Learning maps the problem of building the two DTs in a typical regression problem (Vapnik 1998;Shawe-Taylor and Cristianini 2004). In fact, ML techniques aim at estimating the unknown relationship μ between input and output through a learning algorithm A H which exploits some historical data to learn h and where H is a set of hyperparameters which characterises the generalisation performance of A (Oneto 2020). The historical data consist on a series of n examples of the input/output relation μ and are defined as D n = {(x 1 , y 1 ), . . . , (x n , y n )} where x ∈ X and y ∈ Y. For simplicity, we will indicate with y one of the elements in y, since predicting a series of targets is equivalent to make a model for each one of the targets (Shalev-Shwartz and Ben-David 2014).
In this paper, we will leverage on a machine learning model coming from the Kernel methods family called Kernel regularised least squares (KRLS) (Vovk 2013). The idea behind KRLS can be summarised as follows. During the training phase, the quality of the learned function h(x) is measured according to a loss function (h(x), y) (Rosasco et al. 2004) with the empirical error A simple criterion for selecting the final model during the training phase could then consist in simply choosing the approximating function that minimises the empirical error L n (h). This approach is known as empirical risk minimization (ERM) (Vapnik 1998). However, ERM is usually avoided in machine learning as it leads to severe overfitting of the model on the training dataset. As a matter of fact, in this case, the training process could choose a model, complicated enough to perfectly describe all the training samples (including the noise, which afflicts them). In other words, ERM implies memorisation of data rather than learning from them. A more effective approach is to minimise a cost function where the trade-off between accuracy on the training data and a measure of the complexity of the selected model is achieved (Tikhonov and Arsenin 1979), implementing the Occam's razor principle In other words, the best approximating function h * is chosen as the one that is complicated enough to learn from data without overfitting them. In particular, C(·) is a complexity measure: depending on the exploited machine learning approach, different measures are realised. Instead, λ ∈ [0, ∞) is a hyperparameter, that must be set a-priori and is not obtained as an output of the optimisation procedure: it regulates the trade-off between the overfitting tendency, related to the minimisation of the empirical error, and the underfitting tendency, related to the minimisation of C(·). The optimal value for λ is problem-dependent, and tuning this hyperparameter is a non-trivial task, as will be discussed later in this section. In KRLS, models are defined as where ϕ is an a-priori defined Feature Mapping (FM) (Shalev-Shwartz and Ben-David 2014) allowing to keep the structure of h(x) linear. The complexity of the models, in KRLS, is measured as i.e., the Euclidean norm of the set of weights describing the regressor, which is a standard complexity measure in ML (Shalev-Shwartz and Ben-David 2014; Vovk 2013).
Regarding the loss function, the square loss is typically adopted because of its convexity, smoothness, and statistical properties (Rosasco et al. 2004) Consequently, Problem (2) can be reformulated as By exploiting the Representer Theorem , the solution h * of the Problem (6) can be expressed as a linear combination of the samples projected in the space defined by ϕ It is worth underlining that, according to the kernel trick, it is possible to reformulate h * (x) without an explicit knowledge of ϕ, and consequently avoiding the curse of dimensionality of computing ϕ, using a proper kernel function K ( Several kernel functions can be retrieved in the literature (Scholkopf 2001; Cristianini and Shawe-Taylor 2000), each one with a particular property that can be exploited based on the problem under exam. Usually, the Gaussian kernel is chosen because of the theoretical reasons described in Keerthi and Lin (2003); Oneto et al. (2015) and because of its effectiveness (Fernández-Delgado et al. 2014;Wainberg et al. 2016). γ is another hyperparameter, which regulates the nonlinearity of the solution that must be tuned as explained later. Basically, the Gaussian kernel is able to implicitly create an infinite-dimensional ϕ, and thanks to this, the KRLS are able to learn any possible function (Keerthi and Lin 2003). The KRLS problem of Eq. (6) can be reformulated by exploiting kernels as where y = [y 1 , . . . , y n ] T , α = [α 1 , . . . , α n ] T , the matrix Q such that Q i, j = K (x j , x i ), and the identity matrix I ∈ R n×n . By setting the gradient equal to 0 w.r.t. α, it is possible to state that which is a linear system for which effective solvers have been developed over the years, allowing it to cope with even very large sets of training data (Young 2003). Note that the computational complexity of making the prediction for a new point, see Eq. (8), increases with the number of samples in the training dataset. Nevertheless, in this paper, we did not focus on this specific issue which can be addressed in many ways, but we focus our attention on the ability of making accurate predictions. In fact, the simplest way to reduce the computational requirements of the prediction phase is to sub-sample the training set (as one can observe also from the results of Sect. 4, we do not need so much data to achieve high accuracy) with some advanced method as reported in (Aupetit 2009). Another method is to use a different loss, such as the epsilon-insensitive one (Shawe-Taylor and Cristianini 2004), to reduce the number of points actually used by the training model.
The problems we still have to face is how to tune the hyperparameters for this approach (λ, γ , and − for the second DT) and how to estimate the performance of the final model. Model selection (MS) and error estimation (EE) deal exactly with these problems . Resampling techniques are commonly used by researchers and practitioners, since they work well in most situations, and this is why, we will exploit them in this work . Other alternatives exist, based on the statistical learning theory, but they tend to underperform resampling techniques in practice . Resampling techniques are based on a simple idea: the original dataset D n is resampled once or many (n r ) times, with or without replacement, to build three independent datasets called learning, validation, and test sets, respectively L r l , V r v , and T r t , with r ∈ {1, . . . , n r }, such that where h = A H (L r l ) is a model built with the algorithm A with its set of hyperparameters H and with the data L r l , and where M(h, V r v ) is a desired metric. Since the data in L r l are independent from the data in V r v , H * should be the set of hyperparameters which allows achieving a small error on a data set that is independent from the training set. Then, to evaluate the performance of the optimal model which is h * A = A H * (D n ) or, in other words, to perform the EE phase, the following procedure has to be applied: Since the data in L r l ∪ V r v are independent from the ones in T r t , M( f * A ) is an unbiased estimator of the true performance, measured with the metric M, of the final model . In this work, we will rely on Complete k-fold crossvalidation which means setting n r ≤ n k n− n k k , l = (k −2) n k , v = n k , and t = n k and the resampling must be done without replacement . Note that, in our application, we have a further constraint in terms of dependence in time between the samples. For this reason, when resampling the data form D n , we actually keep data of different periods in L r l , V r v , and T r t (Hamilton 1994). For what concerns the metric M that we will use in our paper, we will rely on the mean absolute error (MAE) which computes the average absolute distance between the prediction and the actual value to predict (Willmott and Matsuura 2005). Since in regression, it is quite hard to synthesise the quality of a predictor in a single metric, we will also rely on visualisation techniques like the scatter plot (Shao et al. 2017) (see Sect. 4).

Experimental results
This section is devoted to the presentation of the results of applying and testing the methodology described in Sect. 3 leveraging on the data presented in Sect. 2.
As first step, we have to report the hyperparameters ranges of the MS phase which is common to all experiments. The set of hyperparameters tuned during the MS phase is H = {γ, λ} chosen in H = {10 −4.0 , 10 −3.8 , · · · , 10 +4.0 } × {10 −4.0 , 10 −3.8 , · · · , 10 +4.0 }. For the second DT, we also investigated the effect of + ∈ {1, 2, 4, 8, 16, 32, 64, 128} and − ∈ {1, 2, 5, 10, 20}. For the second DT, − has also been tuned during the MS phase. All of the tests have been repeated 30 times, and the average results are reported together with their t-student 95% confidence interval, to ensure the statistical validity of the results.  . 3 First DT: scatter plot (real versus predicted) the real distribution, and the error distribution for the axial tension of the different MoLs and bridles

First digital twin
The first DT was designed to predict the MoL tension in healthy conditions to monitor the drift between the expected and true behaviour. This DT has to predict the instantaneous MoL tension from the current factors influencing the MoLs (i.e., the motions of the turbine and the environmental conditions) to accurately forecast the MoLs tension. Table 4 reports the MAE of this first DT for the different MoLs and bridles.
Moreover, Fig. 3 reports the scatter plot (Real versus Predicted) the real distribution, and the error distribution for the axial tension of the different MoLs and bridles.
From Table 4 and Fig. 3, it is possible to observe that the proposed DT is able to accurately predict the MoL tension in healthy conditions from the current factors influencing Fig. 4 Second DT: scatter plot (real versus predicted), real distribution, error distribution, and a piece of trend in time (real versus predicted) of the axial tension of the different MoLs and bridles with + fixed at 64 (more the 1 min) as it is a sufficient time horizon for providing a reasonable warning to avoid accidents and − fixed to the optimal value according to Table 7 Fig. 5 Second DT: scatter plot (real versus predicted), real distribution, error distribution, and a piece of trend in time (real versus predicted) of the axial tension of the ML 1 B 1 varing + , with − fixed to the optimal value according to Table 7 Table 4). At this point, it is important to remember that the dataset is composed of sporadic periods throughout the year (as described in Table 3) and such captures a wide range of operating con-ditions (which is clear from the metocean features reported in the same table). Despite this wide range of operating conditions, the DTs exhibit small error variance, indicating that the models performed consistently well across all conditions. Figure 3 allows us to better understand the quality of the developed DTs, further supporting the discussion raised from the tabular results, showing the behaviour of the predictions  6 Second DT: MAE of the different MoLs and bridles when fixing − = 20 (the value of − which most commonly is the optimal one according to Table 7) and varying + , when + is fixed at 64 (the most reasonable time horizon to be able to avoid accidents) and varying − near future MoL tension for safety purposes (e.g., control and maintenance activities). Table 5 shows how the MAE of the proposed model varies changing − for different values of + . Vice versa, Table 6 shows the MAE's variation of the proposed model as a function of + for different values of − . It is worth mentioning that Tables 5 and 6 actually contain the same information but presented in different ways for the sake of readability. Finally, Tables 7 reports the MAE when − is optimised during the MS phase for different values of + .
Analogously to the results shown in Sect. 4.1, we have reported a series of pictures to improve and enrich the ability of the reader to understand the real performance of the developed models, as in regression, it is harder to synthesise the quality of the model just with a single metric (as we did in Tables 5, 6, and 7). For this reason, Fig. 4 reports the scatter plot (Real versus Predicted), real distribution, error distribution, and a snapshot of the trend in time (real and predicted) of the axial tension for the different MoLs and bridles with + fixed at 64 s (more the 1 min). In fact, this represents a sufficient time horizon for providing a reasonable warning and fixing − to the optimal value according to Table 7. Then, Fig. 5 reports scatter plot (real versus predicted), real distribution, error distribution, and a snapshot of the trend in time (real and predicted) of the axial tension of the ML 1 B 1 varying + fixing − to the optimal value according to Table 7. Finally, Fig. 6 reports how the MAE varies in predicting the axial tension of the different MoLs and bridles when fixing − = 20 (the value of − which most commonly is the optimal one according to Table 7) and varying + and when fixing + = 64 (the most reasonable time horizon to be able to avoid accidents) and varying − .
From Table 5, 6, and 7 and from Figs. 4 and 5, it is possible to observe that • as expected from the theory, observing from Tables 5  and 6, an optimal value for − exists to provide enough information to learn the current behaviour of the MoLs and forecast the tension accurately; • as for the prediction horizon ( + ), the more we increase it, the higher is the prediction error (see Tables 5 and 6) which is again inline with one would expect. The performance's decrease caused by the increase in the prediction horizon is reported in Fig. 6. For small + , as one may expect, the errors are very low (MAE ≈ 5/10 kN). The model performance for the horizon we are interested for practical applications (≈ 1/2 min) is good enough (MAE ≈ 15/20 kN) as reported in Table 7; • there is a very good agreement between the real and predicted MoL tension as reported in Fig. 4. This demonstrates that the models were effective at learning the short-term behaviour of the MoLs.
In conclusion, results indicate that this second family of DTs proved to be an effective tool for operators to utilise where a time-sensitive response to the state of the MoLs is required.

Conclusions
The number of installed FOWTs has exponentially grown in the last decade, quadrupling the total installed capacity. This growth is expected to continue in the next decade leading to an increasing need to address main challenges for FOWT projects: monitoring the system's integrity, extending the lifespan of the components, and maintaining FOWTs safely at scale. Effectively and efficiently addressing these challenges would unlock the wider-scale deployment of FOWTs. For these reasons, in this work, we developed two DTs to predict the MoL tension of an FOWT exploiting state-of-theart data-driven methods and leveraging the data coming from the Hywind Pilot Park to test our proposals. The first DT was able to predict the MoL tension under healthy conditions to monitor the drift between the expected and true behaviour. This DT represents an effective solution to detect long-term drifts in the mechanical response of the MoLs, and accurately predict the behaviour of the healthy system to compare with the actual one. The second DT utilised past, present, and forecasted data to predict the near future MoL tension for safety purposes. In particular, we were able to achieve good results (with an error of around 15 [kN]) for a forecast horizon of approximately 1-2 min, which is enough to better operated the FOWTs generate early safety-related warnings. By changing the FOWT operational settings, according to this second DT prediction, operators can increase the lifespan of the MoLs by reducing the actual mechanical stress and, additionally, in the case where FOWT operational maintenance is in progress, the prediction from the DT can serve as early safety warning to operators.
As a concluding remark, it is important to stress that this work is just a preliminary step forward validated on a limited amount of data. The commercial use of FOWTs is still in its infancy, and the available dataset only captures sporadic periods throughout the year. For this reason, more tests with a more extensive set of data need to be conducted. Moreover, these DTs need to be integrated into real monitoring, control, and safety systems to validate the potentiality of the models.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.