# Best Practices on High Frequency Radar Deployment and Operation for Ocean Current Measurement

Carlo Mantovani, Lorenzo Corgnati, Jochen Horstmann, Anna Rubio, Emma

Reyes, Céline Quentin, Simone Cosoli, Jose Luis Asensio, Julien Mader,

Annalisa Griffa

► To cite this version:

Carlo Mantovani, Lorenzo Corgnati, Jochen Horstmann, Anna Rubio, Emma Reyes, et al. Best Practices on High Frequency Radar Deployment and Operation for Ocean Current Measurement. Frontiers in Marine Science, 2020, 7, 10.3389/fmars.2020.00210. hal-02559129

HAL Id: hal-02559129

https://univ-tln.hal.science/hal-02559129

Submitted on 15 Apr 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Best Practices on High Frequency Radar Deployment and Operation for Ocean Current Measurement

Carlo Mantovani¹*, Lorenzo Corgnati¹, Jochen Horstmann², Anna Rubio³, Emma Reyes⁴, Céline Quentin⁵, Simone Cosoli⁶, Jose Luis Asensio³, Julien Mader³ and Annalisa Griffa¹

OPEN ACCESS

**Edited by:**
Juliet Hermes,
South African Environmental
Observation Network (SAEON),
South Africa

**Reviewed by:**
Jeffrey Dean Paduan,
Naval Postgraduate School,
United States
Brian Emery,
University of California, Santa
Barbara, Unite States

\*Correspondence:
Carlo Mantovani
carlo.mantovani@cnr.it;
carlo.mantovani@sp.ismar.cnr.it

**Specialty section:**
This article was submitted to
Ocean Observation,
a section of the journal
Frontiers in Marine Science

**Received:** 15 November 2019
**Accepted:** 17 March 2020
**Published:** 09 April 2020

**Citation:**
Mantovani C, Corgnati L, Horstmann J, Rubio A, Reyes E, Quentin C, Cosoli S, Asensio JL, Mader J and Griffa A (2020) Best Practices on High Frequency Radar Deployment and Operation for Ocean Current Measurement. *Front. Mar. Sci.* 7:210. doi: 10.3389/fmars.2020.00210

¹Institute of Marine Sciences, National Research Council of Italy, Pozzuolo di Lerici, Italy, ²Helmholtz-Zentrum Geesthacht, Institute of Coastal Research, Geesthacht, Germany, ³AZTI, Marine Research, Basque Research and Technology Alliance (BRTA), Pasaia, Spain, ⁴Coastal HF Radar Facility, SOCIB - Balearic Islands Coastal Ocean Observing and Forecasting System, Palma, Spain, ⁵Université de Toulon, Aix-Marseille University, CNRS, IRD, MIO UM110, Toulon, France, ⁶Ocean Graduate School and the UWA Oceans Institute, The University of Western Australia, Crawley, WA, Australia

High Frequency Radar (HFR) technology refers to land based remote sensing instruments capable of measuring surface currents and ocean waves at ranges up to 200 km or more. HFR technology is widely acknowledged as a cost-efficient tool to monitor coastal regions and has potential use to monitor coastal regions all over the world. Globally, the number of HFR stations is steadily increasing. Regional networks provide real-time data in support of operational activities such as search and rescue operations, fast response in case of maritime accidents and spill of pollutants, and resource management. Each operator needs a general understanding of the working principles in order to ensure that instruments are managed properly. A set of harmonized quality assurance and quality control procedures is recommended, along with an effective approach to HFR data discovery and dissemination, to provide high quality measurements to the end users. Different documents providing best practices for operation and maintenance have emerged in the past years. They are oriented either to Direction Finding (DF) or Beam Forming (BF) systems, or to specific manufacturer's radar systems. The main objective of this paper is to offer a comprehensive "Best Practices" document in an effort of ensuring consistency between different deployments and harmonized operations of HFR systems. This, regardless of system manufacturer, antenna design and setup. A homogeneous approach is given when possible, general concepts and definitions are introduced in order to provide a framework for both data processing and management steps. Examples are also given from the European HFR users with focus on Near Real Time data flow suitable for operational services.

**Keywords:** high frequency radar, best practices, operational oceanography, remote sensing, surface currents, data management, quality control
# INTRODUCTION

High Frequency ocean Radar (HFR) systems are cost-efficient tools to monitor coastal regions at a range of up to 200 km or more, with potential to monitor coastal regions all over the world. Oceanographic HFRs are mainly utilized to measure ocean surface current fields (Paduan and Rosenfeld, 1996; Gurgel et al., 1999) for applications such as search and rescue (Ullman et al., 2006), oil spill monitoring (Abascal et al., 2009), marine traffic information (Breivik and Saetra, 2001), calibration and validation of numerical circulation models, and data assimilation (Paduan and Shulman, 2004; Barth et al., 2008). Further applications of the HFR systems include mapping of sea state parameters (significant wave height, wave period and direction; Wyatt, 1990; Gurgel et al., 2006), wind direction mapping (Heron and Rose, 1986; Shen et al., 2012; Kirincich, 2016), tsunami early warning systems (Lipa et al., 2006; Gurgel et al., 2011), and ship detection and tracking (Ponsford et al., 2001; Maresca et al., 2014).

HFR technology relies on the Bragg resonant mechanism, that is a resonant backscatter resulting from coherent reflection of a transmitted radio wave by a particular train of ocean waves with half the wavelength of the transmitted signal. Since ocean waves travel with a non-zero speed, the frequency of the backscattered radio wave is Doppler shifted. For a radar transmitting at a central frequency $f_0$, two narrow peaks (Bragg peaks) are theoretically expected in the backscattered spectrum, shifted $\pm \Delta f$ from $f_0$, and associated with the waves traveling toward (right peak) and away (left peak) from the radar, along the radial direction (Crombie, 1955). When ocean waves propagate over a current field, an additional Doppler shift adds up to the one previously described.

Such frequency offset corresponds to the velocity of the component of the current approaching to, or receding from, the receiver, and is referred to as radial velocity vector in the standard operational HFR jargon. Further in-depth analysis of the full spectra of the backscattered signals can also provide information on the sea state, wind, tsunami and discrete targets (e.g., vessels).

However, extracting information other than surface currents presents a much greater challenge since these are obtained from much weaker or partial parts of the signal, which are more likely to be corrupted by noise and interference (Barrick, 1977; Wyatt et al., 2006) and require integration times longer than those used for currents (Wyatt et al., 2013). To date, the great majority of HFRs are used to map ocean surface currents; as such, the best practices contained in this document apply to currents if not explicitly mentioned otherwise.

HFR systems need to resolve in range and azimuth their targets, i.e., the ocean patches contributing to the backscattered radio signal, in order to assign them a velocity value. While range gating is commonly used across different platforms to resolve for range, the determination of the direction of arrival (DOA) presents greater challenges. Commonly, HFR systems are classified based on the technique they use for DOA inversion: beam forming (BF) and direction finding (DF) systems.

BF systems, such as for instance the University of Hamburg Wellen Radar (WERA) (Gurgel et al., 1999) and the University of Hawaii LERA (Flament et al., 2016), require a receiving antenna array. They electronically steer the receiving array to the direction of each selected patch of the ocean surface and resolve the Doppler spectrum on a predefined grid (within the angular range and resolution allowed by the array configuration). Radial current velocities are therefore directly associated with each grid cell. Direction-finding (DF) systems acquire the backscattered spectrum as superposition of the echoes generated by all the ocean patches falling within a range ring and resolve the DOA on a polar grid by analyzing the covariance matrix of the complex voltages received by the directional antennas (Emery, 2018). In the latest developments, such as for instance in the SeaSonde systems (Barrick and Lipa, 1997), a tailored version of the Multiple Signal Classification (MUSIC) algorithm (Schmidt, 1986) is applied as DF algorithm.

HFR systems differ also in the transmit approach, as they can use pulse continuous (CW), frequency modulated continuous (FMCW) or interrupted (FMICW) waveforms. Different transmit approaches dictate the use of separate or co-located receive and transmit antennas and the method for resolving the range (Table 1).

**Abbreviations:** ACORN, IMOS Ocean Radar Facility; ADSL, Asymmetric digital subscriber line; AIS, Automatic Identification System; AODN, Australian Ocean Data Network; APM, Antenna Pattern Measurement; BF, Beam Forming; CC, Creative Commons; CF, Climate and Forecast; CMEMS, Copernicus Marine Environment Monitoring Service; COS, Codar Ocean Sensors; DF, Direction Finding; DOA, Direction of Arrival; DOI, Digital Object Identifier; EIRP, Effective Isotropic Radiated Power; FAIR, Findable, Accessible, Interoperable, and Reusable; FMCW, Frequency-Modulated Continuous-Wave; FMICW, Frequency-Modulated Interrupted Continuous-Wave; FOLs, First-Order Limits; GDOP, Geometrical Dilution Of Precision; GOOS, Global Ocean Observing System; GPS, Global Positioning System; HFR, High-Frequency Radar; HTTP, Hypertext Transfer Protocol; IEEE, Institute of Electrical and Electronics Engineers; IMOS, Integrated Marine Observing System; IOOS, Integrated Ocean Observing System; IT, Information Technology; ITU, International Telecommunication Union; JCOMM, Joint Technical WMO-IOC Commission for Oceanography and Marine Meteorology; LAN, Local Area Network; LEMP, Lightning electromagnetic pulse protection; MAC, Media Access Control; NetCDF, Network Common Data Form; QA, Quality Assessment; QARTOD, Quality Assurance/Quality Control of Real-Time Oceanographic Data; QC, Quality Control; RF, Radio Frequency; SDN, SeaDataNet; SIM, Subscriber Identity Module; SIM, Subscriber Identity Module; SMS, Short Message Service; SSH, Secure Shell; THREDDS, Thematic Real-time Environmental Distributed Data Services; UCSB, University of California-Santa Barbara; ULSD, Unweighted Least Squares Fitting; UPS, Uninterruptible Power Supply; VNC, Virtual Network Computing; (V)SWR, (Voltage) Standing Wave Ratio.

<table><caption>**TABLE 1** | List of HF radar systems found in literature and their techniques for obtaining the range and azimuth of the sea echo.</caption><thead><tr><th rowspan="2">HF radar system name</th><th rowspan="2">Range resolution method</th><th colspan="2">Azimuth resolution method</th></tr><tr><th>Direction Finding</th><th>Beam Forming</th></tr></thead><tbody><tr><td>CODAR/NOAA (Barrick et al., 1977)</td><td>Pulse CW</td><td>+</td><td>-</td></tr><tr><td>COSRAD (Heron et al., 1985)</td><td>Pulse CW</td><td>-</td><td>+</td></tr><tr><td>PISCES (Shearman and Moorhead, 1988)</td><td>FMICW</td><td>-</td><td>+</td></tr><tr><td>OSCR (Prandle et al., 1992)</td><td>Pulse CW</td><td>-</td><td>+</td></tr><tr><td>C-CORE (Hickey et al., 1995)</td><td>FMICW</td><td>-</td><td>+</td></tr><tr><td>SeaSonde (Barrick and Lipa, 1997)</td><td>FMICW</td><td>+</td><td>-</td></tr><tr><td>WERA (Gurgel et al., 1999)</td><td>FMCW</td><td>+</td><td>+</td></tr><tr><td>LERA (Flament et al., 2016)</td><td>FMCW</td><td>+</td><td>+</td></tr></tbody></table>
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Concept of surface current derivation from a two HFR site network. Individual sites measure a radial current component, which can be mapped to a surface current vector by a combination algorithm.</figcaption></figure>

Two HFR stations are required at least to resolve the two-dimensional current field from radial velocity data. The two sites must be adequately spaced and overlook at the same ocean area from two different angles. Separation distance is controlled by the operational range of each HFR system, which in turn depends mainly on its operating frequency. Separation also controls the shape and the size of the HFR domain where currents can be resolved through the geometry of the intersection angles of the radial vectors (geometrical dilution of precision, GDOP; **Figure 1**).

Other factors controlling range include water salinity, sea state, interference, and integration time. Range resolution typically varies between a few hundred meters to 12 km, while the azimuthal resolution is both platform and algorithm-dependent, with typical values between 5° and 18°.

Worldwide, approximately 400 HFR systems have been installed, for the most diverse range of applications (Roarty et al., 2016, 2019). In Europe, the number of systems is growing with 64 HFR systems implemented and several others in the planning stage. From the last available survey performed (Mader et al., 2016; Roarty et al., 2019), the most commonly deployed commercial HFRs are SeaSonde systems for DF and WERA systems for BF groups (see **Table 1**). Therefore, recommendations and examples provided here mostly refer to these specific commercial brands.

Nowadays, HF radar systems are integrated in many European coastal observatories with proven potential for monitoring coastal currents and providing inputs for operational data

assimilation and assessment of numerical ocean forecasting models, especially near the coast (Barth et al., 2008, 2011; Marmain et al., 2014; Stanev et al., 2015; Iermano et al., 2016). The growing number of HFR systems, their optimization from installation to maintenance and operation, and the need for standardized data processing and analysis approaches, have driven the HFR community to work at European and international levels toward a coordinated development of coastal HFR technology and its products (Rubio et al., 2017, 2018; Roarty et al., 2019).

In response to the need for optimizing operation performance, different documents providing best practices for operation and maintenance have emerged in the past years. A list of guidelines and best practice documents have been used as a basis for the elaboration of this paper (Cook et al., 2008; SeaSonde® Remote Unit Operator's Manual, 2010; Philip, 2012; International Oceanographic Data and Information Exchange, 2016; Corgnati et al., 2018a; Wera Best Practices, 2018; Horstmann et al., 2019). They refer either to DF or BF systems, or to specific manufacturer's radar systems. In this paper, the authors build on the existing documents and rely on direct experience with the aim of offering a comprehensive review of the existing best practices, and complement them with novel recommendations. This paper thus covers planning, setup, operations and data management aspects in order to ensure a broader approach to optimal operation of HFR systems independently from manufacturer or antenna design and setup.
Usually, a network of HFR systems is established in response to specific scientific questions or operational needs. One of the most important criteria in the planning phase of a HFR network is the selection of the desired spatial and temporal resolution, along with maximum range and spatial coverage, that optimize ocean surface current sampling. These parameters will determine the operating frequency, type of system to be used, number of systems and their relative location. Other practical aspects need to be considered, such as available space, infrastructure (power supply, communication, access roads), noise and interference sources (power lines, industry, nearby antennas, metal fences). All these aspects are discussed in the section "New Deployment Requirements and Planning".

A number of factors affect HFR performances during setup and operation with impact on data accuracy or with an interruption of data flow. Environmental factors, such as low sea state or extreme wave conditions, can introduce spatial and temporal data gaps due either to the low-energy environment or the saturation of the spectrum region delimiting the first order Bragg peaks (1st order region); radio interference may contaminate the 1st order region and bias the radial current maps; failures in the network connectivity can disrupt data flow.

While all the HFRs share the same principles of operation, differences in signal-processing suggest different quality assessment procedures and quality control metrics, along with different metadata or descriptors. Keeping in mind such characteristics, a homogeneous approach, capable of identifying and discussing comparable features, is proposed in the section "Setup and Maintenance," dedicated to recommendations for site setup and maintenance to enable continuous operation.

Lastly, recommendations for data management and the software tools currently available for data pre-processing/post-processing, are provided in the "Data Management and Quality Control" section. In addition to data management at local level, a specific example is provided to demonstrate the near real-time Quality Control (QC) application to surface current data, following Corgnati et al. (2017, 2018a), and their operational ingestion in the European HF radar node (Corgnati et al., 2018b) and dissemination to main European and Global marine data portals.

## NEW DEPLOYMENT REQUIREMENTS AND PLANNING

It is important, from both budget and planning perspectives, to consider the efforts required to source site permissions, site access, transmit licenses, and last but not least, intended use of the data.

Complex regulations at the local, national and international level can delay the deployment of the HFR stations. In the worst case scenario, non-compliance with one of the several required permissions may lead to a failure in proceeding with the installation stage. A further important aspect in the planning phase is the determination of the location of each

individual site as well as the site-to-site distance. The optimal location and site-to-site distance depend on the utilized system, frequency, salinity and shape of the coastline. This task becomes particularly difficult when the coastline does not offer easy access or any suitable infrastructures (e.g., buildings, roads, electrical power line). In case parameters such as sea state are required, the distance between sites has to be reduced accordingly (Wyatt et al., 2006, 2007).

## Site Requirements

The optimal candidate site should match the following characteristics:

* located as close as possible to the shoreline but safe from waves and flooding;

* protected from unauthorized human access and from damage caused by animals;

* located in a flat or slightly sloping area allowing human access without hazards;

* accessible by vehicles;

* have enough space to accommodate antennas, electronics, and cables;

* free of electrically conductive objects (e.g., metallic fences, poles, and containers) in the antenna near-field;

* free of radio interference at the operating frequency band;

* free of obstacles limiting the field of view toward the ocean;

* have nearby access to the electrical grid;

* have stable and broadband internet connectivity, either wired or wireless.

Except for compact-type systems that allow co-located receive (Rx) and transmit (Tx) elements, Rx, and Tx antennas should be separated by a minimum distance depending on the operating frequency and the signal processing technique. For instance, optimal separation distance d between Rx and Tx elements for low-frequency SeaSonde systems is one wavelength ($\lambda_{HF}$) at the system's center frequency: $\lambda_{HF} = c/f_{HF}$, with c the speed of light in vacuum (m/s), $f_{HF}$ the center frequency (Hz). For instance, d = 67, 57, 32 m at 4.463, 5.2625, and 9.33 MHz, respectively. Typical setups of BF and DF are depicted in Figure 2.

In all cases, antennas should be installed as close to water as possible (see manufacturer's recommendations in Table 2), to avoid unnecessary signal loss for propagation over land.

For DF systems, electrically conductive objects such as metal structures or power lines within the antenna near field should be avoided, as they would couple in an unpredictable way with the receive elements thus compromising data quality. If these conditions cannot be satisfied, the site should not be considered as suitable for the installation.

Cliffs and steeply sloped ground can also affect HFR measurements, by acting as a reflector of the transmitted radio signal or introducing unwanted vertical propagation modes. Preferences should be given to gently sloping or flat terrain. "Many operating HF radar sites do not meet every one of these criteria and still produce acceptable measurements. When an ideal site is not available, consideration should be given to mitigate existing obstructions" (Cook et al., 2008).
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Upper panel shows a typical setup of a beam forming system. The photos show the antennas placement of a WERA system on the island of Wangerooge at the German coast of the North Sea. The receive antenna array (Rx) is shown on the left-hand side and the transmit array (Tx) on the right hand side, respectively. The bottom panel shows a typical direction finding HF radar system setup for the low frequency (4, 5, and 9 MHz) bands (left hand side). To the right, the SeaSonde installed at Matxitxako Cape (northern coast of Spain).</figcaption></figure>

## Power

Grid power is required at the installation site to feed the control electronics, communication devices, enclosure, and services such as air-conditioning. If grid power is not available, suitable off-grid sources should be used, such as solar panels, wind turbines, diesel generators, or their combination (Statscewich et al., 2011). “Assuming continuous operation, the cost of creating off-grid electricity to power a HF radar site is in the tens of thousands of Euros” (Cook et al., 2008). In continuous operation, the power consumption of most commercial HF radar systems excluding air conditioning is between 300 and 500 W, based on technical specs and direct measurements statistics of SeaSonde and WERA

systems. In addition, an air conditioning unit draining between 500 and 2,000 W of power, depending on the climate of the HF radar site, the location of the electronics and the type of enclosure, is required in most cases.

Power outages are a very common downtime cause for HFRs. As such, it is strongly suggested to ensure that the power line setup comply with the state of the art and that properly designed and dimensioned backup power sources are available.

## Communication

Remote network connection is considered a prerequisite for any HFR installation. Besides being essential for near real time
<table><caption>**TABLE 2** | Recommended maximum distance to the water for SeaSonde and WERA Rx and Tx, based on operating frequency.</caption><thead><tr><th rowspan="2">Frequency (MHz)</th><th colspan="2">Maximum distance to water (m)</th></tr><tr><th>SeaSonde</th><th>WERA</th></tr><tr><th></th><th></th><th>Rx</th><th>TX</th></tr></thead><tbody><tr><td>4-6</td><td>250</td><td>800</td><td>500</td></tr><tr><td>12-14</td><td>150</td><td>500</td><td>300</td></tr><tr><td>24-27</td><td>150</td><td>400</td><td>200</td></tr><tr><td>47-50</td><td>150</td><td>200</td><td>100</td></tr></tbody></table>

operational usage, remote communication with the HFR site via a broadband internet connection allows for redundant data backup and for a series of monitoring and management operations that reduce significantly the need of on-site maintenance.

Best results are achieved with a dedicated ADSL connection. If this option is not available, a mobile network data connection should be considered as a second option. With the fast development of broadband cellular networks, data transfer and remote management of HFR stations are now easy tasks at almost no cost. Two or more SIM data cards from different mobile phone companies can be used simultaneously with specific 4G modem-router devices, ensuring backup link and improved data rate. Industrial grade modem-router are strongly suggested as they provide wider operational range with temperature, better protection against humidity and dust, and some extremely useful software features, for instance the continuous check of the connection status, a watchdog timer for automatic reboot in case of prolonged network disconnection, remote reboot via SMS, and MAC address filtering, amongst other options.

Wireless outdoor bridges can be used to link the remote site to a hardwired network connection if this is located over a distance of kilometers, in case of poor or unavailable mobile network connection at the site. Wireless outdoor bridges are implemented in several ways following the IEEE 802.11 recommendations and in most cases they rely on a point-to-point communication that requires free line-of-sight between two directional outdoor antennas.

Satellite internet should be considered as a potential alternative option for remote areas, although its performance is sensitive to weather conditions. Satellite internet companies should be contacted in advance to see if an intended HF radar site falls within their coverage area. Common satellite internet plans offer enough bandwidth and data volume at reasonable costs, allowing remote management and transfer of the most important data.

At minimum, an internet connection for the HF radar site needs to be able to transfer approximately 300 KB hourly (radial velocities files). In case of extremely slow connection, some HF radar systems offer the option of remote management using command line through SSH and/or control panels over HTTP, both requiring less bandwidth than screen sharing programs or graphical remote desktop access (e.g., VNC, Teamviewer).

Power surge protectors on Ethernet data lines is strongly suggested to protect from lightning strikes and power surges.

A protection should be placed also on the coaxial cable of a 3G/4G modem-router if an outdoor antenna is used.

## Permissions

A significant amount of time should be allocated to obtain the required permits for operating the HF radar site. As a minimum, this implies authorizations from the local Authorities (e.g., City Councils) and the transmit permits from National Agency for the selected frequency with a sufficient bandwidth. The International Telecommunication Union (ITU) Resolution 612 (2012) allocated dedicated frequency bands to oceanographic radars for radiodetermination purposes, with secondary-type licenses on a basis of non-interferences. Each operator is requested to liaise with the national authorities in order to be granted with the operating frequency, the bandwidth, and a unique call sign that identifies each site. Resolution 612 also defines the separation distances between primary and secondary services so to avoid harmful interferences, and requires international agreements to be in place between neighboring Countries if the distance of the HFR site from the border is shorter than a frequency-dependent limit. In this case, in order to prevent big delays on getting an installation approval, coordination with the neighboring Countries should be undertaken well in advance. As an example, an agreement between Italy, France and Spain is in effect that allows to use the full bandwidth (100 kHz) in the North West Mediterranean area for all the SeaSondes or WERAs, under the requisite that the former operate within the 13 MHz band and the latter operate within the 16 MHz band. Mutual cross interference within the same frequency band is avoided thanks to synchronization techniques.

In addition to the general radio transmission license, some Countries enforce a strict limit on the maximum radiated power. In Europe this threshold is set to 10 W (EIRP, effective isotropically radiated power), and an additional permission is required above this value. For experimental purposes (e.g., testing new deployment sites) HFRs can be operated below this level; potential losses in HFR performances can be partially compensated by optimizing the processing parameters.

Based on the general principles of prevention many Countries define minimum health and safety requirements regarding the exposure of people to electromagnetic fields (European Union, 2013). Although the radiated power of an HF radar is very low and the effect on a person is not well known, it is demonstrated by field measurements that within the distance of few meters from the HFR antennas the electric and magnetic fields intensities may exceed such thresholds. For this reason, the presence of humans in the vicinity of the transmitting antennas should be avoided at least before assessing the level of electromagnetic emissions.

In some countries specific constraints may exist in order to preserve the architectural and natural heritage landscapes, as well as to protect sensitive coastal environments. Even limitations due to archeological and ethnographic concerns and constraints have been faced by the authors and should be foreseen in the general deployment plan. Obtaining permissions often involves different governmental offices and requires professional advice in order to get the final approval.
Other possible limitations imposed by authorities are due to the presence of other facilities using nearby radio frequency bands that could be disturbed by the HF radar.

# Additional Equipment

## HF Radar Enclosure and Air Conditioning

The electronics enclosure can be of different nature depending on the available space at the site, cooling, heating and dehumidifying requirements, the number of devices to be hosted, the need of protection against the sun, water, dust.

If a building is available, the enclosure can be located inside a room and a standard rack cabinet can be used. If required, an air conditioning system can be installed in the room.

If a building is not available, the following options can be chosen:

* weatherproof, climate-controlled shelter or trailer. This solution allows the operator to work in a small but still comfortable environment, and provides a robust protection against natural hazards (weather, animals) and vandalisms. A trailer has also the advantage that it can be relocated with less effort.

* sealed, insulated, air-conditioned, enclosure with minimum fitting size for the electronics. Such compact solutions are less protected and may require specific and tailored air conditioning methods but are very flexible, e.g., can be deployed with very little space needed and relocated.

## Data Acquisition

In most commercial HFRs the control of the electronics and data acquisition are performed by x86-based computers ranging from consumer PCs to entry level servers and running Mac OS X or Linux operating systems. They are typically provided in bundle with the HF radar system and already pre-configured with all the needed control and processing software. They only need to be configured for a few parameters such as network settings, site-specific information, processing options. As the lifetime of the computer is typically much shorter than the electronics of the HF radar system, care should be taken with respect to compatibility between the manufacturer's software and newer operating systems.

A redundant external data storage system is recommended on-site as a backup for the data acquired and saved on the computer's disk.

## Power Line Accessories and Uninterruptible Power Supply

Once a suitable electrical power source is established, a remote HFR station may require additional solutions in order to minimize the need of maintenance on site due to power-related issues. They may include:

* a dedicated electrical panel and line, bypassing any pre-existing potentially problematic panels or electrical lines.

* a dedicated electrical grounding if not already existing or if not reliable (a test is highly recommended), and a separate grounding line for the lightning protection system.

* a circuit breaker with automatic reclosing capability, able to restore the power supply if the cause that triggers the breaker is only temporary.

* a smart power strip that can be switched on/off on schedule or by remote control e.g., if a hardware power reset is needed.

An uninterruptible power supply (UPS) should be used to provide near-instantaneous protection and power backup to the HFR system components. UPS acts as surge suppressor and ensures within certain limits stable sine wave (pure or simulated depending on the model) power through over-voltages and brownouts. UPS minimum requirements should include:

* an adequate output rating, that in most cases can be equal to 1 KW;

* an Ethernet card and a website interface for remote configuration and monitoring; enough battery capacity to ensure 15-20 min of autonomy considering the maximum load, in order to properly shut down the sensitive equipment;

* the possibility to expand the battery pack if upgrade is needed;

* two or more outlet groups that can be managed separately.

UPS systems should be considered emergency power backup solutions meant to protect the most sensitive electronics components only.

## Lightning Electromagnetic Pulse Protection

Lightning is a frequent cause of damage in some geographical regions. The need for lightning electromagnetic pulse protection (LEMP) is also mentioned in Cook et al. (2008): "LEMP should be installed inline on any antenna (i.e., receiver and transmitter channels, GPS, communications) as a safety precaution for personnel and radar electronics. Lightning arrestors provide an alternate path to ground during a high voltage surge from lightning strike. There are a variety of designs, but typically the inline gas discharge types are used for RF communications, including HF radar." At least two levels of lightning protection are recommended for any system:

* At the antenna pole, as a protection of the transmission line;

* At the container/room cable inlet, as a protection to the electronic components.

Furthermore, Cook et al. (2008) indicates that different devices may require different specifications for lightning arrestors, for instance "the transmitter requires a lightning arrestor with a higher sparkover voltage than the receiver. Typically, common lightning arrestors (such as the Altelicon AL-NFNFB) come with gas tubes rated for 90 V sparkover voltage. In this case, replacement gas tubes with 350 V sparkover voltage can be purchased."

To mitigate this risk, special attention should be paid also to the design of the cables path. Wrong cables placement could invalidate indeed the LEMP when distances of protected and unprotected lines are too short, as the
overvoltage is transmitted by electromagnetic induction bypassing the lightning arrestor, both on coaxial cables and power line.

## Deployment Budget

While the total cost of the initial deployment is mostly affected by the cost of the HFR hardware, a series of non-negligible additional costs should be allocated in the operational budget.

The total cost of a HFR system and its planning and setup is determined by the following items:

* cost of the initial equipment purchase (only HFR components, e.g., Rx, Tx, computer, antenna, cables);

* cost of additional equipment (e.g., UPS, data storage);

* cost of submitting requests for authorization (design of the installation and related work);

* cost of electromagnetic impact evaluation if required by authorities;

* tax for the license for radio band usage;

* cost for rental of land or building hosting the HRF system;

* purchasing of the enclosure or the shed or both;

* contract with electricity supplier;

* contract for internet connection;

* cost of working hours of professionals and for site preparation work (electrician, bricklayer);

* training courses by the HFR system manufacturer;

* cost of the third-party insurance;

* cost of insurance for HFR station (e.g., damages from natural events or other causes);

* cost of personnel and travels for site survey, installation planning and execution (either if done by the purchasing subject or by the manufacturer or a third company);

* cost of IT infrastructure for data processing and archiving.

After the initial deployment, some of the items listed above can be easily identified as regular operating costs.

## SETUP AND MAINTENANCE

Following (Pearlman et al., 2019), Quality Assurance (QA) is a series of processes and protocols aimed at ensuring that the best product is delivered. Significant examples of QA guidelines applied to the production of HFR data are provided in the present and in the next section, including not only best practices ensuring the optimal condition for data acquisition (e.g., setup and maintenance of the hardware) but also methods for measuring, examining and testing if the quality standards are satisfied, commonly referred to as Quality Control (QC).

### Setup

The setup of the HFR site involves both hardware and software elements. The site selection and description of HFR components are covered in the previous section, together with an overview of the required permissions. Once the design of the HFR station is drafted and all the elements are in place, the setup of the system has to be tailored to the environment and to the intended application.

### Long, Medium, and Short-Range Configuration

The HF radar settings and performances (range and radial resolution) depend upon the allocated frequency and bandwidth, as specified by ITU 612 Resolution (Rev.WRC-12) of the Radio Regulations. The number of allowed frequency bands can be grouped in three blocks with effect on HF radar range: long, medium and short range. Typical range values for surface current measurements according to the manufactures of SeaSonde and WERA systems are listed in Table 3. Table 4 summarizes the theoretical maximum range for surface current measurements with respect to salinity. Values are estimates in% of the optimum range utilizing a propagation model (Gurgel et al., 1999).

### Antenna Mount and Cables

HF radar systems may have different antenna sizes, numbers, designs and setups depending on the operating frequency, implementation concept and signal processing technique (DF or BF). Antennas can be implemented with fiberglass whips, metal poles, wires or combination of these elements held by rigid structures (Figure 2). In general, the following antenna systems are available:

* compact antennas for SeaSonde DF systems;

* small antennas used in array configurations in particular for BF systems;

* the very small active antenna for receive arrays from BF systems.

Already stated in Cook et al. (2008): structural “stability is maximized when the Rx and Tx antennas are mounted in level concrete pads constructed at the HF radar site (without metal rebar to distort the antenna pattern in case of DF systems). Anchors for the Rx and Tx antenna non-conductive guy wires can be incorporated into the concrete pad as well. Cableways should be trenched from the pads to the electronics enclosure to eliminate exposed above ground wiring or placed in protective tubes. The construction permits, soil disturbance, and additional labor this mounting entails may limit its applicability to many HF radar sites.”

HF radar systems deployment should be carefully planned with respect to the length and paths of the coaxial cables to the antennas. In case that extensions are needed, additional cables could be purchased on the free market, however, they should match the electrical characteristics of the original cables and should be selected with great care. In case of long paths (e.g., >100 m) from antennas to Tx and Rx units, a better cable should be adopted in terms of signal attenuation at the given frequency. A continuity test should be always performed after connector installation to ensure proper insulation between shield and central pin.

The cable and connector role is often underestimated, while they represent a crucial component of the system and may compromise the quality of the signal. Since they have a low impact on the total cost of HF radar systems, they should match with high quality standards. Moisture penetration inside the cables is the main reason for efficiency loss in time, and should be prevented under all circumstances by using specialized cables (some plastic sheaths are more effective than others).
<table><caption>**TABLE 3** | HF radar optimum range performance for surface current measurements with respect to the operating frequency for radio bands allocated in Europe (Region 1) for oceanographic radars.</caption><thead><tr><th></th><th>ITU frequency (kHz)</th><th>Optimum range SeaSonde (km)</th><th>Optimum range WERA (km)</th><th>Bandwidth (kHz)</th><th>Range resolution (km)</th><th>Bragg wavelength (m)</th><th>Bragg frequency peak (Hz)</th><th>Doppler Offset at 0.05 m/s (mHz)</th></tr></thead><tbody><tr><td rowspan="3">Long range</td><td>4438</td><td>220</td><td>500</td><td>25</td><td>6</td><td>34</td><td>0.22</td><td>1.5</td></tr><tr><td>4488</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>5250</td><td>175</td><td>400</td><td>25</td><td>6</td><td>28</td><td>0.23</td><td>1.8</td></tr><tr><td rowspan="5">Medium range</td><td>5275</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>9305</td><td>80</td><td>200</td><td>50</td><td>3</td><td>16</td><td>0.31</td><td>3.1</td></tr><tr><td>9355</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>13,450</td><td>60</td><td>110</td><td>100</td><td>1.5</td><td>11</td><td>0.37</td><td>4.5</td></tr><tr><td>13,550</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="5">Short range</td><td>16,100</td><td>60</td><td>90</td><td>100</td><td>1.5</td><td>9</td><td>0.41</td><td>5.4</td></tr><tr><td>16,200</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>24,450</td><td>30</td><td>60</td><td>150</td><td>1</td><td>6</td><td>0.51</td><td>8.2</td></tr><tr><td>24,600</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>26,200</td><td>30</td><td>55</td><td>150</td><td>1</td><td>6</td><td>0.52</td><td>8.8</td></tr><tr><td rowspan="5"></td><td>26,350</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>39,000</td><td>20</td><td>30</td><td>250</td><td>0.6</td><td>4</td><td>0.64</td><td>13.1</td></tr><tr><td>39,500</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>42,000</td><td>15</td><td>25</td><td>325</td><td>0.46</td><td>4</td><td>0.66</td><td>14.1</td></tr><tr><td>42,500</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

For wave retrieval the optimum range reduces by approximately 30%.

Furthermore, great care should be taken by protecting the connectors from direct exposure, using specific greases and self-fusing insulation tape when connecting them. In general, all cables should be handled with care, not being stressed, twisted or narrow bent.

Cables already used in the field should be avoided for new installations. If they are the only option, they must be carefully checked with visual inspection and specific instruments as multimeters and cable fault finders for first insight, and with more advanced instruments for complete testing of electrical characteristics.

## Antenna Tuning

Within the theoretical limits given by the chosen operating frequency, the HF radar performances can be significantly reduced in case of a bad tuning of the transmit antenna, i.e., when significant fraction of the power is reflected back to the amplifier due to a mismatch of the impedance (usually 50 ohms) between the antenna and the amplifier at the given frequency. For this reason, it is always recommended to check each antenna after installation in the field.

<table><caption>**TABLE 4** | Variation of the HF radar range on salinity in percent, for typical transmit frequencies.</caption><thead><tr><th>Salinity</th><th>8</th><th>16</th><th>35</th></tr></thead><tbody><tr><td>8 MHz</td><td>49</td><td>73</td><td>100</td></tr><tr><td>12 MHz</td><td>45</td><td>70</td><td>100</td></tr><tr><td>16 MHz</td><td>44</td><td>69</td><td>100</td></tr><tr><td>20 MHz</td><td>44</td><td>68</td><td>100</td></tr><tr><td>25 MHz</td><td>43</td><td>67</td><td>100</td></tr><tr><td>30 MHz</td><td>43</td><td>67</td><td>100</td></tr></tbody></table>

The amount of forward and reflected power is expressed through the Voltage Standing Wave Ratio (VSWR, or alternatively SWR) and can be measured with a battery driven antenna analyzer without amplifier, or via the HF radar's software through the internal circuit board when this feature is available. A VSWR value below 2 is recommended, while a greater value means that reflected power is coming back to the electronic devices which, besides resulting in lower efficiency, can lead to damage of the electronic equipment. Refer to the manufacturer's manual respect to the tuning requirements and possibilities of the transmit antenna. If no tuning possibilities are available, an antenna tuner can be used to minimize the reflected power.

Note that the SeaSonde transmit antenna is tuned by the manufacturer prior to delivery and no tuning by operators is foreseen. However, if required, an approximate and limited tuning can be performed by changing the length of the antenna whips by trimming, and by adding or removing the capacity top-hat in the low-frequency band, or optimally matching the length of the transmission line.

WERA systems are pre-tuned by the manufacturer but can be further fine-tuned in the field. After the tuning of each antenna element of the transmit and receive array, an internal calibration of the entire systems is needed, which takes into account the antenna tuning and cable length for each element.

For adjusting the level of the power amplifier, the output power of the transmit antenna should fulfill the following requirements:

* to be below the transmit power allowed by the agreement of the local frequency agency.

* Not to saturate the receiver (voltage rms values < 7 V along the receive array)
For FMiCW systems, the blank delay should be used to adjust the saturation in the receiver.

"Regardless of tuning method, it is important to monitor the transmitted and reflected power to diagnose transmitter health and functionality" (Cook et al., 2008).

## Receiving Antenna Calibration

The analysis of the radar signals to resolve the azimuth needs a good knowledge of the receiving antenna specificities. Any HFR system is subject to the influence of the nearby environment, and should be always calibrated once in place (Kohut and Glenn, 2003; Yang et al., 2018), as an important part of the system setup and maintenance, in order to verify the theoretical antenna response and to introduce correction factors if needed. Calibration can be performed with different techniques, among which two are the most common: internal calibration method, commonly applied to BF systems, and far-field calibration method, also known within the HFR community as antenna pattern measurement (APM), usually applied to DF configurations. "While the assumption of an ideal antenna pattern allows to generate maps of radial currents, it does not account for distortions in the real antenna pattern," which often causes inaccuracies in the measurements (Cook et al., 2008). This is particularly important for compact DF antenna like SeaSonde, where the electromagnetic environment can affect phase and gain more significantly than in case of large antenna array.

For BF configurations, the analytical antenna array response function or antenna manifold is used. It is computed by solving the electromagnetic equations for the waves propagating to the array using the precise positions (less than 0.2 m accuracy) of each antenna element. Phase and gain calibration is required and is performed internally. Furthermore, dedicated software for the BF systems can compensate variations of the antenna characteristic caused by environmental or technical conditions (Helzel and Kniephoff, 2010). For WERA systems, all antenna parameters are automatically monitored, at least once per hour and a warning is automatically generated if any parameter reaches a critical value, which might cause the automatic calibration procedure to fail. In such a case a preventive maintenance should be carried out. In some extreme conditions an APM can be helpful to improve the quality of the BF, in particular when being utilized for wave measurement applications.

For DF systems APM is mandatory and can be performed, following a well established procedure, using a boat equipped with a transponder provided by HFR system manufacturer. The boat moves over a circle around the HF radar antennas while the HFR system acquires the transponder signal and the boat's position is tracked with a GPS device.

Alternative methods recently developed, such as the Automated APM via Automatic Identification Systems (AIS) or APM obtained by drones, can be used. The former uses Doppler echoes from vessels cruising in the area, which in many cases lead to a return signal similar to that of a transponder (Emery et al., 2014), and combining it with the position of the vessels obtained, when available, from the information contained in the automatic identification system (AIS). The latter employs small aerial

drones (Washburn et al., 2017) or RC-boats carrying a light radio emitter. Flying drones speed up the APM procedure but in other hands are subject to strict regulations in order to be operated.

## Network Setup

Correct network setup is crucial for data retrieval and remote management, which in turn make the HF radar operations more sustainable as the operators do not need to travel often to the sites for data backup and system monitoring. With any kind of network connection at the HF radar site (Ethernet, mobile, satellite etc.) the recommended network configuration should rely on the usage of a router under which a Local Area Network (LAN) can be established.

Proper setup of a LAN is out of the scope of this paper, information on internet services running at the HF radar sites or suggestions on Remote Desktop Management Software can be found in the HFRs operation manuals.

## First Order Limits Detection

When setting up a new HF radar system, operators are requested to insert specific configuration parameters such as the identification of the system, the geographical information, the operating frequency and bandwidth, the processing options. One important step in HF radars signal processing chain is the detection, inside the power spectra maps, of the regions containing the first order bragg peaks, i.e., the portion of the spectra containing information on surface currents. Such detection, better known as first order limits (FOLs) determination, is usually performed by the HF radar acquisition and processing software with automatic algorithms that look for the Doppler peaks at one channel (usually the omnidirectional vertical loop for SeaSonde systems) and must be tuned by the operator during the initial period of acquisition through the selection of some parameters and visual inspection of the results.

Improper FOLs identification may introduce large errors or missing values of sea current velocity. Depending on sea state, interferences, presence of other targets (e.g., ships) in the field of view of the radar, the first order region may be partially masked and difficult to detect. This is particularly true in those areas where a high variability in currents and waves regimes leads to an overlap of first and second order. Bragg peaks with high variability in time and along the HF radar range. In these circumstances, operators should carefully monitor the performance of their initial setup and optimize it for different sea conditions. When this is not possible, different subset of data should be reprocessed each one with a specific set of processing parameters.

These considerations apply especially for DF systems, where the power spectra contain the contributions from all the azimuthal directions. In some cases the application of advanced methods for FOLs detection can improve the quality of the data (Kirincich, 2017).

## Maintenance

As all in situ instrumentation, also HF radar suffers from the normal deterioration, in particular - but not only - of
the outdoor components. Besides that, major damages have been experienced due to severe weather events inducing storm surge or lightning, which often lead to electrical damages or antennae breakage. Another not to be underestimated source of damage results from animals (in particular to cables) or human vandalism.

On the other hand, issues may arise from specific hardware or software failures without any external cause, such as data corruption due to storage failure, or software modules and applications unexpectedly freezing. For all these reasons, even if HF radar systems can be considered fully automated acquisition platforms, regular monitoring is necessary to ensure continuous and correct operation.

The first, and more frequent, kind of maintenance can be performed remotely. With an internet connection enabled, the system status as well as its data acquisition process and data themselves can be easily diagnosed and verified by means of web tools, email alerts, screen sharing applications, even on a daily basis.

As a complement, on-site inspection is unavoidable and is recommended on a regular biannual basis, in order to confirm the good status of the equipment and prevent issues, but also to perform scheduled actions such as data backup on external disk, UPS battery replacement. Additional specific on-site inspections are recommended after a severe weather event (e.g., very high winds, hail, ice storms, floods, lightning).

According to Cook et al. (2008), the following actions should be performed regularly as typical on-site maintenance:

* for DF systems and directional antennas, ensure that RX antenna compass bearing is reported correctly and no changes occurred;

* check if antenna mounting and masts are properly secured and placed, and guy wires are intact and in tension;

* inspect electronics, connectors and other metal components for signs of corrosion or humidity;

* check if conduit or even cables are damaged;

* ensure that air conditioning system is operating properly, filters are clean, all the cooling fans are rotating correctly (no friction or vibrations), components inside the enclosure are not overheated;

* perform tests on UPS for battery status and for correct shutdown and reboot sequence;

* verify the correspondence between on-site system status and remote diagnostics;

* execute an incremental data backup on a portable device.

Other specific items to be checked will be dictated by the diagnostic tools output.

Results of remote and on-site checkings should be included in periodic reports, to help the operator to keep track of maintenance history.

## Diagnostic Reporting

Once the HF radar station is deployed, its proper functioning can be assessed by the diagnostic reporting on its performance. HF radar systems, as discussed in the previous sections in this paper,

can be implemented with different hardware and software design, and different signal processing techniques. For this reason, they can provide different performance indicators. However, all HF radar systems allow recording common operating parameters like temperature and voltage of subcomponents (e.g., chassis, amplifier), forward and reflected power, and indicators related to Signal to Noise ratio, and provide alert messages when measured values exceed given thresholds. This aspect must be considered carefully, thresholds must be set and software properly configured according to manufacturer's manuals, and every change in those indicators, especially if permanent, should be investigated.

As an example, any significant permanent variation of the reflected power is most likely related to a hardware failure (e.g., antennae, cables, or amplifiers). Signal to Noise Ratio is usually not stable because of environmental conditions and periodically fluctuates with repeated daily pattern, but any significant permanent change in this pattern or in its variation range should be investigated.

HF radar management software often provides diagnostic tools and warning messages.

The SeaSonde Radial Suite generates hardware and software diagnostics that are saved in so called DIAG files. An extension with \*.hdt refers to hardware diagnostics and an extension with \*.rdt refers to radial diagnostics. Much of this information is reported within the radial file itself, but diagnostic files help to aggregate and show them, through DiagDisplay application, in a bigger picture. Plots of diagnostic parameters over custom time windows can be also exposed in a web page served by the RadialWebServer. Another application, the RadialSiteReporter, is able to perform a scheduled detailed check of all the software and hardware components and to produce and send to the operator detailed email alerts with highly customizable rules. These alerts are also shown by the RadialWebServer, plus are added to an alert log.

The WERA software provides hardware and software diagnostics on all levels of data acquisition or data processing. Automatic status and warning messages (\*.status) will automatically be sent out to defined recipient groups at different intervals depending on defined priority levels. High priority messages will be sent out immediately, lower priority messages only once per day or per week. An anti-spam filter makes sure that no messages within some hours are send out more than once. In case of a detected problem, the message will contain possible solutions for the problem. Plots of time series files for different kinds of parameters, like internal voltages, temperatures, forward and reverse transmit power, hard disk usage, antenna performance, can be accessed by the operator for further analysis and troubleshooting. The described status monitoring system can simply be extended by creating additional \*.status files in the defined data format.

In Tables 5, 6 the recommended ranges of diagnostic parameters indicating good health of the system are listed for Codar and WERA HFRs respectively. Alerts should be thrown when these ranges are exceeded because of risk of hardware or operation failures.
<table><caption>TABLE 5 | Recommended ranges for diagnostic parameters for Codar SeaSonde HFR systems (Emery and Washburn, 2007).</caption><thead><tr><th>Parameter</th><th>Code</th><th>Mean or (accepted values)</th><th>Standard deviation</th></tr></thead><tbody><tr><td colspan="4"><strong>Standard seasonde hardware diagnostic statistics</strong></td></tr><tr><td>Receiver chassis temp (°C)</td><td>RTMP</td><td>&lt;40<sup>a</sup> 26.2<sup>b</sup></td><td>6.0<sup>a</sup> 5.9<sup>b</sup></td></tr><tr><td>AWG board temp (°C)</td><td>MTMP</td><td>&lt;50<sup>a</sup> 36.5<sup>b</sup></td><td>7.0<sup>a</sup> 6.6<sup>b</sup></td></tr><tr><td>MTMP-RTMP (°C)</td><td></td><td>&lt;12<sup>b</sup></td><td></td></tr><tr><td>Receiver + VDC supply</td><td>SP05</td><td>min ≥ 4.5; max ≤ 6<sup>a</sup> 5.2<sup>b</sup></td><td>0.1<sup>b</sup></td></tr><tr><td>Receiver – VDC supply</td><td>SN05</td><td>max ≤ –4.5; min ≥ –6<sup>a</sup> –5.1<sup>b</sup></td><td>0.2<sup>b</sup></td></tr><tr><td>Receiver + 12 VDC supply</td><td>SP12</td><td>min ≥ 10; max ≤ = 14<sup>a</sup> 12.3<sup>b</sup></td><td>0.1<sup>b</sup></td></tr><tr><td>Transmitter + VDC supply</td><td>XP05</td><td>min ≥ 4.5; max ≤ 6<sup>a</sup></td><td></td></tr><tr><td>Transmitter + 28 V Supply</td><td>XP28</td><td>min ≥ 23; max ≤ 35<sup>a</sup></td><td></td></tr><tr><td>Transmitter front panel chassis temp (°C)</td><td>XPHT</td><td>max ≤ 40<sup>a</sup> 28.9<sup>b</sup></td><td>5.2<sup>b</sup></td></tr><tr><td>Transmitter amplifier temp (°C)</td><td>XAHT</td><td>max ≤ 50<sup>a</sup> 34.5<sup>b</sup></td><td>5.0<sup>b</sup></td></tr><tr><td>AWG module temperature (°C)</td><td>AHOT</td><td>max ≤ 50<sup>a</sup></td><td></td></tr><tr><td>Receiver front panel temperature (°C)</td><td>RHOT</td><td>max ≤ 40<sup>a</sup></td><td></td></tr><tr><td>XAHT-XPHT (°C)</td><td></td><td>&gt;10<sup>b</sup></td><td></td></tr><tr><td>Transmitter forward power (W)</td><td>XAFW</td><td>53.0<sup>a</sup></td><td>13.0<sup>a</sup></td></tr><tr><td>Transmitter reflected power (W)</td><td>XARW</td><td>5.0<sup>a</sup></td><td>5.0<sup>a</sup></td></tr><tr><td>Channel 1 signal-to-noise ratio (dB)</td><td>SSN1</td><td>min ≥ 20<sup>a</sup></td><td></td></tr><tr><td>Channel 2 signal-to-noise ratio (dB)</td><td>SSN2</td><td>min ≥ 20<sup>a</sup></td><td></td></tr><tr><td>Channel 3 signal-to-noise ratio (dB)</td><td>SSN3</td><td>min ≥ 20<sup>a</sup></td><td></td></tr><tr><td>Number of GPS satellites<sup>c</sup></td><td>NSAT</td><td>min ≥ 4<sup>a</sup></td><td></td></tr><tr><td>Phase lock loop loss</td><td>PLLL</td><td>0.2<sup>b</sup></td><td>2.42</td></tr><tr><td>Run time (hours)</td><td>RUNT</td><td>190<sup>a</sup> 188<sup>b</sup></td><td>400<sup>a</sup> 402<sup>b</sup></td></tr><tr><td>Number of radial vector –13 MHz band</td><td>RADV</td><td>520<sup>a</sup> 523<sup>b</sup></td><td>310<sup>a</sup> 309<sup>b</sup></td></tr><tr><td>Number of radial vector –25 MHz band</td><td>RADV</td><td>280<sup>a</sup> 276<sup>b</sup></td><td>110<sup>a</sup> 106<sup>b</sup></td></tr><tr><td>Number of radial vector –40 MHz band</td><td>RADV</td><td>960<sup>a</sup> 960<sup>b</sup></td><td>470<sup>a</sup> 469<sup>b</sup></td></tr><tr><td>Average # Solns per range cell –13 MHz</td><td>RAPR</td><td>15<sup>a</sup></td><td>5<sup>a</sup></td></tr><tr><td>Average # Solns per range cell –25 MHz</td><td>RAPR</td><td>10<sup>a</sup> 9<sup>b</sup></td><td>5<sup>a</sup> 4<sup>b</sup></td></tr><tr><td>Average # Solns per range cell –40 MHz</td><td>RAPR</td><td>30<sup>a</sup> 33<sup>b</sup></td><td>15<sup>a</sup> 15<sup>b</sup></td></tr><tr><td>Maximum radial range (km) –13 MHz band</td><td>RADR</td><td>75.7<sup>b</sup></td><td>18.7<sup>b</sup></td></tr><tr><td>Maximum radial range (km) –25 MHz band</td><td>RADR</td><td>29.4<sup>b</sup></td><td>5.2<sup>b</sup></td></tr><tr><td>Maximum radial range (km) 40 MHz band</td><td>RADR</td><td>10.6<sup>b</sup></td><td>2.1<sup>b</sup></td></tr><tr><td>Average bearing of all radials –13 MHz band</td><td>RABA</td><td>–88.4<sup>b</sup></td><td>80.0<sup>a</sup> 79.0<sup>b</sup></td></tr><tr><td>Average bearing of all radials –25 MHz band</td><td>RABA</td><td>–91.8<sup>b</sup></td><td>80.0<sup>a</sup> 84.8<sup>b</sup></td></tr><tr><td>Average bearing of all radials –40 MHz band</td><td>RABA</td><td>–5.7<sup>b</sup></td><td>70.0<sup>a</sup> 72.5<sup>b</sup></td></tr></tbody></table>

## DATA MANAGEMENT AND QUALITY CONTROL

The main objective of this section is to provide general concepts and definitions in order to formalize the data processing and management steps regardless of the HF Radar system adopted, with focus on Near Real Time data flow suitable for operational services. HF radar data is *in situ* gridded data sampled at high-frequency that has to be managed according to its peculiarity and complexity, which derives from the fact that it includes diverse data types (radials and totals) and with different kinds of native formats from different HF radar systems.

Different steps in the data life cycle (as schematized in Figure 3) have been defined following the processing levels specified in Table 7.

### Data Acquisition and Pre-processing

The first step consists of raw electromagnetic data acquisition and extraction of radial components of the ocean currents performed by each HFR system.

The output, the so-called radial files, must be transferred to a processing station responsible for combination of overlapping radial vectors measured by two or more HF radar systems. The processing station is a generic name for the IT infrastructure that will perform the combination, either at the data-provider level

<footnote><sup>a</sup>Recommendations given by remote site monitoring scripts rs_WARN.pl. This script compares the parameters of the hardware against accepted values as listed on the tables. <sup>b</sup>Recommendations given by the documentation (http://hfradar.msi.ucsb.edu/brian_emery/files/reports/2008_diagnostics_evaluation.pdf). <sup>c</sup>Warning: the new GPS modules (MINI-JLT REV v1.0) do not provide information about the number of GPS satellites. The rs_WARN.pl must be modified to include this exception.</footnote>
<table><caption>**TABLE 6** | Recommended ranges for diagnostic parameters for WERA HFR systems as reported by the manufacturer manual.</caption><thead><tr><th>Parameter</th><th>Nominal value/Limit</th></tr></thead><tbody><tr><td colspan="2"><strong>Standard WERA IV hardware diagnostic parameters</strong></td></tr><tr><td>Frequency Control Rack +3.3 V</td><td>3.1-3.5 V</td></tr><tr><td>Frequency control rack +5 V</td><td>4.9-5.15 V</td></tr><tr><td>Frequency control rack +12 V</td><td>11.2-12.8 V</td></tr><tr><td>Frequency control rack -12 V</td><td>-12.8 to -11.2 V</td></tr><tr><td>Frequency control rack +27 V</td><td>26-28 V</td></tr><tr><td>Frequency control rack +5 V</td><td>4.9-5.15 V</td></tr><tr><td>Frequency control rack +12 V</td><td>11.2-12.8 V</td></tr><tr><td>Frequency control rack +18 V</td><td>17.2-18.8 V</td></tr><tr><td>Frequency control rack -21 V</td><td>-22 to -20 V</td></tr><tr><td>Frequency control rack temperature</td><td>Adjustable (40/45/50 degrees)</td></tr><tr><td>Receiver rack 1/2 +20 V</td><td>19-21 V</td></tr><tr><td>Receiver rack 1/2 +21 V</td><td>20-22 V</td></tr><tr><td>Receiver rack 1/2 -21 V</td><td>-22 to -20 V</td></tr><tr><td>Receiver rack 1/2 temperature</td><td>Adjustable (40/45/50 degrees)</td></tr><tr><td>Power amplifier +26.5 V</td><td>25.5-27.5 V</td></tr><tr><td>Power amplifier +15 V</td><td>14.2-15.8 V</td></tr><tr><td>Power amplifier +15 V</td><td>14.2-15.8 V</td></tr><tr><td>Power amplifier forward power</td><td>(Depends on power amplifier gain setting)</td></tr><tr><td>Power amplifier reverse power</td><td>&lt;4 watts</td></tr><tr><td>Power amplifier temperature</td><td>Adjustable (~70°)</td></tr><tr><td>Power amplifier self protection (-6 dB gain)</td><td>Yes/No</td></tr><tr><td>User PC hard disk usage</td><td>85%/90%/95% (warning levels)</td></tr><tr><td>Receive antenna gain value</td><td>Gain factor &gt; 10</td></tr><tr><td>DPT<sup>a</sup> – receive channel (1…16) overdriven</td><td>&gt;7 V RMS</td></tr><tr><td>DPT<sup>a</sup> – receive channel (1…16) calibration value I/Q</td><td>&lt;10</td></tr><tr><td>DPT<sup>a</sup> – receive channel (1…16) calibration value balance</td><td>0.9 &lt; x &lt; 1.1</td></tr><tr><td>DPT<sup>a</sup> – receive channel (1…16) I/Q balance</td><td>90 ± 3°</td></tr><tr><td>DPT<sup>a</sup> – receive channel (1…16) time series of phase value</td><td>Almost constant over time</td></tr><tr><td>DPT<sup>a</sup> – receive channel (1…16) time series of RMS voltage values</td><td>(Similar for all receive channels)</td></tr><tr><td>Internal calibration measurement values</td><td>All within limits</td></tr><tr><td>Low level WERA hardware monitoring (Many automatic checks during boot and for each measurement.)</td><td>All basic tests OK</td></tr></tbody></table>

or as third-party competence center (e.g., manufacturer servers, National nodes and data centers, EU HF radar Node). According to the definition of data levels, this flow is summarized in the following scheme:

**Input product** Level 0 (Signal received by the antenna before the processing stage)

**Output product** Level 2A (HF radar radial velocity data without QC defined)

**Data source** HF radar remote stations

**Data target** Central or Processing Station for combination into totals

HF radar manufacturers provide software suites for data acquisition and radial velocities production. The two most common commercial software packages for HF Radar acquisition and pre-processing are:

* the suite of WERA/Helzel (available for Linux OS only, suggested distribution open SUSE), based, for the signal processing algorithms, on Klaus Werner Gurgel works (Gurgel et al., 1999), and supporting the two HF radar implementations and spectral information processing methods: Beam Forming processed on a regular cartesian grid (or optionally with polar coordinates), and Direction Finding by Least Mean Square developed for compact antenna system.

* the Seasonde radial suite (available for Mac OS X only) from Codar Ocean Sensors, based, for the signal processing algorithms, on the works of Lipa and Barrick (1983) and Barrick and Lipa (1997).

Third party software for data acquisition and preprocessing is available under commercial license for beam forming HF radar systems, see **Table 8** for more details.

Operators are invited to investigate the license conditions, especially the policy for software updates, which may differ significantly, in order to allocate the correct budget for maintaining the systems.

## Data Processing

The combination of radial velocities from two or more remote stations for obtaining total surface current velocity files is covered in this paragraph. According to the definition of data levels, this flow is summarized in the following scheme:

**Input product** Level 2A (HF radar radial velocity data without QC defined)

**Output product** Level 3A (HF radar total velocity data without QC defined)

**Data source** Central or Processing Station

**Data target** Post-processing Station

The most commonly adopted combination algorithm of radial vectors into total vectors is the unweighted least squares fitting (UWLS) algorithm. The UWLS approach (Lipa and Barrick, 1983; Graber et al., 1997) assumes that, for each grid point, the radial velocities within the search radius are produced by a uniform velocity vector, i.e., the correlation of the current vector is assumed to be one everywhere within the search radius and zero outside.

HF radar manufacturers provide combining suite software under commercial license (Table 9) often as separate product with respect to the acquisition and pre-processing software described in the previous section. Besides those solutions, the

<footnote><sup>a</sup>DPT, Direct Path Test measurement.</footnote>
<figure><img src="image_5.png" /><figcaption>FIGURE 3 | Data life cycle of HF radar data from acquisition to distribution.</figcaption></figure>

international community developed and freely distributes the Matlab® tool “HFR_p Programs”<fnref n="1" /> (Kaplan et al., 2005).

## Data Post-processing

Once both radial and total files are produced, it may be convenient or required to post-process them in order to achieve a common level of Quality Control (QC), or a specific file format, or a required variable/attribute naming convention.

Best practices suggest that the NetCDF should be adopted as file format either for radial and total vectors files. This standard is indeed adopted by the US IOOS network, EU community and the Integrated Marine Observing System (IMOS) Ocean Radar Facility Ocean radar Network (ACORN), at least for total vectors.

Mandatory, suggested or ancillary variables and attributes and their naming conventions are a matter of discussion. Different conventions/directives are adopted by different communities, however a coordination work is carried on under the umbrella of the Global HF Radar network (Roarty et al., 2019).

As a data post-processing example, the authors want to report here the case of the application of QC tests and their description, according to the definition of data levels and to the European common data and metadata model for surface currents

(Corgnati et al., 2017, 2018a,b). The flow is summarized in the following scheme:

<table><thead><tr><th></th><th></th></tr></thead><tbody><tr><td>Input product</td><td>Level 2A (HF radar radial velocity data without QC defined) and Level 3A (HF radar total velocity data without QC defined)</td></tr><tr><td>Output product</td><td>Level 2B (HF radar radial velocity data with a minimum set of QC) and Level 3B (HF radar total velocity data with a minimum set of QC)</td></tr><tr><td>Data source</td><td>Post-processing Station</td></tr><tr><td>Data target</td><td>Distribution centers</td></tr></tbody></table>

As described by Chapman et al. (1997); Barrick (2002), Kohut and Glenn (2003); Cosoli and Bolzon (2010), Lipa (2013); Forget (2015) and Cosoli and Grcic (2019), even if we assume that the radar hardware is operating correctly, many different sources of uncertainty in the radial velocities can be identified. It is thus highly recommended to provide information about the validity and correctness of HFR measurements via Quality Control (QC) procedures.

The authors selected from the Quality Assurance/Quality Control of Real-Time Oceanographic Data (QARTOD) manual (U.S. Integrated Ocean Observing System, 2016) a battery of QC

<footnote n="1">¹https://github.com/rowg/hfrprogs</footnote>
<table><caption>**TABLE 7** | Processing levels for HF radar data.</caption><thead><tr><th>Processing level</th><th>Definition</th><th>Products</th></tr></thead><tbody><tr><td>Level 0</td><td>Reconstructed, unprocessed instrument/payload data at full resolution; any and all communications artifacts, e.g., synchronization frames, communications headers, duplicate data removed</td><td>Signal received by the antenna before the processing stage. (No access to these data in Codar systems)</td></tr><tr><td>Level 1A</td><td>Reconstructed, unprocessed instrument data at full resolution, time-referenced and annotated with ancillary information, including radiometric and geometric calibration coefficients and georeferencing</td><td>Spectra by antenna channel</td></tr><tr><td>Level 1B</td><td>Level 1A data that have been processed to sensor units for next processing steps. Not all instruments will have data equivalent to Level 1B</td><td>Spectra by beam direction</td></tr><tr><td>Level 2A</td><td>Derived geophysical variables at the same resolution and locations as the Level 1 source data</td><td>HF radar radial velocity data</td></tr><tr><td>Level 2B</td><td>Level 2A data that have been (re)processed with a minimum set of QC</td><td>HF radar radial velocity data</td></tr><tr><td>Level 2C</td><td>Level 2A data that have been reprocessed for advanced QC</td><td>Reprocessed HF radar radial velocity data</td></tr><tr><td>Level 3A</td><td>Variables mapped on uniform space-time grid scales, usually with some completeness and consistency</td><td>HF radar total velocity data</td></tr><tr><td>Level 3B</td><td>Level 3A data that have been (re)processed with a minimum set of QC</td><td>HF radar total velocity data</td></tr><tr><td>Level 3C</td><td>Level 3A data that have been reprocessed for advanced QC</td><td>Reprocessed HF radar total velocity data</td></tr><tr><td>Level 4</td><td>Model output or results from analyses of lower level data, e.g., variables derived from multiple measurements</td><td>Energy density maps, residence times, etc.</td></tr></tbody></table>

tests to be always applied to HFR radial and total data. These tests, listed as mandatory in the European common QC model for near real-time HF radar current data and textually reported here from Corgnati et al. (2018a), are:

## RADIAL DATA TESTS

* **Syntax check:** this test will ensure the proper formatting and the existence of all the necessary fields within the radial NetCDF file.

* **Over-water test:** this test labels radial vectors that lie on land with a “bad data” flag and radial vectors that lie on water with a “good data” flag.

* **Velocity Threshold:** this test labels radial velocity vectors whose module is bigger than a maximum velocity threshold with a “bad data” flag and radial vectors whose module is smaller than the threshold with a “good data” flag.

* **Variance Threshold:** this test labels radial vectors whose temporal variance is bigger than a maximum threshold with a “bad data” flag and radial vectors whose temporal variance is smaller than the threshold with a “good data” flag. The 2013 Codar Current Newsletter suggests not to use variance data for real-time QC, due to the fact that the CODAR parameter defining the variance is

<table><caption>**TABLE 8** | Software for data acquisition and pre-processing.</caption><thead><tr><th>Software/HF radar type</th><th>Features</th><th>Requirements, license and links</th></tr></thead><tbody><tr><td>Standard WERA® Data Manager and Viewer/WERA HF radar either in Beam Forming or Direction Finding version</td><td>Data processing, real time quality flagging, data archive, plotting tools, system monitoring</td><td>Running on Linux operating system under commercial license, distributed in bundle with the purchase of each WERA HF radar system. Further information on application software packages: https://helzel-mestechnik.de/de/11190-application-software</td></tr><tr><td>Codar SeaSonde Radial Suite/Codar SeaSonde HF radar (direction finding)</td><td>Real time processing, quality flagging, data archive, plotting tools, system monitoring</td><td>Running on Mac Os X operating system under commercial license, distributed in bundle with the purchase of each Codar SeaSonde HF radar system. Further information: http://www.codar.com/Manuals/Suite_Radial_Software.pdf</td></tr><tr><td>Seaview real-time software for WERA (Beam Forming only), PISCES, OSCR HFR systems</td><td>An integrated wave, current and wind measurement package, with data viewer, web access facilities, radar auto-diagnostics, and environmental warning systems</td><td>Running on Solaris-Sparc and linux-x86 operating systems under commercial license. Further information: http://www.seaviewsensing.com/software.html</td></tr></tbody></table>

<table><caption>**TABLE 9** | Standard software for data processing.</caption><thead><tr><th>Software/HF radar type</th><th>Features</th><th>Requirements, license and links</th></tr></thead><tbody><tr><td>Codar SeaSonde Combine Suite/Codar SeaSonde HF radar (direction finding)</td><td>Real-time processing applications and tools for radial vector combination, data visualization and archiving</td><td>Running on Mac Os X operating system under commercial license. Further information: http://www.codar.com/Manuals/Suite_Combine_Software.pdf</td></tr><tr><td>HFR_Progs-2.1.2</td><td>Set of Matlab® functions and scripts performing total currents generation, Open-Boundary Modal Analysis (OMA), interpolation, filtering, tides, EOFs</td><td>Os independent, Matlab® license is required. The code depends also on several freely available software packages: m_map, openMA, t_tide, arrow, mexnc, gshhs shoreline database. https://github.com/rowg/hfrprocs</td></tr><tr><td>Standard WERA® Data Manager and Viewer/WERA HF radar either in Beam Forming or Direction Finding configuration</td><td>Data processing, real time quality flagging, data archive, plotting tools, system monitoring</td><td>Running on Linux operating system under commercial license. Further information on application software packages: https://helzel-mestechnik.de/de/11190-application-software</td></tr></tbody></table>

computed at each time step, and therefore considered not statistically solid. Therefore, this test is applicable only to Beam Forming (BF) systems. Data files from Direction Finding (DF) systems will apply instead the “Temporal Derivative” test reporting the explanation
"Test not applicable to Direction Finding systems. The Temporal Derivative test is applied" in the comment attribute.

* **Temporal Derivative:** for each radial bin, the current hour velocity vector is compared with the previous and next hour ones. If the differences are bigger than a threshold (specific for each radial bin and evaluated on the basis of the analysis of one-year-long time series), the present vector is flagged as bad_data, otherwise it is labeled with a good_data flag. Since this method implies a one-hour delay in the data provision, the current hour file should have the related QC flag set to 0 (no QC performed) until it is updated to the proper values when the next hour file is generated.

* **Median Filter** (text reported from U.S. Integrated Ocean Observing System, 2016): for each source vector, the median of all velocities within a radius of <RCLim> and whose vector bearing (angle of arrival at site) is also within an angular distance of <AngLim> degrees from the source vector's bearing is evaluated. If the difference between the vector's velocity and the median velocity is greater than a threshold, then the vector is labeled with a "bad_data" flag, otherwise it is labeled with a "good_data" flag.

* **Average Radial Bearing:** this test labels the entire datafile with a 'good_data' flag if the average radial bearing of all the vectors contained in the data file lies within a specified margin around the expected value of normal operation. Otherwise, the data file is labeled with a "bad_data" flag. This test is applicable only to DF systems. Of course, data files from BF systems will have this variable filled with "good_data" flags (1) and the explanation "Test not applicable to Beam Forming systems" in the comment attribute.

* **Radial count:** test labeling radial data having a number of velocity vectors bigger than the threshold with a "good data" flag and radial data having a number of velocity vectors smaller than the threshold with a "bad data" flag.

## TOTAL DATA TESTS

* **Syntax check:** this test will ensure the proper formatting and the existence of all the necessary fields within the total NetCDF file.

* **Data density threshold:** this test labels total velocity vectors with a number of contributing radials bigger than the threshold with a "good data" flag and total velocity vectors with a number of contributing radials smaller than the threshold with a "bad data" flag.

* **Velocity threshold:** this test labels total velocity vectors whose module is bigger than a maximum velocity threshold with a "bad data" flag and total vectors whose module is smaller than the threshold with a "good data" flag.

* **Variance threshold:** this test labels total vectors whose temporal variance is bigger than a maximum threshold with a "bad data" flag and total vectors whose temporal variance is smaller than the threshold with a "good data" flag. As stated for radial data, this test is applicable only to Beam Forming (BF) systems. Data files from Direction Finding (DF) systems will apply instead the "Temporal Derivative"

test reporting the explanation "Test not applicable to Direction Finding systems. The Temporal Derivative test is applied." in the comment attribute.

* **Temporal derivative:** for each total bin, the current hour velocity vector is compared with the previous and next hour ones. If the differences are bigger than a threshold (specific for each grid cell and evaluated on the basis of the analysis of one-year-long time series), the present vector is flagged as "bad_data," otherwise it is labeled with a "good_data" flag. Since this method implies a one-hour delay in the data provision, the current hour file should

<table><caption>TABLE 10 | Software for data post-processing for European common HFR data and metadata model.</caption><thead><tr><th>Tool name</th><th>Matlab® (HFR_Combiner)</th><th>Java™ (JRadar)</th></tr></thead><tbody><tr><td>Tool description</td><td>The HFR_Combiner tool is a bundle of scripts developed with Matlab® to perform the tasks of transforming Codar and WERA total current into the European HFR Standard. The tool also allows for combining Codar radial current files into total ones and producing both radial and total data into the European HFR Standard.</td><td>The J(ava)Radar tool is an executable jar file already compiled and packaged ready to work.</td></tr><tr><td>Tool available on GitHub online</td><td>YES<br>https://github.com/LorenzoCorgnati/HFR Node_tools</td><td>YES<br>https://github.com/Ilasensio/JRadar</td></tr><tr><td>Additional required software and packages (must be pre-installed by attendees)</td><td>Valid and licensed Matlab® version, later than R2016b release. Previous releases can be used, but it is recommended to have at least the R2016b version.<ul><li>HFR_Progs-2.1.2</li><li>M_Map</li><li>GSHHS</li><li>Nctoolbox-1.1.3 mysql-connector-java-5.1.17 JDBC driver for connecting to the MySQL database</li><li>Rdir</li></ul></td><td>Java™ 8 or higher runtime environment installed with administrator permissions on the computer.</td></tr><tr><td>Graphical user interface</td><td>NO</td><td>YES</td></tr><tr><td>Converts CODAR TUV files (Totals) to NetCDF European HFR Standard</td><td>YES</td><td>YES</td></tr><tr><td>Performs radial combination</td><td>YES</td><td>NO</td></tr><tr><td>Performs QC</td><td>YES</td><td>YES</td></tr><tr><td>Can run in real time</td><td>YES</td><td>YES</td></tr><tr><td>Requires scripting skills</td><td>YES</td><td>NO</td></tr><tr><td>Allows batch processing</td><td>YES</td><td>YES</td></tr></tbody></table>
have the related QC flag set to 0 (no QC performed) until it is updated to the proper values when the next hour file is generated.

* **GDOP threshold:** this test labels total velocity vectors whose GDOP is bigger than a maximum threshold with a "bad data" flag and the vectors whose GDOP is smaller than the threshold with a "good data" flag.

These mandatory QC tests are manufacturer-independent, i.e., they do not rely on particular variables or information provided only by a specific device.

Each QC test will result in a flag related to each data vector which will be inserted in the specific test variable. These variables can be matrices with the same dimensions of the data variable, containing, for each cell, the flag related to the vector lying in that cell, in case the QC test evaluates each cell of the gridded data, or a scalar, in case the QC test assesses an overall property of the data.

The ARGO QC flag scale is applied as flagging scheme (Corgnati et al., 2018a).

For some of these tests, HF radar operators will need to select the best thresholds. Since a successful QC effort is highly dependent upon selection of the proper thresholds, this choice is not straightforward, and may require trial and error before final selections are made. These thresholds should not be determined arbitrarily but based on historical knowledge or statistics derived from historical data.

The software tools for HF radar Real Time data post-processing into the EU standard, described in Table 10, have been developed by the EU HFR Node and they are freely shared (provided under the license: Attribution-Non-Commercial-ShareAlike 4.0 International -CC BY-NC-SA 4.0). More information about the license in this link: http://creativecommons.org/licenses/by-nc-sa/4.0/) and kept updated on a GitHub repository.

## Data Archival and Preservation

On-site archiving on a separate storage system is recommended as the first step for data preservation. Both the Codar SeaSonde and WERA softwares include scripts and applications for scheduled data archiving.

For Codar systems, at the very least, it is advisable to save all Range Series (.rs) data files, from which radial diagnostic information, Cross Spectra (.cs) files, and all the derived oceanographic variables can be recalculated, allowing also for different reprocessing settings. From a practical point of view, in order to ensure quicker checks of processing settings, it is suggested to save also the Cross Spectra files (at least CSS-type averaged files). Another bunch of deployment-specific files must be archived for SeaSonde systems, they are: in the RadialConfigs folder, Header.txt file and MeasPattern.txt files; gps Track data and Time Series files generated during APM, needed for recalculating Antenna Pattern in case of need; the folder/Codar/SeaSonde/Logs/or any equivalent log file containing the operator comments on all the changes performed on the radial site, including periods of validity of one or more Antenna Pattern files; alert and diagnostic files for later troubleshooting of the station behavior.

The WERA toolbox contains a script named "BackupDataFiles.sh," which is used to automatically create a backup copy of raw data and other data products onto an external backup hard disk. If two hard disks are used for redundancy, the script should be duplicated for optimum performance.

For WERA systems using beam forming software it is not required to save Cross Spectra files. The *.URFI or *.RFI files containing information about radio frequency interference during the measurement and optionally the frequency pre-scan *.RAW files may be archived beside the raw data files *.USORT or *.SORT. The results (*.*CAL) of the automatic direct path test measurements will contain information about antenna health.

<table><caption>TABLE 11 | European and global marine data distribution portals served by the EU HFR Node.</caption><thead><tr><td>Distribution portal</td><td>Description</td><td>Portal link</td><td>Data access</td></tr></thead><tbody><tr><td>EMODNet-Physics</td><td>European Marine Observation and Data Network: Physics</td><td>http://www.emodnet-physics.eu/Map/</td><td>http://thredds.emodnet-physics.eu/thredds/HFRADARCatalog.html</td></tr><tr><td>Sea Data Net</td><td>Pan-European infrastructure for ocean and marine data management. The SeaDataNet portal offers different services: discovering, visualization, access and data downloading</td><td>http://seadatanet.maris2.nl/v_cdi_v3/search.asp</td><td></td></tr><tr><td>Copernicus Marine Service INSTAC</td><td>The In Situ TAC is the component of the Copernicus Marine Service which ensures a consistent and reliable access to a range of in situ data for the purpose of service production and validation. On the in situ TAC dashboard you can explore and download all the multi-source, multi-platform and heterogenous data collected in near-real-time (within 24 h) and delivered to the Copernicus Marine Service.</td><td>http://www.marineinsitu.eu/dashboard/</td><td>ftp://nrt.cmems-du.eu/Core/INSITU_GLO_NRT_ OBSERVATIONS_013_030\penalty-\@M (radar_total dataset to be available from April 2019) ftp://nrt.cmems-du.eu/Core/INSITU_GLO_UV_ NRT_ OBSERVATIONS_013_048/radar_total/(new product to be available from April 2019)</td></tr><tr><td>GEO Global High Frequency Radar Network</td><td>Vision for a global operational system measuring ocean surface currents to support monitoring of marine and coastal ecosystems</td><td>https://rucool.marine.rutgers.edu/geohfr/index.html</td><td></td></tr></tbody></table>
In addition to the on-site archiving, the HF radar remote stations should backup the most important data in real time to a remote repository through a network connection, as much as the bandwidth allows. The full backup of the data can be executed periodically on-site with a portable storage device. The minimum estimated amount of data to be archived per year and per station is in the order of 100 GB for CODAR systems, and 1 TB for WERA systems.

## Data Dissemination

Homogenization of file format, data description and Quality Control procedures are prerequisites for data sharing among data centers and data aggregators through the concepts of data interoperability and machine to machine interaction.

The authors strongly encourage HFR data providers to investigate the state of progress of regional and global HFR standards and infrastructures, and to take advantage of the service provided for enabling a data flow in standard interoperable format toward the community.

The European HF radar Node is an example of a recommended way for channeling HFR standardized data toward European and global data portals (Table 11).

The Node is responsible for the NRT HFR current data collection from the HF radar operators in Europe, for the application of the standard QC model, for the conversion of the QCed data files into the European standard data and metadata model. It is operational since April 2019 for the Copernicus Marine Service (CMEMS-INSTAC), SeaDataNet (SDN/SDC), and EMODnet Physics data delivery. The connection with US network for reciprocal data exchange is under development.

In US, radial files are collected in near real-time at regional level and then pushed into national nodes for total vector maps production and distribution through a Thematic Real-time Environmental Distributed Data Services (THREDDS) server<fnref n="2" />.

The Australian Ocean Data Network provides open access to real time quality controlled and delayed-mode quality-controlled HF radar derived surface currents, winds and waves through their THREDDS catalog<fnref n="3" /> as well as through the data portal<fnref n="4" />. One advantage of the data portal is that it allows for data to be aggregated over a user-defined time window, platform or region. Standard NetCDF format that adhere to the CF-1.6 and IMOS-1.4 conventions are available, either unaggregated or aggregated.

The Global HF Radar network web site<fnref n="5" /> shows a map with all of the locations of the HF Radar sites all over the world automatically connected.

In 2017 “the Global HFR Network was recognized by the Joint Technical WMO-IOC Commission for Oceanography and Marine Meteorology (JCOMM) as an observing network of the Global Ocean Observing System (GOOS)” (Roarty et al., 2019). This is expected to help coordinate the effort at international level and hopefully to define a unique global standard within the next few years.

It is worth mentioning that the data management plan for HF radar data should follow other few very general best practices as:

* follow the FAIR (Wilkinson et al., 2016) paradigm as part of making research data findable, accessible, interoperable and re-usable.

* include, inside metadata and as part of the data management plan, information about data policy, and possibly share unrestricted data, to make them freely available, thus facilitating access to a wider public. Creative Commons license Attribution 4.0 International (CC BY 4.0)<fnref n="6" /> is an example adopted by some HF radar data provider.

* associate a Digital Object Identification (DOI) to each dataset, and thus helping to make the data discoverable, accessible and citable.

Other general best practices on How to develop a Data Management Plan (International Oceanographic Data and Information Exchange, 2016) are available at the Ocean Best Practices repository.

## CONCLUSION

In response to the need for optimizing the operation performance of HF radars, different documents providing best practices for radar systems operation and maintenance have emerged in the past years. Most of them are either oriented to DF or BF systems, or to specific manufacturer's radar systems. In this document we compiled and completed existing documentation with the aim of offering a broad "Best Practices" manual for optimal operation of HF radar systems with independence from manufacturer or antenna design/setup.

This "Best Practices" paper, fed with direct experience of the authors on HF radar systems planning, setup, operations and data management and on the recent literature, provides some advancement with respect to the existing documents in several aspects. In particular:

* a more general approach to the technology is provided, including other examples of products and suggesting best practices referring when possible only to the antenna design and setup, and to how the spectral information is processed in order to determine the direction of arrival of the received signal (DF/BF techniques).

* the latest development in the field of HF radar set-up, maintenance and operation, are included. Among them, communication and remote management of the HF radar station, and new calibration (antenna pattern measurement) methods for DF configurations.

* a section dedicated to HF radar real time harmonized data management has been developed, from acquisition to post-processing, archiving, preservation and dissemination, with a look on interoperability at a global level.

* the need of distributing standardized and harmonized HF radar files has been highlighted, together with further

<footnote n="2">²https://hfrnet-tds.ucsd.edu/thredds/catalog.html</footnote>

<footnote n="3">³http://thredds.aodn.org.au/thredds/catalog/IMOS/ACORN/catalog.html</footnote>

<footnote n="4">⁴https://portal.aodn.org.au/</footnote>

<footnote n="5">⁵http://global-hfradar.org/</footnote>

<footnote n="6">⁶https://creativecommons.org/licenses/by/4.0/</footnote>
details on the European common data and metadata model for real time HF radar surface currents and Quality Control Tests and Flags as adopted by CMEMS *in situ* Thematic Assembly Center for HF radar data distribution since April 2019.

Many details, impossible to list in this paper, can be found in the reference manuals and in the literature. Although this document provides significant improvements and additional information to what is available in the literature, many aspects of HF radar management best practices are under discussion or at an early stage.

Future developments relate to: the description of other types of HF radars currently in use or emerging; other kinds of data measured by HFR systems (like tsunami detection, winds or waves); advanced methods for Quality Control and gap filling (either in real time or in delayed mode); new algorithms for signal processing optimization; a global standard for file format, data and metadata description.

## AUTHOR CONTRIBUTIONS

CM coordinated the work and contributed to all the sections. AR, JM, and AG contributed to the introduction and conclusion sections. CQ, JH, and SC contributed to introduction, setup and maintenance sections for beam forming HFR systems. LC and JA contributed to the data management and Quality Control section. ER contributed to Antenna Pattern Measurement, Maintenance and Data Management sections, editing of **Figure 3** and acronyms list. All authors contributed to the conception, writing and revising of the work for the final submission.

## FUNDING

Research was supported by the following projects: JERICO-NEXT (EU-H2020 GA No. 654410); RITMARE (Italian flagship project); INCREASE (CMEMS Service Evolution Call

for Tenders 21-SE-CALL1), SeaDataCloud (EU-H2020 GA No. 730960), Copernicus Marine Environment Monitoring Service (CMEMS) through IBISAR (CMEMS User Uptake Call for Tenders 67-UU-DO-CMEMS-DEM4_LOT7). Publication fees are funded by: JERICO-NEXT project (Joint European Research Infrastructure network for Coastal Observatory – Novel European eXpertise for coastal observaTories), funded by the European Union's Horizon 2020 Research and Innovation Program under Grant Agreement No. 654410.

## ACKNOWLEDGMENTS

This paper is based on work carried on within the JERICO-NEXT project (Joint European Research Infrastructure network for Coastal Observatory – Novel European eXpertise for coastal observaTories), funded by the European Union's Horizon 2020 Research and Innovation Program under grant agreement no. 654410. We would like to thank the projects RITMARE (Italian national funds), INCREASE (CMEMS Service Evolution Call for Tenders 21-SE-CALL1), SeaDataCloud (EU-H2020 GA n. 730960), in whose framework the activities of data harmonization and standardization have been carried on. We acknowledge financial support from Copernicus Marine Environment Monitoring Service (CMEMS) through IBISAR (CMEMS User Uptake Call for Tenders 67-UU-DO-CMEMS-DEM4_LOT7) and INSTAC (http://www.marineinsitu.eu/). CMEMS is implemented by Mercator Ocean in the framework of a delegation agreement with the European Union. IMOS was supported by the Australian Government through the National Collaborative Research Infrastructure Strategy and the Super Science Initiative. Thanks to Dr. Charles-Antoine Guerin (Université de Toulon, Laboratoire MIO) for useful discussion on signal processing techniques, and Dr. Mark Otero (Scripps Institution of Oceanography) for reviewing the final draft of the paper. Thanks to Integrated Ocean Observing System and Radiowave Operators Working Group for sharing best practices documents.

## REFERENCES

Abascal, A. J., Castanedo, S., Medina, R., Losada, I. J., and Álvarez-Fanjul, E. (2009). Application of HF radar currents to oil spill modelling. *Mar. Pollut. Bull.* 58, 238-248. doi: 10.1016/j.marpolbul.2008.09.020

Barrick, D. E. (1977). The ocean waveheight nondirectional spectrum from inversion of the HF sea-echo doppler spectrum. *Remote Sens. Environ.* 6, 201-227. doi: 10.1016/0034-4257(77)90004-9

Barrick, D. E. (2002). *Geometrical Dilution of Statistical Accuracy (GDOSA) in Multi-Static HF Radar Networks*. Mountain View, CA: Codar ocean sensors.

Barrick, D. E., Evans, M. W., and Weber, B. L. (1977). Ocean surface currents mapped by radar. *Science* 198, 138-144. doi: 10.1126/science.198.4313.138

Barrick, D. E., and Lipa, B. J. (1997). Evolution of bearing determination in HF current mapping radars. *Oceanography* 10, 72-75. doi: 10.5670/oceanog.1997

Barth, A., Alvera-Azcárate, A., Beckers, J. M., Staneva, J., Stanev, E. V., and Schulz-Stellenfleth, J. (2011). Correcting surface winds by assimilating HFR surface currents in the german bight. *Ocean Dynam.* 61, 599-610. doi: 10.1007/s10236-010-0369-0

Barth, A., Alvera-Azcárate, A., and Weisberg, R. H. (2008). Assimilation of high-frequency radar currents in a nested model of the West Florida shelf. *J. Geophys. Res.* 103:1311-1308. doi: 10.1029/2007JC004585

Breivik, O., and Saetra, O. (2001). Real time assimilation of HF radar currents into a coastal ocean model. *J. Mar. Syst.* 28, 161-182. doi: 10.1016/s0924-7963(01)00002-1

Chapman, R. D., Shay, L., Graber, H., Edson, J., Karachintsev, A., Trump, C., et al. (1997). On the accuracy of HF radar surface current measurements: intercomparisons with ship-based sensors. *J. Geophys. Res. Oceans* 102, 737-718.

Cook, T., Hazard, L., Otero, M., and Zelenke, B. (2008). *Deployment and Maintenance of a High-frequency Radar (HF radar) for Ocean Surface Current Mapping: Best Practices*. La Jolla, CA: University of California.

Corgnati, L., Mantovani, C., Novellino, A., Rubio, A., Mader, J., Reyes, E., et al. (2017). Recommendation Report 1 for HFR data implementation in European marine data infrastructures, Deliverable D5.13 of the EU project Jerico-NEXT. Available at: http://www.jerico-ri.eu/download/jerico-next-deliverables/JERICO-NEXT-Deliverable-5.13_V1.pdf

Corgnati, L., Mantovani, C., Novellino, A., Rubio, A., Mader, J., Reyes, E., et al. (2018a). Recommendation Report 2 on improved common procedures for HFR
QC analysis, Deliverable D5.14 of the project H2020 JERICO-RI, in submission on OceanBestPractice repository. Available at: http://www.jerico-ri.eu/download/jerico-next-deliverables/JERICO-NEXT-Deliverable_5.14_V1.pdf

Corgnati, L., Mantovani, C., Rubio, A., Igoa, J. L., Reyes, E., Novellino, A., et al. (2018b). Building strong foundations towards the pan-European high frequency radar network. *Boll. Geof. Suppl.* 59, 246–247.

Cosoli, S., and Bolzon, G. (2010). A real-time and offline quality control methodology for seasonde high-frequency radar currents. *J. Atmos. Ocean. Tech.* 29, 1313–1328. doi: 10.1175/JTECH-D-11-00217.1

Cosoli, S., and Grcic, B. (2019). *Quality Control Procedures for IMOS Ocean Radar Manual Version 2.1*. Hobart, TAS: Integrated Marine Observing System, 113. doi: 10.26198/5c89b59a931cb

Crombie, D. D. (1955). Doppler Spectrum of Sea Echo at 13.56 Mc./s. *Nature* 175, 681–682. doi: 10.1038/175681a0

Emery, B. M. (2018). Evaluation of alternative direction-of-arrival methods for oceanographic HF radars. *IEEE J. Ocean. Eng.* 99, 1–14. doi: 10.1109/JOE.2019.2914537

Emery, B. M., and Washburn, L. (2007). *Evaluation of SeaSonde Hardware Diagnostic Parameters as Performance Metrics*. Silver Spring, MA: NOAA.

Emery, B. M., Washburn, L., Whelan, C., Barrick, D., and Harlan, J. (2014). Measuring antenna patterns for ocean surface current HF radars with ships of opportunity. *J. Atmos. Ocean. Technol.* 31, 1564–1582. doi: 10.1175/jtech-d-13-00181.1

European Union (2013). *Official Journal, Directive /35/EU of the European Parliament and of the Council of 26 June 2013 on the Minimum Health and Safety Requirements Regarding the Exposure of Workers to the Risks Arising from Physical Agents (Electromagnetic Fields) (20th Individual Directive within the Meaning of Article 16(1) of Directive 89/391/EEC) and Repealing Directive 2004/40/EC*. Brussels: European Union.

Flament, P., Harris, D., Flament, M., Fernandez, I. Q., Hlivak, M., Flores, X., et al. (2016). *A Compact High Frequency Doppler Radio Scatterometer for Coastal Oceanography*. Washington, DC: American Geophysical Union.

Forget, P. (2015). Noise properties of HF radar measurement of ocean surface currents. *Radio Sci.* 50, 1–14. doi: 10.1002/2015RS005681

Graber, H. C., Haus, B. K., Chapman, R. D., and Shay, L. K. (1997). HF radar comparisons with moored estimates of current speed and direction: expected differences and implications. *J. Geophys Res.* 102, 18749–18766. doi: 10.1029/97JC01190

Gurgel, K.-W., Antonischki, G., Essen, H.-H., and Schlick, T. (1999). Wellen radar WERA: a new ground-wave HF radar for ocean remote sensing. *Coast. Eng.* 37, 219–234. doi: 10.1016/S0378-3839(99)00027-7

Gurgel, K.-W., Dzvonkovskaya, A., Pohlmann, T., Schlick, T., and Gill, E. (2011). Simulation and detection of Tsunami signatures in ocean surface currents measured by HF radar. *Ocean Dyn.* 61, 1495–1507. doi: 10.1007/s10236-011-0420-9

Gurgel, K.-W., Essen, H.-H., and Schlick, T. (2006). An empirical method to derive ocean waves from second-order bragg scattering: prospects and limitations. *IEEE J. Ocean. Eng.* 31, 804–811. doi: 10.1109/JOE.2006.886225

Helzel, T., and Kniephoff, M. (2010). “Software beam forming for ocean radar WERA,” in *Proceedings of the IEEE Oceans Conference*, Seattle.

Heron, M. L., Dexter, P. E., and McGann, B. T. (1985). Parameters of the air-sea interface by high-frequency ground-wave Doppler radar. *Aust. J. Mar. Freshw. Res.* 36, 655–670. doi: 10.1071/MF9850655

Heron, M., and Rose, R. (1986). On the application of HF ocean radar to the observation of temporal and spatial changes in wind direction. *J. Oceanic Eng.* 11, 210–218. doi: 10.1109/JOE.1986.1145173

Hickey, K., Khan, R. H., and Walsh, J. (1995). Parametric estimation of ocean surface currents with HF radar. *J. Oceanic Eng.* 20, 139–144. doi: 10.1109/48.376678

Horstmann, J., Corgnati, L., Mantovani, C., Quentin, C., Reyes, E., Rubio, A., et al. (2019). Report on Best Practice in the implementation and use of new systems in JERICO-RI. Part 1: HF-radar systems. Deliverable D2.4 of the project H2020 JERICO-RI. In submission on OceanBestPractice repository. Available at: http://www.jerico-ri.eu/download/jerico-next-deliverables/JERICO_NEXT_Deliverable_2.4_final.pdf

Iermano, I., Moore, A. M., and Zambianchi, E. (2016). Impacts of a 4-dimensional variational data assimilation in a coastal ocean model of Southern Tyrhenian Sea. *J. Mar. Sys.* 154, 157–171. doi: 10.1016/j.jmarsys.2015.09.006

International Oceanographic Data and Information Exchange (2016). *Guidelines for a Data Management Plan*. Paris: UNESCO.

Kaplan, D. M., Largier, J., and Botsford, L. W. (2005). HF radar observations of surface circulation off Bodega Bay (northern California, USA). *J. Geophys. Res.* 110:C10020. doi: 10.1029/2005JC002959

Kirincich, A. (2016). Remote sensing of the surface wind field over the coastal ocean via direct calibration of HF radar backscatter power. *J. Atmos. Oceanic Technol.* 33, 1377–1392. doi: 10.1175/JTECH-D-15-0242.1

Kirincich, A. (2017). Improved detection of the first-order region for direction-finding HF radars using image processing techniques. *J. Atmos. Oceanic Technol.* 34, 1679–1691. doi: 10.1175/JTECH-D-16-0162.1

Kohut, J. T., and Glenn, S. M. (2003). Improving HF radar surface current measurements with measured antenna beam patterns. *J. Atmos. Ocean. Tech.* 20, 1303–1316. doi: 10.1175/1520-0426(2003)020<1303:ihrscm>2.0.co;2

Lipa, B. J. (2013). “Uncertainties in SeaSonde current velocities,” in *Proceedings of the IEEE/OES Seventh Working Conference on Current Measurement Technology*, San Diego, CA.

Lipa, B. J., and Barrick, D. E. (1983). Least-squares methods for the extraction of surface currents from CODAR crossed-loop data: application at ARSLOE. *IEEE J. Oceanic Eng.* 8, 1–28.

Lipa, B. J., Barrick, D. E., Bourg, J., and Nyden, B. B. (2006). HF radar detection of Tsunamis. *J. Oceanogr.* 62, 705–716. doi: 10.1007/s10872-006-0088-9

Mader, J., Rubio, A., Asensio, J. L., Novellino, A., Alba, M., Corgnati, L., et al. (2016). *The European HF Radar Inventory*. Southampton: EuroGOOS.

Maresca, S., Braca, P., Horstmann, J., and Grasso, R. (2014). Maritime surveillance using multiple high-frequency surface-wave radars. *IEEE Trans. Geosci. Remote Sens.* 52, 5056–5071. doi: 10.1109/tgrs.2013.2286741

Marmain, J., Molcard, A., Forget, P., and Barth, A. (2014). Assimilation of HF radar surface currents to optimize forcing in the North Western Mediterranean Sea. *Nonlinear Process. Geophys.* 21, 659–675. doi: 10.5194/npg-21-659-2014

Paduan, J. D., and Rosenfeld, L. K. (1996). Remotely sensed surface currents in monterey bay from shore-based HF radar (coastal ocean dynamics application radar). *J. Geophys. Res. Oceans* 101, 20669–20686. doi: 10.1029/96JC01663

Paduan, J. D., and Shulman, I. (2004). HFR Data Assimilation in the Monterey Bay Area. *J. Geophys. Res.* 109:C07S09. doi: 10.1029/2003JC001949

Pearlman, J., Bushnell, M., Coppola, L., Karstensen, J., Buttigie, P. L., Pearlman, F., et al. (2019). Evolving and sustaining ocean best practices and standards for the next decade. *Front. Mar. Sci.* 6:277. doi: 10.3389/fmars.2019.00277

Philip, A. (2012). *Oceanographic Applications of Coastal Radar*. Report Oceanography No 41. Stockholm: Environmental Protection Agency.

Ponsford, A. M., Sevgi, L., and Chan, H. C. (2001). An integrated maritime surveillance system based on high-frequency surface-wave radars. operational status and system performance. *Antennas Propag. Mag.* 43, 52–63. doi: 10.1109/74.979367

Prandle, D., Loch, S. G., and Player, R. (1992). Tidal flow through the straits of Dover. *J. Phys. Oceanogr.* 23, 23–37. doi: 10.1175/1520-0485(1993)023<0023: tftso>2.0.co;2

Roarty, H., Cook, T., Hazard, L., George, D., Harlan, J., Cosoli, S., et al. (2019). The global high frequency radar network. *Front. Mar. Sci.* 6:164. doi: 10.3389/fmars. 2019.00164

Roarty, H., Hazard, L., and Fanjul, E. A. (2016). Growing network of radar systems monitors ocean surface currents. *Eos* 97. doi: 10.1029/2016EO049243

Rubio, A., Mader, J., Corgnati, L., Mantovani, C., Griffa, A., Novellino, A., et al. (2017). HF radar activity in european coastal seas: next steps toward a pan-European HF radar network. *Front. Mar. Sci.* 4:8. doi: 10.3389/fmars.2017.00008

Rubio, A., Mader, J., Corgnati, L., Mantovani, C., Griffa, A., Novellino, A., et al. (2018). “Present and future of the European HF radar network: outcomes of the INCREASE project,” in *Proceedings of the 4th Orca Meeting*, Okinawa.

Schmidt, R. O. (1986). Multiple emitter location and signal parameter estimation. *Trans. Antennas Propag.* 34, 276–280. doi: 10.1109/TAP.1986.1143830

SeaSonde® Remote Unit Operator's Manual (2010). *SSRS-100 Product Series October 19*. Mountain View, CA: Codar Ocean Sensors.

Shearman, E. D. R., and Moorhead, M. D. (1988). PISCES: a coastal ground-wave HF radar for current, wind and wave mapping to 200-km ranges. *Proc. IGARSS* 88, 773–776.
Shen, W., Gurgel, K. W., Voulgaris, G., Schlick, T., and Stammer, D. (2012). Wind-speed inversion from HF radar first-order backscatter signal. *Ocean Dyn.* 62, 105–121. doi: 10.1007/s10236-011-0465-9

Stanev, E. V., Ziemer, F., Schulz-Stellenfleth, J., Seemann, J., Staneva, J., and Gurgel, K. W. (2015). Blending surface currents from HF radar observations and numerical modeling: tidal hindcasts and forecasts. *J. Atmos. Oceanic Technol.* 32, 256–281. doi: 10.1175/jtech-d-13-00164.1

Statscewich, H., Weingartner, T., Danielsen, S., Grunau, B., Egan, G., and Timm, J. (2011). A high-latitude modular autonomous power, control, and communication system for application to high-frequency surface current mapping radars. *Mar. Technol. Soc. J.* 45, 59–68. doi: 10.4031/MTSJ.45.3.7

U.S. Integrated Ocean Observing System, (2016). *Manual for Real-Time Quality Control of High Frequency Radar Surface Current Data: A Guide to Quality Control and Quality Assurance for High Frequency Radar Surface Current Observations. Version 1.0.* Silver Spring, MD: U.S. Department of Commerce.

Ullman, D. S., O'Donnell, J., Kohut, J., Fake, T., and Allen, A. (2006). Trajectory prediction using HF radar surface currents: monte carlo simulations of prediction uncertainties. *J. Geophys. Res.* 111:C12005. doi: 10.1029/2006jc003715

Washburn, L., Romero, E., Johnson, C., Emery, B., and Gotschalk, C. (2017). Measurement of antenna patterns for oceanographic radars using aerial drones. *J. Atmos. Oceanic Technol.* 34, 971–981. doi: 10.1175/jtech-d-16-0180.1

Wera Best Practices (2018). *WERA Remote Ocean Sensing*. WERA Best Practice. Version 1. Kaltenkirchen: Helzel Messtechnik GmbH.

Wilkinson, M. D., Dumontier, M., Aalbersberg, I. J., Appleton, G., Axton, M., Baak, A., et al. (2016). The FAIR Guiding Principles for scientific data management and stewardship. *Sci. Data* 3:160018. doi: 10.1038/sdata.2016.18

Wyatt, L. R. (1990). A relaxation method for integral inversion applied to HF radar measurement of the ocean wave directional spectrum. *Int. J. Remote Sens.* 11, 1481–1494. doi: 10.1080/01431169008955106

Wyatt, L. R., Green, J. J., and Middleditch, A. (2007). *Data Quality and Sampling Requirements for Reliable Wave Measurement with HF Radar*. Vancouver, BC: IEEE.

Wyatt, L. R., Green, J. J., Middleditch, A., Moorhead, M. D., Howarth, J., Holt, M., et al. (2006). Operational wave, current and wind measurements with the pisces HFR. *J. Ocean Eng.* 31, 819–834. doi: 10.1109/JOE.2006.888378

Wyatt, L. R., Jaffrés, J. B., and Heron, M. L. (2013). Spatial averaging of HF radar data for wave measurement applications. *J. Atmos. Oceanic Technol.* 30, 2216–2224. doi: 10.1175/JTECH-D-12-00206.1

Yang, L., Dang, R., Li, M., Zhao, K., Song, C., and Xu, Z. (2018). A fast calibration method for phased arrays by using the graph coloring theory. *Sensors* 18:4315. doi: 10.3390/s18124315

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2020 Mantovani, Cognati, Horstmann, Rubio, Reyes, Quentin, Cosoli, Asensio, Mader and Griffa. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.