# Return Period of a Sea Storm with at Least Two Waves Higher than a Fixed Threshold

Authors: Felice Arena; Giuseppe Barbaro; Alessandra Romolo
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

Practical applications in ocean engineering require the long-term analysis for prediction of extreme waves, that identify design conditions. If extreme individual waves are investigated, we need to combine long-term statistical analysis of ocean waves with short-term statistics. The former considers the distribution of standard deviation of free surface displacement in the considered location in a long-time span, of order of 10 years or more. The latter analyzes the distribution of individual wave heights in a sea state, which is a Gaussian process in time domain. Recent advanced approaches enable the combination of the two analyses. In the paper the analytical solution is obtained for the return period of a sea storm with at least two individual waves higher than a fixed level. This solution is based on the application of the Equivalent Triangular Storm model for the representation of actual storms. One of the corollaries of the solution gives the exact expression for the probability that at least two waves higher than fixed level are produced during the lifetime of a structure. The previous solution of return period and the relative probability of exceedance may be effectively applied for the risk analysis of ocean structures.


---

## Body

## 1. Introduction

The concept of Equivalent Triangular Storm (ETS) was introduced by Boccotti in the eighties [1, 2]. A complete analysis was then given in Boccotti [3], where the ETS model was applied for solving in a closed form of some long-term problems, which are very useful for the design of coastal and offshore structures.In other  words, the ETS model represents each actual storm as a triangle, where the triangle height is equal to the maximum significant wave height during the actual storm and the triangle base is such that the maximum expected wave height in the actual storm is equal to the maximum expected wave height in the triangular storm. Then, the triangle height is achieved from the significant wave height time series of the storm, while the triangle base is calculated by means of an iterative procedure.If the equivalent sea is considered, which is given by substituting a sequence of triangles to the sequence of actual storms, some analytical solutions were found for the return period of severe storms. Boccotti [3] found the analytical solution for the return period R(H) of a sea storm in which the highest wave has height exceeding a fixed threshold H.The solution for the return periodR(H) was then given, with a different formal derivation, by Arena and Pavone [4]. They considering the ETS approach combined the long-term and the second-order crest height distribution [5, 6] to achieve the return period R(C) of a sea storm in which the highest nonlinear crest height exceeds a fixed threshold C.Recently, Arena and Pavone [7] proposed a generalized approach for the long-term modelling of extreme waves, by giving the general expression of the return period RN(H) of a sea storm in which exactly N waves higher than H occur. Following this approach the solution is also given for the return period R≥N(H) of a sea storm in which at least N waves larger than a fixed threshold H occur, with N=1, 2, 3, and so forth.Different approaches to combine short- and long-term statistics were proposed by [8], who gave a generalization of Borgman’s results [9] and by Tromans and Vanderschuren [10] and Forristall [11], who introduced the analysis of extreme waves during storms and considered the loads and response on structures for engineering applications.In this paper, by following a different logic, with respect to that one adopted by Arena and Pavone [7], a new solution is found for the return period R≥2′(H) of a sea storm in which at least two waves higher than a fixed threshold H occur. The new solution for R≥2′(H) is achieved from direct combination of short-term and long-term wave statistics, while the Arena and Pavone solution gave the R≥2(H) solution as a function of both R≥1(H) and R1(H) by means of the following expression:
(1)1R≥2(H)=1R≥1(H)-1R1(H),
where R≥1(H) and R1(H) are the return periods of a sea storm in which, respectively, (i) at least a wave higher than H occurs; (ii) a single wave higher than H occurs.The final comparison between the two solutions, which is proposed starting from NOAA-NODC (USA) buoys data and from Italian buoys network (RON), from ISPRA, Institute for Environmental Protection and Research (Italy), shows a full agreement. This result confirms that the Equivalent Triangular Storm (as well as the Equivalent Power Storm by Fedele and Arena [12]) model is very powerful to determine the long-term statistics of extreme waves during storms (see also Arena et al. [13]).The results are of interest for engineering applications, because they enable to perform a complete analysis of extreme waves that we may expect in the lifetime of a structure, which may be, for example, either an offshore structure [14–19] or on an upright breakwater (see, e.g., Boccotti et al. [14, 15] and Romolo and Arena [20, 21]). This may be done in terms of return value of the extreme individual crest-to-trough wave height, as well as by achieving the second wave (the third and so on), in terms of height, occurring in a fixed time span on the considered structure.


## 2. Return PeriodR≥2′(H) of a Sea Storm in Which at Least Two Waves Higher than a Fixed Threshold H Occur

The return periodR≥2′(H) of a sea storm in which at least two waves higher than a fixed threshold H occur may be defined, in alternative, as the return period of a sea storm in which the second wave height (if waves are ordered by decreasing crest-to-trough wave heights) is larger than H.If a large time spanτ is considered, we may define
(2)R≥2′(H)=τNII(H,τ),
with NII(H,τ) the number of waves, during τ, which are (i) higher than H; (ii) the second waves (if waves occurring in the storm are ordered by decreasing height) in their own storm. In the following, the solution is given for NII(H,τ). Let us consider the Triangular Sea, which is given by the sequence of equivalent triangular storms.The number of waves with height betweenx and x+dx, occurring in sea states with significant wave height Hs in (h,h+dh) in the triangular storms with height between a and a+da and base between b and b+db, during τ, is
(3)[pA(a)daN(τ)pB(b∣a)db]Δt(h,dh,a,b)·p(x;Hs=h)dxT¯(h),
where N(τ)is the number of triangles during τ, PA(a) is the probability density function of the triangle heights, PB(b∣a) is the probability density function of the triangle base, with given height a, and Δt(h,dh,a,b) is the time in which the significant wave height Hs is in the range (h, h+dh) in those triangles with height a and base b. Finally, p(x;Hs=h) is the probability density function of the wave height in a sea state with Hs=h, and T-(h) is Rice’s mean period [22, 23].It follows that the numberNII(x,h;a,b) of waves with height between x and x+dx, occurring in sea states with significant wave height Hs in (h,h+dh) in the triangular storms with height between a and a+da and base between b and b+db, during τ, which are the second in order of height in their own storm, is given by
(4)NII(x,h;a,b)=[pA(a)daN(τ)pB(b∣a)db]Δt(h,dh,a,b)·p(x;Hs=h)dxT¯(h)P2(x,h;a,b),
where, in a triangular storm with height a and base b, given that a wave with height between x and x+dx occurs in a sea state with significant wave height between h and h+dh, P2(x,h;a,b) is the probability that:(i)
just a wave higher thanx will occur;(ii)
all the other waves of the sea storm will have height smaller thanx.More in general, if an actual storm is considered,P2 is the probability that a wave with height between x and x+dx, occurring in a sea state with significant wave height between h and h+dh, is the second wave in height in its own storm. It is then given by
(5)P2=11-P(x;Hs=h)×[-P(x;Hs=h)1-P(x;Hs=h)+∑i=1NP(x;Hs=hi)1-P(x;Hs=hi)Ni]×∏i=1N[1-P(x;Hs=hi)]Ni,
where N is the number of the sea states in the storm. Equation (5), in integral form, may be rewritten as [9, 24]
(6)P2=11-P(x;Hs=h)×[-P(x;Hs=h)1-P(x;Hs=h)+∫0DP[x;Hs=h(t)]1-P[x;Hs=h(t)]1T¯[h(t)]dt]·exp{∫0Dln[1-P(x;Hs=h(t))]1T¯[h(t)]dt},
where D is the storm duration.If the triangular sea is considered, for a sea storm with heighta and base b, the probability P2(x,h;a,b) is given by
(7)P2(x,h;a,b)=11-P(x;Hs=h)×[-P(x;Hs=h)1-P(x;Hs=h)+b-(a)ahjhh×∫0aP(x;Hs=h′′)1-P(x;Hs=h′′)1T-(h′′)dh′′].exp[b-(a)a∫0a1n[1-P(x;Hs=h′)]T-(h′)dh′].The number of waves, duringτ, which are higher than H and the second in their own storm, is
(8)NII(H,τ)=∫H∞∫0∞∫h∞∫0∞NII(x,h;a,b)dbdadhdx.In conclusion, combining (2), (4), (7), and (8), the expression is given for the return period R≥2′(H) of a sea storm in which at least two waves higher than a fixed threshold H occur:
(9)R≥2′(H)={ln[1-P(x;Hs=h′)]T-(h′)∫H∞∫0∞1T¯(h)p(x;Hs=h)11-P(x;Hs=h)·∫h∞-dp(Hs=a)da×exp[b-(a)a∫0aln[1-P(x;Hs=h′)]T-(h′)dh′]·[b-(a)a∫0a1T-(h′′)P(x;Hs=h′′)1-P(x;Hs=h′′)dh′′+-P(x;Hs=h)1-P(x;Hs=h)]dadhdxln[1-P(x;Hs=h′)]T-(h′)}-1,
where, according to conclusions of Arena and Pavone [4], the bases of the triangular storms may be considered constant (b¯) with respect to a, because a and b are stochastically independent to each other;  alternatively the Dirac distribution may be considered [3], because it is slightly conservative: in this case the function b¯(a) may be represented by means of a regression, as given in the Appendix.Figure1 shows the return period R≥2′(H), which has been calculated by considering functions defined in the Appendix. Values of parameters used for calculation are defined in Table 1.Table 1
Parameters of the Weibull significant wave height distribution (A4) and of the base height regression (A5) for some locations (from [25]).
BuoyLocation
C
1 (hour)
C
2   (m−1)
u
w (m)
h
l (m)
NDBC 46004 51.0°N–136.0°W
110.25
6.43
·
1
0
-
2
1.484
2.489
0.65
NDBC 44008 40.5°N–69.4°W
75.6
5.57
·
1
0
-
2
1.100
1.146
0.55
RON Crotone 39.0°N–17.2°E
91.84
3.48
·
1
0
-
2
0.956
0.590
0.08Figure 1
Return periodsR≥2′(H), R≥1(H), and R1(H), calculated from data of RON buoy of Crotone (Italy), with the parameters defined in the Appendix.


## 3. An Alternative Approach for Calculation ofR≥2(H)

A different approach for calculation of the return periodR≥2(H) was given by Arena and Pavone [7], by means of (1). In that equation, R≥1(H) represents the return period of a sea storm in which the maximum wave height exceeds the threshold H, and R1(H) represents the return period of a sea storm during which one wave only with crest-to-trough height larger than a fixed threshold H occurs.For the  calculation ofR≥1(H), the solution given by Arena and Pavone [4] has been applied:
(10)R≥1(H)=-{∫0∞dp(Hs=a)daab¯×[1-exp(b¯a∫0a1T-(h′)×ln[1-P(H;Hs=h′)]dh′1-exp(b¯a∫0a1T-(h′)b¯a)]da}-1.
The return period R1(H) is calculated by considering the following expression [7]:
(11)R1(H)=-{∫0∞dp(Hs=a)daab¯×[∫0aP(H;Hs=h)T-(h)[1-P(H;Hs=h)]dh]·exp[b¯a∫0aln[1-P(H;Hs=h′)]T-(h′)dh′]da}-1.


## 4. The Encounter Probability

In general, the occurrences of severe storms are assumed to form a homogeneous Poisson process. Then, the probability that in the given time intervalL (which, in the design of ocean structures, may be considered equal to the lifetime of the structure) at least a sea storm with given properties occurs, may be written as
(12)P(L,R)=1-exp(-LR),
where R is the return period of the considered storm.Then, if we assume that(a)
the occurrences of the sea storms in which the highest wave is larger thanH,(b)
the occurrences of the sea storms in which just a wave higher thanH occurs,(c)
the occurrences of the sea storms in which at least two waves with height larger thanH occurwill represent the Poisson processes, we have, for example, that the probability that during L at least a storm will occur, with the maximum wave height larger than H is
(13)Pa[L,R≥1(H)]=1-exp[-LR≥1(H)].
Equation (13) gives the probability that the maximum wave height in the lifetime (time) L will be greater than H.More in general, if we define(i)
Πa(L,N) as the probability that N occurrences of the process (a) will occur in the time span L,(ii)
Πb(L,N) as the probability that N occurrences of the process (b) will occur in the time span L,(iii)
Πc(L,N) as the probability that N occurrences of the process (c) will occur in the time span L,(the processes a, b, and c being defined above), it follows that
(14)PI=1-Πa(L,0)-Πb(L,1)Πc(L,0)
represents the probability that the second wave, in order of height, during the time L will be higher than H;
(15)PII>Πc(L,1)Πb(L,0)
represents the probability that the second wave, in order of height, higher than H, during the time L, will occur in the same storm in which the highest wave happens;
(16)PIII≡PIIPI>Πc(L,1)Πb(L,0)1-Πa(L,0)-Πb(L,1)Πc(L,0)
represents the probability that, given that the second wave in order of height during L is higher than H, it will belong to the same storm of the highest wave. In other words, PIII(≡PII/PI) is the probability that the two highest waves in the time (lifetime) L will occur during the same sea storm.Note that bothPII and PIII have been defined by considering the stochastic independence of the processes (b) and (c).Finally, if the Poisson processes are considered, the probabilities defined in this section may be calculated as(17)Πa(L,0)=exp[-LR≥1(H)],Πb(L,0)=exp[-LR1(H)],Πc(L,0)=exp[-LR≥2(H)],Πb(L,1)=LR1(H)exp[-LR1(H)],Πc(L,1)=LR≥2(H)exp[-LR≥2(H)].
Figures 2 and 3 show the probabilities of occurrence Pa[L,R≥1(H)], PI, PlII, and PlIII, which are calculated from data of RON buoy of Crotone (Italy) and of NOAA-NODC 46004 buoy, for L=50 years. Note that PlII and PlIII are defined by right hand side of (15) and (16), respectively: they represent the lower bound of probabilities PII and PIII, respectively.Probabilities of occurrencePa[L,R≥1(H)] (11), PI (12), and the lower bound of probabilities PII (13) and PIII (14), which are indicated as PlII and PlIII, respectively, for L=50 years. (a) RON buoy of Crotone (Italy); (b) NOAA-NODC 46004 buoy.
(a)
(b)Figure 3
See caption of Figure2. Data of NOAA-NODC 44008 buoy.It is interesting to note that if, for a fixed value ofL, a large value of the probability is considered (0.8-0.9, e.g.), the height value achieved from PI curve is 3-4% smaller than wave height given by Pa[L,R≥1(H)] probability. This difference increases as smaller values are considered of probability.For example, if data of Figure2 are considered, we have that(1)
for a value of probabilityP=0.8, we find a wave height equal to 25.3 m from Pa[L,R≥1(H)] and to 24.4 m from PI (the ratio is equal to 0.964); in other words, for a lifetime L=50 years and a probability equal to 0.8, we find a wave height, which is maximum in its own storm, equal to 25.3 m; the height of the wave that in 50 years will be the second one in height, for P=0.8, will be equal to 24.4 m;(2)
for a value of probabilityP=0.2, we find a wave height equal to 28.6 m from Pa[L,R≥1(H)] and to 26.8 m from PI (the ratio is equal to 0.939);(3)
for a value of probabilityP=0.05, we find a wave height equal to 30.9 m from Pa[L,R≥1(H)] and to 28.4 m from PI (the ratio is equal to 0.919).The probabilityPIII, in conditions (1)–(3), will be greater than 0.05, 0.29, and 0.49, respectively.



---

*Source: 416212-2013-03-21.xml*