1. Introduction
Ocean currents are driven by several forces, including tides, density gradients, wind stress and geostrophic deflections [
1]. Understanding ocean currents is fundamental to a diverse range of marine activities [
2]. Therefore, acquiring ocean current data with fine temporal and spatial resolution over a long time to analyze the mechanisms of its movement is important. However, it is difficult to obtain continuous data over a long time relying only on conventional equipment [
3].
Conventional ocean current measuring instruments, such as current recorders and acoustic doppler current profilers (ADCP), measure data at single points. High-frequency radars (HFR) can be developed as an alternative for measuring surface currents at higher temporal and spatial resolutions; HFR is now widely applied [
4,
5,
6]. HFR can provide highly resolved temporal (tens of minutes) and spatial data vital to understanding marine mechanisms over relatively large domains [
7]. HFR-derived ocean surface currents are utilized at multiple spatiotemporal scales along coastal regions to understand the spatial variability of tidal currents [
8,
9,
10] and associated mechanisms [
4,
7], to investigate mesoscale and sub-mesoscale eddy dynamics [
6,
10,
11,
12], to understand coastal upwelling features [
13,
14], and to unveil seasonal surface circulation variability [
5,
15,
16,
17]. Several radar systems are now widely applied, such as Seasonde CODAR, WERA, and OSMAR. They can be categorized into beam-forming and direction-finding types. Beam-forming HFR systems electronically steer a linear phased array of receiving antennas toward a sector of ocean surface, such as OSCR and WERA [
18,
19], while direction-finding HFR systems, such as CODAR, exploit directional properties of loop antennas to determine bearing [
20].
Regular semidiurnal tides are important, and interactions between tide and bathymetry occur over different ranges of depth. Shen et al. [
21] conducted a study on the spatial distribution of surface tidal currents in southwestern Taiwan based on the HFR dataset using the method of quasi-harmonic analysis. They found that tidal currents were affected by semi-diurnal and shallow tidal water constituents. Shallow water tidal constituents are also highly significant off the west coast of Ireland and are investigated herein. Many processes can benefit from knowledge of the maximum currents in a region; these can be assessed in a number of ways, such as probable maximum current velocity (PMCV) and measured maximum current velocity (MMCV). These currents are particularly relevant to offshore engineering projects, such as the Green Atlantic Project proposed by the Irish Electricity Supply Board (
https://esb.ie, accessed on 1 January 2021). In this research, remote sensing HFR was used to determine the surface circulation of a proposed large offshore floating wind farm site. Surface currents measured using the HFR network during the period from December 2015 to December 2016 were analyzed. In order to study the spatiotemporal distribution characteristics and variations of ocean current flows off the western coast of Ireland, we have looked at seasonal variations, the significance of wind impacts, tidal ellipses, rotational spectral analysis, and maximum current analysis. This is important research, as it provides new insights into a region of huge societal potential at early stages of sustainable economic exploitation; when this installation is completed, 1400 MW of power will be generated from hundreds of floating wind turbines. The electricity generated will be used to power the Irish national grid through the development of the largest synchronous compensator in the world.
The structure of this paper is as follows.
Section 2 presents methodologies, including the study area, HFR system, and tide. Results are presented in
Section 3. Discussions and main conclusions are given in
Section 4.
2. Methodologies
2.1. Study Area
Located on the northwest European Continental Shelf, the waters off the Irish west coast span the transitional zone of oceanic to coastal waters, exhibiting characteristics of both coastal and Atlantic behavior. The width of the continental shelf (defined by the 200 m contour) approaches within 100 km of the coast in the region. In general, the seabed slopes rapidly to 80–100 m within 20 km of the coast. From here, the bathymetry plateaus out to the shelf break. Thus, shoaling is significant and shallow water tidal constituents contribute to the hydrodynamics of the region. The Atlantic weather systems strongly influence the region; average wind speeds in January are on the order of 11 m/s, and during June, the average is around 7 m/s [
22]. Consequently, wind plays a significant role in determining residual current circulation in the region [
23,
24]. The meteorological conditions exert significant influence on magnitudes and directions of surface currents locally; surface currents generally contain contributions from wind stress, whereas winter surface currents may actually be dominated by wind forcing in the study area [
25]. Oceanographic data for west coast Ireland is sparse. McMahon et al. [
22] described some oceanographic features of the region. A major feature of Ireland’s western continental shelf is a thermohaline front, which separates Ireland’s coastal waters from offshore waters.
This study focuses on analyses of sea surface currents observed by an HFR system from a full year’s data; the system covered an offshore area with maximum meridional and latitudinal lengths of approximately 100 km and 90 km, (see
Figure 1b). To ensure reliability and integrity of the dataset, sea surface currents monitored for more than 80% of the observation period were considered monitoring points of high density (in total 162) (see
Figure 1b). This data selection criterion is the same as that used by Liu, Zhou, and Wen [
6]. The maximum latitudinal length and meridional lengths of the high-density monitoring area are approximately 40 km and 50 km, respectively.
Figure 1b shows that the majority of the high-density points are distributed in the intersection area between two HFR stations.
Observational datasets include surface currents obtained from the HFR system and tidal water elevations from a tidal gauge close to the study area. Detailed descriptions of each dataset are presented below.
2.2. HFR System
Despite the potential for economic development off the Irish west coast, little commercial exploitation has taken place. Large oil and gas deposits are known to exist on the continental shelf. The region is rich in sea life and has large economic potential for fishing and aquaculture. Studies undertaken in recent years have identified the west coast of Ireland as being in one of the regions with the highest offshore wind and wave energy potentials in the world [
26]. In April 2021, the Irish National Electricity Supply Board (ESB) announced detailed plans to immediately begin the development of a 1.4 GW offshore floating wind farm in the region. This development, Green Atlantic, will include a major hydrogen production facility to store excess wind energy. Post-completion, Ireland will become an exporter of clean energy and will benefit significantly in economics terms from this development.
The study area is located on the west coast of Ireland (9.26°–10.79°W, 52.33°–53.26°N), as is shown for the Atlantic in
Figure 1. It is located in the westerlies and ocean current movements are affected by the North Atlantic Oscillation and low-pressure system [
27,
28].
Figure 1b shows that the bathymetry in the study area ranges from 0–140 m, and the isobaths nearshore are roughly parallel to the coastline. To the west of the study area, the water depth increases suddenly up to approximately 1000 m deep ocean waters. The region with depths of 80–100 m consists of alternating ridges and troughs. Ridges are scattered throughout the northwest of the study area between 80 m and 100 m isobaths and to the west of study area with 120 m isobaths. This complex bathymetry with ridges, troughs, and basins leads to rapid changes in water depth in the study area.
Figure 1b also shows the deployment location of HFR stations (R1–R4) and their area of coverage.
HFR is a shore-based sea state monitoring system. Compared to traditional current measuring instruments, it has several advantages: its high spatial and temporal resolution and wide area of coverage, and that it is on an all-weather system. HFR senses fields of sea surface currents, waves, and sea surface winds using Bragg scattering and Doppler frequency shift effects produced by the interactions between high-frequency radio waves and the sea surface [
27,
29].
The University of Galway deployed two sets of Seasonde HFR stations on Galway Bay in 2011. The four observation stations were connected at the end of 2014 and had monitoring parameters within approximately 200 km of the west coast of Ireland. Characteristics of the radar stations are presented in
Table 1. Surface currents monitored by the HFR system over a complete year (December 2015 to December 2016) were selected for this research. Herein, the authors focused on analyzing spatial and temporal characteristics of the surface flow fields with high coverage density.
Inter-comparison between surface currents observed using HFR and ADCP at Galway Bay was undertaken by O’Donncha et al. [
30]. The accuracy of HFR data, with root mean squared error (RMSE) by 10–12 cm/s, was comparable to datasets in other coastal areas, such as surface currents to the southeast of China and currents on Monterey Bay, located on the west coast of the USA [
6,
31].
2.3. Tidal Current Analysis
2.3.1. Tidal Current Patterns
To investigate tidal hydrodynamics, the major semi-diurnal constituents (M2, S2, and N2), major diurnal constituents (K1 and O1), and shallow water constituents (M4, M6, S4, and MS4) are considered in this study. The major and minor axes of tidal ellipses are determined by the tidal current harmonic constant, and the major axis points in the direction of the maximum velocity of the tidal constituents [
32].
Coefficient F is used to describe characteristics of tidal ellipses and tidal currents. Tidal current patterns were assessed using [
33]:
where
,
, and
are the lengths of the semi-major axes of the tidal ellipses of
O1,
K1, and
M2, respectively.
When
F ≤ 0.5, the area is mainly affected by regular semi diurnal tidal currents; when 0.5 <
F ≤ 2, the area is mainly affected by irregular semi diurnal tidal currents; when 2 <
F ≤ 4, the area is mainly affected by irregular diurnal tides; and when
F > 4, the area is affected by regular diurnal tides [
33].
Rotational characteristics of tidal currents can be expressed by the rotation rate
R [
33]:
where
is the magnitude of the short axis of its tidal ellipse and
is the magnitude of the major axis of the tidal ellipse.
When the absolute value of the rotation rate
R is greater than 0.25, the tidal current has strong rotational flows; when the absolute value of
R is less than 0.25, it indicates that the tidal current is characteristic of reciprocating flows. The positive or negative value of R indicates the direction of tidal ellipse rotation. When
R is positive, the tidal ellipse moves counterclockwise; when
R is negative, the tidal ellipse moves clockwise [
33].
To first explore the spatial characteristics of tidal currents, surface current data were spatially averaged at each monitoring time [
21]. Secondly, harmonic analysis of surface currents was carried out using the t-tide package in Matlab 2014b.
2.3.2. Rotary Spectral Analysis
The Coriolis force is an important force to include when analysing ocean currents. Traditional analysis splits ocean current vectors into zonal (u) and meridional components (v). Power spectra of both zonal and meridional components are computed respectively; the direction of ocean current motion is determined from both components. However, traditional methods cannot reflect the action of the Coriolis force or the influence of the earth’s rotation on sea water while analyzing the movement of sea surface currents [
34,
35,
36].
In order to study the frequency characteristics of ocean current movements, the methods of spectrum analysis were used in characterizing features of different high frequency dependent data. Rotary spectra, which consider movement of ocean currents as a function of u(t) and v(t) and reflect effect of rotation, have been widely used to analyze the characteristics of currents movements [
35,
37]. The rotary spectral analysis method is an effective way to determine the strength of different periodic signals in ocean current motion and may be used to determine the significance of each tidal constituent.
Here, the rotary spectral analysis method is used to analyze ocean current energy [
34,
35,
36]. Firstly, we define current vectors in complex form as [
34]:
The counterclockwise rotation energy spectrum is defined as:
The clockwise rotation energy spectrum is defined as:
The total energy spectrum
S is given by:
where
and
indicate counterclockwise and clockwise rotary spectra, respectively;
indicates the motion of current vector with time;
represents the motion process of the east component of the surface current vector;
indicates the motion process of the northern component of the surface current vector; and
and
represent the auto spectra of
and
, respectively.
indicates the cross spectrum of
and
, Im indicates imaginary part of variable,
W(f) represents the Fourier–Stepian transformation of
,
E represents expectation, and the increment
dW(
f) represents the complex amplitude associated with the differential of the circular frequency (f) [
34,
35,
36].
The rotary spectral method categorizes the total density spectrum into a clockwise spectrum and a counterclockwise spectrum. The difference spectrum, computed by differentiating the clockwise spectrum from the counterclockwise spectrum, can be used as an important index to evaluate patterns of ocean flow vectors. If the value of the difference spectrum is 0, it indicates absolute reciprocating motion of the current on a line segment and that the rotation coefficient of the current vector is 0 [
34,
35,
36]. The absolute value of the rotation rate is the ratio of the differential spectral energy density to the total spectral energy density. If the rotation rate of the current vectors is 1, then the current vector is in complete rotational motion [
35].
In this research, hourly surface current data over the year 2016 obtained from HFR observation system were analyzed as follows. Firstly, since few gaps exist among the observations, missing data for observational points with high coverage density were filled using linear interpolation, a commonly used approach in data filling. Secondly, surface currents were averaged in space. Thirdly, the spatially averaged surface current dataset was used for rotary spectral analysis.
2.3.3. Effects of Bathymetry Variation on Currents
Marine bathymetry influences shallow water constituents of tidal currents. The coefficient of shallow tidal water constituent, G, is used to describe effects of bathymetry in tidal wave propagation and the movements of tidal currents; G is calculated as follows [
33]:
where
is the major axis of the M4 tidal ellipse,
is the major axis of MS4 tidal ellipse, and
is the major axis of the M2 tidal ellipse.
G describes the influence of bathymetry on the propagation of a tidal wave. When
G is greater than 0.04, the effect of the shallow tidal water constituent is considered significant. The higher the value of
G, the more significant the effect of coefficient of shallow tidal water constituents [
33].
2.4. Maximum Ocean Currents
In order to study the main driving force of surface current movement at some particular time, the time of maximum ocean current was selected. PMCVs and MMCVs are widely used in the field of marine engineering; the comparison between PMCVs and MMCVs is often used to study whether wind or tidal forces dominate the movements of ocean currents at the time of the maximum current speeds. PMCVs are estimated vectors of the maximum current speeds based on the outcome of harmonic tidal analysis. MMCVs are a measured vector of the maximum current speeds, in this case based on HFR observations. MMCVs represent extreme current speed at each spatial point during the analysis period [
33]. Values of PMCVs can be obtained from:
where
,
,
, and
represent the semi-major axis vectors of the tidal ellipses of
M2,
S2,
O1, and
K1, respectively.
3. Results
3.1. Mean Surface Flow Fields
To explore the general spatial distribution characteristics of sea surface flow fields in the study area at a seasonal scale, HFR surface currents were averaged at each observation point over a whole year and four seasons separately.
Figure 2a–d show the mean vector maps in spring, summer, autumn and winter, respectively.
Figure 2e shows the mean surface vector fields over the full analysis year.
Figure 2 shows that flow directions and magnitudes of mean surface vectors are significantly different across the four seasons. Flow directions of mean surface vectors are predominantly northern or western. In spring, the distribution of the mean surface flow field is quite disparate, dominant flow directions are northwestern and western (see
Figure 2a). However, flow directions of mean surface vectors in summer, autumn, winter and the whole analysis year are more consistent and are mainly northeast and northwest. In summer, the mean surface flows had northwestern directions only near the north of the study area due to the blocking of the islands (see
Figure 2b). In autumn, the trend of the mean surface flow field is mainly northwest. Northeastern flow trends only occurred at a few points close to the West Atlantic Ocean.
Figure 2d shows that mean surface flows in winter dominantly had a northeastern flow trend. Only in the northern part of the observation area does a northwest flow exist due to blocking of the islands.
Figure 2e shows that northward flow is the dominant trend of annual mean surface flow fields. The annual mean surface flow field presents northward flow in the southern part of study area, while a northwestern flow trend exists in the northern part. This difference may result from the blocking effect of the islands (Aran Islands).
To further analyze characteristics of mean flow fields at a seasonal scale, sea surface currents data (at 162 points) in the study area were firstly spatially averaged, and then the maximum, minimum, and mean values in space were computed; these are presented in
Table 2.
Table 2 shows that the maximum and minimum values of the seasonally averaged surface current velocity occur in winter and spring; they are 54.08 cm/s and 0.26 cm/s respectively. The mean current speed of winter is considerably greater than that of the other seasons. The mean value of averaged surface velocity for the full year is less than the values during autumn and winter.
3.2. Tidal Elevation Analysis
Tidal water levels from the nearest tidal gauge (Galway port, as is shown in
Figure 1b) were used for harmonic analysis using MATLAB’s t-tide package [
38]. Based on the original HFR current data, 67 tidal constituents (with 95% confident interval) were extracted using the method proposed by Halverson and Pawlowicz [
39]. Note that only those tidal constituents with signal to noise ratios (SNR) greater than 2 are considered significant and quantified in this study. The diurnal tidal constituents (K1, O1, P1, Q1), semidiurnal tidal constituents (M2, N2, S2, K2), and shallow tidal water constituents (M4, MS4, M6) are presented in
Table 3.
Table 3 shows that M2 and K1 constituents have the maximum amplitudes of the semidiurnal and diurnal constituents, respectively. The M2 amplitude is approximately 3 times greater than the S2 constituent, 5 times greater than the N2 constituent, and 10 times greater than the K2 constituent. The amplitude of the K1 tidal constituent is approximately 2 times greater than the O1 constituent and 5 times greater than the Q1 constituent. The averaged amplitudes of the semi-diurnal constituents are approximately 11 times greater than the diurnal constituents and 25 times greater than the shallow tidal water constituents. The values of the amplitudes of different shallow tidal water constituents are similar, while the amplitudes of diurnal tidal constituents and semi-diurnal tidal constituents are significantly different.
3.3. Characteristics and Patterns of Tidal Currents
Based on the spatially averaged sea surface current data, the value of the coefficient F calculated, Equation (1), is 0.3009; this indicates that the tidal currents exhibit strong semi diurnal characteristics. To further study the characteristics of tidal currents, harmonic, spectral, and tidal ellipse analyses were performed on HFR data and are presented in the following sections.
3.3.1. Harmonic Analysis on Spatially Averaged Surface Currents
The dominant astronomical and shallow water tidal constituents are presented in descending order of amplitude in
Table 4. Meteorological tidal constituents SA (solar annual constituent) and SSA (solar semiannual constituent) are determined using solar radiation and other meteorological factors on an annual scale [
40]. Since the length of the surface current data analyzed in this research is only one year, it is not possible to undertake an effective analysis of the characteristics of these constituents. Thus, meteorological tidal constituents were not included in this study. In order to show trends in surface current velocity components, time series at six representative points (see
Figure 1b) are shown in
Figure 3. The value of u component ranges from −47 cm/s to 45 cm/s at 6 points and v component ranges from −16 cm/s to 62 cm/s. During the analysis period, there are significant semidiurnal characteristics at the movements of sea surface currents in the study area.
Table 4 shows that amplitude of the M2 tidal constituent is approximately three times and four times greater than those of the S2 and N2 constituents, respectively. The amplitude of the K1 tidal constituent is approximately equal to the S1 and P1 constituents and about twice as large as the O1 constituent. The average amplitudes of semi-diurnal tides are approximately three times that of diurnal tidal constituents, and 16 times that of shallow tidal water constituents. The values of the amplitudes of different shallow tidal water constituents are similar to each other, while the amplitudes of diurnal tidal constituents and semi-diurnal constituents are significantly different from each other. Thus, the energy of the semi-diurnal signals is dominant in the study area compared to other periodic signals. The value of F here is 0.30 (<0.50), indicating that the tidal surface currents had significant contributions from semidiurnal constituents which is consistent with above.
In addition,
Table 4 shows that ocean current movement is significantly affected by the semidiurnal tidal constituents, including M2, S2, and N2; M2 plays a major role. To investigate rotational characteristics of tidal currents, dominant constituents M2 and S2 are selected to calculate rotation rate using Equation (2). The results show that the rotation rates (R) of M2 and S2 are 0.56 and 0.32, respectively. The absolute values of rotation rates are greater than 0.25, and their rotation rates are also positive; this indicates that the tidal ellipses of M2 and S2 exhibit rotational flow and that their direction of rotation is predominantly counterclockwise. Additionally, results indicate that the rotationality of the M2 tidal ellipse is more significant than that of the S2 tidal ellipse.
3.3.2. Rotary Spectral Analysis
To further examine hydrodynamics in the study area, rotational spectral analysis was conducted on spatially averaged surface tidal currents, see
Figure 4. Nowhere in the study area did flows exhibit absolute rotational flows.
Figure 4a,b show clockwise spectra and counter clockwise spectra of currents, and
Figure 4c shows the total spectrum of currents.
The magnitude of the rotary spectrum is significant in the frequency range 0.03 cph to 0.1 cph; thus, rotary spectra in the range 0.03 cph to 0.1 cph are chosen for further analysis,
Figure 4b. Through rotary spectral analysis, dominant tidal constituents are chosen to describe the characteristics of tidal currents. From analysis of the spectra in
Figure 4b, it can be shown that this counterclockwise spectrum of surface currents accounts for 52% (the counterclockwise spectrum at all frequencies/the total spectrum at all frequencies) of the total spectral energy density, while the clockwise spectrum accounts for the other 48% (the clockwise spectrum at all frequencies/the total spectrum at all frequencies). Furthermore, 45% of the frequency of the difference between spectra behaves as a counterclockwise spectrum, and 55% of the frequency behaves as a clockwise spectrum. The counterclockwise spectrum accounts for more of the total spectral energy density and less of the of the difference spectrum. This indicates that the counterclockwise density spectrum is dominant in the total energy density spectrum. However, the frequency distribution of the counterclockwise density spectrum is relatively concentrated; the spectrum peaks are concentrated in a relatively narrow frequency range, while the clockwise density spectrum is dispersed over a broader frequency range. The energy density of the clockwise spectrum is higher than that of the counterclockwise spectrum in most frequencies. The energy density of the semidiurnal tidal constituents, dominated by M2 and S2, accounts for 16% of the total spectral energy. The counterclockwise energy density of the frequency corresponding to the M2 tidal constituent accounts for 13%. In combination with analysis of
Figure 4c, counterclockwise energy spectrum of semidiurnal constituents, especially M2, contributes the most to the total energy spectrum. Spectral energy presents a cusp around 0.067 cph (cycle per hour) in
Figure 4, which corresponds to the period of inertial flow in the study area.
Analysis of tidal patterns and rotary spectra indicates that the pattern of tidal currents in the study area is mainly driven by semi-diurnal tidal forcing. The range of low frequency corresponds to the period of one year, half a year, and several months. The spectral peak density may be affected by wind stress, showing strong monthly and seasonal characteristics.
3.3.3. Spatiotemporal Rotation of Surface Currents
To further investigate the characteristics of tidal currents movement, the spatial distribution of tidal ellipses is analyzed as shown in
Figure 5. Based on the harmonic analysis of the HFR surface currents and tidal water elevations at Galway port, M2 and S2 (astronomical constituents) tidal ellipses are shown in
Figure 5. According to the rotary spectral analysis, the energy density of M2 tidal constituent accounts for the highest proportion of all constituents. In the high-density coverage area, the major and minor axes of the M2 tidal ellipses are generally larger than those of the S2 tidal constituent. The average lengths of the major and minor axes of M2 tidal ellipses are 9.51 cm/s and 4.47 cm/s, respectively; the average lengths of the major and minor axes of S2 tidal ellipses are 3.70 cm/s and 1.14 cm/s, respectively. The mean lengths of the major and minor axes of M2 tidal ellipses are approximately 3 to 4 times greater than those of S2 tidal ellipses. Thus, the effect of the M2 constituent on tidal currents was more significant than that of the S2 constituent. For the semidiurnal tidal signals, M2 and S2, most of the current ellipses mainly exhibit counterclockwise flow. Only one tidal ellipse of M2 is in a clockwise direction; this ellipse is located in the southeast region of the study area. Tidal ellipses of the S2 constituent exhibit counterclockwise flow over the entire study area.
M2 and S2 tidal ellipses mainly exhibit rotating flow; the rotation of M2 ellipses is more pronounced than that of S2 ellipses. As shown in
Figure 5a, the M2 tidal ellipses in the northwest of the study area tend to exhibit rotating flow, whereas the tidal ellipses in the southeast of the study area exhibit reciprocating flow tendency. The spatial distribution characteristics of the tidal ellipses of both S2 and M2 are similar. They show the characteristics of rotating flow in the northwest of the study area and more reciprocating flow characteristics in the southeast direction due to topographic steering. Although S2 ellipses generally exhibit rotating flow, the rotational characteristics are not as pronounced as the M2 ellipses. Considering the two dominant tidal ellipses, their flow structures are influenced by both coastal topography and submarine bathymetry. The tidal ellipses in the northwest of the study area exhibit rotational flow tendencies since they are less affected by the islands, while the tidal ellipses in the east of the study area are inclined to show reciprocating trends because they are more affected by coastal topography.
However, the tidal ellipses in the south of the study area show that whether reciprocating flows are far from or close to the coast, the patterns are different from the findings of Leaman et al. [
41]. This is because the movements of sea surface currents are not only influenced by the coastal topography but also by the submarine topography. Sea surface currents in the area of seabed ridges tend to show reciprocating characteristics, while the currents in the area of seabed slopes show rotational characteristics; this explains why currents in the north of the study show rotational characteristics despite proximity to the islands and why sea surface currents in the south show reciprocating characteristics in the area of seabed ridges regardless of distance from the coast.
In the middle of the study area, bathymetry slope changes slowly compared to the northern and southern parts of the study area. Offshore tidal ellipses (west of the study area) exhibit reciprocating flow tendencies, while tidal ellipses close to shore (east of the study area) exhibit rotating flow tendencies. However, tidal ellipses exhibit reciprocating flow tendencies in the northern and southern part of study area because of energy propagating caused by significant bathymetric slope changes.
Shallow tidal water constituents are caused by abrupt changes in water depth when a tidal wave propagates from the deep sea to shallow water. The shallow tidal water constituent coefficient based on M4, MS4, and M2 tidal constituents can be used to characterize the influence of how the marine topography in an area affects propagation of a tidal wave. Comparison between water depths and the distribution of the coefficients of shallow tidal water constituents at high-density points (see
Figure 6) shows that areas with small coefficients of shallow tidal water constituents correspond to seabed areas with shallow slopes. Conversely, areas with large coefficients of shallow tidal water constituents are mainly located in areas with significant bathymetric changes. However, contours of coefficients of shallow water tidal constituents are not parallel to the coastline or local isobaths. This indicates that there is an abrupt change in water depth at points with high coefficient of shallow tidal water constituents in the south of study area. Two other areas of abrupt change in bathymetry occur in the middle and south of the study area, they are between 60 m and 80 m isobaths (see
Figure 6a).
The coefficients of the shallow tidal water constituent in this study area are on the order of 0.08 (>0.04), this suggests that the shallow tidal water constituent effect is significant. This is important, as it illustrates that hydrodynamically, this region is in the transition between deep sea and shallower coastal waters and is influenced by this feature. Furthermore, there are a few areas, such as in the south of the study area, with large shallow tidal water constituent coefficients where values are greater than 0.10. These values occur in areas where marine bathymetry was varying more abruptly.
3.4. Characteristics of Maximum Surface Currents Maxima
In order to provide useful information for coastal engineering operations and to further study the mechanisms of surface current movements, the maximum surface current velocities, both PMVCs and MMCVs, were computed using Equation (11); summary results are presented in
Figure 7.
Figure 7 shows range and distribution of PMCVs and MMCVs throughout the study area. The ranges of PMCVs and MMCVs are 7 cm/s to 13 cm/s and 60 cm/s to 130 cm/s, respectively; very large differences exist between the two velocity fields. Large values of PMCV prevail in the south and north of the region, while large values of MMCVs prevail in the center of the region. The spatial distribution of PMCVs does not agree well with that of MMCVs, primarily due to wind influences. To investigate effects of wind on MMCVs, statistics of spatially annually averaged ECMWF wind speeds, and maximum, minimum, and mean values, are presented in
Table 5.
Table 5 shows that the maximum annual wind speed occurs in winter, with a value of 21.96 m/s. Average wind speeds are the highest (9.03 m/s) and lowest (6.38 m/s) in winter and summer, respectively. At times of maximum wind speeds, more energy is transferred from to the top surface layer of the water column; surface current speeds are highly influenced by strong winds [
42]. Wind stresses cause the large differences between PMCVs and MMCVs. PMCVs are based on the computation of tidal constituents alone and only reflect tidally driven currents; however, in this region strong prevailing winds significantly increase speeds of surface currents. PMCVs do not correlate well with MMCVs in the study area; our analysis has shown that PMCVs should not be used for estimating maximum surface currents in the region; doing so will result in grossly under-estimated velocities.
The above approach presents a most useful methodology to determine significant characteristics of complex oceanographic features. This is the first time these techniques have been applied to a highly wind-dominated environment.
4. Discussion
As expected, tidal currents in the study area are shown to be dominated by semi-diurnal constituents, especially M2 and S2 tidal constituents. Shallow water tidal constituents are generated by abrupt changes in bathymetry. These constituents give rise to a complex hydrodynamics regime, including changes to the slopes of tidal ellipses.
Previously, studies on mechanics of current movements have been undertaken using spectral analysis. Li et al. [
43] studied the characteristics of tidal currents in the northern South China Sea using rotational spectral analysis based on ADCP data and found that the power spectral density is stronger at the diurnal frequency band than at the semidiurnal frequency band. However, the tidal currents in this study are dominated by semidiurnal tides. Li et al. [
43] studied the characteristics of tidal currents over water columns, and they found that the clockwise spectrum was stronger than the counterclockwise spectrum and that the counterclockwise spectrum dominated the spectral power density of the total rotary spectra. Rotary spectral analysis on HFR surface currents in our study area indicates that the counterclockwise spectrum dominated current movements. Yang et al. [
44] found that tidal ellipses around Tangshan International Tourism Island exhibited reciprocating characteristics, whereas tidal ellipses exhibited rotational characteristics in our study area. Additionally, Ren et al. [
42] studied residual current circulation and its response mechanism to wind at a seasonal scale for the same area; this study mainly focused on characterizing the patterns of tidal currents. Detailed analysis herein of mean surface flow fields, spatially averaged surface currents, and maximum surface currents provide a new perspective when assessing ocean currents.
Shen et al. [
33] found that the PMCV was in accordance with the MMCVs around the Zhoushan Islands; PMCVs and MMCVs exhibit very significant differences in our study area; wind is one of the main driving forces of surface currents at the time of MMCVs rather than tidal forces alone. Because only a one-year dataset was used for analysis in this research, effects of SA and SSA constituents were not considered in detail. Because the movements of currents at the time of the maximum current are mainly driven by wind stress rather than tidal force, PMCVs are different from MMCVs. Thus, in this area, it is important not to use the PMCV method when estimating maximum currents. In locations in which strong prevailing winds exist, it is very important that MMCVs are determined, and PMCVs should not be relied upon to provide estimates of maximum surface currents.
5. Conclusions
Surface currents were obtained from HFR systems that were used to study characteristics of tidal and non-tidal currents of the west coast of Ireland island in this research using a rigorous methodology. Characteristics and patterns of tidal currents were analyzed and discussed above. The main conclusions from our research are as follows:
(1) The direction of temporally averaged flow in the four seasons and over the whole year in the study area shows significant northly flow. The flow direction of currents in spring and autumn is primarily westerly; the directions of temporally averaged currents throughout the four seasons and over the entire year are strongly influenced by wind directions. The surface currents are affected not only by the wind stress but are also steered by bathymetry and topography.
(2) Surface flow fields in the study area are mainly affected by the tidal dynamics of the M2 constituent, showing significant regular semidiurnal tide characteristics. The tidal ellipses mainly exhibit rotating flow patterns.
(3) Shallow water tidal constituents are significantly influenced by changes in bathymetry. Values of the coefficient of shallow tidal water constituents are relatively large in areas with abrupt change of bathymetry.
(4) There is a significant difference between distribution of PMCVs and MMCVs. High values of PMCVs exist in the south of the study area, while high values of MMCVs exist in the middle of the study area. This is because the sea surface currents are mainly driven by the wind stress at the time of MMCVs, while PMCVs are mainly determined by tidal currents.
Forecasting of surface currents in the study area will be carried out in further studies using alternative methods, such as machine learning methods. Through the analysis of surface currents in the study area, the significant driven factors of the currents will be chosen as inputs to develop prediction models. New information and approaches have been presented in this paper that will be significant for the commercial exploitation of marine assets of the west coast of Ireland; additionally, tracking pollutant slicks, such as oil spills, will require highly accurate estimates of surface currents. More ocean current data at depth, which were not available during this study, should be measured and used in the analysis of the vertical structure of currents in this hydrodynamically complex region. Additionally, the relevance of wind stress with PMCVs or MMCVs and effects of wave-current interaction on PMCVs and MMCVs will be explored in future studies also.
Author Contributions
Conceptualization, L.R. and G.P.; methodology, G.P.; validation, L.R. and G.P.; formal analysis, G.P., Y.W. and P.Y.; investigation, L.Y. and P.Y.; resources, L.R. and M.H.; data curation, L.R. and M.H.; writing—original draft preparation, L.Y., Y.W., G.Z. and Q.Z. and Z.Z.; writing—review and editing, L.R., Z.Z. and M.H.; visualization, G.P., L.Y. and Y.W.; supervision, L.R.; project administration, L.R., G.Z. and Z.Z.; funding acquisition, L.R. and Z.Z. All authors have read and agreed to the published version of the manuscript.
Funding
This work was supported by the Belt and Road Special Foundation of the National Key Laboratory of Water Disaster Prevention (No. 2022491711); the State Key Laboratory of Satellite Ocean Environment Dynamics, Second Institute of Oceanography, MNR and the Zhejiang Provincial Natural Science Foundation of China under Grant (No. LR21D060002).
Data Availability Statement
The data presented in this study are available on request from the corresponding author.
Acknowledgments
We would like to thank ECMWF for providing wind data and the Irish Marine Institute for providing the tidal data.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Capodici, F.; Cosoli, S.; Ciraolo, G.; Nasello, C.; Maltese, A.; Poulain, P.-M.; Drago, A.; Azzopardi, J.; Gauci, A. Validation of HF radar sea surface currents in the Malta-Sicily Channel. Remote Sens. Environ. 2019, 225, 65–76. [Google Scholar] [CrossRef]
- Hays, G.C. Ocean currents and marine life. Curr. Biol. 2017, 27, R470–R473. [Google Scholar] [CrossRef] [PubMed]
- Mantovani, C.; Corgnati, L.; Horstmann, J.; Rubio, A.; Reyes, E.; Quentin, C.; Cosoli, S.; Asensio, J.L.; Mader, J.; Griffa, A. Best practices on high frequency radar deployment and operation for ocean current measurement. Front. Mar. Sci. 2020, 7, 210. [Google Scholar] [CrossRef]
- Tran, M.C.; Sentchev, A.; Nguyen, K.C. Multi-scale variability of circulation in the Gulf of Tonkin from remote sensing of surface currents by high-frequency radars. Ocean. Dyn. 2021, 71, 175–194. [Google Scholar] [CrossRef]
- Ren, L.; Chu, N.; Hu, Z.; Hartnett, M. Investigations into synoptic spatiotemporal characteristics of coastal upper ocean circulation using high frequency radar data and model output. Remote Sens. 2020, 12, 2841. [Google Scholar] [CrossRef]
- Liu, F.; Zhou, H.; Wen, B. DEDNet: Offshore eddy detection and location with HF radar by deep learning. Sensors 2020, 21, 126. [Google Scholar] [CrossRef] [PubMed]
- Mandal, S.; Sil, S.; Gangopadhyay, A. Tide-current-eddy interaction: A seasonal study using high frequency radar observations along the western Bay of Bengal near 16° N. Estuar. Coast. Shelf Sci. 2020, 232, 106523. [Google Scholar] [CrossRef]
- Chavanne, C.; Flament, P.; Carter, G.; Merrifield, M.; Luther, D.; Zaron, E.; Gurgel, K. The surface expression of semidiurnal internal tides near a strong source at Hawaii. Part I: Observations and numerical predictions. J. Phys. Oceanogr. 2010, 40, 1155–1179. [Google Scholar] [CrossRef]
- Gough, M.K.; Garfield, N.; McPhee-Shaw, E. An analysis of HF radar measured surface currents to determine tidal, wind-forced, and seasonal circulation in the Gulf of the Farallones, California, United States. J. Geophys. Res. Ocean. 2010, 115, C04019. [Google Scholar] [CrossRef]
- Mandal, S.; Sil, S.; Gangopadhyay, A.; Murty, T.; Swain, D. On extracting high-frequency tidal variability from HF radar data in the northwestern Bay of Bengal. J. Oper. Oceanogr. 2018, 11, 65–81. [Google Scholar] [CrossRef]
- Mandal, S.; Sil, S.; Pramanik, S.; Arunraj, K.S.; Jena, B.K. Characteristics and evolution of a coastal mesoscale eddy in the Western Bay of Bengal monitored by high-frequency radars. Dyn. Atmos. Ocean. 2019, 88, 101107. [Google Scholar] [CrossRef]
- Schaeffer, A.; Gramoulle, A.; Roughan, M.; Mantovanelli, A. Characterizing frontal eddies along the E ast A ustralian C urrent from HF radar observations. J. Geophys. Res. Ocean. 2017, 122, 3964–3980. [Google Scholar] [CrossRef]
- Lorente, P.; Piedracoba, S.; Montero, P.; Sotillo, M.G.; Ruiz, M.I.; Álvarez-Fanjul, E. Comparative Analysis of Summer Upwelling and Downwelling Events in NW Spain: A Model-Observations Approach. Remote Sens. 2020, 12, 2762. [Google Scholar] [CrossRef]
- Paduan, J.D.; Cook, M.S.; Tapia, V.M. Patterns of upwelling and relaxation around Monterey Bay based on long-term observations of surface currents from high frequency radar. Deep. Sea Res. Part II Top. Stud. Oceanogr. 2018, 151, 129–136. [Google Scholar] [CrossRef]
- Kohut, J.T.; Glenn, S.M.; Chant, R.J. Seasonal current variability on the New Jersey inner shelf. J. Geophys. Res. Ocean. 2004, 109, C07S07. [Google Scholar] [CrossRef]
- Mandal, S.; Sil, S.; Gangopadhyay, A.; Jena, B.K.; Venkatesan, R.; Gawarkiewicz, G. Seasonal and tidal variability of surface currents in the Western Andaman Sea using HF radars and buoy observations during 2016–2017. IEEE Trans. Geosci. Remote Sens. 2020, 59, 7235–7244. [Google Scholar] [CrossRef]
- Paduan, J.; Cook, M.; Fernandez, D. Two-dimensional diurnal to monthly period surface currents in Monterey Bay from CODAR-type HF radar. In Proceedings of the IGARSS’97. 1997 IEEE International Geoscience and Remote Sensing Symposium Proceedings. Remote Sensing—A Scientific Vision for Sustainable Development, Singapore, 3–8 August 1997; pp. 1814–1816. [Google Scholar]
- Gurgel, K.-W.; Antonischki, G.; Essen, H.-H.; Schlick, T. Wellen Radar (WERA): A new ground-wave HF radar for ocean remote sensing. Coast. Eng. 1999, 37, 219–234. [Google Scholar] [CrossRef]
- Hammond, T.; Pattiaratchi, C.; Eccles, D.; Osborne, M.; Nash, L.; Collins, M.B. Ocean surface current radar (OSCR) vector measurements on the inner continental shelf. Cont. Shelf Res. 1987, 7, 411–431. [Google Scholar] [CrossRef]
- Emery, B.M.; Washburn, L.; Harlan, J.A. Evaluating radial current measurements from CODAR high-frequency radars with moored current meters. J. Atmos. Ocean. Technol. 2004, 21, 1259–1271. [Google Scholar] [CrossRef]
- Shen, Z.; Wu, X.; Lin, H.; Chen, X.; Xu, X.a.; Li, L. Spatial distribution characteristics of surface tidal currents in the southwest of Taiwan Strait. J. Ocean. Univ. China 2014, 13, 971–978. [Google Scholar] [CrossRef]
- McMahon, T.; Raine, R.; Titov, O.; Boychuk, S. Some oceanographic features of northeastern Atlantic waters west of Ireland. ICES J. Mar. Sci. 1995, 52, 221–232. [Google Scholar] [CrossRef]
- Pingree, R.; Griffiths, D. Currents driven by a steady uniform wind stress on the shelf seas around the British-Isles. Oceanol. Acta 1980, 3, 227–236. [Google Scholar]
- Pingree, R.; Le Cann, B. Three anticyclonic Slope Water Oceanic eDDIES (SWODDIES) in the southern Bay of Biscay in 1990. Deep Sea Res. Part A Oceanogr. Res. Pap. 1992, 39, 1147–1175. [Google Scholar] [CrossRef]
- Polonsky, A.; Sukhonos, P. Variability of the wind stress, the field of currents, wind stress curl and vorticity of surface currents in the North Atlantic. Izv. Ross. Akad. Nauk. Seriya Geogr. 2017, 62–73. [Google Scholar] [CrossRef]
- Saviano, S.; Esposito, G.; Di Lemma, R.; de Ruggiero, P.; Zambianchi, E.; Pierini, S.; Falco, P.; Buonocore, B.; Cianelli, D.; Uttieri, M. Wind Direction Data from a Coastal HF Radar System in the Gulf of Naples (Central Mediterranean Sea). Remote Sens. 2021, 13, 1333. [Google Scholar] [CrossRef]
- Dobrynin, M.; Kleine, T.; Düsterhus, A.; Baehr, J. Skilful seasonal prediction of ocean surface waves in the Atlantic Ocean. Geophys. Res. Lett. 2019, 46, 1731–1739. [Google Scholar] [CrossRef]
- Gonzalez, P.L.; Brayshaw, D.J.; Zappa, G. The contribution of North Atlantic atmospheric circulation shifts to future wind speed projections for wind power over Europe. Clim. Dyn. 2019, 53, 4095–4113. [Google Scholar] [CrossRef]
- Barrick, D. First-order theory and analysis of MF/HF/VHF scatter from the sea. IEEE Trans. Antennas Propag. 1972, 20, 2–10. [Google Scholar] [CrossRef]
- O’Donncha, F.; Hartnett, M.; Nash, S.; Ren, L.; Ragnoli, E. Characterizing observed circulation patterns within a bay using HF radar and numerical model simulations. J. Mar. Syst. 2015, 142, 96–110. [Google Scholar] [CrossRef]
- Paduan, J.D.; Rosenfeld, L.K. Remotely sensed surface currents in Monterey Bay from shore-based HF radar (Coastal Ocean Dynamics Application Radar). J. Geophys. Res. Ocean. 1996, 101, 20669–20686. [Google Scholar] [CrossRef]
- Kundu, P.K.; Blanton, J.O.; Janopaul, M.M. Analysis of current observations on the Georgia shelf. J. Phys. Oceanogr. 1981, 11, 1139–1149. [Google Scholar] [CrossRef]
- Shen, Z.; Wu, X.; Fei, Y.; Xu, X.a.; Chen, X. Surface tidal currents in the open sea area to the east of the Zhoushan Islands measured with high frequency surface wave radar. Acta Oceanol. Sin. 2013, 32, 5–10. [Google Scholar] [CrossRef]
- Mooers, C.N.K. A technique for the cross spectrum analysis of pairs of complex-valued time series, with emphasis on properties of polarized components and rotational invariants. Deep Sea Res. Oceanogr. Abstr. 1973, 20, 1129–1141. [Google Scholar] [CrossRef]
- Gonella, J. A Rotary-component method for analysing meteorological and oceanographic vector time series. Deep Sea Res. Oceanogr. Abstr. 1972, 19, 833–846. [Google Scholar] [CrossRef]
- Perkin, H. Inertial oscillations in the Mediterranear. Deep Sea Res. Oceanogr. Abstr. 1972, 19, 289–296. [Google Scholar] [CrossRef]
- Calman, J. On the interpretation of ocean current spectra, part II: Testing dynamical hypotheses. J. Phys. Oceanogr. 1978, 8, 644–652. [Google Scholar] [CrossRef]
- Pawlowicz, R.; Beardsley, B.; Lentz, S. Classical tidal harmonic analysis including error estimates in MATLAB using T_TIDE. Comput. Geosci. 2002, 28, 929–937. [Google Scholar] [CrossRef]
- Halverson, M.; Pawlowicz, R. Tide, wind, and river forcing of the surface currents in the Fraser River plume. Atmos. Ocean. 2016, 54, 131–152. [Google Scholar] [CrossRef]
- Arabelos, D.; Asteriadis, G.; Contadakis, M.; Spatalas, S.; Sachsamanoglou, H. Atmospheric tides in the area of Thessaloniki. J. Geodyn. 1997, 23, 65–75. [Google Scholar] [CrossRef]
- Leaman, K.D.; Sanford, T.B. Vertical energy propagation of inertial waves: A vector spectral analysis of velocity profiles. J. Geophys. Res. Ocean. 1975, 80, 1975–1978. [Google Scholar] [CrossRef]
- Ren, L.; Yang, L.; Pan, G.; Zheng, G.; Zhu, Q.; Wang, Y.; Zhu, Z.; Hartnett, M. Characterizing Residual Current Circulation and Its Response Mechanism to Wind at a Seasonal Scale Based on High-Frequency Radar Data. Remote Sens. 2022, 14, 4510. [Google Scholar] [CrossRef]
- Li, R.; Chen, C.; Xia, H.; Beardsley, R.C.; Shi, M.; Lai, Z.; Lin, H.; Feng, Y.; Liu, C.; Xu, Q.; et al. Observed wintertime tidal and subtidal currents over the continental shelf in the northern South China Sea. J. Geophys. Res. Ocean. 2014, 119, 5289–5310. [Google Scholar] [CrossRef]
- Yang, J.; Ding, W.; Cui, J.; Guo, S. Characteristical analysis of tidal and residual currents in the sea area around Tangshan international tourism island. IOP Conf. Ser. Earth Environ. Sci. 2021, 632, 022044. [Google Scholar] [CrossRef]
| Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).