Edinburgh Research Explorer Simulating Extreme Directional Wave Conditions

: Wave tank tests often involve simulating extreme wave conditions as they enable the maximum expected loads to be inferred: a vital parameter for structural design. The deﬁnition, and simulation of, extreme conditions are often fairly simplistic, which can result in conditions and associated loads that are not representative of those that would be observed at the deployment location. Here we present a method of deﬁning, simulating at scale, and validating realistic site-speciﬁc extreme wave conditions for survival testing of wave energy converters. Bivariate inverse-ﬁrst order reliability method (I-FORM) environmental contours deﬁne extreme pairs of signiﬁcant wave height and energy period ( H m 0 – T E ), while observed extreme conditions are used to deﬁne realistic frequency and directional distributions. These sea states are scaled, simulated and validated at the FloWave Ocean Energy Research Facility to demonstrate that the site-speciﬁc extreme wave conditions can be re-created with accuracy. The presented approach enables greater realism to be incorporated into tank testing with survival sea states. The techniques outlined and explored here can provide further and more realistic insight into the response of offshore structures and devices, and can help make important design decisions prior to full-scale deployment.


Introduction
Offshore structures and devices, including wave energy converters (WECs), are tested in wave tanks to assess behaviour, performance and loads at scale prior to deployment at sea. For survivability, it is paramount to identify the expected extreme conditions, which often initially focusses on the identification of extreme significant wave heights. This is achieved using extreme value distributions to extrapolate the observed data to the desired return period and is suggested by Goda [1] to be the first step in coastal structure design. Extreme sea states are hence often used in such test programmes and enable the largest expected loads to be inferred, which can feed back into design and help validate numerical prediction.
When defining extreme wave conditions, it is common to extrapolate extreme significant wave height, H m0 , values as detailed in Coles [2], Caires [3], Teena et al. [4]. To define and simulate sea states, however, appropriate wave period values are also required. These can be obtained using joint probability methods such as the inverse-first order reliability method (I-FORM) technique (implemented in Haver and Nyhus [5], Winterstein et al. [6], Berg [7]) recommended by the EquiMar Protocols [8] for defining extreme conditions for WECs. These methods provide environmental contours relating to certain probabilities or return periods. The I-FORM method has been implemented in Doherty and Folley [9] to define extreme significant wave height-period (H-T) combinations for Aquamarine's Oyster WEC, which is used a basis to create extreme wave spectra for tank testing. In this case, Bretschneider spectra were used in all cases, and it is assumed that they have been 1.
Defining extreme combinations of significant wave height and period relating to specified return periods 2.
Calculating, and subsequently incorporating, observed spectral and directional complexity in the defined extreme sea states 3.
Sea state correction and validation in the wave basin An approach is presented in this paper which, with the aid of an example, deals with all of these considerations: from site data to tank validation. The method is demonstrated using data from EMEC's Billia Croo wave site, and the resulting sea states are re-produced at the FloWave Ocean Energy Research Facility (FloWave). Section 2 covers points 1 and 2 above: the definition of realistic extreme sea states. In Section 2.1, I-FORM contours are used to define extreme combinations of significant wave height and energy period, H m0 and T E . Section 2.2 details the method used to calculate sea state directionality, and an approach is presented to use this information in the extreme sea state definition. Simulation of the conditions at scale is dealt with in Section 3, with Section 3.1 discussing the practical issues associated with scaling extreme wave conditions. Sea state re-creation, correction and validation, along with final results are presented in Section 3.2. Issues with scaling extreme seas in wave tanks are further discussed in Section 4, as well as assessing the validity of the sea state creation methodology described in this paper. Concluding remarks are made in Section 5.

The European Marine Energy Centre: Site and Data Overview
Based in the Orkney Isles, UK, EMEC is a full-scale grid connected test facility for wave and tidal energy devices. Over the last 14 years, EMEC has gained vast experience in both device installation and resource measurement at these sites. Located on the west coast of the Orkney Isles, EMEC's Billia Croo wave site is exposed to distant North Atlantic storms from the west/north-west, in addition to more locally generated wind-seas.
Half-hourly sea states collected using Datawell's Directional Waverider buoys have been made available for use, providing good quality directional wave data from January 2010 to January 2014.
Frequency spectra and directional Fourier coefficients have been provided, with frequency spectra used directly to compute the statistical parameters, H m0 and T E , used throughout this work, via Equation (1).
where m n are spectral moments of the frequency spectrum, S, i.e., m n = f n S( f )δ f . The frequency spectra, in combination with directional Fourier coefficients, have been used to estimate the directional distributions required in this analysis. This process is described in Section 2.2. The calculated mean frequency and directional spectra for the site, along with the mean weighted directional spreading function (DSF mean ) is shown in Figure 1. Mean values of H m0 and T E are 2.05 m and 8.33 s respectively, with the most energetic recorded sea state having a H m0 of around 11.5 m, and T E of 15 s.

The FloWave Ocean Energy Research Facility
All experiments presented in this paper were carried out at the FloWave Ocean Energy Research Facility, Edinburgh, UK. The 2 m depth facility consists of a 25 m diameter wave basin encircled by 168 active-absorbing wavemakers, enabling waves to be generated from all angles. Sea states can therefore be created with arbitrary directional characteristics, and FloWave has been shown in to be able to effectively re-create complex multi-directional multi-modal spectra [17,18]. This allows the directional complexity of observed ocean conditions such as those defined in Section 2 to be effectively simulated.

Defining Extreme Directional Wave Conditions
This section deals with the extreme sea state definition, and aims to define a set of realistic directional spectra that are representative of extremes at the site. Section 2.1 shows how statistical parameters, H m0 and T E , can be appropriately associated with extreme events at specified return periods. In Section 2.2, these values are used to create realistic extreme sea states, which are based off observed directional spectra associated with storms at the site.

Significant Wave Height
There are a variety of approaches for inferring the largest expected significant wave heights in a given time frame. For the initial peak identification a variety of methods can be employed, including the commonly used block maxima (BM) or peaks over threshold (POT) approaches [3], while an r-largest or average conditional exceedance rate (ACER) [19,20] methodology could equally be applied. It is recommended by Caires [3] that if the dataset is relatively small (a few years) then a POT approach should be taken. Additionally, if multiple storm events occur within relatively small time-frames, it is recommended in Teena et al. [4] to use this POT method to include all the extreme events in the analysis. As this is the case for the dataset available a POT approach will be used in this work. Figure 2 shows the points identified as POT-valid points for different thresholds. These points must be independent [3] and are ensured here by only considering the maximum points over the threshold, in the 24 h preceding and after the considered event.   BM derived data (such as annual maxima) require fitting with a generalised extreme value (GEV) distribution [3], which consists of well known different types [21]: Gumbel (type I), Fréchet (II) or reversed Weibull (type III) distributions. For threshold derived extreme data e.g., POT, Leadbetter [22] showed that a generalised pareto distribution (GPD) should be used, and as such, a POT-GPD analysis is deemed most appropriate for the dataset available. Similarly to GEV distributions there are different types of GPD depending on the value of the shape parameter ξ: exponential (type I), Pareto (type II) and a form of beta-distribution (type III). The equations for the GPD is shown in Equations (2) and (3).
where: ξ = shape parameter, σ u = scale parameter, λ u = rate, u = threshold. Type I: ξ = 0, Type II: ξ > 0, Type III: ξ < 0 Choosing the appropriate threshold for the POT analysis is a trade-off between bias and variance and must be small enough to ensure there are sufficient points to determine the GPD parameters, but also must be large enough so that the data is truly capturing extreme events, and as such converges to a GPD [23]. Based on the stability parameter defined in Caires [3] a threshold of 4 m was eventually chosen, considered large enough as to only consider extremes but include a sufficient number of points, while being associated with a 'stable' region, and with a high associated r 2 value. The final GPD distribution with the 4 m threshold is shown in Figure 3, with the observed data overlaid.

Energy Period
As mentioned in Section 1, to define extreme sea states effectively, appropriate wave period values are required. Log-normal distributions have been shown to effectively describe wave period distributions by Burrows and Salih [24], and therefore a parametrised conditional log-normal distribution has been fitted to values of energy period similar to Berg [7]. The parameters have only been fitted for H m0 values roughly corresponding to those above the POT threshold, under the assumption that there may be a different relationship between wave period and extreme wave heights than for the total dataset. Bin sizes for the parametrisation are also increased for large values of H m0 , to ensure there is enough data in each bin to get meaningful standard deviations. The resulting fit is shown in Figure 4, with the log-normal parametrisation described by Equations (4)- (7).

Combined Conditions
Using the I-FORM method the two probability distributions can now be used to create contours relating to specified return periods. This is achieved by defining β to be the exceedance value of a standard normal variable relating to the specified return period, m. Two standard normal variables U 1 and U 2 are then defined along a circle so that their radius is equal to β. The standard normal probability of exceedance corresponding to the values of U 1 and U 2 can then be used to find the equivalent values of H m0 and T E via the obtained distributions. As in Berg [7], this is described formally in Equations (8)-(10): It is typical to inflate β to account for approximating true stochastic response by its median value.
From Winterstein et al. [6] and Berg [7] it is suggested to inflate β using β = β/ 1 − α 2 0 , with α 0 values of between 0.1 and 0.2. A value of 0.15 has therefore been chosen for this analysis, with the resulting contours for return periods of 1, 10, 50 and 100 years shown in Figure 5. The choice of extreme H m0 -T E combinations for tank testing, or numerical modelling, is dependent on the nature of the resulting probability distributions and the sensitivities of the device being tested (e.g., [9])). As wave energy converters are typically very sensitive to wave frequency (see [16,[25][26][27]) it seems advantageous to choose a range of T E values, enabling exploration of system response to extreme conditions at various dominant forcing frequencies. Device and mooring design will influence the decision significantly. Typically, however, conditions of large wave height but relatively low wave period values may typically correspond to resonant conditions of the device and thus maximum motions, or may be likely to cause steepness induced snatch loads in mooring lines (e.g., [28]). Maximum mooring extension may occur from large surge loads within high T E sea states, while the largest impact loads may occur when state power or significant wave height is at a maximum.
Various points have been chosen for 1:10 and 1:50 year extremes, to serve as an example cases for further analysis. The 1:10 year extremes are used as example extreme conditions which prototype devices may be subject to, with short deployments of a couple years. Extreme 1:50 year events will be more representative of full deployment lifetimes in the region of 20 years or more, as recommended by the International Electrotechnical Commission: IEC 61400-3 for offshore wind turbines [29], IEC 62600-2 for marine energy converters [30], and used in Ruiz [31] and Valamanesh et al. [32]. Points have been chosen corresponding to locations on each contour for maximum energy period, maximum power (normally maximum H m0 as well), and minimum T E above a H m0 threshold (10 m for 1:10, and 11 m for 1:50). The T E conditions specified, in combination with the extreme contours, define the combined extreme H m0 -T E values. These are shown in Table 1 and are highlighted in Figure 6. In this figure the closest unique observed data points are also highlighted, whose spectra are scaled to the desired H m0 -T E values in Section 2.2.2.    (Table 1), and arrows indicate which observed sea state is chosen to match each specified extreme condition.

Reconstructing Directional Spectra from Buoy Data
The buoy data available for the Billia Croo wave site does not explicitly include directional spectra, yet has directional Fourier coefficients which can be used as the basis for reconstruction. There are variety of methods available for the reconstruction of directional spectra, all of which have different assumptions and formulations, which results in DSFs which differ from each other significantly.
It is suggested in Nwogu et al. [33], Benoit [34], Benoit et al. [35], Kim et al. [36] that the maximum entropy principle (MEM2/MEP) is the most reliable for directional spreading estimates from single-point measuring systems like buoys. Despite this, the solution to the MEP is non-linear, computationally challenging and does not always converge. This led Hashimoto [37] to propose using Newton's method of local linearisation to help the convergence and Kim et al. [36] to suggest approximate solutions to the problem. To avoid the convergence issues experienced with both the full MEP solution and the linearised approach, the aforementioned approximate solutions [36] have been used for this work, which have been found to give very similar directional spectra outputs to full solution equivalents.

Applying Scaling-Shifting Procedure to Observed Extreme Spectra
To define the extreme directional spectra, the 'closest' observed sea state to the defined extreme H m0 -T E points Table 1 are identified, using the M number defined in Equations (11) and (12). If a given observed sea state is the closest to multiple desired extreme points, distinct points are chosen that minimises total M number. Due to the relatively short dataset and low abundance of very high power sea states, there are conflicts for extreme conditions 1 and 2, along with 4 and 5. Forcing these sea states to be based on different observations allows the inclusion of differing conditions, highlighted in Figure 6.
The selected sea states do not match the specified H m0 and T E values, and as such go through a scaling procedure. This required an amplification and a shifting in frequency. Noting that (13) can be used to get the correct H m0 by means of applying a single spectrum amplification factor. To obtain the correct T E , the whole spectrum is first interpolated to give a large number of frequency values. The spectrum is then shifted on the frequency axis until the desired T E is obtained. An iterative, rather than analytic, method was used here as T E is a function of , to find the shift, the desired energy distribution would need to be known, which remains a function of the unknown frequency shift itself and the known spectral shape. The results from this scaling-shifting procedure for both the frequency and directional spectra are shown in Figures 7 and 8 respectively. From Figure 8, it is evident that the resulting (observed) directional distributions are highly frequency dependent and would not have been obtained from a purely parametric approach to the sea state definition.  Figure 7. Frequency spectra for sea states identified in Figure 6, scaled to desired extreme conditions (Table 1). M value overlaid corresponding to Equation (11).  Figure 8. Directional spectra for identified extreme sea states, scaled to desired extreme conditions. 1) to 6) correspond to sea states 1 to 6 defined in Table 1.

Simulating Extreme Directional Wave Conditions
This section covers the experimental aspect of this work: scaling, generating and validating the sea states. The defined sea states are scaled in Section 3.1, where difficulties are highlighted relating to a combination of depth and wavemaker constraints. Section 3.2 details the generation, measurement, and correction procedure and shows the final validated extreme spectra.

Scaling Sea States for Tank Tests
Froude scaling is used to scale sea states in order to obtain the correct ratio between inertial and gravitational forces, which are dominant for ocean waves. Ideally, all environmental conditions of a wave test are therefore Froude scaled, but there are constraints to consider, related to the sea state in question along with the characteristics and limitations of test facilities. Considering frequency dependent wave height limits of the FloWave facility, as well as absolute frequency limits of wave generation, it was found that the largest scale that the defined extreme sea states can be generated at effectively are 1:61. The desired scale factor is 1:26 (52 m site and 2 m tank), and as such, the sea state needs to be created at a non-depth-ratio.
In simple parametric cases, H m0 and T E /T p values are typically used for scaling, which effectively results in Froude scaled frequency spectra, S( f ). However, if these are scaled and reproduced at any ratio other than the depth ratio, there are wavelength and power discrepancies to consider, differing from the desired Froude scaled values. For all but deepwater conditions, non-depth ratio scaling makes it is impossible to Froude scale all elements of a sea state [38]. Table 2 details how these errors can be calculated. In this table the subscript s refers to Froude scaled values, f s to full-scale values, and o to obtained values i.e., those that are not directly scaled and where errors arise. Alternative scaling approaches are also presented, which allow preservation of energy distributions across wavenumber, S(k), and wave power across frequency, p( f ).
The 1:61 Froude scaled extreme frequency spectra are shown in Figure 9, with the power and wavelength errors incurred by non-depth-ratio scaling. It is interesting to note the changing sign of the frequency dependent wave power error at around 0.4 Hz: low frequency components have an over-representation of power while higher frequency components display lower power than desired. Above around 1.1 Hz (tank scale) wave frequencies that were in deep water at full scale remain to be in deep water at tank scale and as such power and wavelength are Froude scaled effectively. Importantly, total sea state power, P = p( f )d f , does not Froude scale in this case, with the total discrepancy dependent on the spectral shape. For this study, to demonstrate the ability to re-create the extreme sea states the S( f ) scaled spectra have been used. The difference between the different scaling approaches are highlighted further in Appendix A: Figure A1, whereby the six extreme spectra are scaled using the three approaches described. This highlights the difference between the desired spectra at that scale and those obtained using the different scaling methodologies. Table 2. Effect of various scaling approaches on other parameters (guaranteed Froude scaled parameters shown in blue).

S( f ) Scaling S(k)
Scaling  Figure 9. Extreme sea frequency spectra Froude scaled at 1:61. Resulting relative frequency dependent wavelength and power errors at this scale shown above. Dashed blue line indicated the frequency where power errors change sign (at this scale and depth ratio). 1) to 6) correspond to sea states 1 to 6 defined in Table 1.

Generating, Measuring and Validating Extreme Directional Sea States
The Froude scaled directional spectra, E( f , θ), have been created in the FloWave basin, prior to measurement and correction based on the incident spectrum. To generate the directional sea states, the single-summation method [39] has been used, minimising phase-locking between frequency components and thus increasing the ergodicity of the wave field. Sea states have been created with a repeat time, T, of 2048 s and a run-time of 2100 s. This test length provides roughly a thousand waves, suggested in [40] to be the typical length of a storm, and should provide a good representation of individual extreme wave heights. A directional array of eight resistance type wave gauges (shown in [18]) has been deployed in the centre of the FloWave basin.
The Single-summation PTPD Approach with In-line Reflection (SPAIR) method [18] has been used to isolate the incident and reflected directional spectra. The method has been applied over the final 2048 s of the 2100 s measurement, corresponding to the defined sea state repeat time. This ensures that high frequency waves have had sufficient time to propagate to the measurement area, and that the assumption of stationarity is more valid. In addition, taking the measurement over exactly one repeat time ensures that every generation frequency component has undergone an even number of cycles and spectral leakage should be minimised. Due to slight under-generation, the incident wave spectrum is corrected using frequency dependent amplitude based correction factors as defined in Equation (14). This is repeated until mean spectral errors, (Equation (15)), are below 5%, which happens for all six sea states after a single iteration.
The final incident frequency spectra are shown in Figure 10, along with the desired spectrum and the calculated frequency dependent reflection coefficients. Corresponding weighted DSFs are shown in Figure 11, while the total directional spectra are displayed in Figure 12. It is evident that the measured sea states are very close to the desired and the directional and frequency distributions have been created to a high degree of accuracy. Discrepancies in frequency spectra are extremely low, as demonstrated by Figure 10, while the shape of the incident directional spreading match very closely with the desired distributions (Figures 11 and 12). Reflection coefficients are generally low (below about 1 Hz), and discrepancies are well defined. This demonstrates that the complex defined energy distributions have been well produced, and that realistic directional extreme sea states are possible to simulate.  Figure 10. SPAIR incident frequency spectrum outputs for the extreme sea states, compared with desired. Shown with calculated frequency dependent reflection coefficients. 1) to 6) correspond to sea states 1 to 6 defined in Table 1.  Figure 11. SPAIR method weighted DSF outputs for the extreme sea states, compared with desired. 1) to 6) correspond to sea states 1 to 6 defined in Table 1. Figure 12. SPAIR method directional spectrum outputs for the example extreme spectra. Total output directional spectra shown, including reflections. Energy densities are shown by the colour scale and are normalised by the maximum of the desired spectrum for each sea state. 1) to 6) correspond to sea states 1 to 6 defined in Table 1.

Discussion
This paper demonstrates a process of simulating realistic extreme directional sea states using wave buoy data and observed extremes: from site data to validation at tank scale. There is some interesting discussion to be had over the validity of the presented process, with particular issues evident when applied to the current dataset and re-created in the FloWave facility. This includes the H m0 -T E contour creation, the scaling-shifting procedure at full scale, and the Froude scaling trade-offs highlighted in Section 3.1. These issues are expanded upon below.

H m0 -T E Contour Creation
In Section 2, four years of buoy data was used to define extreme distributions of H m0 and combined H m0 -T E . This is not enough data to effectively capture inter-annual variability, and ideally a larger dataset would be used, particularly for gaining confidence in extrapolating 1:50 year events. This would reduce the large uncertainty in the return values (Figure 3). To check the obtained values were reasonable, comparison has been made to Lawrence et al. [41] where an extreme H m0 distribution is obtained for the site using a 20-year hindcast model. 1:10 and 1:100 year values were found to be around 10.7 m and 12 m respectively, which compares favourably to the 11 m and 12.2 m found in this analysis. This provides increased confidence in the presented analysis technique as these hindcast H m0 values should be relatively accurate.
I-FORM contours were created in order to obtain combined H m0 -T E values suitable for extreme sea state definition. In this process a parametrisation of the conditional probability of T E on H m0 is carried out, with the fit of log-normal parameters shown in Figure 4. Figure 5 shows the resulting contours, and it is evident that some observed values are outwith the 1:100 year contour. This is due to the requirement to parametrise the conditional distributions of T E , and again is somewhat a result of having a small dataset. A different (and likely more accurate for low H m0 values) parametrisation would have been found if more data was available.
It is worth re-visiting the fact that the H m0 and T E values these contours are derived from have been based on 30 min data samples. If these were calculated using longer time periods, such as the three hours sometimes used (and typically used in numerical models), then the observed extreme values would be slightly reduced. This would consequently provide smaller estimates of the expected extreme wave conditions. Assessing the current dataset, using three hour moving averages causes the maximum H m0 to reduce from 11.36 m to 10.77 m, and maximum T E down from 17.39 s to 16.79 s (for 3 hour blocks this becomes 10.34 m and 16.27 s respectively). This demonstrates the potential discrepancy caused by the sample length choice, whilst also highlighting one of the inherent difficulties in analysing ocean wave data. A compromise must be made so that the assumption of stationarity is valid, while simultaneously obtaining enough samples to have reasonably small uncertainty in spectral and statistical estimates. Whether 30 min or longer samples are more appropriate for extreme value analysis is therefore a matter of debate. For this work, however, the 30 min sample length was chosen because it is widely suggested (e.g., Euan et al. [42] and Alvarez-Esteban et al. [43]) to ensure that both the assumption of stationarity and ergodicity are reasonable.

Magnitude of Transformation
Observed spectra were used to define the test conditions. In order to ensure they matched the defined H m0 -T E values they were amplified (H m0 ) and shifted (T E ). The resulting alteration to observed spectra are shown in Figure 7. This process is clearly only reasonable for small shifts in H m0 and T E . It is noted, however, that variation in observed spectral shapes over short time-frames (30 min) are many times larger than the ones proposed here, suggesting that the magnitude of the applied transformations may be reasonable, if not physically representative of wave energy transfer mechanisms. This change in spectral shapes for the 3 h prior to each identified extreme point is shown in Figure 13. From this figure, it is evident that there are quite significant changes in both H m0 and T E over relatively short time-frames. This aids somewhat in putting the proposed changes to the observed spectra into perspective.
This extreme value analysis aims to identify and create improbable conditions, yet care must be taken to ensure the test spectra are not impossible. Despite the example transformations being reasonable, it must be ensured that changes are not too large and do not breach any physical laws. In particular, if increasing H m0 significantly and/or reducing T E , then the new steepness must be calculated to ensure that the proposed spectrum is still physically maintainable and wave steepness limits are not breached. Equally, shifting a spectrum to have a much larger wave period may also represent an impossible energy distribution given the fetch length available. Therefore, in addition to checking steepness, it seems that a limit should be placed on the overall transformation.
From observation of the transformed spectra in Figure 7, it is suggested that the required scaling for the 1:10 year example events are very minor, well within typical spectral uncertainty, and hence represent a valid change. Relative to the large scale variations over short time frames, it also appears that scaling sea states for the 1:50 year events may also be reasonable. However, some of the large changes in the energy period required for some of these sea states, particularly Extreme 6, make this more questionable. Despite these concerns, the applied frequency shift seems to still create plausible conditions, due to the occurrence of many sea states with larger T p values, up to 25 s (0.04 Hz) in the dataset. 0 hours before 0.5 hours before 1 hours before 1.5 hours before 2 hours before 2.5 hours before Figure 13. Frequency spectra change for the three hours prior to the extreme sea state identified. 1) to 6) correspond to sea states 1 to 6 defined in Table 1.
On the basis of these observations, and noting the associated M values from Equation (11), it seems reasonable to suggest that any M value of less than 0.05 would be a very minor required change and hence clearly acceptable. It is suggested, based on the magnitude of the changes relative to observed variations over half-hour periods, that M values of less than 0.1 are also potentially acceptable, as long as both general and site-specific physical laws are not broken.

Spectral Variability
The methodology applied in this paper gives emphasis on preserving the form of observed directional spectra. The idea behind this approach is that it ensures resulting spectral definitions of extreme events are both possible and characteristic of those in the defined H m0 -T E region. For certain values of H m0 and T E the possible and likely range of spectral shapes, both in frequency and direction, may vary significantly i.e. there may be a wide range of underlying ocean conditions that can provide the same statistics. For other values, particularly those which represent highly energetic sea states, the range of possible conditions will be smaller. Inherent limits on steepness and the necessity to have large fetch lengths will impose limits on plausible predominant propagation directions, directional spreading and peak frequencies of the wave system. It is both important, and interesting, to explore the likely spectral variation, as to put into context the work presented.
The expected variation for very extreme events at the site is difficult to explore in this instance as the dataset is relatively small and observations of highly extreme conditions are limited. To explore the expected spectral variation, two examples are presented where there a number of observation points, one for low T E conditions and one for high T E . These examples allow for identification of the expected spectral variability for different wave conditions, but also helps put into perspective the emphasis placed on individual spectral forms throughout this work.
All points with M < 0.05 to the specified conditions are scaled to the desired values and compared. The two points identified are shown in Figure 14, along with the observed points within the M < 0.05 region. These desired points were chosen as the differing T E values will relate to different underlying physical wave conditions (and there are a reasonable number of observation points to analyse).  The spectra and frequency-averaged DSF's for the high and low T E examples are shown in Figures 15 and 16 respectively. It is evident from these figures that, for the high T E case, observed frequency spectra are reasonably similar, with low peak frequencies and narrow spectral bandwidths. Also suggested from the selected data is that large period storms are likely to originate from a particular direction (north-west). The necessity to have a large fetch length also means that there is essentially a spatial filter applied which limits, and practically standardises, the observed directional spreading. For the low T E example this is not the case. Frequency spectra are similar to each other, with increased peak frequencies and spectral bandwidths, yet much greater variation is observed in the directional spreading than observed in the high T E case. The lower energy periods mean that these sea states do not have the same fetch requirements, and comprise of more locally driven, and very steep, 'wind-driven' seas which can originate from a larger variety of directions. The results of this analysis demonstrate that, although the method presented ensures likely and realistic directional wave conditions to be used, it does not guarantee that all of the possible 'types' of spectral conditions are identified. This is potentially an issue that could be resolved with a larger dataset. Increased extreme observations would allow for the identification of different spectral types related to specified H m0 -T E definitions. With a more extensive dataset, or by choosing less extreme H m0 -T E combinations which are more abundant, the different associated spectral types can incorporated into a tank testing program. This approach, when achievable, will improve the presented methodology and further aid de-risking WEC's in realistic extreme wave conditions.

Froude Scaling Issues
As mentioned in Section 3.1, it was not possible to Froude scale all elements of the extreme sea states in FloWave. This is because Froude scaling is never fully achievable for intermediate water conditions when depth-ratio scales are not attainable. Froude scaled S( f ) & E( f , θ) were used to demonstrate the process of re-creating, correcting and validating sea states in Section 3.1, however, this was essentially an arbitrary choice. If the device being tested is very sensitive to wavelength or wave steepness (e.g., attenuator type WEC) then it is probably more appropriate to scale wavenumber spectra, S(k), ensuring the correct distribution of energy across the different wavelengths. A point absorber WEC may be mostly responsive to different frequencies of waves and as such S( f ) scaling is favourable, while p( f ) scaling may provide the best representation for power production tests rather than testing in extreme conditions. The outcomes of extreme tank tests feed back into design of offshore structures and renewable energy devices, and can affect the design and total cost significantly. It is, therefore, important to maximise both the relevance and accuracy of tests and results so that design decisions are well informed. Where possible it is recommended to use facilities where depth-ratio scaling can be achieved. This, however, is often problematic for the creation of extreme sea states due to wavemaker limits and limited ranges of tank depths. In such scenarios, as presented here, where depth-ratio scaling is not achievable, it is suggested that both S( f ) and S(k) scaled sea states are used in testing and the results compared. The worst case (e.g., maximum load case) should then be considered for design.

Conclusions
Extreme wave conditions are simulated in physical model testing to assess survivability of offshore structures and devices, and can have a large influence over design and eventual cost. Tank tests that simulate extremes should therefore ensure the conditions generated are representative of the deployment site, and are accurately re-produced. In this paper we present a method of achieving this, using site data, inverse-first order reliability method (I-FORM) contours, and observed extreme events to define realistic directional sea states for the Billia Croo wave site. Six differing extreme conditions are chosen relating to 1:10 and 1:50 year events, and are simulated at scale in the FloWave basin. The whole process from site data to tank scale validation is presented. Careful consideration must be given at each of the detailed stages, allowing scale reproduction of realistic directional extreme sea states to a high degree of accuracy. The inclusion of such conditions in tank testing programmes allows for increased, and more relevant, learning to occur prior to full-scale development and deployment.