Interoperability of Direction-Finding and BeamForming High-Frequency Radar Systems : An Example from the Australian High-Frequency Ocean Radar Network

Direction-finding SeaSonde (4.463 MHz; 5.2625 MHz) and phased-array WEllen RAdar WERA (9.33 MHz; 13.5 MHz) High-frequency radar (HFR) systems are routinely operated in Australia for scientific research, operational modeling, coastal monitoring, fisheries, and other applications. Coverage of WERA and SeaSonde HFRs in Western Australia overlap. Comparisons with subsurface currents show that both HFR types agree well with current meter records. Correlation (R), root-mean-squares differences (RMSDs), and mean bias (bias) for hourly-averaged radial currents range between R = (−0.03, 0.78), RMSD = (9.2, 30.3) cm/s, and bias = (−5.2, 5.2) cm/s for WERAs; and R = (0.1, 0.76), RMSD = (17.4, 33.6) cm/s, bias = (0.03, 0.36) cm/s for SeaSonde HFRs. Pointing errors (θ) are in the range θ = (1◦, 21◦) for SeaSonde HFRs, and θ = (3◦, 8◦) for WERA HFRs. For WERA HFR current components, comparison metrics are RU = (−0.12, 0.86), RMSDU = (12.3, 15.7) cm/s, biasU = (−5.1, −0.5) cm/s; and, RV = (0.61, 0.86), RMSDV = (15.4, 21.1) cm/s, and biasV = (−0.5, 9.6) cm/s for the zonal (u) and the meridional (v) components. Magnitude and phase angle for the vector correlation are ρ = (0.58, 0.86), φ = (−10◦, 28◦). Good match was found in a direct comparison of SeaSonde and WERA HFR currents in their overlap (ρ = (0.19, 0.59), φ = (−4◦, +54◦)). Comparison metrics at the mooring slightly decrease when SeaSonde HFR radials are combined with WERA HFR: scalar (vector) correlations for RU, V, (ρ) are in the range RU = (−0.20, 0.83), RV = (0.39, 0.79), ρ = (0.47, 0.72). When directly compared over the same grid, however, vectors from WERA HFR radials and vectors from merged SeaSonde–WERA show RU (RV) exceeding 0.9 (0.7) within the HFR grid. Despite the intrinsic differences between the two types of radars used here, findings show that different HFR genres can be successfully merged, thus increasing current mapping capability of the existing HFR networks, and minimising operational downtime, however at a likely cost of slightly decreased data quality.


Introduction
Shore-based high-frequency radar (HFR) systems, operating in the frequency range between 3-30 MHz, are extensively used to remotely sense ocean currents in coastal areas, and for operational purposes [1].More than 100 HFR systems are presently in operation in the US [2], with increasing numbers of installations in the Asia-Pacific region [3,4] and Europe [5,6].Extensive analyses have proven the general reliability and limitations of HFR systems, including their ability to record oceanographic features at different spatial and temporal scales, their usefulness for the validation of numerical circulation models, and the potential for search and rescue purposes [7][8][9][10][11][12][13][14][15][16][17].The commercial direction-finding (SeaSonde) and phased-array (WEllen RAdar, or WERA HFR systems, provided respectively by Codar Ocean Sensors (COS) and Helzel MessTechnik [18,19], are the most commonly deployed systems, although other HFR systems, such as the PISCES developed by Neptune Radar Ltd, or the Least Expensive RAdar (or LERA) HFRs, are also available (e.g., [3,20]).
High-frequency radar systems rely on the Doppler shift of a radio wave signal transmitted to the ocean and backscattered from resonant Bragg-matching ocean waves.SeaSonde systems transmit a frequency-modulated interrupted continuous wave (FMICW) pulse from an omnidirectional vertical antenna and use a compact antenna system (two orthogonal loops and a monopole element) in combination with a direction-finding algorithm to resolve the azimuth of the incoming signal on a polar grid [21].At frequencies above 13.5 MHz, the transmit and receive elements can be combined in a single element with significant improvement on system compactness.At lower frequencies, separate transmit and receive elements are required, typically spaced one radar wavelength apart.
Phased-array WERA HFR systems transmit a frequency-modulated continuous wave (FMCW) pulse from a phased-array transmitter, usually composed of four vertical elements in a rectangular arrangement.A linear array of up to 16 elements, each typically spaced 0.45 radar wavelengths apart, is used to resolve the bearing of a signal through a beam-forming approach [22].Decoupling between the transmit and receive array is achieved through physical separation of the two arrays.
Direct comparisons between different types of HFR systems, in general, and between SeaSonde and WERAs in particular, are limited in number and temporal duration, with the majority performed about 20 years ago when the technology was less mature.This was also acknowledged in Reference [23], where an extensive comparison between SeaSonde and WERA HFR deployed in the West Florida shelf was conducted, and a comprehensive review of previous inter-comparisons was also given.For the sake of documentation, a summary of the reviews and the main findings from Reference [23] are presented here.Simultaneous current measurements from SeaSonde and beam-forming ocean surface current radars (OSCRs) were collected in Monterey Bay and compared against two days of acoustic Doppler current profiler (ADCP) observations at a 9-m depth collected in May 1995 [24].Multifrequency phased-array radars, OSCRs, and SeaSonde HFRs were deployed to measure the Chesapeake Bay outflow plume over a six-week experiment [25].When compared to near-surface (2-m depth) measurements, multi-frequency phased-array data were found to be accurate within 10 cm/s [26].A direct comparison of direction-finding and beam-forming accuracy was performed in 1996 between CODAR and WERA HFR along the German coasts [27].Additional comparisons between CODAR, WERA, and point-wise current meter data were undertaken in the same region in 2000 [28]; root-mean-square difference (RMSD) values ranged between 14-15 and 11-13 cm/s for CODAR and WERA radial velocity components, respectively.Compared to more recent findings, statistics were generally poorer in the earlier studies (see [1,[24][25][26][27][28][29]).).Reference [23] for instance reported RMSD in the range 5.1-9.2 and 3.8-6.5 cm/s for hourly-averaged SeaSonde and WERA HFR radial velocities over a three-month period in a low-energy shelf.As suggested in Reference [23], this may be due to a number of factors, such as improved reliability of in situ current meter data, and advancements in the HFR processing algorithms [30][31][32].Validation with ADCP data and inter-comparisons between SeaSonde and phased-array systems (not commercial WERA systems) have also been performed in other regions (e.g., [33]).The reported RMSD values between HFR and ADCP data were in the range (6.62, 11.3) cm/s for current velocity, and in the range (5.75, 13.7) cm/s for hourly-averaged radial velocity.
In the present work, we use Reference [23] as a starting point and extend the analyses to the combined SeaSonde-WERA HFR data, with the aim of proving the compatibility and interoperability of the two HFR types, for instance in case of downtime of one or more HFR components in an observing network.While merging radials from different frequencies appears to be a common practice (especially for SeaSonde HFR systems; see for instance Reference [29]), to our knowledge this is the first attempt to combine radial data from different HFR types and estimate the accuracy of the resulting vector maps for operational purposes.
The importance of this analysis is straightforward, considering the increasing number of deployments around the globe and the potential to fill in gaps in cases of operational downtime in an observing network without the need to apply gap-filling or optimal-interpolation techniques.
The paper is organized as follows.The Australian HFR network is introduced in Section 2, along with the data analysis approach and a description of the validation metrics.Subsurface current meter data are also introduced in this section.Results are given in Section 3, followed by sections containing a discussion of the main findings and the conclusions.

Materials and Methods
The Australian Ocean Radar facility is part of the Integrated Marine Observing System (IMOS), a national collaborative research infrastructure tasked with collection and dissemination of ocean data [34].The facility, based at the University of Western Australia, manages both direction-finding (SeaSonde) and phased-array (WERA) HFR systems.Long-range (9.335 MHz) WERA and long-range (4.463 MHz) SeaSonde HFR systems, with overlapping coverage, are in operation in the Rottnest (ROT) Shelf and Turquoise (TURQ) Coast (Western Australia, WA; Figure 1).
Separate intermediate range (13.5 MHz) WERA and long-range (5.2675 MHz) SeaSonde HFR stations are located at Coffs Harbor (COF) and Newcastle (NEWC), New South Wales (NSW).Long-range (9.335 MHz) WERA HFR systems are also operational in the South Australian Gulf (SAG) regions.Additional details on operational settings, such as frequency band and bandwidth, accuracy levels, and quality-control (QC) procedures specific to the WERA HFR systems, and main limitations of the Australian HFR systems, including strong interference within the operating frequency bands, are provided in References [34,35].
The two SeaSonde systems deployed in WA are located at Lancelin (LANC; 31 The standard software provided by the manufacturer is used to invert radial maps from the ocean Doppler spectra; hourly radial velocity data within the radar coverage were obtained by averaging around the cardinal hour a sequence of observations collected within an 80 min interval.The final output was produced every 60 min.At least two radial observations were required at each range and bearing in the final radial map.
To avoid mutual interferences between SeaSonde systems, transmitted chirps are synchronized between the two systems.The nominal angular resolution was set to 2 • ; calibrated antenna patterns were used at both stations, and provided an angular coverage of 180 • and approximately 150 • , respectively.To avoid interference to primary users, transmit power was limited to 8 W; maximum operational range exceeded 200 km offshore which reduced to 50 km or less with a typical day-to-night variability.
Two WERA systems, located at Guilderton (GUI; 31  1).Both sites use a four-element transmitter array, and a receive array of 16 equally-spaced elements.The systems operate at 9.335 MHz with a 33.3 kHz bandwidth, 1024-point fft, and 0.26 s sweep rate, an integration time of 5 min, and collect radial velocity data on a Cartesian grid with 4-km spacing every 10 min in alternate mode.As discussed in Reference [43], the 5-min sampling used here for the phased-array HFR systems was considered equivalent to the 80 min integration for the crossed-loop systems in terms of noise in the receive antennas.At the working frequency, sampling depth was approximately 1.3 m.Two WERA systems, located at Guilderton (GUI; 31°20.4960'S,115°29.3640'E),and Fremantle-Port Beach (FRE; 32°.0198'S, 115°44.7480'E)constitute the Rottnest Shelf (ROT) deployment in WA (Figure 1).Both sites use a four-element transmitter array, and a receive array of 16 equally-spaced elements.The systems operate at 9.335 MHz with a 33.3 kHz bandwidth, 1024-point fft, and 0.26 seconds sweep rate, an integration time of 5 minutes, and collect radial velocity data on a Cartesian grid with 4-km spacing every 10 minutes in alternate mode.As discussed in Reference [43], the 5-min sampling used here for the phased-array HFR systems was considered equivalent to the 80 min integration for the crossed-loop systems in terms of noise in the receive antennas.At the working frequency, sampling depth was approximately 1.3 m.
Hourly-averaged radial currents were generated from the QC 5 min radial velocity data by averaging at least three valid observations within the hour.For real-time operations, hourly averages were computed at minute 30 (i.e., at 00:30, 01:30, and on).In the following comparisons and merging, however, the averaging was performed consistently with SeaSonde HFR or the mooring data, to minimize temporal offsets between the datasets.Hourly-averaged radial currents were generated from the QC 5 min radial velocity data by averaging at least three valid observations within the hour.For real-time operations, hourly averages were computed at minute 30 (i.e., at 00:30, 01:30, and on).In the following comparisons and merging, however, the averaging was performed consistently with SeaSonde HFR or the mooring data, to minimize temporal offsets between the datasets.
The IMOS-Australian National Mooring Network (ANMN) Facility collects subsurface current data in the Rottnest Shelf region (Figure 1; Table 1).Quality-controlled data are available for download through the data portal [44][45][46][47][48].A set of compulsory quality-control tests including: impossible date, impossible depth, impossible horizontal and vertical velocities, correlation magnitude, and depth-correction for mooring knock-down, are performed on the time-series of subsurface velocities, along with other optional tests [49].Instruments deployed at the moorings include Nortek 190 KHz Continental current profilers at the WA Perth Canyon, 200 m (WACA20) mooring, a Teledyne RD Instruments (RDI) Workhorse Sentinel ADCP at the WA Rottnest National Reference Station (NRSROT), a Teledyne RDI Long Ranger ADCP at the WA Two Rocks 490 m (WATR50) mooring, and a Nortek Aquadopp Pro at the WA Two Rocks 103 m mooring (WATR10).Temporal resolution of the mooring records ranges between 10-30 min (Table 1).For consistency with HFR data, hourly-averaged subsurface currents were considered.The HFR and subsurface current data used in the following analyses span a four-month time period between 01 October 2017 and 31 January 2018.Comparison metrics include correlation coefficient (R), root-mean-square difference (RMSD), vector correlation magnitude and phase [50], as defined through Equations ( 1)-( 5).
In Equations ( 4) and ( 5), (u r,m , v r,m ) stand for the zonal and meridional components of the HFR (mooring) velocity components, respectively, i identifies the imaginary unit, and < > represents the ensemble average operator.
The analysis method described in Reference [9] is used to quantify bearing offsets in the radial currents for both SeaSonde and WERA HFR systems: the components of moored current velocity are first translated into radial current components in the direction of the HFR site, then R and RMSD values are calculated using HFR radial data at fixed range and at all bearings.In the presence of bearing offsets, the angle θ (min) that maximizes (minimizes) R (RMSD) does not correspond to the angular sector where the mooring is located (θ moor ).The difference between the two angles is a quantitative measure of the bearing errors.
The HFR current components are derived from radial velocity data on a Cartesian grid following Equations ( 6)-( 8), in which urad 1,2 > and θ 1,2 stand for the radial speed and direction from the two radar stations at each grid point: Additionally, mean bias between velocity components is calculated as follows: In the following analyses, vectors are mapped onto the WERA HFR grid and constraints are applied on the intersecting beam geometry in order to reduce errors of geometrical dilution of precision.SeaSonde HFR radials are converted from their polar coordinate systems and mapped on the same Cartesian grid before the radial-to-vector conversion step by averaging radial data within a 1 km distance from each WERA grid point.In order to avoid temporal mismatches between WERA and SeaSonde HFR radials, WERA radial data are averaged consistently with SeaSonde HFR data.
Standard statistical tests are used to derive significance levels for scalar correlation values.Methods described in [51][52][53][54] are used to derive confidence levels on amplitude of the vector correlation from the percentiles of the sample distribution obtained following a bootstrap approach, and assuming variables to be independent and identically distributed (IID) variables.First, the horizontal components of surface and subsurface currents are resampled to produce new current time series.This procedure was repeated 1000 times and amplitudes and phase angles were computed for each resampled pair to generate a distribution of correlation amplitudes.Then, the 95% confidence intervals for vector correlation were obtained from the 2.5% and 97.5% percentiles of the resulting distribution.
Finally, the Welch modified periodogram method [55] was used to extract a smooth estimate of the power density distribution of surface HFR and subsurface currents at selected locations.

Accuracy of WERA and SeaSonde Radial Data
Comparisons between HFR radial and subsurface currents in the direction of the radar are summarized in Table 2 for both WERA and SeaSonde HFR radars.Results shown here refer to the radar grid point closest to the mooring locations.For WERA HFR systems, comparison metrics are in relatively good agreement with previously reported findings in different ocean regions (for instance [1,2]), but poorer compared to a previous validation in the same area for a different time period [34].This is particularly true for FRE when compared to all the mooring data, with the exception of mooring WATR10.A similar trend is observed for the GUI HFR station, but in this comparison an improvement in comparison metrics can be observed for moorings WACA20 and WATR10 that in Reference [34] had the lowest agreement.One possible explanation for this pattern may be related to a different distribution of the noise level and the corresponding radial signal-to-noise ratio (SNR).
When the error analysis described in Reference [9] was performed, bearing offsets were observed (Table 3) for both beam-forming and direction-finding systems.For FRE radar, bearing offsets ranged between θ = [2 • , 9 • ], and between θ = [0 • , 8 • ] for the GUI station.For a 16-element receive array with 0.45 λ HF spacing as the ROT HFR WERA systems, the theoretical bearing resolution (λ HF /D, where D is the length of the antenna array and λ HF is the wavelength of the transmit signal) was approximately 8.5 • .As such, the bearing offsets are within the system's nominal angular resolution.
Comparison metrics between SeaSonde HFR radials and subsurface currents were lower than corresponding values for the WERA HFR systems: R was in the range R = (0.1, 0.42) for both GHED and LANC, and RMSD = (26.2,37.2) cm/s.Bearing errors for SeaSonde were also typically larger than beam-forming systems, in the range θ = (2 • , 21 • ) for LANC (Figure 2) and in the range θ = (8 • , 13 • ) for GHED radar, however, they were consistent with values reported at different installations.
The results of the comparisons for GHED HFR were likely biased due to the significant lower data coverage for this station at the moorings (Tables 2 and 3), located at the far edges of the range and angular coverages (Figure 1).Table 2. Comparison of surface HFR and subsurface radial data for the WERA and SeaSonde HFR stations.Radar look angle in the direction of the mooring (counterclockwise north) is also provided along with percent data available for the comparisons.R is the correlation coefficient between radial velocities as defined in Equation (1).Units for RMSD (bias) are cm/s.

Comparison of WERA and SeaSonde Vector Maps
For both SeaSonde and WERA HFR vector maps, comparison was limited to a marginal sector of the total grid (less than 25 points), located at the boundaries of the geometrically stable areas, and affected by the available data, as well as by the different grid resolution.
The WERA and SeaSonde vector grids have different spatial resolutions; as such, comparison is performed at the coarser grid cell points after remapping current data onto a common grid.A minimum of 50% data availability was required at each point to perform the statistical comparison.
Correlation values for the U (V) velocity components ranged between RU = (-0.02,0.45), RV = (0.29, 0.73) (p-values < 0.001), while RMSDU, V values ranges between (22,49) cm/s and (19,59) cm/s for the U, V components respectively.Typical values for amplitude and phase angle of the vector correlation are in the range  = (0.19, 0.59), and  = (−4°, 54°).The bearing offset reported from the LANC-NRSROT mooring (Figure 2) can originate both from distortions in the antenna beam pattern, and from incorrect settings in the processing software, specifically a missing land mask which is typically used to exclude specific directions or islands.It is worth noting that low correlation values do not necessarily suggest the presence of errors in the data, but may originate instead from the predominant current direction, especially if the low R is associated with low RMSD values as was the case of NRSROT (WACA20) and FRE in Table 2. Also, all correlation values are statistically significant (p-values < 0.001).
Finally, Tables 2 and 3 show that SeaSonde HFR perform better than WERAs in terms of mean bias, with values in the range (0.03, 0.33) cm/s and (−5.2, +5.2) cm/s, respectively.

Comparison of WERA and SeaSonde Vector Maps
For both SeaSonde and WERA HFR vector maps, comparison was limited to a marginal sector of the total grid (less than 25 points), located at the boundaries of the geometrically stable areas, and affected by the available data, as well as by the different grid resolution.
The WERA and SeaSonde vector grids have different spatial resolutions; as such, comparison is performed at the coarser grid cell points after remapping current data onto a common grid.A minimum of 50% data availability was required at each point to perform the statistical comparison.

Comparison of Merged WERA-SeaSonde Merged Vector Maps
The WA HFR network configuration (Figure 1) presents several possible pairs of WERA-SeaSonde combinations (FRE-LANC; FRE-GHED; GUI-LANC; GUI-GHED) (Figure 3); however, only one of them (FRE-LANC; Figure 3b) provides vector maps over a spatial domain comparable with the FRE-GUI configuration and with enough temporal coverage for reliable statistical analyses, thus offering the possibility of a direct comparison with subsurface currents at the mooring locations.Others provide either more limited spatial coverage and no overlap with the moorings, in spite of the large temporal data availability (GUI-GHED; Figure 3c), or too low data return (FRE-GHED; Figure 3d) at the moorings.

Comparison of Merged WERA-SeaSonde Merged Vector Maps
The WA HFR network configuration (Figure 1) presents several possible pairs of WERA-SeaSonde combinations (FRE-LANC; FRE-GHED; GUI-LANC; GUI-GHED) (Figure 3); however, only one of them (FRE-LANC; Figure 3b) provides vector maps over a spatial domain comparable with the FRE-GUI configuration and with enough temporal coverage for reliable statistical analyses, thus offering the possibility of a direct comparison with subsurface currents at the mooring locations.Others provide either more limited spatial coverage and no overlap with the moorings, in spite of the large temporal data availability (GUI-GHED; Figure 3c), or too low data return (FRE-GHED; Figure 3d) at the moorings.Statistics of the comparison between the FRE-GUI HFR pair and the subsurface currents at the mooring locations are provided in Table 4; for the sake of comparison, results of the cross-validation of unmerged and merged vector data at the same mooring locations are provided in Table 5.
At the mooring locations (Table 4), correlation values for unmerged data range between (-0.12, 0.83) for the U component, lower than corresponding values for V, which is found to range between Statistics of the comparison between the FRE-GUI HFR pair and the subsurface currents at the mooring locations are provided in Table 4; for the sake of comparison, results of the cross-validation of unmerged and merged vector data at the same mooring locations are provided in Table 5.At the mooring locations (Table 4), correlation values for unmerged data range between (-0.12, 0.83) for the U component, lower than corresponding values for V, which is found to range between (0.6, 0.86).Amplitude of the vector correlation generally exceeds 0.58, with the highest value (0.86) for the deep-water mooring (WATR50), which also shows the lowest angular offset (-6 • ) with HFR currents in spite of the distance from surface.Relatively poor comparison metrics are observed at the NRSROT mooring, and in particular for the U component (R U = (-0.12,0.31)), the RMSD V and the veering angles, which exceed 20 • (absolute values).
It is interesting that at several deep-water locations (for instance: WATR50; WATR10; WACA20), the cross-validation statistics suggest good matches between time-series at the mooring locations, with ρ in the range (0.65, 0.86), angular offset in the range ϕ = (-5 • , 28 • ), RMSD U = (13.2,15.7) cm/s, RMSD V = (15.4, 17.3) cm/s, and low bias values (bias U = (-5.1,-1.9) cm/s; bias V = (6.5, 9.6) cm/s).If a time-lagged cross-correlation is performed between surface HFR currents and subsurface data, results (not shown here) suggest that there is no lag between the currents at the two different depths.A spectral analysis of HFR and subsurface currents for the WATR50 mooring suggests that the majority of the differences are found within the diurnal-inertial frequency bands while the low-and high-frequency bands have comparable energy levels between mooring and HFR currents (Figure 4).Similar differences in this frequency band are commonly observed when HFR and subsurface currents are compared.This feature was also documented by Reference [56], who reported a surface-intensified sea breeze response in their early SeaSonde observations.high-frequency bands have comparable energy levels between mooring and HFR currents (Figure 4).Similar differences in this frequency band are commonly observed when HFR and subsurface currents are compared.This feature was also documented by Reference [57], who reported a surfaceintensified sea breeze response in their early SeaSonde observations.Comparison metrics for the merged SeaSonde−WERA and subsurface currents at the mooring locations (Table 6) match the trend reported in Table 4 for the unmerged data, although with lower correlation and higher RMSD values: R U = (−0.2,0.83), RMSD U = (14.3,24.2) cm/s, bias U = (−9.5,1.1); R V = (0.39, 0.79), RMSD V = (19.9,26.2) cm/s, bias V = (−1.3,9.6), ρ = (0.47, 0.72), and ϕ = (3 • , 51 • ).The spatial distributions of R U , R V , and ρ between merged and unmerged data are provided in Figure 5, along with the mean current pattern for the analysis period.R U values typically exceeded 0.95 within the radar footprint (Figure 5a), while R V values also exceeded 0.9 (Figure 5b), suggesting good agreement between datasets.Good agreement was confirmed in the general circulation pattern: the offshore southward current and the weaker northward current mapped by the WERA HFR radar (black arrows), associated, respectively, with the Leeuwin Current (LC) and the Capes Current (CC) [40,42], they are also well reconstructed in the merged WERA-SeaSonde data (red arrows).Typical bias for the zonal component (Figure 6a) was below 0.15 m/s across the radar domain, and increases for the meridional component, particularly in areas of weak currents (Figure 6b).Lower R U , R V , and ρ were typically observed offshore; this poorer agreement was reflected in the anomalous and inconsistent mean current vectors at the outer edges of the domain (Figure 5d), where bias values also become significant.
inconsistent mean current vectors at the outer edges of the domain (Figure 4d), where bias values also become significant.

Discussion
Developed initially as a research-based tool, HFR technology has reached a maturity stage such that national and international operational programs and networks make extensive use of HFR capabilities in order to provide reliable sea surface maps over a variety of temporal and spatial scales [1][2][3][4].Several genres of HFR systems are available [18,22], classified on the basis of the approach used to resolve the direction of the sea echo.Extensive validation exercises have proven the general reliability and provided accuracy estimates of both types of HFR systems over a wide variety of deployments [1].With the exception of a limited number of test cases or cross-validation exercises [22][23][24][25][26][27][28], limited attempts have been made to compare vectors from different systems in overlapping areas.Even fewer attempts have been made to combine surface radial current data from different systems in their areas of common overlap, and assess the accuracy of the resulting current maps.In this study, we used Reference [23] as a starting point and further explored the possibility of combining radial velocity data from the two major commercial ocean HFR systems (direction finding SeaSonde and beam forming WERA HFR) that are operationally managed along the coast of Western Australia (WA; Figure 1; [34]).
First, we use independent subsurface current data measured at various locations within the HFR footprint to provide estimates of radial data accuracy for each of the two HFR systems.Results show that in general, both HFRs provide reliable and accurate information.Beam-forming WERA HFRs seem to perform slightly better than SeaSonde in terms of correlation, RMSD values, and bearing offsets.Bearing offsets are commonly observed in direction-finding HFR systems [9], and are not necessarily uniform within range from the receive antennas.In the case of the LANC SeaSonde HFR system, offsets can be as high as 30°.Large offsets are also found for GHED SeaSonde HFR systems; however, they should be interpreted in terms of the low percent return rates at the radar far ranges at the mooring locations.If large bearing offsets occur at far ranges from the HFR receiver, they translate into significant spatial offsets with potential to bias the accuracy of the radial maps, thus decreasing the quality of the ocean current maps.Bearing offsets are also present in the WERA HFR systems; however, they are typically lower than SeaSonde and are within the angular accuracy of the beam-forming systems [43].While antenna pattern distortions for direction-finding HFRs, such as the SeaSonde, manifest themselves in pointing errors, distortions of the antenna pattern leading to

Discussion
Developed initially as a research-based tool, HFR technology has reached a maturity stage such that national and international operational programs and networks make extensive use of HFR capabilities in order to provide reliable sea surface maps over a variety of temporal and spatial scales [1][2][3][4].Several genres of HFR systems are available [18,22], classified on the basis of the approach used to resolve the direction of the sea echo.Extensive validation exercises have proven the general reliability and provided accuracy estimates of both types of HFR systems over a wide variety of deployments [1].With the exception of a limited number of test cases or cross-validation exercises [22][23][24][25][26][27][28], limited attempts have been made to compare vectors from different systems in overlapping areas.Even fewer attempts have been made to combine surface radial current data from different systems in their areas of common overlap, and assess the accuracy of the resulting current maps.In this study, we used Reference [23] as a starting point and further explored the possibility of combining radial velocity data from the two major commercial ocean HFR systems (direction finding SeaSonde and beam forming WERA HFR) that are operationally managed along the coast of Western Australia (WA; Figure 1; [34]).
First, we use independent subsurface current data measured at various locations within the HFR footprint to provide estimates of radial data accuracy for each of the two HFR systems.Results show that in general, both HFRs provide reliable and accurate information.Beam-forming WERA HFRs seem to perform slightly better than SeaSonde in terms of correlation, RMSD values, and bearing offsets.Bearing offsets are commonly observed in direction-finding HFR systems [9], and are not necessarily uniform within range from the receive antennas.In the case of the LANC SeaSonde HFR system, offsets can be as high as 30 • .Large offsets are also found for GHED SeaSonde HFR systems; however, they should be interpreted in terms of the low percent return rates at the radar far ranges at the mooring locations.If large bearing offsets occur at far ranges from the HFR receiver, they translate into significant spatial offsets with potential to bias the accuracy of the radial maps, thus decreasing the quality of the ocean current maps.Bearing offsets are also present in the WERA HFR systems; however, they are typically lower than SeaSonde and are within the angular accuracy of the beam-forming systems [43].While antenna pattern distortions for direction-finding HFRs, such as the SeaSonde, manifest themselves in pointing errors, distortions of the antenna pattern leading to significant side lobe weighting in a phased-array system, such as the WERA system, may cause biases in the amplitude estimates of the radial currents.
Effects of the vertical shear associated with different operating frequencies (4.463 MHz and 9.335 MHz for the SeaSonde and WERA HFR systems, respectively) are neglected at this stage, due to the fact that most of the comparisons with subsurface current data are performed well below the surface.For instance, the nominal distance from surface for the first usable current meter bin exceed 50 m in case of the deep-water WATR50 mooring, yet R and RMSD values for both direction finding and beam forming systems are consistent.
Second, we directly compare vector currents from the beam forming (FRE and GUI) HFRs to those derived from pairs of direction finding (LANC and GHED) HFR systems.The overlap area is limited to a small sector of their coverages, and is controlled by the geometry of the intersecting radar-look angles.Comparisons in this common area show a relatively good match despite different grid resolutions and vector mapping methods.
Differences may originate from a number of causes.For WERA HFRs, vectors are derived as described in Equations ( 6)-( 8) using a 1-to-1 matching of radial velocities from the two stations at each grid point.Hourly WERA radial data are averaged from a sequence of 5-min integration samples.SeaSonde vector maps are derived through a similar least-squares fit, which uses radial data within a 10-km distance from each grid point.This second approach biases direction and speed of the resulting current maps, especially at the far edges of the radar grid, both in the case of simplified and more complex current fields [57], and propagates radial velocity errors across adjacent vector grid points.
Other biases may arise from the temporal sampling scheme used in the hourly SeaSonde radial maps; however, this is accounted for properly in the radar configurations.Typically, hourly radial maps are derived collecting several intermediate short-term intermediate products ("short term radials" following the standard operational SeaSonde jargon), that are then "merged" to define a final product through a median filtering method.The WERA radials are, on the other hand, derived by averaging a sequence of radial maps collected every 10 min, and requires at least three valid observations at each grid point.To avoid biases from different merging approaches, the Australian SeaSonde systems use an averaging merging approach consistent with that used for WERA HFR systems.
Third, we compared vector maps from WERA systems with merged SeaSonde-WERA radial data on the same grid.We used Equations ( 6)-(8) to compute vectors from LANC-FRE HFR systems for a period of 4 months between October 2017 and January 2018.The period was chosen on purpose: due to hardware problems with one of the WERA systems, no real-time (RT) vectors were available over the Rottnest Shelf (Figure 1); and both SeaSonde systems had been recently calibrated and optimized for real-time operations.
To set a confidence limit to WERA vector accuracy, comparisons with subsurface data at the mooring locations were first performed.Then, comparison metrics were derived across the entire radar domain.Comparison metrics for unmerged and merged data against subsurface data were consistent, although merged SeaSonde-WERA had slightly poorer performances.
This may be related to a combination of factors, including the re-mapping of LANC radar radial velocities from a polar coordinate system to a higher resolution Cartesian grid where FRE radials are sampled; or the effects of the vertical shear associated with the different operating frequencies of the two HFR systems.As this was a preliminary test of system interoperability, the latter factor was neglected in first approximation, also due to the lack of high-resolution surface current data in the region from either drifters or glider observations that would provide quantitative estimates of the velocity shear effects.Further investigation is planned to quantify the contribution of similar possible error sources through the deployment of near surface drifters within the radar coverage.
When comparing the time-averaged current pattern, results show good agreement between the two data sets and further support the hypothesis of interoperability of the different HFR systems.There is good agreement between the main offshore meandering structure as well as the current reversal in the coastal shelf area (Figure 5d), that are associated respectively with the Leeuwin Current (LC) and Capes Current (CC) system [40,42].
There are regions where comparison is qualitatively and quantitatively poor, but they are located either at far ranges from LANC HFR, where radial velocity errors are commonly observed due to the poor SNR constraints, or at the edges of the FRE station where side lobes may bias amplitude estimates of the radial currents.

Conclusions
Along the coast of Western Australia, the IMOS Ocean Radar Facility operates the two major commercial-type HFR systems, the direction-finding SeaSonde and the beam-forming WERA ocean radars.
In spite of their overlapping radial coverage, they routinely operate separately collecting ocean surface maps representative of the shelf areas of the Turquoise Coast (SeaSonde systems) and the Rottnest Shelf (WERA systems).The relative short distance between individual HFR stations, combined with the extensive availability of subsurface currents from different moorings within the HFR coverage, provide a unique opportunity to determine accuracy levels of radial and current velocities, and investigate the possibility of combining radial velocity data from different platforms.This latter aspect is of particular interest in light of the fact that it can potentially fill in the gaps in case of downtime of one component of an HFR observing network.
We tested this hypothesis using a four-month data period between October 2017 and January 2018.This time period was chosen on purpose: one of the WERA HFR systems was shut down for a prolonged period of time due to hardware issues, so it was necessary to investigate the interoperability of the two HFR in order to fill the gap (although in delayed mode only).
In spite of the intrinsic differences between the direction-finding SeaSonde HFR systems (polar coordinate grid, 2 • nominal angle resolution, 6 km range resolution, measurement depth of 2.7 m, integration time 80 minutes with 1-h output rate, non-uniform bearing errors over the operational range), and the beam-forming WERA (cartesian coordinate system with 4.5 km resolution, effective measurement depth 1.3 m, integration time of 5 minutes with 1-h output rate, low bearing errors), surface vector maps can be derived with reasonable accuracy.
Combining radials from different HFR genres is not commonly performed, however it is shown here that this is possible, thus increasing the spatial and temporal coverage of the HFR observing network, especially for instance in case of prolonged downtime of an HFR component.

22 Figure 1 .
Figure 1.Coverage of the high-frequency radar (HFR) systems deployed along the Western Australia (WA) coast, with the dominant bathymetric features.Red (black) markers are used for the WEllen RAdar WERA (SeaSonde) Fremantle (FRE) and Guilderton (GUI) stations (Green Head (GHED) and Lancelin (LANC)); red (black) areas show the typical coverage of each HFR station.Mooring locations are marked as black diamonds inside the HFR footprints.NRSROT, WATR50, WATR10 and WACA20 stand for Rottnest Reference Station, WA Two Rocks 490 m, WA Two Rocks 103 m, and WA Perth Canyon 200 m moorings, respectively.

Figure 1 .
Figure 1.Coverage of the high-frequency radar (HFR) systems deployed along the Western Australia (WA) coast, with the dominant bathymetric features.Red (black) markers are used for the WEllen RAdar WERA (SeaSonde) Fremantle (FRE) and Guilderton (GUI) stations (Green Head (GHED) and Lancelin (LANC)); red (black) areas show the typical coverage of each HFR station.Mooring locations are marked as black diamonds inside the HFR footprints.NRSROT, WATR50, WATR10 and WACA20 stand for Rottnest Reference Station, WA Two Rocks 490 m, WA Two Rocks 103 m, and WA Perth Canyon 200 m moorings, respectively.

Figure 2 .
Figure 2. Bearing error analysis for LANC HFR SeaSonde site, showing the color-coded magnitude of correlation between HFR radial velocity over bearing and subsurface currents at the mooring locations.Black markers are used for the mooring locations, while red markers identify the locations that maximize the correlation at each specific range.

Figure 2 .
Figure 2.Bearing error analysis for LANC HFR SeaSonde site, showing the color-coded magnitude of correlation between HFR radial velocity over bearing and subsurface currents at the mooring locations.Black markers are used for the mooring locations, while red markers identify the locations that maximize the correlation at each specific range.

Figure 3 .
Figure 3. Data return maps for different HFR network configuration: (a) refers to FRE-GUI stations (default ROT WERA configuration); (b) shows the coverage obtained merging FRE and LANC radial data; (c) same as (b) but for the GUI-LANC pair; and (d) coverage obtained with the FRE-GHED configuration.The colormap documents the percent data availability for the different configurations.Locations of the moorings are marked as red squares.Refer to Figure 1 for the codes of the HFR systems and the moorings.

Figure 3 .
Figure 3. Data return maps for different HFR network configuration: (a) refers to FRE-GUI stations (default ROT WERA configuration); (b) shows the coverage obtained merging FRE and LANC radial data; (c) same as (b) but for the GUI-LANC pair; and (d) coverage obtained with the FRE-GHED configuration.The colormap documents the percent data availability for the different configurations.Locations of the moorings are marked as red squares.Refer to Figure 1 for the codes of the HFR systems and the moorings.

Figure 4 .
Figure 4. Rotary spectra for the WATR50 data at the 56 m level and corresponding rotary spectra for HFR surface currents at the radar grid cell closest to the mooring location.Units for frequency and power spectral density are cycles per hour (cph) and (m/s) 2 /cph, respectively.

Figure 5 .Figure 5 .
Figure 5. Spatial distribution of correlation between WERA and merged SeaSonde-WERA radial current maps, and corresponding time-averaged current patterns: (a) refers to the zonal component; (b) refers to the meridional component; (c) is the amplitude of the vector correlation; and (d) presents the time-averaged current patterns for the WERA (red arrows) and the merged SeaSonde-WERA Figure 5. Spatial distribution of correlation between WERA and merged SeaSonde-WERA radial current maps, and corresponding time-averaged current patterns: (a) refers to the zonal component; (b) refers to the meridional component; (c) is the amplitude of the vector correlation; and (d) presents the time-averaged current patterns for the WERA (red arrows) and the merged SeaSonde-WERA (black arrows).Vectors are plotted every two grid points and offset by one grid point to improve readability.Locations of the moorings are marked as red squares.

Figure 6 .
Figure 6.Spatial distribution of the mean bias (m/s) for the two velocity components between WERA and merged SeaSonde-WERA: (a) refers to the zonal component; (b) refers to the meridional component.Locations of the moorings are marked as white squares.

Figure 6 .
Figure 6.Spatial distribution of the mean bias (m/s) for the two velocity components between WERA and merged SeaSonde-WERA: (a) refers to the zonal component; (b) refers to the meridional component.Locations of the moorings are marked as white squares.

Table 1 .
Details of the mooring in the Rottnest Shelf region.See Figure1for their locations within the HFR coverage area.NRSROT, WATR50, WATR10 and WACA20 stand for Rottnest Reference Station, WA Two Rocks 490 m, WA Two Rocks 103 m, and WA Perth Canyon 200 m moorings, respectively.

Table 3 .
Same as Table1but for the radar angle at which correlation with mooring currents is maximized.Results are displayed for both WERA and SeaSonde HFR stations.Units for RMSD (bias) are cm/s.

Table 4 .
Comparison between WERA surface HFR currents and subsurface current data at the mooring locations.Units for RMSD and phase angle are cm/s and degrees, respectively.Positive (negative) values for phase angle indicate that surface HFR currents are rotated counterclockwise (clockwise) from subsurface currents.

Table 5 .
Comparison between WERA HFR currents and SeaSonde−WERA HFR currents at the mooring locations.Units for RMSD and phase angle are cm/s and degrees, respectively.Positive (negative) values for phase angle indicate that WERA HFR currents are rotated counterclockwise (clockwise) from SeaSonde−WERA HFR currents.

Table 6 .
Same as Table3abut refers to the combined SeaSonde−WERA HFR data.

Table 6 .
Same as Table3abut refers to the combined SeaSonde−WERA HFR data.