1. Introduction
In numerical and physical experiments, nonlinear focusing has been identified as a possible mechanism for rogue wave generation. Its role for the formation of rogue waves in the real ocean, however, remains unclear. Nonlinear focusing was first described by Benjamin and Feir (1967): due to a modulational instability, a uniform wave train in deep water may dissolve into groups, which subsequently produce one large wave that is growing at the expense of the surrounding waves. Ten years later, Lake et al. (1977) demonstrated this so-called Benjamin–Feir instability to work in wave tank experiments. Onorato et al. (2004) showed—also in a wave flume—that as the instability develops, the rogue wave occurrence frequency increases. The Benjamin–Feir instability has therefore been proposed as a possible explanation for the formation of rogue waves (Janssen 2003).
Alber (1978) showed analytically that the requirement for a Benjamin–Feir-type instability to occur in a random wave field in deep water, is a sufficiently low bandwidth and small directional spreading of the wave spectrum. This has been substantiated numerically (Yuen and Ferguson 1978) and experimentally (Stansberg 1995). Waseda et al. (2009) showed in wave tank experiments that narrow-banded conditions favor an increased rogue wave occurrence.
A narrow and peaked frequency spectrum indicates regular wave conditions and more pronounced wave groups than a broad-banded spectrum. The bandwidth parameters described above are thus also a measure of wave groupiness (Holthuijsen 2007). High Qp values and low ε and ν values are therefore expected in swell-dominated sea states or wave fields including wave groups.
Besides a narrow bandwidth, another prerequisite for the modulational instability to occur is a narrow directional spreading. Numerical calculations have shown that unidirectional seas may result in much higher waves than predicted by second-order theory (Gibson and Swan 2006). This has been interpreted as an increase in the probability for rogue wave occurrences due to the modulational instability in unidirectional waves (Gramstad and Trulsen 2007). A broadening of the spectrum on the contrary leads to a reduction in rogue wave occurrence probability (Onorato et al. 2002; Waseda 2006; Janssen and Bidlot 2009). Latheef and Swan (2013) reported from a laboratory study that higher-order nonlinear effects can be important in directionally spread seas as well, provided the waves are sufficiently steep and not too short crested. In Onorato et al.’s (2009) laboratory experiments in a 3D wave basin, waves with a directional spreading larger than σΘ = 15° showed an almost Gaussian distribution, thus, a rogue wave occurrence frequency close to the expectations of second-order theory. In the real ocean, this would correspond, for example, to prevailing swell conditions. In pure wind-wave conditions, a high BFI does not necessarily indicate the presence of a modulational instability (Waseda 2006).
Gramstad and Trulsen (2007) stated that the probability of large waves depends on both the spectral bandwidth and the directional spreading, and they recommended including both parameters in the rogue wave probability prediction.
Meanwhile, the validity of the modulational instability has been extended to intermediate water depths (0.8 < kh < 1.36) (e.g., Toffoli et al. 2013) It was confirmed numerically that the destabilization of a wave train due to oblique perturbations is possible in intermediate water (Fernandez et al. 2014). Toffoli et al. (2009) found that unlike in deep water (kh > 1.36), the deviation from Gaussian statistics in intermediate water is independent of the directional spread. Thus, in intermediate water, modulational instabilities may also occur for realistic directional spreading (Karmpadakis et al. 2019).
Dysthe et al. (2008) argued that the deep-water preconditions for a modulational instability are unlikely to occur in wind waves conditions. Based on hindcast data from the Sea of Japan, which include typhoon conditions, several authors found that narrow wave spectra that favor a modulational instability are possible in some parts of the ocean (Waseda et al. 2009; Tamura et al. 2009). A number of authors found, based on measurement studies, that the occurrence frequency of rogue waves was dependent on the spectral bandwidth in the underlying wave field. Karmpadakis et al. (2020) described, based on radar measurements in the North Sea, a reduction in wave heights with an increase in spectral bandwidth. Cattrell et al. (2018) investigated buoy data from the U.S. coast and found that the spectral bandwidth parameters of rogue seas displayed different probability distributions than those estimated from nonrogue seas. Christou and Ewans (2014) stated, based on radar and laser data from the North Sea and other locations, that the spectral bandwidth might be an indicator for distinguishing rogue waves from high nonrogue waves. Most recent findings by Häfner et al. (2021b), based on machine learning algorithms applied to buoy data from the U.S. coast, showed that the spectral bandwidth was much more informative about rogue wave probability than the BFI.
On the other hand, some authors have found that the rogue wave occurrence frequency was not or only weakly dependent on the spectral bandwidth (e.g., Stansell 2004). Goda (1970) wrote, supported by numerical experiments, that wave heights defined by the zero-upcrossing method practically followed the Rayleigh distribution, independently of the spectral bandwidth. When examining the broadness parameter ε, Christou and Ewans (2014) found little difference between rogue wave samples and the highest nonrogue samples.
Also concerning the directional spreading during rogue wave occurrence, investigations have come to different conclusions in different parts of the ocean. While Waseda et al. (2011) confirmed, based on radar measurement data from a platform in the North Sea, that on days with a high occurrence of rogue waves, the directional spreading of the wave spectrum was narrower than on other days, Christou and Ewans (2014) found no significant differences between nonrogue samples and rogue wave samples and concluded that the environmental conditions generating normal waves, were also able to form rogue waves.
As a result, a number of authors have come to the conclusion that the modulational instability was not the main reason for the formation of rogue waves in their measurement data (e.g., Cattrell et al. 2018). Fedele et al. (2016) found, based on a large collection of field data from various locations in Europe, that the main generation mechanism of rogue waves is the constructive interference of elementary waves (dispersive and directional focusing), enhanced by second-order bound nonlinearities and not the modulational instability. They concluded that rogue waves are likely to be rare occurrences of weakly nonlinear random seas.
In a previous study, we investigated a measurement dataset from the southern North Sea (Teutsch et al. 2020). We found an enhanced occurrence frequency of rogue waves compared to the Forristall distribution in the overall dataset, but especially when recorded by radar devices. The Forristall distribution was developed based on waves in the Gulf of Mexico during hurricane conditions (Forristall 1978). Although it was found to be a good fit to measurement data in other ocean regions (e.g., Casas-Prat et al. 2009; Waseda et al. 2011), it should be kept in mind that this empirical distribution accounts for nonlinearities present in the original buoy dataset. Nonlinearities might be different in the southern North Sea, which is especially relevant when investigating rogue waves. In the present study, we investigate the possibility that those rogue waves exceeding the expectations of random superposition, were generated by modulational instabilities. Even at stations showing a good agreement with the Forristall distribution, modulational instabilities might be present. These are not necessarily reflected by the wave statistics (Fernandez et al. 2016). Based on a field measurement dataset of approximately 123 000 samples from radar stations and 63 000 samples from wave buoys in the southern North Sea, our aim is to investigate whether the requirements for a modulational instability, that is, narrow bandwidth in both frequency and angular direction, is given during rogue wave occurrence. This is done in terms of the bandwidth parameters [Eqs. (2)–(5)], the directional spreading [Eq. (9)] and the combined parameter R for directional spreading and spectral bandwidth [Eq. (8)]. Finally, the BFI [Eq. (1)] in our dataset is evaluated and compared for time series with and without rogue wave occurrence.
The measurement area and the dataset, as well as the estimation method of the directional spectrum, are described in section 2. The results of the evaluation of the bandwidth parameters and the BFI are described in section 3. In section 4, our results are related to previous experimental studies, and differences in the findings from the two different measurement devices are discussed. From the results and discussions, we draw our conclusions in section 5.
2. Data and methods
a. Data
Surface elevation measurement data from 2011 to 2016 at eight measurement stations in the southern North Sea were investigated (Fig. 1). The data are part of a larger set that was introduced in an earlier study (Teutsch et al. 2020). At four of the stations, the data were provided by fixed radar devices, measuring the air gap to the sea surface at a sampling frequency of either fs = 2 Hz or fs = 4 Hz. The remaining four stations are equipped with surface-following buoys of type MkIII, measuring at a sampling rate of fs = 1.28 Hz. The buoy data were delivered in samples of 30-min duration. The radar data, which were available as continuous time series, were split into 30-min samples accordingly. All samples were quality controlled according to the procedure described in Teutsch et al. (2020).

Measurement sites considered in the study. Red circles: radar stations; blue squares: wave buoys.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Measurement sites considered in the study. Red circles: radar stations; blue squares: wave buoys.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Measurement sites considered in the study. Red circles: radar stations; blue squares: wave buoys.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
To exclude low-energy sea states, only samples with a significant wave height Hs above the long-term 70th percentile of the significant wave height at each station Hs,70 were included in the analysis. The significant wave height Hs is here defined as the mean of the highest 30 % of the wave heights in a 30-min sample. The Hs,70 was calculated from the significant wave heights Hs of all 30-min samples during the six years of available measurement data. On the one hand, this excludes possible measurement uncertainties caused by small waves that are only described by a few points, and on the other hand, it includes only rogue waves of heights relevant for offshore activities. The Hs,70 is presented for each station in Table 1.
Water depth h and long-term 70th percentile of the significant wave height Hs,70 at all stations included in the analysis.


The 99th percentile of the peak period Tp,99, in all samples above Hs,70, compared to the peak period threshold Tp, below which samples are analyzed with respect to preconditions for a modulational instability.


Number of available 30-min measurement samples with Hs above the long-term 70th percentile of the significant wave height Hs,70. The data are classified according to the station and the sampling frequency fs of the recording.


b. Methods
For all 30-min samples in Table 3, the spectral bandwidth was estimated from the FFT spectrum in terms of the broadness parameter ε [Eq. (2)], the narrowness parameter ν [Eq. (4)] and Goda’s peakedness parameter Qp [Eq. (5)]. In addition, the BFI was calculated for each sample, according to Eq. (7). As described in Teutsch et al. (2020), each sample was assigned to one of the following categories:
-
“Nonrogue samples”—measurement samples that did not contain a rogue wave.
-
“Height rogue samples”—measurement samples that contained a rogue wave according to the height criterion (Haver and Andersen 2000)
with H denoting the height of the rogue wave from crest to trough.
-
“Extreme rogue samples”—measurement samples that contained a rogue wave according to a more strict height criterion of
-
“Crest rogue samples”—measurement samples that contained a rogue wave according to the crest criterion (Haver and Andersen 2000)
with C denoting the crest height of the rogue wave above still water level.
-
“Double rogue samples”—measurement samples that contained a rogue wave according to both the criteria defined in Eqs. (11) and (13). Double rogue samples were excluded in the groups of height and crest rogue samples.
Table 4 shows the number and percentage of available samples in each category.
Number of samples in each category for the groups of radar and buoy stations.


The Datawell Waverider buoys provide information on their position on the water surface in heave, north and west directions. Based on the three-dimensional information, it is possible to estimate a directional wave spectrum for each 30-min sample. We calculated the wave spectrum from buoy data according to Huntley et al. (1977), who applies the iterated maximum likelihood method (IMLM) developed by Pawka (1983). According to Young (1994), a maximum likelihood method yields a broader spectrum with a lower peak than the real spectrum. However, as the author points out, it is “the analysis technique with the highest directional resolving capabilities” for wave buoy data (Young 1994). Based on the directional spectrum, the directional spreading of each sample, and subsequently the parameter R, which relates the directional spreading to the spectral bandwidth, were calculated [Eqs. (8) and (9)].
3. Results
a. Spectral bandwidth
The measurement data analyzed in this study, were recorded by different measurement instruments and at different sampling rates (Table 3). When comparing spectral bandwidth parameters and BFI at all stations, we observed that the parameter range differed with the measurement frequency fs (Figs. 2 and 3). This became most obvious for the broadness parameter ε, whose median value showed to be approximately 30% higher in radar samples measured at fs = 4 Hz than in buoy samples with fs = 1.28 Hz (Fig. 2). To test the hypothesis that the investigated spectral parameters are sensitive to the measurement frequency, we subsampled the 4-Hz data at station L9 at 2 and at 1 Hz, respectively. The result demonstrates that the same time series yield different spectral parameters when sampled at different frequencies (Fig. 4). The changes are most pronounced for the broadness parameter ε. The reason for this dependency of the parameters on the sampling frequency is a change in spectral shape with the measurement frequency (Fig. 5). This results in a change in spectral moments (Table 5). A change in sampling frequency strongly affects those bandwidth parameters that are dependent on the higher moments of the spectrum [Eqs. (2) and (4)]. This issue has already been raised by Goda (1988a,b), who introduced a peakedness parameter that is independent of the higher moments of the frequency spectrum [Eq. (5)]. Indeed, the peakedness parameter Qp is least of all bandwidth parameters affected by the measurement frequency (Fig. 2).

Bandwidth parameters ε, ν, and Qp, calculated from the available data at all measurement stations.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Bandwidth parameters ε, ν, and Qp, calculated from the available data at all measurement stations.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Bandwidth parameters ε, ν, and Qp, calculated from the available data at all measurement stations.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Benjamin–Feir index, calculated from the available data at all measurement stations.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Benjamin–Feir index, calculated from the available data at all measurement stations.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Benjamin–Feir index, calculated from the available data at all measurement stations.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Data from station L9, originally sampled at 4 Hz, and the same data subsampled at 2 and 1 Hz.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Data from station L9, originally sampled at 4 Hz, and the same data subsampled at 2 and 1 Hz.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Data from station L9, originally sampled at 4 Hz, and the same data subsampled at 2 and 1 Hz.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Change in spectral shape with a change in sampling frequency.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Change in spectral shape with a change in sampling frequency.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Change in spectral shape with a change in sampling frequency.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Spectral moments mn, bandwidth parameters, and BFI, calculated from the same time series, measured starting at 0000 LT 1 Jan 2016. It was originally sampled at fs = 4 Hz and subsamples at 2 and 1 Hz.


Since our dataset consists of buoy and radar measurement data, which have been sampled at different frequencies, buoy and radar data will be treated separately throughout this study. Within the set of radar measurements, we subsampled all time series recorded at 4 Hz with fs = 2 Hz. This way, we still retain measurement data at a high frequency, while enlarging the set of 2-Hz samples and obtaining bandwidth parameters in a range that is more comparable to the buoy data than the 4-Hz data. Figure 6 is a repetition of Fig. 2, but with all 4-Hz radar data subsampled at 2 Hz. The BFI for 2-Hz radar data and 1.28-Hz buoy data is shown in Fig. 7.

Bandwidth parameters ε, ν, and Qp, calculated from the available data at all measurement stations. The 4-Hz time series measured at radar stations have been subsampled at fs = 2 Hz for comparability.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Bandwidth parameters ε, ν, and Qp, calculated from the available data at all measurement stations. The 4-Hz time series measured at radar stations have been subsampled at fs = 2 Hz for comparability.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Bandwidth parameters ε, ν, and Qp, calculated from the available data at all measurement stations. The 4-Hz time series measured at radar stations have been subsampled at fs = 2 Hz for comparability.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Benjamin–Feir index, calculated from the available data at all measurement stations. The 4-Hz time series measured at radar stations have been subsampled at fs = 2 Hz for comparability.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Benjamin–Feir index, calculated from the available data at all measurement stations. The 4-Hz time series measured at radar stations have been subsampled at fs = 2 Hz for comparability.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Benjamin–Feir index, calculated from the available data at all measurement stations. The 4-Hz time series measured at radar stations have been subsampled at fs = 2 Hz for comparability.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
The differences between the respective stations are insignificant. Within the subset of radar measurement samples, one could infer an increased broadness of the frequency spectrum in terms of both ε and ν with a decrease in water depth, with L9 and Clipper situated in water depths much shallower than Leman (see Table 1). This trend, however, is not seen in the peakedness parameter Qp, and it is not confirmed by the buoy stations, with the deepest stations being FN1 and LTH with a water depth of h = 30 m (Fig. 6). The BFI does not display any relation to the water depth in which the samples were recorded (Fig. 7).
Serio et al. (2005) suggested to include information on the water depth in the formulation of the BFI. We have recalculated the BFI according to their proposed equations and found no significant differences compared to the BFI calculated from Eq. (1). However, an extension of the BFI to directional effects, as proposed by Mori et al. (2011) [

The 2D Benjamin–Feir index according to Mori et al. (2011), compared with the BFI according to Eq. (1) at the buoy stations.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

The 2D Benjamin–Feir index according to Mori et al. (2011), compared with the BFI according to Eq. (1) at the buoy stations.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
The 2D Benjamin–Feir index according to Mori et al. (2011), compared with the BFI according to Eq. (1) at the buoy stations.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Figure 9 shows the broadness parameter ε, combined for each instrument category, but separated into the different sample categories defined in section 2. Figure 9 and Table 6 indicate that, compared to nonrogue samples, spectral bandwidth is slightly enhanced in crest rogue samples and slightly reduced in height and extreme rogue samples. The differences are somewhat more apparent when displayed in a histogram, comparing nonrogue to height and extreme rogue samples, and nonrogue to crest rogue samples (Figs. 10 and 11). Double rogue samples display a lower bandwidth than usual (Table 6), but they will not be treated further, as they belong to both height and crest rogue categories.

Broadness parameter ε, calculated from the different time series categories. Distributions are shown as box-and-whisker plots (box: interquartile range; whiskers: 1.5 times the interquartile range; horizontal line inside the box: median; dots: data outside the whiskers).
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Broadness parameter ε, calculated from the different time series categories. Distributions are shown as box-and-whisker plots (box: interquartile range; whiskers: 1.5 times the interquartile range; horizontal line inside the box: median; dots: data outside the whiskers).
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Broadness parameter ε, calculated from the different time series categories. Distributions are shown as box-and-whisker plots (box: interquartile range; whiskers: 1.5 times the interquartile range; horizontal line inside the box: median; dots: data outside the whiskers).
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Broadness parameter ε, calculated at the radar stations, comparing nonrogue to (a) height and extreme rogue samples and (b) crest rogue samples.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Broadness parameter ε, calculated at the radar stations, comparing nonrogue to (a) height and extreme rogue samples and (b) crest rogue samples.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Broadness parameter ε, calculated at the radar stations, comparing nonrogue to (a) height and extreme rogue samples and (b) crest rogue samples.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

As in Fig. 10, but calculated at the buoy stations.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

As in Fig. 10, but calculated at the buoy stations.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
As in Fig. 10, but calculated at the buoy stations.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
To test whether the difference between the ε distributions in the different sample categories is statistically significant, a Monte Carlo technique was applied to ε, and subsequently for all bandwidth parameters and the BFI, comparing the parameters calculated from nonrogue samples to height and extreme rogue samples on the one hand and to crest rogue samples on the other hand. An overview of the Monte Carlo calculations that were performed, is given in Table 7. The applied Monte Carlo technique is described in the following. In each execution, 10 000 samples of the length of the comparison population (height and extreme or crest rogue time series) were drawn with replacement from the population of nonrogue time series. For these samples, the respective parameters in Table 7 were calculated. The mean and standard deviation of the 10 000 parameter distributions is displayed and compared to the comparison population in Figs. 12–15. A significant difference between nonrogue time series and the comparison population is given if the distribution of the comparison population is located outside of two standard deviations of the distribution calculated from the nonrogue time series.

Result of the Monte Carlo technique for the broadness parameter ε. Results from (left) height and extreme rogue samples and (right) crest rogue samples are compared to results from 10 000 samples, drawn from the population of nonrogue time series. For these samples, the mean and the standard deviations of the distributions are plotted (Campbell 2022).
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Result of the Monte Carlo technique for the broadness parameter ε. Results from (left) height and extreme rogue samples and (right) crest rogue samples are compared to results from 10 000 samples, drawn from the population of nonrogue time series. For these samples, the mean and the standard deviations of the distributions are plotted (Campbell 2022).
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Result of the Monte Carlo technique for the broadness parameter ε. Results from (left) height and extreme rogue samples and (right) crest rogue samples are compared to results from 10 000 samples, drawn from the population of nonrogue time series. For these samples, the mean and the standard deviations of the distributions are plotted (Campbell 2022).
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

As in Fig. 12, but for the narrowness parameter ν.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

As in Fig. 12, but for the narrowness parameter ν.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
As in Fig. 12, but for the narrowness parameter ν.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

As in Fig. 12, but for the peakedness parameter Qp.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

As in Fig. 12, but for the peakedness parameter Qp.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
As in Fig. 12, but for the peakedness parameter Qp.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

As in Fig. 12, but for the BFI.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

As in Fig. 12, but for the BFI.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
As in Fig. 12, but for the BFI.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Parameters and stations for which a Monte Carlo technique was performed: nonrogue samples were compared with height/extreme rogue samples and with crest rogue samples, respectively.


From the result of the Monte Carlo technique for the broadness parameter ε, it can be inferred that the frequency spectra calculated from height and extreme rogue samples at the investigated radar and buoy stations were significantly more narrow than the frequency spectra calculated from time series without rogue waves (Fig. 12). It can be concluded that in our case, height and extreme rogue waves typically occurred in slightly more narrow sea states than observed under nonrogue conditions. For crest rogue samples, Fig. 12 shows a different result: these did not differ significantly from nonrogue time series.
The standard deviation in the 10 000 samples drawn from radar data, is very small, due to a narrow ε distribution in radar samples, as seen in Fig. 10.
The result of the Monte Carlo technique for the narrowness parameter ν is slightly different. A significantly more narrow frequency spectrum for height and extreme rogue samples can only be identified at the buoy stations (Fig. 13). The distribution of ν for height and extreme rogue samples at radar stations is situated within two standard deviations of the ν distribution calculated from 10 000 nonrogue sample realizations. For crest rogue samples at both radar and buoy stations, the result remains that the distribution of ν is not significantly different from the ν values in nonrogue samples.
In terms of the peakedness parameter Qp, for which a mean value slightly above Qp = 2 is typical at all stations, no significant difference between nonrogue and rogue wave samples was identified (Fig. 14). Although the peakedness/groupiness of the wave spectrum appears slightly lower for crest rogue samples than usual and slightly higher for extreme rogue samples than usual, the curves of the comparison populations are within two standard deviations of the distribution from nonrogue samples. The same result applies to the BFI. Under the hypothesis that the modulational instability caused the rogue waves measured at our stations, a higher BFI would be expected in rogue wave samples. However, this does not seem to be the case (Fig. 15). Although within two standard deviations, low BFI values (BFI ≈ 0.2) seem to be more unusual in height and extreme rogue samples than in nonrogue samples.
b. Directional spreading
The directional spreading of the wave energy in each sample is only available at the buoy stations, which provide three-dimensional information on their location, but not at the radar stations with only one-dimensional information on the air gap. Since a narrow directional spreading represents a requirement for the process of the modulational instability to occur, we investigated the directional spreading at all buoy stations and compared it for time series with and without rogue waves. Following the reasoning in the previous section, this was done by means of a Monte Carlo technique, in which a sample of 10 000 time series was randomly drawn from the population of nonrogue time series. Figure 16 shows that the directional spreading in the measured rogue wave samples was not significantly different from the directional spreading in samples without rogue waves. Thus, the condition of a specifically low directional spreading in rogue wave samples as a prerequisite for the modulational instability to operate, was not present in our buoy data.

As in Fig. 12, but for the directional spreading.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

As in Fig. 12, but for the directional spreading.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
As in Fig. 12, but for the directional spreading.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Finally, the combination of spectral bandwidth and directional spreading was examined. This was done because the favoring condition for the formation of a modulational instability has been formulated as a narrow spectrum in both direction and frequency (e.g., Alber 1978; Stansberg 1995). Figure 17 shows bandwidth-directional spreading pairs for all buoy samples. It shows that nonrogue samples, as well as the plotted rogue wave samples, cluster around medium bandwidth and directional spreading values. No specific accumulation of rogue wave samples at low bandwidth and directional spreading is seen, which would point toward favorable conditions for the modulational instability.

Comparison of the bandwidth-directional spreading pairs calculated from nonrogue and from rogue wave samples. The bandwidth is quantified in terms of the parameters ε and ν. Color coding: density of nonrogue samples. Black dots: rogue wave samples, as defined in the legend to each panel.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Comparison of the bandwidth-directional spreading pairs calculated from nonrogue and from rogue wave samples. The bandwidth is quantified in terms of the parameters ε and ν. Color coding: density of nonrogue samples. Black dots: rogue wave samples, as defined in the legend to each panel.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Comparison of the bandwidth-directional spreading pairs calculated from nonrogue and from rogue wave samples. The bandwidth is quantified in terms of the parameters ε and ν. Color coding: density of nonrogue samples. Black dots: rogue wave samples, as defined in the legend to each panel.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Additionally, the Monte Carlo technique was repeated for parameter R [Eq. (8)] in buoy measurement samples, which quantifies the importance of the directional spreading with respect to the spectral bandwidth. Although it appears that the maximum in the rogue wave distributions is shifted toward a higher R value, compared to nonrogue samples, the comparison with the standard deviation reveals that there are no significant differences apparent between rogue wave samples and nonrogue samples (Fig. 18). We conclude that rogue waves in our buoy dataset did not occur in sea states that are narrow-banded in frequency and angular direction. This indicates that the modulational instability is an unlikely mechanism for the formation of these rogue waves.

As in Fig. 12, but for the combination parameter R.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

As in Fig. 12, but for the combination parameter R.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
As in Fig. 12, but for the combination parameter R.
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
4. Discussion
Wave measurement time series from the southern North Sea that were evaluated in the present study, generally showed a wide spectral bandwidth, which is an unlikely condition for nonlinear focusing to occur in deep water (Alber 1978; Fedele et al. 2016). A similar observation of unfavorable conditions for the modulational instability was made by Fedele et al. (2016), who presented sea state parameters for three historical rogue wave occurrences in the North Sea. Due to a high bandwidth parameter ν and a large directional spreading σΘ, they came to the conclusion that these rogue waves were probably not the result of nonlinear focusing. This conclusion is transferable to our dataset, which displays high spectral bandwidth and directional spreading values for a substantially larger number of rogue waves in the North Sea.
Theoretical considerations and experimental results clearly confirm an increased rogue wave occurrence in narrow bandwidth and directional spreading (Waseda et al. 2009). Observations of field measurement data, however, have come to contradictory conclusions. In our data, the spectral bandwidth was significantly lower than usual in height and extreme rogue wave samples, when quantified by the broadness parameter ε [Eq. (2)]. This result is basically consistent with a finding by Christou and Ewans (2014), who discovered, mainly based on radar data from the North Sea, recorded at fs = 2 or 4 Hz, that rogue wave samples are more narrow-banded than the highest samples of nonrogue waves, and concluded that the spectral bandwidth might be an indicator for distinguishing time series with rogue waves from time series with high nonrogue waves. Christou and Ewans (2014) did not distinguish between height and crest rogue waves, but given that in our data, the sample of height rogue samples is approximately 5 times larger than that of crest rogue samples (Table 4), our results are still comparable to theirs. One major difference in the result, however, is that Christou and Ewans (2014) based their observation on the narrowness parameter ν. In terms of the broadness parameter ε, where we detected the strongest deviations, Christou and Ewans (2014) found little difference between rogue wave samples and the highest nonrogue samples. The dataset of Christou and Ewans (2014) included data from the same radar stations as we used in our study, however, for a different time period. Another difference to our dataset is the use of additional data from other regions from a large range of water depths (7–1311 m). The differences in the investigated data could be the reason for our discrepancies from their study. At the buoy stations, we also recognized differences between nonrogue and rogue wave samples in the narrowness parameter ν. Also, Karmpadakis et al. (2020) described a clear influence of the spectral bandwidth ν on the wave height distribution in radar measurements in the North Sea. They found wave heights to decrease with increasing bandwidth. Most of their time series were recorded at fs = 2–4 Hz. Also, the dataset of Karmpadakis et al. (2020) includes the radar stations we used in the present study, with the time series starting somewhat earlier and continuing until February 2017. In addition to the difference in measurement period, it should be noted that Karmpadakis et al. (2020) did not exclude any Hs values, which makes our study not exactly comparable to theirs. However, Karmpadakis et al. (2020) distinguish between deep and shallow water in their results. Most recent findings by Häfner et al. (2021b) revealed, based on machine learning algorithms applied to Coastal Data Information Program (CDIP) buoy data from the U.S. coast with Hs > 1 m, that the spectral bandwidth (ν and Qp) provides some information on the probability of rogue waves according to the height criterion. Furthermore, their study revealed that bandwidth effects are not relevant for wave crests. Based on our measurement data, we can confirm these results: although we observed slightly higher bandwidths than usual in crest rogue wave time series, these differences showed to be statistically insignificant. The CDIP dataset includes data from water depths between approximately 10 and 4500 m. Häfner et al. (2021b) supposed that the influence of second-order nonlinearities on wave crests might be higher than estimated on the basis of buoy measurement data. Wave buoys tend to cancel second-order nonlinearities by their own Lagrangian movement and thus to overestimate the mean water level, which in turn leads to an underestimation of crest heights (Forristall 2000). In addition, wave buoys may be dragged through or slide away from short wave crests, which might result in missing the maximum amplitudes. Especially for strongly nonlinear waves, this leads to significant differences compared to radar devices (Longuet-Higgins 1986). These in turn may overestimate crests by misinterpreting spray, breaking waves, or even fog (Grønlie 2006). Forristall (2005) noted that there is no standard way to calibrate measurement instruments and that it is not possible to decide which instrument yields the “most correct” results. With our dataset consisting of both buoy and radar time series, we had the possibility to compare the results for wave crests based on both instruments and came to the conclusion that the influence of spectral bandwidth on crest rogue waves is similarly negligible in both types of time series. Also, Cattrell et al. (2018) stated an influence of the spectral bandwidth on rogue wave occurrence. However, they came to different conclusions. In their investigation of 80 of the 161 CDIP buoys that were available to Häfner et al. (2021b), the distributions of the spectral bandwidth parameters ε and ν indicated that the probability of observing a height rogue wave increased in seas with a higher bandwidth, while that of observing a crest rogue wave increased in seas with a narrow spectral bandwidth. As opposed to Häfner et al. (2021b), Cattrell et al. (2018) investigated the entire range of Hs values, which could lead to the difference in findings. Our data from the southern North Sea seem to resemble the results of Häfner et al. (2021b). Finally, Stansell (2004) claimed, based on a statistical analysis of storm waves at the North Alwyn platform in the North Sea, measured by laser altimeters (Hs > 3 m, h = 130 m, fs = 5 Hz), that the occurrence probability of rogue waves is only weakly dependent on the spectral bandwidth ν, which confirms our results regarding the narrowness parameter ν at the radar stations.
It is possible that the results on the influence of spectral bandwidth on rogue wave occurrence are dependent on the sampling frequency of the measurement device: while buoy measurements at fs = 1.28 Hz identified strong dependencies of at least height rogue wave occurrence on the spectral bandwidth (Häfner et al. 2021b; Cattrell et al. 2018), and these results were still valid (Karmpadakis et al. 2020) or less pronounced (Christou and Ewans 2014) in 2–4-Hz radar measurements, the influence of the spectral bandwidth was described as insignificant in a study based on 5-Hz laser measurements (Stansell 2004). We can confirm that in our dataset, the influence of spectral bandwidth on height rogue waves was more pronounced in the time series recorded at fs = 1.25 Hz than in the 2 Hz and subsampled at 2 Hz, radar data. This leads to the question to what extent the influence of the spectral bandwidth on rogue wave occurrence is actually a physical mechanism and how much of the effect should be attributed to the measurement instrument. Measurement instruments with low sampling frequencies are known to underestimate the occurrence frequency and the heights of the largest waves (Stansell et al. 2002). A false estimation of the spectral bandwidth could be a consequence. For further investigations on this issue, it would be valuable to have a radar device and a wave buoy installed not only in the same area, but at exactly the same place and compare measurement results.
Other reasons for differences in measurement results could be the restriction to specific Hs ranges applied by some authors (e.g., Stansell 2004; Häfner et al. 2021b), or a different behavior of waves at shallow water stations, which were not considered separately in all studies. Although rogue waves may occur in all sea states (Teutsch et al. 2020), we have focused our study on rogue waves in high sea states with a great relevance for offshore activities. This choice may have an influence on the results and conclusions of our study. An alternative option would be to group time series by the metocean conditions during the recordings or by sea state characteristics like spectral shape or steepness, while including sea states of all heights. However, including low sea states introduces measurement uncertainties. Subdividing the dataset of samples above Hs,70 into sea state groups reduces the size of the different groups in a way that results cannot be drawn on a reliable basis. Therefore, our results and conclusions are presented for all samples above Hs,70.
In our data, the BFI as the commonly accepted indicator for the modulational instability, was not characteristic during rogue wave occurrence. This again confirms a result of Häfner et al. (2021b), who showed for their data that the BFI as a measure of nonlinear effects does not play a role for the prediction of rogue waves. Furthermore, it should be kept in mind that the formulation of the BFI was derived under the assumption of a narrow-banded Gaussian spectrum (Janssen and Bidlot 2009). With this in mind, it should be questioned whether the BFI is a valid indicator for the application to real ocean data. Waseda (2006) showed in tank experiments that the BFI was only informative on the non-Gaussianity of a sea state in very narrow-banded waves.
Also in terms of the directional spreading, we could not identify any characteristic tendency during rogue wave occurrence, which agrees with findings by Christou and Ewans (2014) and Häfner et al. (2021b). We agree with Christou and Ewans (2014), who conclude from this observation that “the environmental conditions generating normal waves, are also able to form rogue waves.” In a recent study that evaluated numerical models in comparison with in situ data, it was confirmed that both linear and nonlinear effects are important for the formation of rogue waves and thus, rogue waves may also be expected in sea states with large directional spreading (Kirezci et al. 2021). On the other hand, Waseda et al. (2011) observed in a hindcast that on days with a high occurrence of height rogue waves, as identified by radar measurements at a platform in the North Sea, the directional spreading of the wave spectrum was narrower than on other days. In our combined analysis of directional spreading and spectral bandwidth, we found that our rogue waves did not occur in sea states that were narrow-banded in frequency and angular direction.
We ensured for all parameters in rogue wave samples, by reducing the height of each identified rogue wave by 50%, that the bandwidth parameters were not influenced by the presence of the rogue wave itself. In fact, this test showed that while the significant wave height of the samples decreased as a result of the reduction of the rogue wave height, the bandwidth parameters and the BFI did not change. Thus, the influence of one rogue wave on the spectral parameters of a 30-min time series was small, as opposed to the 20-min windows investigated by Stansell (2004). He found that the difference in spectral bandwidth disappeared when the rogue wave was removed. In our 30-min time series, it appeared that the information in the bandwidth on the possible presence of a rogue wave was contained within the wave spectrum, not in the rogue wave itself. Analyzing measurement data and comparing the results to previous measurement studies raises the issue of sampling variability: statistical results become less reliable due to limited data availability. In 30-min time series, the effect of sampling variability may dominate over nonlinear effects, which are then difficult to assess (Bitner-Gregersen et al. 2021). This effect is certainly important when studying rogue waves, which are located at the tail of the wave height distribution and therefore naturally represented by only few measurements. The uncertainty imposed by the limited number of data, may be counteracted by spatial measurements (Bitner-Gregersen et al. 2021). This possibility was not available to us, a problem of many measurement studies in the ocean. The analysis of longer records could be another solution, but samples longer than 30 min cannot be assumed to be stationary, which again poses problems to the statistical evaluation (Holthuijsen 2007). The use of 30-min samples is therefore a compromise between long enough records (to produce reliable statistical results) and short enough records (to reasonably assume stationary wave conditions). The measurement data from the eight stations in the southern North Sea have been part of an earlier study (Teutsch et al. 2020). In that study, the rogue wave frequency was documented with respect to the total number of measured waves at each station (although without the Hs,70 restriction). When comparing the bandwidth parameters ε and ν, as presented in Fig. 6, with the rogue wave frequencies from the previous study, it is seen that a high rogue wave occurrence corresponds to a narrow frequency spectrum (Fig. 19).

Comparison of the bandwidth parameters ε and ν with rogue wave occurrence frequencies from Teutsch et al. (2020).
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1

Comparison of the bandwidth parameters ε and ν with rogue wave occurrence frequencies from Teutsch et al. (2020).
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
Comparison of the bandwidth parameters ε and ν with rogue wave occurrence frequencies from Teutsch et al. (2020).
Citation: Journal of Physical Oceanography 53, 1; 10.1175/JPO-D-22-0059.1
The parameters Qp and BFI did not show any correlation with the rogue wave frequency. The directional spreading could not be compared at all stations, since the radar measurements are one-dimensional. From Teutsch et al. (2020), it is noted that all buoy stations showed rogue wave frequencies below the expectation from second-order theory. Here, Janssen and Bidlot (2009) should be cited, who state that the sea state defocusses above a certain threshold of σΘ, which yields even fewer rogue waves than expected according to second-order theory. Having the thresholds of 14°–15° in wave tank experiments (Waseda 2006; Onorato et al. 2009) and 30° in the ocean (Waseda et al. 2009) in mind, this could indeed be the case in our measurement data.
5. Conclusions
In a dataset consisting of radar and buoy measurements from high sea states in the southern North Sea, we identified lower values of spectral bandwidth than usual during height rogue wave occurrence. Samples with rogue waves according to the crest criterion in turn could not be attributed to specific bandwidth conditions. The directional spreading did not give any indication on the occurrence of rogue waves of any kind, neither did the BFI as a commonly applied indicator for the modulational instability. We conclude that the majority of rogue waves in our dataset were probably not generated by a modulational instability, since high spectral bandwidth and directional spreading are unfavorable conditions for nonlinear focusing.
Acknowledgments.
Ina Teutsch received funding for this work from the Federal Maritime and Hydrographic Agency (BSH). The buoy data were kindly provided by BSH. The authors are grateful to Graham Feld and Shell for providing the radar data. The authors are grateful to Jens Möller and Christian Senet for their valuable input. The article processing charges for this publication were covered by a Research Centre of the Helmholtz Association.
Data availability statement.
The underlying wave buoy and radar data are the property of and were made available by the Federal Maritime and Hydrographic Agency, Germany, and Shell, United Kingdom, respectively. They can be obtained upon request from these organizations.
REFERENCES
Alber, I. E., 1978: The effects of randomness on the stability of two-dimensional surface wave trains. Proc. Roy. Soc., 363A, 525–546, https://doi.org/10.1098/rspa.1978.0181.
Alber, I. E., and P. G. Saffman, 1978: Stability of random nonlinear deepwater waves with finite bandwidth spectra. TRW Rep. 31326-6035-RU-OO, 50 pp.
Benjamin, T. B., 1967: Instability of periodic wavetrains in nonlinear dispersive systems. Proc. Roy. Soc., 299A, 59–76, https://doi.org/10.1098/rspa.1967.0123.
Benjamin, T. B., and J. E. Feir, 1967: The disintegration of wave trains on deep water. J. Fluid Mech., 27, 417–430, https://doi.org/10.1017/S002211206700045X.
Bitner-Gregersen, E. M., O. Gramstad, A. K. Magnusson, and M. P. Malila, 2021: Extreme wave events and sampling variability. Ocean Dyn., 71, 81–95, https://doi.org/10.1007/s10236-020-01422-z.
Campbell, R., 2022: raacampbell/shadederrorbar. GitHub, https://github.com/raacampbell/shadedErrorBar.
Cartwright, D. E., and M. S. Longuet-Higgins, 1956: The statistical distribution of the maxima of a random function. Proc. Roy. Soc., 237A, 212–232, https://doi.org/10.1098/rspa.1956.0173.
Casas-Prat, M., L. Holthuijsen, and P. Gelder, 2009: Short-term statistics of 10,000,000 waves observed by buoys. Coastal Engineering 2008, Vol. 1, J. McKee Smith, Ed., World Scientific, 560–572, https://doi.org/10.1142/9789814277426_0047.
Cattrell, A. D., M. Srokosz, B. I. Moat, and R. Marsh, 2018: Can rogue waves be predicted using characteristic wave parameters? J. Geophys. Res. Oceans, 123, 5624–5636, https://doi.org/10.1029/2018JC013958.
Christou, M., and K. Ewans, 2014: Field measurements of rogue water waves. J. Phys. Oceanogr., 44, 2317–2335, https://doi.org/10.1175/JPO-D-13-0199.1.
Dysthe, K., H. E. Krogstad, and P. Müller, 2008: Oceanic rogue waves. Annu. Rev. Fluid Mech., 40, 287–310, https://doi.org/10.1146/annurev.fluid.40.111406.102203.
Fedele, F., 2015: On the kurtosis of deep-water gravity waves. J. Fluid Mech., 782, 25–36, https://doi.org/10.1017/jfm.2015.538.
Fedele, F., J. Brennan, S. P. de León, J. Dudley, and F. Dias, 2016: Real world ocean rogue waves explained without the modulational instability. Sci. Rep., 6, 27715, https://doi.org/10.1038/srep27715.
Fernandez, L., M. Onorato, J. Monbaliu, and A. Toffoli, 2014: Modulational instability and wave amplification in finite water depth. Nat. Hazards Earth Syst. Sci., 14, 705–711, https://doi.org/10.5194/nhess-14-705-2014.
Fernandez, L., M. Onorato, J. Monbaliu, and A. Toffoli, 2016: Occurrence of extreme waves in finite water depth. Extreme Ocean Waves, E. Pelinovsky and C. Kharif, Eds., Springer, 45–62, https://doi.org/10.1007/978-3-319-21575-4_3.
Forristall, G. Z., 1978: On the statistical distribution of wave heights in a storm. J. Geophys. Res., 83, 2353–2358, https://doi.org/10.1029/JC083iC05p02353.
Forristall, G. Z., 2000: Wave crest distributions: Observations and second-order theory. J. Phys. Oceanogr., 30, 1931–1943, https://doi.org/10.1175/1520-0485(2000)030<1931:WCDOAS>2.0.CO;2.
Forristall, G. Z., 2005: Understanding rogue waves: Are new physics really necessary? Rogue Waves: Proc. ‘Aha Huliko‘a Hawaiian Winter Workshop, Honolulu, HI, University of Hawai‘i at Mānoa, 29–35, http://www.soest.hawaii.edu/PubServices/2005pdfs/Forristall.pdf.
Gibson, R., and C. Swan, 2006: The evolution of large ocean waves: The role of local and rapid spectral changes. Proc. Roy. Soc., 463A, 21–48, https://doi.org/10.1098/rspa.2006.1729.
Goda, Y., 1970: Numerical experiments on wave statistics with spectral simulation. Port and Harbour Research Institute Rep. 009-03-01, 3–57, https://www.pari.go.jp/search-pdf/vol009-no03-01.pdf.
Goda, Y., 1988a: Numerical investigation on plotting formulas and confidence intervals of return values in extreme statistics (in Japanese). Port and Harbour Research Institute Rep. 027-01-02, 31–92, https://www.pari.go.jp/search-pdf/vol027-no01-02.pdf.
Goda, Y., 1988b: On the methodology of selecting design wave height. Coastal Eng. Proc., 1, 67, https://doi.org/10.9753/icce.v21.67.
Gramstad, O., and K. Trulsen, 2007: Influence of crest and group length on the occurrence of freak waves. J. Fluid Mech., 582, 463–472, https://doi.org/10.1017/S0022112007006507.
Grønlie, O., 2006: Wave radars: Techniques and technologies. Sea Technol., 47, 39–43.
Häfner, D., J. Gemmrich, and M. Jochum, 2021a: FOWD: A free ocean wave dataset for data mining and machine learning. J. Atmos. Oceanic Technol., 38, 1305–1322 ,https://doi.org/10.1175/JTECH-D-20-0185.1.
Häfner, D., J. Gemmrich, and M. Jochum, 2021b: Real-world rogue wave probabilities. Sci. Rep., 11, 10084, https://doi.org/10.1038/s41598-021-89359-1.
Haver, S., and O. J. Andersen, 2000: Freak waves: Rare realizations of a typical population or typical realizations of a rare population? 10th Int. Offshore and Polar Engineering Conf., Seattle, WA, International Society of Offshore and Polar Engineers, ISOPE-I-00-233, https://onepetro.org/ISOPEIOPEC/proceedings-abstract/ISOPE00/All-ISOPE00/ISOPE-I-00-233/7014.
Holthuijsen, L. H., 2007: Waves in Oceanic and Coastal Waters. Cambridge University Press, 404 pp., https://doi.org/10.1017/CBO9780511618536.
Huntley, D. A., R. T. Guza, and A. J. Bowen, 1977: A universal form for shoreline run-up spectra? J. Geophys. Res., 82, 2577–2581, https://doi.org/10.1029/JC082i018p02577.
Janssen, P. A. E. M., 2003: Nonlinear four-wave interactions and freak waves. J. Phys. Oceanogr., 33, 863–884, https://doi.org/10.1175/1520-0485(2003)33<863:NFIAFW>2.0.CO;2.
Janssen, P. A. E. M., and J.-R. Bidlot, 2009: On the extension of the freak wave warning system and its verification. ECMWF Tech Memo. 588, 42 pp., https://www.ecmwf.int/file/22682/download?token=1Rv5H5Z_.
Karmpadakis, I., C. Swan, and M. Christou, 2019: Laboratory investigation of crest height statistics in intermediate water depths. Proc. Roy. Soc., 475A, 20190183, https://doi.org/10.1098/rspa.2019.0183.
Karmpadakis, I., C. Swan, and M. Christou, 2020: Assessment of wave height distributions using an extensive field database. Coastal Eng., 157, 103630, https://doi.org/10.1016/j.coastaleng.2019.103630.
Kirezci, C., A. V. Babanin, and D. Chalikov, 2021: Probabilistic assessment of rogue wave occurrence in directional wave fields. Ocean Dyn., 71, 1141–1166, https://doi.org/10.1007/s10236-021-01487-4.
Lake, B. M., H. C. Yuen, H. Rungaldier, and W. E. Ferguson, 1977: Nonlinear deep-water waves: Theory and experiment. Part 2. Evolution of a continuous wave train. J. Fluid Mech., 83, 49–74, https://doi.org/10.1017/S0022112077001037.
Latheef, M., and C. Swan, 2013: A laboratory study of wave crest statistics and the role of directional spreading. Proc. Roy. Soc., 469A, 20120696, https://doi.org/10.1098/rspa.2012.0696.
Longuet-Higgins, M. S., 1975: On the joint distribution of the periods and amplitudes of sea waves. J. Geophys. Res., 80, 2688–2694, https://doi.org/10.1029/JC080i018p02688.
Longuet-Higgins, M. S., 1986: Eulerian and Lagrangian aspects of surface waves. J.Fluid Mech., 173, 683–707, https://doi.org/10.1017/S0022112086001325.
Mori, N., M. Onorato, and P. A. E. M. Janssen, 2011: On the estimation of the kurtosis in directional sea states for freak wave forecasting. J. Phys. Oceanogr., 41, 1484–1497, https://doi.org/10.1175/2011JPO4542.1.
Naess, A., 1985: On the distribution of crest to trough wave heights. Ocean Eng., 12, 221–234, https://doi.org/10.1016/0029-8018(85)90014-9.
Onorato, M., A. R. Osborne, M. Serio, and S. Bertone, 2001: Freak waves in random oceanic sea states. Phys. Rev. Lett., 86, 5831–5834, https://doi.org/10.1103/PhysRevLett.86.5831.
Onorato, M., A. R. Osborne, and M. Serio, 2002: Extreme wave events in directional, random oceanic sea states. Phys. Fluids, 14, 25–28, https://doi.org/10.1063/1.1453466.
Onorato, M., A. R. Osborne, M. Serio, L. Cavaleri, C. Brandini, and C. T. Stansberg, 2004: Observation of strongly non-Gaussian statistics for random sea surface gravity waves in wave flume experiments. Phys. Rev., 70E, 067302, https://doi.org/10.1103/PhysRevE.70.067302.
Onorato, M., and Coauthors, 2009: Statistical properties of mechanically generated surface gravity waves: A laboratory experiment in a three-dimensional wave basin. J. Fluid Mech., 627, 235–257, https://doi.org/10.1017/S002211200900603X.
Pawka, S. S., 1983: Island shadows in wave directional spectra. J. Geophys. Res., 88, 2579–2591, https://doi.org/10.1029/JC088iC04p02579.
Rye, H., 1977: The stability of some currently used wave parameters. Coastal Eng., 1, 17–30, https://doi.org/10.1016/0378-3839(77)90004-7.
Saulnier, J.-B., A. Clément, A. F. de O. Falcão, T. Pontes, M. Prevosto, and P. Ricci, 2011: Wave groupiness and spectral bandwidth as relevant parameters for the performance assessment of wave energy converters. Ocean Eng., 38, 130–147, https://doi.org/10.1016/j.oceaneng.2010.10.002.
Serio, M., M. Onorato, A. Osborne, and P. Janssen, 2005: On the computation of the Benjamin-Feir Index. Nuovo Cimento C, 28, 893–903, https://doi.org/10.1393/ncc/i2005-10134-1.
Stansberg, C. T., 1995: Effects from directionality and spectral bandwidth on non-linear spatial modulations of deep-water surface gravity wave trains. 24th Int. Conf. on Coastal Engineering 1994, Kobe, Japan, American Society of Civil Engineers, https://doi.org/10.1061/9780784400890.044.
Stansell, P., 2004: Distributions of freak wave heights measured in the North Sea. Appl. Ocean Res., 26, 35–48, https://doi.org/10.1016/j.apor.2004.01.004.
Stansell, P., J. Wolfram, and B. Linfoot, 2002: Effect of sampling rate on wave height statistics. Ocean Eng., 29, 1023–1047, https://doi.org/10.1016/S0029-8018(01)00066-X.
Tamura, H., T. Waseda, and Y. Miyazawa, 2009: Freakish sea state and swell-windsea coupling: Numerical study of the Suwa-Maru incident. Geophys. Res. Lett., 36, L01607, https://doi.org/10.1029/2008GL036280.
Tayfun, M. A., 1983: Effects of spectrum band width on the distribution of wave heights and periods. Ocean Eng., 10, 107–118, https://doi.org/10.1016/0029-8018(83)90017-3.
Teutsch, I., R. Weisse, J. Moeller, and O. Krueger, 2020: A statistical analysis of rogue waves in the southern North Sea. Nat. Hazards Earth Syst. Sci., 20, 2665–2680, https://doi.org/10.5194/nhess-20-2665-2020.
Toffoli, A., M. Benoit, M. Onorato, and E. M. Bitner-Gregersen, 2009: The effect of third-order nonlinearity on statistical properties of random directional waves in finite depth. Nonlinear Processes Geophys., 16, 131–139, https://doi.org/10.5194/npg-16-131-2009.
Toffoli, A., and Coauthors, 2013: Experimental evidence of the modulation of a plane wave to oblique perturbations and generation of rogue waves in finite water depth. Phys. Fluids, 25, 091701, https://doi.org/10.1063/1.4821810.
Vandever, J. P., E. M. Siegel, J. M. Brubaker, and C. T. Friedrichs, 2008: Influence of spectral width on wave height parameter estimates in coastal environments. J. Waterw. Port Coastal Ocean Eng., 134, 187–194, https://doi.org/10.1061/(ASCE)0733-950X(2008)134:3(187).
Waseda, T., 2006: Impact of directionality on the extreme wave occurrence in a discrete random wave system. Ninth Int. Workshop on Wave Hindcasting and Forecasting, Victoria, BC, Canada, American Society of Civil Engineering, http://www.waveworkshop.org/9thWaves/Papers/Waseda.pdf.
Waseda, T., T. Kinoshita, and H. Tamura, 2009: Evolution of a random directional wave and freak wave occurrence. J. Phys. Oceanogr., 39, 621–639, https://doi.org/10.1175/2008JPO4031.1.
Waseda, T., M. Hallerstig, K. Ozaki, and H. Tomita, 2011: Enhanced freak wave occurrence with narrow directional spectrum in the North Sea. Geophys. Res. Lett., 38, L13605, https://doi.org/10.1029/2011GL047779.
Young, I. R., 1994: On the measurement of directional wave spectra. Appl. Ocean Res., 16, 283–294, https://doi.org/10.1016/0141-1187(94)90017-5.
Yuen, H. C., and W. E. Ferguson, 1978: Relationship between Benjamin–Feir instability and recurrence in the nonlinear Schroedinger equation. Phys. Fluids, 21, 1275, https://doi.org/10.1063/1.862394.