Intrinsic modes of variability have a significant role in driving the climatic oscillations in the oceanic processes. In this paper, we investigate the influence of an inter-annual mode of variability, the Indian Ocean Dipole (IOD), on the wave climate of the eastern Arabian Sea (AS). Using measured, modeled and reanalysis wave data and reanalysis wind data, we show that the IOD plays a major role in the variability of wave climate of the study region. Due to the IOD-induced changes in equatorial sea surface temperature and sea level pressure, the winds from the northern AS gets modified and cause inter-annual variability in the wave climate over the eastern AS. The changes in wind field over the AS due to the IOD influence the generation or dissipation of the wave field and hence cause a decrease in northwest short-period waves during positive IOD and an increase during negative IOD.
The northern Indian Ocean (IO) is a unique ocean as compared to northern Atlantic and Pacific oceans. Because of the land-locked northern boundary, the wind pattern in this region shows semiannual reversal and causes boreal summer (June–September) and winter (December–February) monsoons. Strong westerly wind in the equatorial IO is limited to short transition periods between the monsoons – both southwest (SW) and northeast (NE). The equatorial zonal wind reaches a maximum in April–May and in October–November (Hastenrath and Polzin, 2004). The wind in the eastern Arabian Sea (AS) shows a decreasing trend during October and November (Zieger et al., 2014, supplementary information).
Wind waves are a prominent feature of the ocean surface and play a major
role in planning activities in the open ocean and in coastal zones (Anoop et
al., 2015). Hence, comprehensive understanding of the properties of the
waves and their potential changes are the major knowledge required for
sustainable management of both the coastal and offshore region. The wave
climate of the eastern AS shows a large response to seasons and has a maximum
wave height during the SW monsoon season (Chempalayil et al., 2012; Glejin et
al., 2013a; Sanil Kumar et al., 2014; Anoop et al., 2015). Glejin et al. (2012) analyzed the wave parameters in three locations of the eastern AS during the
SW season and found that the wave height increased from south to north.
Shanas and Sanil Kumar (2014) studied the changes in wind speed and
significant wave height (SWH) in the eastern AS by analyzing 34 years data. The
average SWH in the eastern AS during pre-monsoon (FMAM), SW monsoon (JJAS)
and post-monsoon (ONDJ) are about 1, 2.7 and 0.7 m respectively with an
annual average value of
The monsoon wind patterns in the northern IO affect the spatial distribution of
sea surface temperature (SST) in the tropical IO. The SST distribution
differs from that observed in the Pacific and Atlantic which are warmer in the
west (Vinayachandran et al., 2009). In the IO, warm water is on the eastern
side and cold water is on the western side (Vinayachandran et al., 2009).
This SST difference varies during the coupled oceanic and atmospheric
phenomena in the equatorial IO known as the Indian Ocean dipole (IOD) (Saji et
al., 1999). Positive IOD (PIOD) is associated with decreases (increases) of
SST and increases (decreases) of sea level pressure (SLP) over the eastern
(western) tropical IO. The negative phase of IOD is intensification of the
normal condition (Vinayachandran et al., 2009). In the equatorial IO, the IOD
appears as a dominant contributor of SST variability during the boreal fall
season (October–December) (Saji and Yamagatta, 2003). About 12% of the
SST variability in the IO is associated with dipole mode events
(Vinayachandran et al., 2009). There is a phase lag in the SST evolution
between the eastern and western tropical IO (Saji et al., 2003) and the
regions of positive anomaly (during PIOD) vary with year (Vinayachandran et
al., 2009). The Dipole Mode Index (DMI) is a quantitative representation
of the strength of the IOD and is a measure of the anomalous zonal SST gradient
across the equatorial IO. It is defined as the difference between SST
anomaly in a western (60–80
The tropical IO displays strong inter-annual climate variability associated with the El Niño–Southern Oscillation (ENSO) and the IOD (Murtugudde et al., 2000; Slingo and Annamalai, 2000). Baquero-Bernal et al. (2002) found that the IOD shows good correlation with the ENSO in the equatorial Pacific Ocean. However, the correlation between the strength of ENSO and IOD is not linear (Shinoda and Han, 2005). IODs co-occurring with ENSO are forced by a zonal wind shift in the descending branch of Walker circulation in the eastern IO and the process that initiates IODs in the absence of ENSO is not clear (Vinayachandran et al., 2009).
The impact of IOD on the wind pattern in the equatorial IO is examined in the following studies. SST and SLP variation produced by the IOD cause easterly zonal wind anomaly especially in its zonal component in equatorial IO (Reverdin, 1985; Murtugudde et al., 2000; Saji et al., 1999; Webster et al., 1999; Sreenivas et al., 2012). The IOD-forced wind anomalies are maximum in the central equatorial IO (Sreenivas et al., 2012) – significant anomalies appear around June, intensify in the following months and peak in October. The anomalous easterlies weaken the eastward Wyrtki jets (Wyrtki, 1971) in the equatorial IO (Reverdin, 1985). The wind anomaly produced during IOD has longer duration, but ENSO has shorter duration (Rao et al., 2002).
Although the influence of IOD on the wind pattern of the IO has been reported (Saji et al., 1999), the role of this event on the wave climate of the IO has not yet been studied. Glejin et al. (2013c) pointed out the possibility of influence of IOD on the wave climate of the southeast coast of India, but further analysis on this topic in this region has not been carried out. Most of the studies in the past have focused on the influence of IOD on the wind pattern of the equatorial IO. In this paper, we examine the impact of IOD on the surface wind field of the AS and its impact on the wave climate of the eastern AS. The analysis is confined to October, when the impact of the IOD is strongest and wind conditions in the AS are calmest. Figure 1a shows the study area. The data sets used in this study and the details of the numerical models are described in Sect. 2. Section 3 describes results and discussion, and the main findings are summarized in Sect. 4.
The major challenge for wave climate study in the eastern AS is the scarcity
of long-term observational data. In this study we use the
available measured data of the Datawell directional waverider buoy off
Ratnagiri (available from 2010 to 2014) and off Honnavar (available from
2008 to 2014) off the central west coast of India. Among these, Ratnagiri is the
northern and Honnavar is the southern location and these locations are
spaced
Due to a lack of sufficient measured data, we used reanalysis products of the
ECMWF (European Centre for Medium-Range Weather Forecasts), (i) ERA-40 (Uppala
et al., 2005) for the period 1958 to 1978 and (ii) ERA-I (Dee et al., 2011)
for the period 1979 to 2014, to derive the wind climatology. The spatial
resolution of ERA-40 is 1.5
In order to simulate the directional wave spectrum at the buoy location, we
have used two third-generation spectral wave models; WAVEWATCH III (WW3)
4.18 and Simulating Waves Nearshore (SWAN) 41.01. WW3 is the wave model
developed by NOAA/NCEP (Tolman, 1991, 2009) and is based on finite
difference solving of the energy balance equation of the spectral wave
action in the approximation of phase averaging. The coastal wave model SWAN
is a third-generation, phase-averaged numerical wave model for the
simulation of waves in waters of deep, intermediate and finite depth (Booij
et al., 1999). The physical parameterization of model physics of WW3 is
described in several works (e.g., Tolman, 1991, 2009) and that for SWAN by
Booij et al. (1999), Ris et al. (1999) and Bunney (2011). We have
implemented a coarse-resolution WW3 model with a resolution of 0.25
Cross-comparison of ERA-40 and ERA-I for October from 1979 to
2001.
Wave frequencies were discretized over 25 bins on a logarithmic scale from
0.04 to 1 Hz; wave direction was binned into 36 intervals of 10
Correlation and partial correlation of SWH and MWP with DMI and NINO3 time series during October, from 1979 to 2014.
Comparison of hindcast SWH values with measured values. Left panel shows the values during October 2009 at deep-water location and the right panel shows the values during September–November 2011 at shallow-water location (Ratnagiri).
During October and November, the waves in the AS are dominated by
low-amplitude (SWH < 1 m) wind seas (Young, 1994; Chempalayil et
al., 2012; Glejin et al., 2013a), and during this time the DMI reaches its
maximum value. For this particular period we examined the role of DMI on the
wave climate of the eastern AS and observed that during October surface waves in
this region show a response to DMI. We selected six locations (Fig. 1a) along
the eastern AS which are at more than 100 km away from the coast of India
(Table 1). The time series plot of SWH and mean wave period (MWP) with DMI
are shown in Fig. 4a. Correlation and partial correlation of the wave
parameters with DMI and SOI (southern oscillation index) are shown in Table 1
and this shows that SWH is negatively related and MWP is positively related
to DMI. The influence of IOD increases towards the south, but after removing
ENSO (SOI) the correlation values in all locations decrease and maximum
correlation is observed in the central AS (location L4). This implies that
if influences of ENSO are removed, then IODs have more impact on the central
eastern AS than further north or south. The impact of ENSO with and without IOD
was checked here and it is identified that without the IOD, the impact of ENSO is
significant only in the southern part of the eastern AS. From this it is certain
that in the eastern AS region IOD has more impact than ENSO and its effect
is dominant off the central west coast of India. The variability of SWH and MWP
with respect to the DMI index from 1979 to 2013 is shown in time series plot at
the same six locations (Fig. 4a). The anomaly in the SWH varies between
Averaged wind pattern of October from 2008 to 2014
Scatter plot of zonal wind within 14 to 20
SST anomaly (
Wind vector (m s
Modeled directional spectrum off Honnavar and Ratnagiri for selected positive and negative IOD years (first and third columns for positive IOD and second and fourth columns for negative IOD). Corresponding years are shown inside the figures and DMI are shown in Fig. 4. The spectral energy is shown in logarithmic scale base 10.
The composite climatology (using ERA-40 and ERA-I) of October wind pattern (from 1958 to 2014) of AS and a part of the equatorial IO is shown in Fig. 1b. The wave measuring locations (Ratnagiri and Honnavar) in central eastern AS are shown as red dots (Fig. 1a). The eastern AS region shows comparatively strong wind compared to the western AS. The wind from the northern AS passes parallel to the Indian west coast to the eastern equatorial IO after merging with equatorial westerly wind, while over the central AS the wind vectors are from the northerly direction and change to northwesterly before merging with the westerly equatorial winds and northwesterly winds nearer the Indian coast. This pattern of wind is due to the low pressure in the eastern equatorial IO due to the warm water in this region compared to the western equatorial IO (Vinayachandran et al., 2009).
Pure positive IOD and combined events are considered following Aparna et al. (2012). A pure IOD event is that which occurs in the absence of an ENSO event (Rao et al., 2002). A positive (negative) IOD event that co-occurs with an El Niño (La Niña) is a combined IOD event. A pure ENSO event is defined similarly. The combined events are strong only if both El Niño and positive IOD events are strong. Figure 1c and d show the wind pattern during pure positive and combined IOD events respectively. As observed from climatology, these figures illustrate that in the eastern AS the northwesterly wind vectors are modified to northeasterly. Instead of going towards the eastern equatorial IO the wind blows towards the western IO. This can be perceived as more strong in combined event of PIOD with a comparatively weak wind vectors along the southwest of India. During the negative IOD (NIOD) events (Fig. 1e and f) wind blows over the central and eastern AS as northwesterly which is in contrast to the wind pattern during PIOD. This modification in the wind pattern over the central AS to either westward or eastward direction from its general climatology pattern during the PIOD or NIOD events is due to the IOD-induced temperature variability in the equatorial IO.
Measured frequency–directional spectra of surface waves from two locations off the central west coast of India are shown in Fig. 5. The wind pattern of corresponding years is shown in Fig. 5a and year wise SST anomaly in west, east and DMI index for October is shown in Fig. 4a and b. For the years 2008 to 2014, it can be seen that maximum positive DMI index is observed during 2011 (Fig. 4a). During this year the wind in the eastern AS has shifted its direction to northeasterly. Its influence is clearly visible in the frequency–directional spectra of waves at Ratnagiri and Honnavar, due to the decrease in short-period waves from the northwest (NW) direction. This difference in higher energy at Ratnagiri than Honnavar is clearly evident from Fig. 5b and c and is caused by the alteration of wind direction to the NE before reaching Honnavar. So the dissipation of NW waves and weak wind in the Honnavar region causes a decrease in the short-period wave energy. Similar patterns are observed in other years such as 2008, 2012 and 2014 which has got comparatively higher DMI. Endo and Tozuka (2015) classified 2008 as IOD Modoki, and 2011 and 2012 as Canonical IOD years, and hence higher DMI is observed in these years.
In contrast to the above, the only year with negative DMI during 2008–2014 is in 2010. During 2010, the wind pattern is in the NW direction and the dominance of short waves from NW is higher. A similar pattern can be seen in 2009 and 2013 as these years also have low DMI index. Unlike PIOD, during negative IOD, the short-period waves are slightly higher at Ratnagiri than at Honnavar and this may be due to the dissipation of the waves from the NW due to the larger distance traveled by the waves to reach Honnavar than Ratnagiri. From Fig. 4b we can see that the maximum positive anomaly in SST of the western equatorial IO since 2008 is observed in 2009. However, the influences of high SST anomalies are not observed either in the wind pattern or wave climate. This shows that the modification of wind vectors not only depends on the SST anomaly in the west or east separately but also on the DMI index.
To confirm the influence of IOD-induced turning wind pattern (Fig. 1b to f)
on the waves off the west coast of India, we spatially averaged the zonal wind
component within the region where wind shows alteration in its direction
(14 to 20
During the negative phase of IOD, the swell height is slightly less than
that during the positive phase (Fig. 5b and c). The turbulent sea state
generated by the dominance of short-period waves during the negative phase of IOD
leads to an increased decay rate of swell (Ardhuin et al., 2009; Young et al.,
2013). The sign reversal of the air–sea momentum flux depends on a parameter
known as inverse wave age (ratio of sea surface wind speed to the speed of
the wave corresponding to peak frequency) and it is very useful to
understand the present sea state (Grachev and Fairall, 2001; Hanley and
Belcher, 2008; Hanley et al., 2010). Here, we calculated the monthly
(October) average inverse wave age off Ratnagiri and the values are 0.32,
0.11, 0.27, 0.35 and 0.17 from 2010 to 2014. From this it is further evident
that during the positive phase of IOD, the eastern AS is a wave-driven wind
region (inverse wave age < 0.15) – i.e., waves transfer momentum to
wind. During the negative IOD period, the region becomes mixed state (0.15 < inverse wave age < 0.83; both wind-driven wave and wave-driven wind are present). The modeled frequency–directional spectra off
Ratnagiri and Honnavar from 2010 to 2014 are shown in Fig. 5d and e.
Comparison of measured and modeled spectra off Honnavar and Ratnagiri shows
that the spectral energy is overestimated by the model since the ERA-I wind
data used in the study are with a coarse resolution of 1
Since IOD has large influence on SST variability in equatorial IO, here we show the SST of AS and western equatorial Indian Ocean in Fig. 7. We considered comparatively strong PIOD and NIOD events from 1979 to 2014. The PIOD years are 1982, 1994, 1997 and 2006 (left panel of Fig. 7). Among these in 1994, IOD occurs in the absence of ENSO. The NIOD years are 1980, 1984, 1996 and 1998 (right panel). This SST variability produces variability in SLP (Fig. 8). The wind pattern of the respective years for AS and part of equatorial IO are shown in the same figure. Since the NIOD is the intensification of the normal condition of October (Vinayachandran et al., 2009), the SLP in the western equatorial IO is higher and lower in the eastern equatorial IO. Hence, the winds from the northern AS will propagate towards the eastern equatorial IO along the eastern AS. This shows the influence of IOD-induced temperature variability on SLP and wind direction.
Since the performance of the numerical model is good we have carried out further analysis using this model for the period where measured data are not available. Simulated directional spectra for the years 1980, 1982, 1984, 1994, 1996, 1997, 1998 and 2006 are shown in Fig. 9. The first and third columns are for PIOD years off Ratnagiri and off Honnavar respectively; similarly the second and fourth columns are for NIOD years. From the figure it is evident that for the positive phase of IOD, the short-period waves from the NW direction are less at both locations, whereas 2006 shows a conflicting result. From Fig. 8, it is clear that in 2006 the wind pattern in the central AS is northeasterly as observed in the other PIOD events, but along the eastern AS the wind vector extends slightly further south compared to the strong events observed during 1994 and 1997. A similar pattern, but comparatively less spatial extension in wind vectors towards the south, is also evident during 1982 (Fig. 8).
In the case of 1996, which is the strongest NIOD event considered in this
study, we expect a more evident signature in SWH at the study locations of
Ratnagiri and Honnavar. However, in the wave spectrum its signature is
weaker than that observed in 1980 and 1984 (Fig. 9). The wind pattern observed
during the period, in 1996, is northerly down to 15
We analyzed the wind pattern over the AS and examined the influence of IOD events on the October wave climate of the eastern Arabian Sea using reanalysis, observation, and model data sets. Analysis of wind pattern influenced by the IOD has been carried out by removing the influence of ENSO from IOD and vice versa. It has been found that IOD has significant influence on the wave climate off the central west coast of India compared to the northern and southern parts. The decreasing of wave height during PIOD is due to the decrease of short-period waves from the northwest direction. Wind blowing from the northern AS is the major determining factor on wave climate. In general, the climatology of wind pattern in central and eastern AS during October is northerly and northwesterly respectively. During PIOD events, the wind vector modifies its direction and blows as northeasterly because of the IOD-induced SST anomaly in the equatorial Indian Ocean. This brings a weakening of wind field over the central and south-of-central west coast of India. This change in direction of wind pattern in the AS causes a decrease in wave height off the central and southwest coast of India during PIOD, whereas during NIOD events the wind vector turns to northwesterly instead of northerly winds in the composite climatology over the region during October. This leads to an increase of short-period waves in the same region. The influence of IOD on the wave climate mainly depends on the modification of wind field caused by the phases of IOD event. If this wind pattern is absent even during a strong IOD event, then the signature of the IOD on the wave climate is also absent. This alteration of wind pattern mainly depends on the IOD-induced SST variability in the eastern and western equatorial IO, and some other unknown factors may also cause slight modification in this wind – more analysis is required to understand this variability. From this study it is clear that IOD has an impact on the wave climate off the west coast of India especially off the central west coast of India due to the decrease of northwesterly short-period waves.
The authors acknowledge the CSIR, New Delhi for funding the wave measurement at Honnavar and ESSO-INCOIS, Ministry of Earth Sciences, Government of India for funding the wave measurement at Ratnagiri. The director of CSIR-NIO, Goa provided encouragement to carry out the study. The deep water buoy data used for validation of the numerical model were collected by the National Institute of Ocean Technology, Chennai and provided by the Indian National Centre for Ocean Information Services (INCOIS) Ministry of Earth Sciences, Hyderabad. T. M. Balakrishnan Nair, A. Nherakkol and J. Singh provided support during data collection. We thank A. Vedpathak and P. Kulkarni, Centre for Coastal & Marine Biodiversity, Dr. Babasaheb Ambedkar Marathwada University, Ratnagiri for providing the logistics required for wave data collection at Ratnagiri. This work forms part of the Ph.D. thesis of the first author. We thank the topic editor J. M. Huthnance and the two reviewers for their critical review and suggestions, which improved the scientific content of the paper. This work is NIO contribution no. 5866. Edited by: J. M. Huthnance