High-Frequency Radar Observations of Surface Circulation Features along the South-Western Australian Coast

: A new merged high-frequency radar (HFR) data set collected using SeaSonde and WERA (WEllen RAdar) systems was used to examine the ocean surface circulation at diurnal, seasonal and inter-annual time scales along the south-west coast of Australia (SWWA), between 29° – 32° S. Merging was performed after resampling WERA data on the coarser SeaSonde HFR grid and averaging data from the two HFR systems in the area of common overlap. Direct comparisons between WERA and SeaSonde vectors in their overlapping areas provided scalar and vector correlation values in the range Ru = [0.24, 0.76]; Rv = [0.39, 0.83]; ρ = [0.44, 0.75], with mean bias between velocity components in the range [ − 0.02, 0.28] ms −1 , [ − 0.16, 0.16] ms −1 for the U, V components, respectively. The lower agreement between vectors was obtained in general at the boundaries of the HFR domains, where the combined effects of the bearing errors, geometrical constraints, and the limited angular field of view were predominant. The combined data set allowed for a novel characterization of the dominant features in the region, such as the warmer poleward-flowing Leeuwin Current (LC), the colder Capes Current (CC) and its northward extensions, the presence of sub-mesoscale to mesoscale eddies and their generation and aggregation areas, along with the extent offshore of the inertial-diurnal signal. The contribution of tides was weak within the entire HFR domain (<10% total variance), whilst signatures of significant inertial- and diurnal-period currents were present due to diurnal – inertial resonance. A clear discontinuity in energy and variance distribution occurred at the shelf break, which separates the continental shelf and deeper offshore regions, and defined the core of the LC. Confined between the LC and the coastline, the narrower and colder CC current was a feature during the summer months. Persistent (lifespan greater than 1 day) sub-mesoscale eddies (Rossby number O (1)) were observed at two main regions, north and south of 31.5° S, offshore of the 200 m depth contour. The majority of these eddies had diameters in the range 10 – 20 km with 50% more counter clockwise rotating (CCW) eddies compared to clockwise (CW) rotating eddies. The northern region was dominated by CCW eddies that were present throughout the year whilst CW eddies were prevalent in the south with lower numbers during the summer months.


Introduction
Surface circulation along the eastern margin the southern Indian Ocean is very different compared to other ocean basins, with the poleward Leeuwin Current (LC) transporting warmer water ( Figure 1) and promoting a downwelling system [1]. The eastern margins of all other ocean basins comprise equatorward cold water currents associated with upwelling. The anomalous ocean currents together with diurnal tides with low range [2], strong sea breezes [3], resonance at the critical latitude [4], high air-sea heat and freshwater fluxes [5] provide a unique environment. In addition, the LC is associated with elevated levels of mesoscale eddy activity compared with other eastern boundary current regions [6,7]. The region is dominated by two major currents ( Figure 1): the warmer Leeuwin Current (LC), flowing southwards along the shelf break, and the colder Capes Current (CC), flowing northwards along the continental shelf inshore the LC [1,8]. The LC is a poleward eastern boundary current driven by an alongshore pressure gradient, generated through the transport of warmer, lower salinity water into the tropical Indian Ocean by the Indonesian Throughflow [9][10][11]. It is augmented with contributions from the central South Indian Counter Current (SICC) and the Holloway Current in the shelf region around 20°-23° S and from southern SICC at around 32° S. LC flows southwards as a relatively narrow (~100 km wide), shallow (upper 250 m) current along the continental shelf break (~200 m) to the southern tip of Tasmania [11,12]. The LC has a seasonal pattern with intensification during Austral autumn/winter (April-September) and a weakening during the Austral summer months (November-March). This seasonal variability is controlled by seasonal changes in the wind stress: (1) during the Australasian monsoon season (January to March) to the north of Australia [13]; and, (2) equatorward wind stress within the study region [3,9]. The LC is also influenced by El Niño and La Niña cycles with the current being stronger (weaker) during La Niña (El Niño) events due to fluctuations in the ITF and equatorial monsoonal winds [10,11,14]. Typical features of the LC include meanders, eddies, and filaments [15] and are highlighted in Figure 1. One of the major features of the LC is the presence of a highly energetic meso-scale eddy field with diameters >200 km and lifetimes of one month or more [15][16][17][18]. Anti-clockwise (warm-core) and clockwise (cold-core) rings are generated by the LC at preferential locations ( Figure 1) that are associated with changes in the coastline orientation, and they migrate westward into the interior of the southern Indian Ocean [6]. The formation of these eddies has been postulated to be from different mechanisms such as instabilities through meanders, shear between the LC and Leeuwin Undercurrent/Capes Current and bathymetric features [16][17][18][19] and are located at ~28° S and ~30.5° S (Figure 1). These eddies are important circulation features, as they influence the heat and momentum balance [7]; regional upwelling and therefore contribute to local biological activity [16,19,20] and, influence the coastal wave climate [21].
The CC is a wind driven current that transport upwelled cooler, salty waters northwards inshore of the LC (Figure 1b) originating at ~34° S and extending as far north as 29° S [8]. It is mainly present during the summer months (November-March) when the upwelling favorable equatorward winds are the strongest. Here, due to Ekman dynamics, the upper layers move offshore in response to the northward wind stress pushing the LC offshore [7,8]. The CC is usually visible as a band of colder water on the continental shelf ( Figure 1b).
Wind action is the primary driving mechanism for surface current variability on the continental shelf. The study region is micro-tidal, dominated by diurnal tides and with a mean spring tidal range of ±0.6 m [2]. The inertial period is around 24 h. The strong land-sea breeze (LSB) regime that characterizes this region [3,4,22] results in strong counterclockwise diurnal currents in the upperocean dynamics through the diurnal-inertial resonance at critical latitudes, that can penetrate below the mixed-layer depth [4]. Extra-tropical storms or tropical cyclones may also contribute to the intensification of upper layer current variability in the region [23].
The study region consists of several distinct bathymetric features (  High-frequency radar (HFR) measurement of surface currents in the coastal ocean has become a standard and cost-effective component for ocean observing systems globally [24]. They were developed more than three decades ago [25,26] and are now widely used with a high level of accuracy [27,28] for a variety of purposes, including ocean current forecasting [29] and long-term trends in surface circulation in coastal areas [30]. HFR systems allow for the synoptic observation of large-scale coastal circulation features with temporal and spatial resolution not readily available with conventional instrumentation a high level of accuracy [27]. As part of the Australian Integrated Marine Observing System (IMOS), HFR systems have been deployed along the south-west Australian coast since 2009.
Previous studies in the region have mainly used satellite imagery [6,8,15], shipborne measurements [1,7,31] and, numerical modelling [17,19]. Mihanovic et al. [4] used oceanographic moorings, satellite sea surface temperature (SST) and HFR data in the southern section of the study region to examine the diurnal-inertial resonance at the critical latitude during the summer months. Studies that have examined the diurnal, seasonal and eddy dynamics in the study region have been limited either by low spatial and temporal resolution (e.g., satellite data), or restricted by small spatial domains (e.g., moorings or single HFR systems).
We build on previous work by merging data from two HFR systems that doubles the spatial coverage whilst retaining the high temporal and spatial information. The aim of this paper is to describe the surface oceanographic features along the coast of south-west Western Australia (SWWA) using the new merged near-surface HFR derived currents. The emphasis is on characterising the: (1) surface currents at diurnal (including inertial) lower frequency (>4 days) and seasonal time scales; and, (2) dynamics of the eddy field. The merged surface current data were collected over a 33-month period (January 2015 to September 2017) and covered a region 180 km offshore and 250 km alongshore at a resolution of 6 km. This paper is arranged as follows: Section 2 describes the data collection and analysis techniques. Results and discussion are presented in Sections 3 and 4. The conclusions are presented in Section 5.

Materials and Methods
Surface current data (top 2 m of the water column) collected by long-range HFR systems ( Figure  2) over the 33-month period January 2015 to September 2017 were used in this study. HFR current data were sourced through commercial direction-finding (SeaSonde) and phased-array (WERA) HFR systems, provided by Codar Ocean Sensors (COS) and Helzel MessTechnik, respectively.
Two SeaSonde systems, located at Lancelin (LANC; 31°01.588′ S, 115°19.708′ E) and Green Head (GHED; 30°04.393′ S, 114°58.003′ E), were used to collect ocean currents at hourly intervals with ranges up to 200 km offshore with a spatial resolution of 6 km. Both stations operated at a central frequency of 4.463 MHz (measurement depth is approximately 2.5 m) with a 25 kHz bandwidth, 512point fft, 1 Hz sweep rate, 20 min Doppler spectra averaging with 15-min output rate. Calibrated antenna patterns are used at both stations. Two WERA systems, located at Guilderton (GUI; 31°20.4960′ S, 115°29.3640′ E) and Port Beach, Fremantle (FRE; 32°.0198′ S, 115°44.7480′ E) provided hourly ocean surface currents at ranges up to 180 km offshore at 4 km spatial resolution. Both sites used 4-element transmit and receive arrays of 16 equally spaced elements. Systems operated at 9.335 MHz (sampling depth approximately 1.3 m) with 33.3 kHz bandwidth, 1024-point fft and 0.26 s sweep rate, and an integration time of 5 min. Radial velocity data were collected on a Cartesian grid with 4 km spacing every 10 min in alternate mode. The same Cartesian grid was used to compute hourly surface current vectors from the hourly-averaged radial data.
Sea-surface temperature (SST) data originate from quality-controlled, merged, day and night, multi-sensor L3S foundation SST product from IMOS level L3S gridded single-sensor composites (IMOS-SRS SATELLITE-SST L3S-01 day composite). SST data for the WA coast were downloaded from the AODN (Australian Ocean Data Network) data portal (https://portal.aodn.org.au). HRPT AVHRR (High Resolution Picture Transmission/Advanced Very High Resolution Radiometer) SST (sea surface temperature) retrievals are produced by the Australian Bureau of Meteorology as a contribution to the Integrated Marine Observing System [32].
The analysis approach used hourly, daily, monthly and seasonally-averaged delayed-mode quality controlled (DMqc) data [33], and the most appropriate statistical tools were used to characterize the dominant circulation features in the study region.
Smoothed estimates of the power density distribution (PSD) of surface HFR currents were derived from hourly HFR data in order to partition current variance over resolvable frequencies at each HFR grid cell. Spectra of hourly HFR observations were derived following the Welch-modified periodogram method [34] using 1024-point segments (approximately 42-days long segments), a Hanning window and 50% overlap between the data segments. Band-averaged energy levels were then extracted by integrating the resulting PSD curves in the selected frequency bands: diurnalinertial frequency [0.038, 0.046] cph (cycles per hour), corresponding to periods in the range [ (1); [35]), were then derived from the smooth PSD estimates of the clockwise (CW)-counterclockwise (CCW) spectra: where ± ( ) are, the CCW and CW spectral estimates at each HFR grid cell at the frequency respectively. CR quantifies the degree of polarization of current observations at each frequency band; for purely rectilinear flow, CR equals to zero, and CR = 1 for purely circular (counterclockwiseclockwise) motion.
The MATLAB t tide package [36] was used to quantify the contribution of tidal forcing to the current variance. Synthetic time series at each grid cell included in the analysis are derived using only constituents with signal-to-noise ratios greater than 2.
For spectral analysis, gaps in the time series were linearly interpolated before the fft calculations, whilst no interpolation was performed before the tidal analyses. Grid points with at least 70% data return were used for both analyses, in order to increase the frequency resolution and reduce biases from missing data.
The merged SeaSonde-WERA data set was obtained as follows. DMqc procedures described in [33,37] were applied separately to SeaSonde and WERA current vector data at each HFR grid point. DMqc flags were applied to the hourly data sets to exclude anomalous currents from the records. We then computed daily-averaged currents fields separately for the SeaSonde and WERA HFR data sets using grid points with at least 50% data return (after removing spikes). Daily SeaSonde and WERA currents were remapped on a common longitude-latitude grid by searching velocity data within 2 km from each grid cell, and averaging separately the U, V components of the velocity vectors. In case of overlap between SeaSonde and WERA currents, U, V from both HFR genres were averaged.
Validation of SeaSonde and WERA HFR data in the region against current meter data for the same data period were performed previously [33,38], showing that both HFR types captured the variability sampled by the current meters consistent with accuracy levels reported in the literature. The interoperability of the two HFR genres in the region was also demonstrated for a shorter time period (1 month), showing that radial data from SeaSonde can be merged successfully with radial data from WERA systems with no significant performance losses [38]. Extending the comparisons to the overlapping areas between the two grids, scalar (vector [39]) correlation (Equation (2); Equation In Equations (4), (u1,2, v1,2) stand for the zonal and meridional components of the SeaSonde (WERA) HFR velocity components respectively, i identifies the imaginary unit, and < > represents the ensemble average operator.
Values for U, V biases (Equation (4)) were in the range biasU = [−0.02, 0.28] ms −1 , biasV = [−0.16, 0.16] ms −1 . Statistics reported here are relative to a sample of 67 grid points, with a lowest population (number of common data) in excess of 500. It must be pointed out that lower agreement between vectors was obtained in general at the boundaries of the SeaSonde domain, where the combined effects of the bearing errors [33,38], the geometrical constraints common to both systems, and the limited angular field of view of the WERA systems were predominant.
Ocean motions can be divided into four classes based on horizontal length scales [40]: (1) large scale motions that are > 300 km; (2) meso-scale motions are related to eddies with diameters 50-300 km; (3) sub-mesoscale motions have horizontal scales 1-50 km; and, (4) small-scale motions have scales <1 km. With HFR coverage, extending 100-170 km offshore with a resolution of 6 km, only submesoscale motions could be resolved. Statistics of eddies (lifespan, dimensions, locations, rotation sense, vorticity, eddy Rossby number) were extracted from the daily-averaged, merged SeaSonde-WERA data set. The vector geometry approach described in [41] was used to identify eddies and their properties in the merged data set. Sub-mesoscale eddies that have time scales of O(days) are commonly observed in the region, particularly in the shear region between the LC and the CC currents, and only persistent eddies (minimum lifespan 1 day) could be detected in the data set used here.
Relative vorticity ( ; Equation (6)), divergence ( ; Equation (7)) and Rossby number (Ro; Equation (8)) were computed from the merged, daily-averaged HFR currents at each grid point by a locally least-squares fitting a velocity plane to each grid point using current measurements from nearby locations [42]. Planetary vorticity (f) is defined in Equation (9), with Ω = 7.2921 × 10 −5 radians/s the rotation rate of Earth and the latitude.

The Mean State
The mean HFR currents over the study period (January 2015 to September 2017) indicated southward currents, offshore the 200 m contour over the whole HFR coverage region associated with warmer SSTs (Figure 2). This is the Leeuwin Current. In the northern region of HFR coverage, the range of the SeaSonde system extended further offshore compared to the southern region where the range of WERA was lower. Inshore of the 100 m depth contour, the currents were variable but the majority contained a northward component and were associated with colder water, particularly in the southern region of HFR coverage ( Figure 2). This is the signature of the wind-driven Capes Current, most prominent during the summer months. There was also a northward component of wind stress, during the winter months that resulted in northward currents.
The mean cross-shelf profile of HFR currents and SST both indicated a Gaussian shape with weaker currents inshore, increasing to a maximum (the LC 'core') and then decreasing further offshore ( Figure 3). The SST profiles indicated colder water inshore, increasing offshore to a maximum and then decreasing further offshore (Figure 3b). The peak in the southward velocity maxima was located further offshore compared to the SST maxima. For example, at 31° S, the peak in velocity (−0.17 ms −1 ) was located at 114.5° E whilst for SST it was at 114.8° E (Figure 3). The LC flows parallel the depth contours ( Figure 2). However, the depth contours (and the shoreline) is oriented 337° with respect to north. Therefore, the 200 m depth contour is shifted to the east as we move south. This is reflected in the cross-shelf profiles of both HFR currents and SST with the maxima migrating eastward (Figure 3). Profiles of HFR currents also indicated northward mean currents at the two southern latitudes (0.07 ms −1 at 32° S and 0.025 ms −1 at 31.7° S) whist at 31° S, although the currents so not indicate a positive component mean currents were close to zero (−0.02 ms −1 ). The maximum velocity of the LC 'core' increased southward, from -0. 17

Seasonal and Inter-Annual Variability
Seasonally averaged HFR and SST data indicated the typical mean circulation structures in the region with the relatively warmer LC as a southward current and the CC as a weaker, colder coastal current (Figures 4-6). The presence of colder water adjacent to the coast during the autumn and winter months (Figures 4c,5c,6c) is due to heat loss [5]. In general, the LC was stronger (mean speeds up to 0.50 ms −1 ), during autumn/winter and weaker during spring/summer.
During the summer months, a band of colder water associated with northward currents is the CC (Figures 4a,5a,6a). Note that the water adjacent the shoreline is warmer due to summer solar heating and the band of colder water, associated with the upwelling CC, is located further offshore but inshore of the warmer LC. The boundary between the colder and warmer water corresponded with a change in current direction from northward to southward with strong horizontal velocity gradients. This was also highlighted by the cross-shelf velocity profiles at different latitudes where northward (positive) current component were recorded inshore (Figure 7a). At 31.0° S, the LC was present as a broader current with the strongest currents west of 114.5° E whilst at 31.7° S and 32.0° S the cross-shelf profiles indicated positive and negative maxima associated with the LC and CC (Figure 7a).    During the winter months, the LC was the most dominant feature mostly associated with the warmer water. Although the HFR coverage captured the 'core' of the current, the larger meso-scale eddies were not captured fully. For example, in 2015, there was a large mesoscale clockwise eddy centered around 31.5° S, 114.0° E during autumn and winter that was only captured partially ( Figure  4b,c) although altimetry data indicated the presence (and the persistence) of this eddy through autumn to winter (http://oceancurrent.imos.org.au/). Similar partial capture of a clockwise mesoscale eddy was also recorded during winter months of 2016 and 2017 (Figures 5c,6c). Cross-shelf velocity profiles at different latitudes indicated strengthening of the southward (negative) current component, compared to those during the summer at all three latitudes (Figure 7). In particular, at 31.0° S, there was a pronounced 'core' in the LC with a maximum at 114.6° E, a significant onshore migration during winter compared to the summer (see Figure 7a,b).
Autumn months indicated the strengthening of the LC between summer and winter (Figures 4b,  5b, 6b). The currents during 2015 autumn were stronger when compared to 2016 and 2017. The spring months represented weakening of the LC from winter to summer (Figures 4d, 5d, 6d)  It is well documented that LC is weaker during El Niño events and stronger during La Niña events [11,14]. HFR data, presented here, were collected during the strong El Niño event of 2015-2016 and thus represented weaker LC conditions. Interestingly, there were differences in the cross-shelf current profiles where at latitude 31° S the annual mean currents recorded in 2017 were almost double the magnitude of 2015 with maximum intensity further offshore (Figure 8a). In contrast, at latitudes 31.7° S and 32.0° S there was not a large difference between the years with the maximum currents recorded in 2016 (Figure 8b,c). This could be due to the inflow from the southern South Indian Counter Current into the LC that occurs around 32° S, as well as other factors such as eddy activity or other larger scale features outside the HFR domain.

Frequency Domain Analyses
Within the HFR domain, rotary spectra of hourly current vectors indicated common features that included energetic low-frequency components, well-defined peaks in the inertial band, (the inertial period ranged between 21 h 45 m and 24 h 17 m within the HFR domain), diurnal frequency band (K1 and S1 tidal constituents or wind), and limited energy within the semidiurnal tidal frequency band (M2, S2). These characteristics are highlighted in the spectra (Figure 9) at the shelf break in the Turquoise coast (location 'A' in Figure 2) and the Rottnest shelf (location 'B' in Figure 2). The points A and B were selected based on their data return (in excess of 90% over the analysis period at both locations), to minimize the effects of gaps and data interpolation.
Both systems yielded similar spectral energy level at lower frequencies. However, there were differences in both the high and low frequency tails of the spectra ( Figure 9): (1) at higher frequencies, SeaSonde currents contained significantly higher noise levels than that for WERA data ( Figure 9a); (2) at the low-frequency bands, both SeaSonde and WERA spectra showed a predominant counterclockwise (CCW) component over the clockwise (CW) spectrum (Figure 9b) although the difference between the CCW and CW peaks were larger for SeaSonde data; and, (3) at the diurnal frequency band, the CCW spectral component was more energetic in the WERA currents and exhibited a stronger polarization (CR = 0.98 at a 24 h period) than SeaSonde data (CR = 0.86 at a 24 h period) over a broader frequency range (Figure 9c). Different variance levels within the CW and CCW spectra and the corresponding rotary coefficients (CR) at point locations (A and B, Figure 2) implied dominance of quasi-circular motions within the diurnal-inertial frequency band (Figure 9b). Similar features were observed along two zonal (East-West), cross-shore transects centered at 30°48′ S ( Figure 10a) and 32° S (Figure 10b), respectively. A different behavior between the two transects can be identified in the semidiurnal tidal frequency band, with CW variance increasing (decreasing) in the offshore direction at 30°48′ S (32° S; Figure 10a,b). Polarization at this frequency band was negative (CW) for both transects but changed to positive (CCW) at higher range (>60 km offshore). Variance levels of the currents increased both southwards and offshore, particularly within the low-and diurnal-frequency bands, with a distinct discontinuity around 30-40 km from the coast, approximately at the location of the 200 m depth contour (e.g., shelf break, Figure 2), indicating spatial distribution of near-inertial and diurnal variance and the low-frequency band ( Figure 11). Here, the near inertial to diurnal band spans [0.038, 0.046] cph, corresponded to periods in the range [21 h 44 m 20.8 s, 26 h 18 m 56.8 s]. The study region experiences strong diurnal-inertial resonance due to strong diurnal sea breezes and the critical latitude that generate near-inertial waves that extend to depths >200 m [4]. Higher variance in deeper water (Figure 11a), therefore, reflect that the diurnalinertial resonance being more effective in deeper water (see also [4]]). The low-frequency band was defined as frequencies below 0.0088 cph (periods longer than 4 days) and the distribution of variance indicated a band of higher values between the 200 m and 3000 m depth contours and is representative of the LC variability (Figure 11b). Tidal currents were generally negligible, consistent with microtidal conditions in this region. Semidiurnal tides and diurnal tides represented 4%-10% of the total variance, with typical amplitudes that barely exceeded 0.02 ms −1 for semidiurnal constituents within the HFR coverage. Amplitudes of diurnal constituents were generally larger but they most likely reflected contamination from the strong wind-driven currents in this frequency band rather than being 'true' tidal oscillations.

Eddy Statistics
The LC is associated with an energetic meso-scale eddy field with maximum eddy diameters > 200 km. However, the HFR coverage was limited to examining the sub-mesoscale (eddies with diameter < 50 km). The automated eddy identification method revealed ~300 eddies with at least 1day persistence over the period 2015-2017. There was a prevalence of counterclockwise-rotating (CCW, 177 cases) over clockwise-rotating eddies (CW, 120 cases). Eddies appeared to persist at the same location for extended periods, with typical drift distance of ~10 km around their center; however, in some cases they were advected with the mean flow and travelled in a southwards direction with speeds up to 8 km/day. Their occurrence was limited to deeper water; offshore 200 m depth contour (Figure 12a). There were two main regions for the aggregation of eddies (Figure 12a): (1) between latitudes 30° S and 31.5° S where larger (50-60 km diameter; Figure 12b), CCW eddies were prevalent. Here, the eddies were spread over a wider region; and, (2) offshore the Rottnest island, associated with the Perth canyon (Figure 2), where both isolated or paired CW and CCW eddies were observed. The Perth canyon region was where the most persistent eddies were found concentrated in a smaller region. This is most likely due to the topographic effect of the canyon and the horizontal shear associated with the LC and the CC.
Overall, there were 50% more CCW eddies compared to CW with eddies with diameter 15-25 km being the most prevalent (Figure 12b). Maximum eddy diameters were up to 60 km. The Rossby numbers (Ro) were up to 1.5 (mean value Ro = 0.46) for CW eddies (Figures 12c), and to Ro = 0.64 for CCW eddies (mean value: Ro = 0.32). Figure 12. (a) Shapes and locations of persistent (lifespan longer than 1 day) eddies along the study region; (b) probability density distribution (pdf) curves for eddy diameters (km); (c) probability density distribution (pdf) curves for Rossby number (Ro; Equation (8)). Blue and red curves represent counterclockwise (CCW) and clockwise (CW) eddies.
Time series of monthly number of days with eddies over the period January 2015 to September 2017 highlighted the seasonal and inter-annual variability ( Figure 13). It should be noted that the yaxis in Figure 13 does not represent the number of unique eddies but rather number of days in each month ('eddy-days') where an eddy was present. For example, October 2015 indicates >20 days of CW eddy activity in the southern region. This was due to a single eddy that persisted for most of this month that also dominated the seasonal currents (see Figure 4d). However, even neglecting October 2015 it is obvious that 2015 had more eddies compared to 2016 and 2017. In the northern region (north of 31.5° S) there were more days with CCW eddies compared to CW eddies (Figure 13b). There was also only a small month-to-month variability in the northern region (except for the second half of 2016 where there only a few 'eddy-days'). In contrast in the southern region, there was a quasiseasonal signal with fewer 'eddy-days' during the summer (Figure 13c).

Discussions
Ocean current measurement using HFR is well established globally over more than three decades and is one of the few techniques that provide near-real time synoptic surface currents in the coastal ocean [24,27]. Currently, two HFR systems are available commercially: SeaSonde by Codar Ocean Sensors (USA) and the WERA by Helzel MessTechnik (Germany). Although they both use HFR signals and Bragg scattering as the fundamental basis for current measurements, they use contrasting methods to estimate the surface currents. The SeaSonde system uses a direction-finding technique whilst the WERA uses a phased-array system. Although these systems are deployed at many locations [24] globally, the West Australian coast is the only region in which the two different systems are co-located and provides data along a 270 km of the coastline (Figure 2). The systems operate at different frequencies and resolutions: the SeaSonde at 4.463 MHz and 6 km resolution with WERA at 9.335 MHz and 4 km. It was previously shown that data from the two major HFR genres do map accurately the observed variability of the ocean currents [33,38], and that data from the two platforms can be successfully merged into a combined product. Here, we show that the combined data sets, especially at the boundaries of the two domains, have the potential to greatly increase the capabilities of coastal observing networks, in general, and along the southwestern Western Australia (SWWA) coastal region, in particular. Previous research used WERA HFR data to focus on the smaller domain in the Rottnest Shelf [4,5], and neglected the area to the North, which has proven to be important for the mesoscale eddy field. The merged SeaSonde-WERA surface current data, mapped onto a 6 km grid, provided robust and consistent current measurements with negligible differences at the diurnal, seasonal and inter-annual time scales. Here we used a three-year record to document unreported features in the SWWA region, and plan to extend the data set to the time period 2010-2019. This data set can be used for a variety of purposes, such as calibration and validation of ocean circulation models, characterization of long-term trends, and in support of fishery management.
There was seamless transition between the two systems as evidenced by the mean (Figure 2), seasonal currents (Figures 4-6) and the variance ( Figure 11). However, there were differences between the two systems due to their local configuration: the maximum horizontal ranges were 170 km and 105 km for the SeaSonde and WERA systems, respectively. This was due to the operating frequencies and grid resolution (see above). There were also differences in both the high-and lowfrequency tails of the spectra from each system with SeaSonde currents having higher noise levels at higher frequencies whilst at the diurnal frequency band, the CCW spectral component was more energetic in the WERA currents ( Figure 9).
The merged SeaSonde-WERA surface current data over the study period (January 2015 to September 2017) reflected the two major current systems in the region: Leeuwin (LC) and Capes (CC) Currents with the latter present mainly during the summer months. The southward flowing LC was associated with warmer water and was located offshore of the 200 m contour (Figure 2). Inshore of the 100 m depth contour, the currents were more variable. The northward flowing CC was present in the summer months, inshore of the LC, and was associated with colder water. The LC and CC systems are well documented along this region of the coast and the HFR data accurately reflected their presence and variability [1,7,8,11,12]. The 'core' of the LC was generally located offshore of the 200 m contour with the maximum velocity increasing from 31° S to 32° S most likely due to inflow from the South Indian Counter Current (SICC). This is consistent with the long-term numerical simulations [3] that indicated inflow from the southern arm of the SICC into the LC at around 32° S (see also [43,44]). The seasonal variability of the LC was well captured by the HFR data with higher currents (mean speeds up to 0.50 ms −1 ), during the winter months whilst weaker currents were present during the summer consistent with previous studies. In general, the LC was stronger during autumn/winter and weaker during spring/summer. The CC was mainly a summer feature with strong horizontal gradients at the interface between the southward flowing LC and the northwards flowing CC (see also [8]).
The study period coincided with the 2015-2016 strong El Niño event and represented conditions when the LC was weaker [11,14]. Nevertheless, there were inter-annual differences with the maximum currents recorded in 2017 almost double that recorded in 2015 (Figure 8a) as the system was adjusting to neutral ENSO conditions in 2017.
The study region has many processes that contribute to diurnal variability that include tides, inertial oscillations and sea breezes. Although the K1 and O1 diurnal tidal constituents dominate the tidal signal [45], analysis of the HFR data indicated that the tidal currents were generally negligible representing <10% of the total variance. The study region experiences one of the strongest and persistent diurnal coastal land sea breeze (LSB) systems globally [3]. The inertial period over the HFR coverage region varies between 21 h 45 m and 24 h 17 m. When the period of local wind forcing is close to the local inertial period as is the case in the study region, a resonance condition occurs. These resonant wind-current responses have been addressed in a variety of field and theoretical studies and have been shown to enhance the upper ocean velocity field and vertical mixing [4,46]. Using HFR and oceanographic mooring data, [4] documented the presence of near-inertial waves in the study region that were generated through diurnal-inertial resonance. During the austral summer, when southerly winds and the LSB system dominated the wind regime, strong counterclockwise diurnal motions with surface amplitudes >0.3 ms -1 , penetrated to 300 m depth with diurnal vertical isotherm fluctuations up to 60 m with the strongest response in water depths >200 m [4]. The HFR data indicated strong CCW oscillations throughout the study region ( Figure 10) with higher variance at the diurnal band in the deeper water (>200 m) consistent with the results presented by [4].
Eastern margins of ocean basins consist of large scale upwelling and equatorward flow except off the West Australian coast where the Leeuwin Current (LC) transports warmer water poleward. The LC is present due to the Indonesian Throughflow that transports water between the Pacific and Indian oceans through the Indonesian Archipelago. The LC is associated with high eddy kinetic energy that is much higher than that of other eastern boundary currents [7]. The highest eddy kinetic energy (>600 cm 2 s −2 ) occurs offshore of the main current flow (~200 m depth contour) between 28.75 and 32° S (Figure 12a). The LC supports a highly energetic eddy field from the sub-mesoscale to mesoscale with eddy diameters ranging from <20 km to >200 m [18]. Fang and Morrow [6], using altimeter data, found that majority of the mesoscale eddies were generated between 28° and 31° S and propagated into the interior of the Indian Ocean. Over a 6 year period, [6] found that 3-9 large eddies were spawned annually. These mesoscale eddies have diameters >100 km and were not wholly identified in the HFR data. However, there was evidence for the presence of these eddies, during the winter months, where part of the eddy was present in the HFR data (Figures 4c, 5c, 6c). HFR data, together with the automated eddy detection algorithm, were able to identify CW and CCW submesoscale eddies with diameters up to 60 km with the majority in the range 10-20 km (Figure 12b). There were 50% more CCW eddies compared to CW eddies. The 177 eddies identified by the automated eddy detection algorithm had Rossby numbers (Ro) of order 1 confirming that these eddy features were part of the sub-mesoscale motions. There were two main regions of eddy presence, separated to the north and south of 31.5° S. The northern region was dominated by CCW eddies and occurred throughout the year. In contrast, the southern region was dominated by CW eddies and was subject to seasonal variability with lower numbers during the summer months ( Figure 13). The differences in the two regions may be explained by the flow patterns associated with the LC and interaction with topography. The LC follows that 200 m depth contour which is similar to the changes in the shoreline orientation. At locations where the shoreline orientation changes, the LC becomes unstable and mesoscale eddies are formed [15]. This is highlighted in Figure 1a where two large CCW eddies are seen to be spinning off the LC at 27.5° S (eddy center at 27.5° S) and 30.5° S (eddy center at 31.5° S) respectively. The circulation around the southern eddy centered at 30.5° S 31.5° S follows an 'S' shaped pattern and are reflected in the mean seasonal HFR currents (Figures 5c and 6c,d) and shown schematically in Figure 14. Here, the LC pathway separates the CCW and CW eddies to the north and south, respectively. This is the reason for the prevalence of CCW (CW) in the north (south).
In the south, the concave shape of the coastline associated with the Perth canyon ( Figure 2) promotes a CW eddy as reflected in HFR currents (e.g., Figures 4d, 5c and 6c,d). However, the LC pathway changes often with the whole system shown in Figure 14 shifting to the north and thus changing the orientation of eddies at different locations. Figure 14. Schematic of the typical Leeuwin Current flow path (black line) and the generation of CCW (blue) and CW (red) vorticity and thus sub-mesoscale eddies.

Conclusions
A multi-year merged HFR data collected using SeaSonde and WERA systems, along the southwest Australia (SWWA) between 29°-32° S were used to examine the ocean surface circulation at diurnal, seasonal and inter-annual time scales. The results indicated that: • SeaSonde and WERA HFR data can be merged successfully with the potential to increase the extent of the coastal ocean under monitoring.

•
The warmer poleward-flowing Leeuwin Current (LC) was the dominant feature of the circulation with a strong seasonal signal: stronger during winter and weaker during summer. The colder Capes Current (CC), located inshore of the LC, was mainly present during the summer months.

•
We observed a zonal migration of the core of the Leeuwin Current and differences in magnitude between years, however more investigation is needed to determine the driving mechanisms and long term variability.

•
The contribution of tides to the circulation was small (<10% total variance). The energy contained in the diurnal period currents was dominant due to diurnal-inertial resonance generated strong sea breezes and inertial period at the critical latitude. • A clear discontinuity in energy and variance distribution occurred at the shelf break that separated the continental shelf and deeper offshore regions for both diurnal and low-frequency bands reflecting reflected the diurnal-inertial resonance and the LC, respectively. • Persistent (lifespan greater than 1 day) sub-mesoscale eddies (Rossby number O(1)) were observed at two main regions, north and south of 31.5° S, offshore of the 200m depth contour. Majority of these eddies had diameters in the range 10-20 km with 50% more counter clockwise rotating (CCW) eddies compared to clockwise (CW) rotating eddies. The northern region was dominated by CCW eddies that were present throughout the year whilst CW eddies were prevalent in the south with lower numbers during the summer months. It is planned to use the merged data set to study the generation mechanisms and the propagation of the eddies within the study area.
Funding: This research was funded by Integrated Marine Observing System IMOS-Ocean Radar-2017-2019 grant number 53000300.