Waves in Shallow Water
Waves in Shallow Water
Waves in shallow water off west coast of India during the onset of
summer monsoon
V. Sanil Kumar1 , C. Sajiv Philip1 , and T. N. Balakrishnan Nair2
1 Ocean Engineering Division, National Institute of Oceanography (Council of Scientific & Industrial Research), Dona Paula,
Goa 403 004, India
2 Indian National Centre for Ocean Information System (Ministry of Earth Sciences), P.O. 21, IDA-Jeedimetla,
Abstract. The wave growth characteristics during the onset normally less than 1.5 m during rest of the period (Kumar and
of summer monsoon in a swell dominated open ocean at a Anand, 2004). The wave energy spectra generally have more
location off the west coast of India at 14 m water depth are than one peak depending on the generation source and also
studied. 67% of the measured waves are due to the swells the waves at any point in the ocean will be sum of local wind
arriving from south and south-west and the balance was due sea and swell components propagating from distant storms
to the seas from south-west to north-west. Wave age of the (Hanson and Phillips, 1999). Kumar et al. (2003) found
measured data indicates that the measured waves are young that along the Indian coast, about 60% of the wave spec-
sea with presence of swells. Even when the wind speed re- tra observed was multi-peaked and they were mainly single
duced to less than 3 m/s, significant wave height more than peaked when the Hm0 was more than 2 m. The double peaked
2 m is present due to the swells in the Arabian Sea. The max- spectra observed were mainly swell dominated with average
imum wave height increased from 2 to 8 m within 60 h. The value of the ratio of the two spectral peaks around 0.6 and
mean wave directions at the high frequencies align with the the average value of the difference between peak frequen-
westerly wind direction and gradually shift to south-westerly cies around 0.09 Hz. Even though the spectra were having
swell direction at low frequencies during the wave growth. two peaks, Kumar et al. (2003) found that Scott and Scott-
The strong westerly winds present between longitude 72◦ Wiegel spectra estimate the maximum spectral energy den-
and 72.5◦ at latitude 12.5◦ has created the high waves (Hm0 sity reasonably well except during the high wave along the
upto 5.65 m) during the measurement period. Indian coast. Off Goa, TMA spectrum (Bouws et al., 1985)
matched the observed sea spectra well when the sea breeze
Keywords. History of geophysics (Ocean sciences)
was active (Neetu et al., 2006) and during this period, Hm0
was less than 2 m. For high waves, it was found that TMA
model, without modifying the model parameters; over pre-
1 Introduction dict the maximum spectral energy, similar to the JONSWAP
model (Hasselmann et al., 1973). Kumar and Kumar (2008)
The wave height along the west coast of India is influenced found that at different water depth (12 to 70 m) around the In-
by the wind conditions in Indian Ocean and Arabian Sea. dian coast, the JONSWAP spectra with the modified spectral
AnGeo Communicates
Wind conditions in the Indian Ocean and Arabian Sea in- parameters represented the measured spectra well for Hm0
crease during mid May and high winds are observed upto more than 2 m.
September which is the period of the summer monsoon.
Along the west coast of India, significant wave height up to A number of studies were carried out to understand wave
6 m is reported during the summer monsoon period (Kumar generation and wave growth (Phillips, 1957; Miles, 1957;
et al., 2006) and the significant wave height (Hm0 ) will be Donelan and Hui, 1990). Young (1997) studied the effect of
fetch length on the growth rate in a finite depth and found that
for longer fetches, the growth rate decreases in finite depth
Correspondence to: V. Sanil Kumar compared to deep water. Wave characteristics during growth
([email protected]) and decay off the Goa coast at 80 m water depth during
6 (a)
Total concept of statistically stationary wave heights was proposed
Sea
Swell
by Longuet-Higgins (1952) and as per this concept; the ratios
4
Hm0 (m)
6 that for shallow water waves the upper limit value for the
4
ratio of Hmax to water depth is 0.55 and is less than that of-
2
143 144 145 146 147 148 149 150 151 152 153 154 155
ten used in engineering practice (0.78). Massel (1996) re-
ported that the ratio of Hmax to water depth is less than 0.6.
300 (c)
The value of 0.57 observed in the present study is similar
Direction (deg)
270
240
to the observation of Massel (1996) and is slightly higher
210 than the value proposed by Nelson (1994). For the present
180 data, significant wave height based
√ on zero crossing analysis
143 144 145 146 147 148 149 150 151 152 153 154 155 (H1/3 ) was found to be √ 3.75 m0. In practice, Hm0 is op-
90 erationally defined as 4 m0. For deep-water narrow band
(d)
70 spectra, H1/3 and Hm0 are equivalent estimates of significant
Percentage
50
wave height (Sarpkaya and Isaacson, 1981). When a nar-
row frequency spectrum is assumed, all characteristic wave
30
heights are theoretically proportional to the standard devia-
10
tion of the water surface. Goda (1979) analysed field data
143 144 145 146 147 148 149 150 151 152 153 154 155
JULIAN DAY 2006
and found that for wind-driven waves in deep water, the ra-
Figure 2. Variation of (A) Significant wave height, sea and swell height, (B) Average mean
tio is approximately 3.8 instead of the narrow-band value of
Fig. 2. Variation
wave period, of (a)
mean wave period Significant
of sea, waveof swell
mean wave period height, seawave
(C) Mean anddirection
swell 4 which gives a reduction of 5%. The reduction in present
height. (b)
of swell and Average
mean mean
wave direction wave
of sea period,ofmean
(D) Percentage sea and wave
swell period of sea, case is 6.2% and the reduction is due to the presence of sea
mean wave period of swell. (c) Mean wave direction of swell and and swells in the measured data which makes the wave spec-
mean wave direction of sea. (d) Percentage of sea and swell. trum relatively broad banded. For wind-driven sea waves,
17 the assumption of a narrow-banded frequency spectrum is no
longer valid. In deep water, the Rayleigh distribution still
The logarithmic approximation is more accurate and univer- holds to a very good approximation for the zero-crossing
sal, but the Prandtl 1/7 approximation (Streeter et al., 1998) wave heights (Longuet-Higgins, 1980), but the ratios of wave
is more convenient to apply and commonly used (Kamphuis, height to standard deviation have to be reduced to account for
2001; Resio et al., 2002). Also the logarithmic approxima- the finite frequency bandwidth. We have found H1/10 /H1/3
tion requires information on the roughness parameter which to be 1.25 and were slightly lower than the Rayleigh value of
is unknown for the study location. 1.275 due to the shallow water effects.
Guedes Soares et al. (2004) showed that in addition to in-
dividual wave parameters, the asymmetry of waves can be
3 Results and discussions described by the statistics of the time series, such as the
skewness and kurtosis. Linear sea states will have no skew-
To identify the swell components from the measured data, the ness and the positive skewness value indicate that the wave
locally generated waves and the swell were separated. Hm0 crests are bigger than the troughs. Present study showed that
upto 5.65 m was recorded during the measurement period the skewness varied from −0.12 to 0.22 with a mean value
(Fig. 2a). Hm0 increased from 2 to 5.65 m in 29 h. The in- of 0.05. The high waves (Hm0 > 4 m) were associated with
crease of Hm0 from 3 to 5.65 m took 15 h and the decay from high skewness value (≈0.2). The kurtosis varied from −0.37
5.65 to 3 m took only 7.5 h. Mean T m02 was 7 s (Table 1) to 0.77. The correlation between the abnormality index (ra-
and the mean period of swell was 10 s (Fig. 2b). The mea- tio of Hmax to Hm0 ) and the kurtosis of the sea surface el-
sured waves were predominantly swells arriving from south evation was 0.6. The dependence of the abnormality index
and south-west and seas were from south-west to north-west on the kurtosis agrees with the results of Guedes Soares et
(Fig. 2c). During the measurement period, 67% of the wave al. (2004). The mean value of the sea surface elevation dur-
height was due to the swells and the remaining was due to ing 30 min interval varied from −0.003 to 0.003 m and the
the seas (Fig. 2d). Good positive correlation (correlation co- average value for all the data is 0.1×10−4 m indicating that
efficient = 0.96) was observed between Hm0 and Hmax . The the height of crest and the trough were nearly equal.
Table 1. Range and average value of wave parameters during the onset of summer monsoon at 14 m water depth.
Me
0.25
Tm 02= 5.4 H m0
2
1 2 3 4 5 6
290
(a)
270
οp
250
230
143 145 147 149 151 153 155
0
143 145 147 149 151 153 155
0.14
fp
0.1
(c)
Figure 4. Variation in sea surface elevation during 27 to 29 May 2006 (for 62 h)
Fig. 4. Variation in sea surface elevation during 27 to 29 May 2006 0.06
143 145 147 149 151 153 155
(for 62 h). 0.6 (d)
nu
0.5
18
a=3.54; b=0.6 for Lake Huron
0.4
a=4.45; b=0.45 for Lake Ontario
143 145 147 149 151 153 155
a=3.94; b=0.38 for North Sea 4
(e)
a=4.04; b=0.47 for Dubai 3
Qp
The variation of T m02 with Hm0 was not following the re- 1
lation proposed in SPM (1984). The high waves (Hm0 =5 143 145 147 149 151 153 155
JULIAN DAY 2006
to 6 m) were associated with mean wave period of 7 to 8 s
and the high wave period (8 to 9 s) was associated with Hm0 Figure 5. Variation of (A) peak wave direction, (B) wind speed, (C) peak frequency, (D)
Fig.
spectral5.narrowness
Variation of and
parameter (a)(E)spectral
peak wave direction,
peakedness parameter (b) wind speed,
varying from 2.5 to 3.5 m (Fig. 3) due to the fact that the
(c) peak frequency, (d) spectral narrowness parameter and (e) spec-
measurements were carried out in open ocean and swells tral peakedness parameter.
that have originated from a distant storm might have decayed
19
with time resulting in low wave height whereas the wave pe-
riod remains high. uninterrupted wind wave measurements and found that the
The variation in sea surface elevation during 27 to 29 May process of shifting toward lower frequency is not persistent
2006 is presented in Fig. 4. The wave height increased from and it is more of an incipient occurrence at the early stage
2 to 8 m during 60 h. The peak wave direction (mean wave of the growth. As the waves continue to grow, the shifts be-
direction corresponding to the maximum spectral energy) de- come less conspicuous and proceeded to coalesce into a uni-
creased from 270◦ to 240◦ during the wave growth period fied predominant peak frequency during the vigorous wave
(Fig. 5a). The peak mean wave directions at the high fre- growth stage.
quencies align with the westerly wind direction and gradu- During the measurements, the spectral narrowness pa-
ally shift to south-westerly swell direction at low frequen- rameter (ν) has a mean value of around 0.5 (Table 1) and
cies. During the study period, the wind speed (U10 ) ranged is marginally high (around 0.55) when bimodal nature of
from 0.1 to 13.7 m/s (Fig. 5b). The value of Hm0 more than the spectrum consisting of low-energy swell from the north
2 m when the local wind speed was less than 3 m/s indicate Indian Ocean and a high frequency local component was
the presence of swell along the west coast of India. present. When the wave height was high, ν was relatively
The shifting of peak frequency towards lower frequen- small with value around 0.44 (Fig. 5d) indicating that the
cies (from 0.12 to 0.08 Hz) is noticed during 15 h on 24 wave spectrum was relatively narrow banded during high
May to 8 h on 29 May 2006, i.e. during the initial stage of waves. The spectral width parameter (ε) varied from 0.74
wave growth (Fig. 5c). The peak frequency is uniform dur- to 0.85. During the high waves, the value of ε was around
ing 8 to 23 h on 29 May 2006. The shifting of peak fre- 0.78. The value of the spectral peakedness parameter (Qp )
quency towards the lower frequency is recognized as one varied from 2.5 to 3 for high waves (Fig. 5e). Large values
of the wave growth characteristics. As a result of the dy- of Qp also indicate the wave spectrum is narrow banded for
namical processes the peak shifts to lower frequencies and high waves.
the spectrum broadens. Liu et al. (2007) presented complete
episode of wave growth process from 10 h of continuous and
1000 where the spectral components below this value are classi-
fied as swell and above as wind sea. In the present case,
only 3% of the waves have inverse wave age more than 0.83
100
(Fig. 6) which indicates that the pure wind sea was very less.
In the open ocean, where most wave systems are close to be-
Non-dimensional energy (E)
During the high waves, the fetch estimated based on Goda, Y.: On the Methodology of Selecting Design Wave Height,
JONSWAP peak-frequency-fetch relationship is from 80 to in: Proceedings 21st Coastal Engineering Conference, ASCE,
100 km and it is found that strong westerly winds present at 899–913, 1988.
between longitude 72◦ and 72.5◦ along the latitude 12.5◦ has Goda, Y.: Numerical experiments on wave statistics with spec-
created the high waves. The propagation speed of the wave tral simulation, in: Report Port and Harbour Research Institute,
Japan, 9, 3–57, 1970.
is between 9 and 10 m/s.
Goda, Y.: A review on statistical interpretation of wave data, in:
Acknowledgements. Authors thank Satish R. Shetye, Director, NIO Report of the Port and Harbour Research Institute, Japan, 18, 5–
and Satheesh C. Shenoi, Director, INCOIS for providing facilities 32, 1979.
and encouragement. Also thank P. Pednekar, K. Ashok Kumar, Guedes Soares, C., Cherneva, Z., and Antao, E. M.: Steepness and
R. Gowthaman, Ganesh N. Naik and Manohar Manohar for the help asymmetry of the largest waves in storm sea states, Ocean Eng.,
during the wave data collection. The data on wind speed and direc- 31, 1147–1167, 2004.
tion recorded by NDBP buoy was provided by National Institute of Hanson, J. L. and Phillips, O. M.: Wind sea growth and dissipation
Ocean Technology, Chennai. This paper is NIO contribution num- in the open ocean, J. Phys. Oceanogr., 29, 1633–1648, 1999.
ber 4709. Hasselmann, K., Barnett, T. P., Bouws, F., Carlson, H., Cartwright,
Topical Editor P. Drobinski thanks H. Branger for his help in D. E., Enke, K., Ewing, J. A., Gienapp, H., Hasselmann, D. E.,
evaluating this paper. Krusemann, P., Meerburg, A., Müller, P., Olbers, D. J., Richter,
K., Sell, W., and Walden, H.: Measurements of wind-wave
growth and swell decay during the Joint North Sea Wave Project
References (JONSWAP), Dtsch. Hydrogr. Z. Suppl., A8(12), 95, 1973.
Kalnay, E., Kanamitsu, M., Kistler, R., Collins, W., Deaven, D.,
Allen, M. and Callaghan, J.: Extreme wave conditions for the Gandin, L., Iredell, S., Saha, G., White, J., Woollen, Y., Zhu,
south Queensland coastal region. Environment technical report A., Leetmaa, B., Reynolds, M., Chelliah, W., Ebisuzaki, W.,
(Queensland. Environmental Protection Agency), 32, 1–7, 1999. Higgins, W., Janowiak, J., Mo, K. C., Ropelewski, C., Wang,
Aparna, M., Shetye, S. R., Shankar, D., Shenoi, S. S. C., Mehra, P., J., Leetmaa, A., Reynolds, R., Jenne, R., and Joseph, D.: The
and Desai, R. G. P.: Estimating the seaward extent of sea breeze NCEP/NCAR 40-year reanalysis project, B. Am. Meteorol. Soc.,
from QuikSCAT scatterometry, Geophys. Res. Lett., 32, L13601, 77(3), 437–471, 1996.
doi:10.1029/2005GL023107, 2005. Kamphuis, W.: Introduction to Coastal Engineering and Manage-
Bala Subrahamanyam, D., Sen Gupta, K., Sudha Ravindran, and ment, World Scientific Publishing Company, 1–472, 2001.
Praveena Krishnan: Study of sea breeze and land breeze along Kudryavtsev, V. and Makin, V.: Impact of swell on the marine atmo-
the west coast of Indian sub-continent over the latitude range spheric boundary layer, J. Phys. Oceanogr., 34, 934–949, 2004.
15◦ N to 8◦ N during INDOEX IFP-99 (SK-141) cruise, Curr. Kumar, V. S., Anand, N. M., Kumar, K. A., and Mandal, S.: Multi-
Sci. India, 80, 85–88, 2001. peakedness and groupiness of shallow water waves along Indian
Barstow, S. B. and Kollstad, T.: Field trials of the directional wa- coast, J. Coastal Res., 19, 1052–1065, 2003.
verider, in: Proceedings of the First International Offshore and Kumar, V. S. and Anand, N. M.: Variations in wave direction es-
Polar Engineering Conference, Edinburgh, III, 55–63, 1991. timated using first and second order Fourier coefficients, Ocean
Bouws, E., Günther, H., Rosenthal, W., and Vincent, C. L.: Similar- Eng., 31, 2105–2119, 2004.
ity of the Wind Wave Spectrum in Finite Depth Water. 1. Spectral Kumar, V. S.: Variation of wave directional spread parameters along
Form, J. Geophys. Res., 90(C1), 975–986, 1985. the Indian coast, Appl. Ocean Res., 28, 98–103, 2006.
Cartwright, D. E. and Longuet-Higgins, M. S.: The statistical dis- Kumar, V. S. and Kumar, K. A.: Spectral representation of high
tribution of the maxima of a random function, Proc. Roy. Soc. shallow water waves, Ocean Eng., 35, 900–911, 2008.
London, A 237, 212–232, 1956. Kumar, V. S., Pathak, K. C., Pednekar, P., Raju, N. S. N., and
Dattatri, J., Raman, H., and Jothi Shankar, N.: Height and period Gowthaman, R.: Coastal processes along the Indian coastline,
distribution for waves off Mangalore Harbour – west coast of Curr. Sci. India, 91(4), 530–536, 2006.
India, J. Geophys. Res., 84(C7), 3767–3772, 1979. Liu, P. C., Chen-Han Tsai, C. C., and Chen, H. S.: On the growth
Donelan, M. A.: The effect of swell on the growth of wind waves, of ocean waves, Ocean Eng., 34, 1472–1480, 2007.
Johns Hopkins APL Technical Digest, 8(1), 18–23, 1987. Long, C. E. and Resio, D. T.: Wind wave spectral observation in
Donelan, M. A., Skafel, M., Graber, H., Liu, P., and Schwab, D.: On Currituk sound, North Carolina, J. Geophys. Res., 112, C05001,
the growth rate of wind-generated waves, Atmos. Ocean, 30(3), doi:10.1029/2006JC003835, 2007.
457–478, 1992. Longuet-Higgins, M. S.: On Statistical Distribution of the Heights
Donelan, M. A. and Hui, W. H.: Mechanics of ocean surface waves, of Sea Waves, J. Mar. Res., 11, 245–266, 1952.
in: Surface Waves and Fluxes, edited by: Geernaert, G. and Longuet-Higgins, M. S.: On the Distribution of the Heights of Sea
Plant, W., Kluwer Academic Publishers, The Netherlands, 1– Waves; Some Effects of Nonlinearity and Finite Band Width, J.
336, 1990. Geophys. Res., 85(C8), 1519–1528, 1980.
Donelan, M., Dobsen, F., Smith, S., and Anderson, R.: On the Mandal, J. C. and Halder, S. R.: Sea breeze like cloud-free zones
dependence of sea surface roughness on wave development, J. during monsoon months, Mausam, 43, 163–168, 1992.
Phys. Oceanogr., 23, 2143–2149, 1993. Massel, S. R.: On the largest wave height in water of constant depth,
Drennan, W., Graber, H., and Donelan, M.: Evidence for the ef- Ocean Eng., 23, 553–573, 1996.
fects of swell and unsteady winds on marine wind stress, J. Phys. Miles, J.: On the generation of surface waves by shear flows, J.
Oceanogr., 29, 1853–1864, 1999.
Fluid Mech., 3, 185–204, 1957. Sarpkaya, T. and Isaacson, M.: Mechanics of wave forces on off-
Neetu, S., Shetye, S. R., and Chandramohan, P.: Impact of sea shore structures, Van Nostrand Reinhold, New York, 491–492,
breeze on wind-seas off Goa, west coast of India, J. Earth Syst. 1981.
Sci., 115, 229–234, 2006. Shore Protection Manual: U.S. Army Coastal Engineering Research
Nelson, R. C.: Depth limited design wave heights in very flat re- Center, Department of the Army, Corps of Engineers, US Govt.
gions, Coastal Eng., 23, 43–59, 1994. Printing Office, Washington, D.C., USA, vols. 1 and 2, 1984.
Phillips, O. M.: On the generation of waves by turbulent wind, J. Streeter, V. L., Wylie, E. B., and Bedford, K. W.: Fluid Mechanics.
Fluid Mech., 2, 417–445, 1957. McGraw-Hill, Singapore, 1–740, 1998.
Portilla, J., Ocampo-Torres, F. J., and Monbaliu, J.: Spectral Par- Thompson, T. S., Nelson, A. R., and Sedivy, D. G.: Wave group
titioning and Identification of Wind Sea and Swell, J. Atmos. anatomy, in: Proceeding of 19th conference on Coastal engineer-
Oceanic Technol., 26, 117–122, 2009. ing, American Society of Civil Engineers, 1, 661–677, 1984.
Rao, C. V. K. P. and Baba, M.: Observed wave characteristics during Violante-Carvalho, N., Ocampo-Torresb, F. J., and Robinson, I. S.:
growth and decay: a case study, Cont. Shelf Res., 16, 1509–1520, Buoy observations of the influence of swell on wind waves in the
1996. open ocean, Appl. Ocean Res., 26, 49–60, 2004.
Resio, D., Bratos, S., and Thompson, E.: Meteorology and wave Young, I. R.: The growth rate of finite depth wind-generated waves,
climate, in: Coastal Engineering Manual, Part II, Hydrodynam- Coastal Eng., 32, 181–195, 1997.
ics, edited by: Vincent, L. and Demirbilek, Z., U.S. Army Corps
of Engineers, Washington, D.C., Chapter II-2, Engineer Manual
1110-2-1100, 2002.