Technical Note: A revised parameterisation of the
Jonswap spectrum
R. M. ISHERWOOD
Orcina Ltd, Plumpton Hail, Ulverston, Cumbria LA12 70N
Equations are derived for estimating the governing parameters of the Jonswap spectrum for a sea
state specified in terms of significant waveheight and average period, Two of the equations are
obtained mathematically from the Jonswap equation and are of general validity. The third is derived
from empirical data published by Houmb and Overvik" and applies only to North Sea conditions.
Key Words: sea state, wave spectrum, Jonswap spectrum, mathematical analysis, computer analysis,
offshore engineering.
INTRODUCTION
The Jonswap spectrum formulation was derived by Hasse
mann et al? The governing parameters were subsequently
defined by Houmb and Overvik! for engineering use in
terms of significant waveheight and average petiod, the
results being presented in tabular form, The purpose of the
work reported here is to revise Houmb and Overvik’s para
meterisation in the light of certain mathematical properties,
of the Jonswap equation and to express the results in the
form of algebraic expressions for convenience in use.
PRELIMINARY OBSERVATIONS
‘The Jonswap spectrum for wave amplitude 7 is defined by
the equation:
Saf) = a8? (2ay FS exp 1.25¢FFmy
wae |
= 0, for ffm
% for f> fm
a
where Sy(f) = spectral density, g = acceleration due to
gravity, "= frequency in Hz and the remaining parameters
are constant for a given sea state, fm isthe frequency corre-
sponding to maximum S,(f).
OF the five parameters in equation (1), %,~0.07 and
= 0.09 are usually taken as absolute constants leaving
«, fin, 7 to be determined in such a way as to give a spec-
‘rum with the required significant waveheight, Hs, and aver-
age period, 72; .2. we have to obtain three unknown from
‘two known parameters. Two of the necessary relationships,
can be obtained mathematically from equation (1) but the
third must be based on empirical data.
In the following sections it is shown that a and fn, non-
dimensionalised with respect to Hs and Tz, can be expressed
as functions of 7 only. Using the empirical data of Houmb
and Overvik it is then shown that 7 is a unique function of
Accepted February 1986, Discussion closes March 1986.
141-1187/87/010047-04 $2.00
(© 1987 Computational Mechanics Publication
a single dimensionless parameter combining Hs and Tz.
‘These expressions avoid the need for interpolation between,
tabulated values and are of great practical convenience in
deriving Jonswap spectra with given waveheight and period.
MATHEMATICAL PROPERTIES OF THE JONSWAP
SPECTRUM
Equation (1) may be rewritten in terms of the relative
frequency, r=ffm, noting that Sq(/)a=S\)e to
maintain constant spectral energy:
oe
(Qnfiny*
els] o
This expression is the product of three terms whose signifi.
cance is as follows:
SO) = fin-5y(N) PSexp(-1.257°4)
og
( © is independent of frequency and is a scale
apne OOP
factor.
(8) FSexp(~1.25r4) isthe spectrum shape term adopted
by the Jonswap. investigators from Pierson and
Moskowitz?
lS]
ii) 20° |is the Jonswap peak enhancement factor.
We adopt the usual definition of spectral moments, viz
i Sy) ar @
and use the standard results that significant waveheight, Hs,
and average zero up crossing period, Tz, are given by
Hs = 4.0Jing @
n -[2 ©
Applied Ocean Research, 1987, Vol. 9,.No, 1 47A revised parameterisation of the Jonswap spectrum: Technical Note: R. M. Isherwood
= (fm Ty and hence
©
Combining equations (2) and (3) we find that
ag?
amy °) aM
where fy is @ function of and ¢ only, This leads immedi-
ately to the result from equation (6) that
m,
f-Tz = 831, 9) @
and from equation (4), with some manipulation, that
eg?
Tepe HOD o
where gi, g2 are aguin functions of and o only. The
dependence of @ on the initially unknown fn can be
removed by substituting for fm from equation (8). Defining
an equivalent wave steepness as
2aHts
ett
(0)
by analogy with monochromatic waves in deep water we
find, after further manipulation, that
ay
Assuming 0, and 0p are constants then g1, f2, £3 Teduce to
functions of 7 only. ot
Equations (8) and (11) then give fim and ¥ in terms of
the known parameters Hs and Tz and the one remaining
unknown parameter, 7.
als? =a37,
EVALUATION OF FUNCTIONS FOR ay fm AND
The functions gy, g3 are complex, difficult and perhaps
impossible to integrate analytically. They can, however, be
evaluated numerically and this has been done for a range of
7 values. The upper limit of integration in evaluating the
moments was ©, achieved by integrating numerically up to
4 predetermined truncation point and adding the contribu-
tions from the tail of the spectrum using analy tically derived
expressions. The following curves have been fitted by
regression analysis to the results so obtained:
Jn Tz = 0.6063 + 0.11647" ~ 0.012247 ay
-85 ——
+83
82
81
-8
73
78
7?
fm. Tz
+76
175 y
+24 2)
+73
+72 :
71
Re 0 o x
Figure 1. fm +Tz against -y— the continuous line represents the theoretical relationship of equation (12). The spots are
Jrom How’ and Overvik for comparison
48. Applied Ocean Research, 1987, Vol. 9, No. 1A revised parameterisation of the Jonswap spectrum: Technical Note: R. M. Isherwood
Alphass2
Figure 2. a/s* against y
Houmb and Overvik for comparison
valid for 0.6 <7 <8,0 with standard error = 0.0001.
as? = 2.964 + 0.4788"? — 0.3430 + 0.042259?
a3)
valid for 0.6 <7 < 8.0 with standard error = 0.0002.
Note that equations (12) and (13) are valid for 6, = 0.07
and 0, = 0,09 only.
‘The functions are shown graphically in Figs. 1 and 2
together with the corresponding values obtained from
Houmb and Overvik. In general the values derived empiri-
cally by Houmb and Overvik agree well with those obtained
hire from purely mathematical considerations.
Houmb and Overvik's results were then used to obtain a
relationship between 7, Hs, and Tz. It was found that
plotting 7 on 5, collapsed the data onto a unigue line and
the following curve fit was obtained by regression analysis
y= 10.54 ~ 1.345" — exp(—19 + 3.7755")
fors>0037 4)
7 = 09 + exp(18.86~ 3.675"), fors <0.037
valid for 0.03