Simulation and detection of tsunami signatures in ocean surface currents measured by HF radar

High-frequency (HF) surface wave radars provide the unique capability to continuously monitor the coastal environment far beyond the range of conventional microwave radars. Bragg-resonant backscattering by ocean waves with half the electromagnetic radar wavelength allows ocean surface currents to be measured at distances up to 200 km. When a tsunami propagates from the deep ocean to shallow water, a specific ocean current signature is generated throughout the water column. Due to the long range of an HF radar, it is possible to detect this current signature at the shelf edge. When the shelf edge is about 100 km in front of the coastline, the radar can detect the tsunami about 45 min before it hits the coast, leaving enough time to issue an early warning. As up to now no HF radar measurements of an approaching tsunami exist, a simulation study has been done to fix parameters like the required spatial resolution or the maximum coherent integration time allowed. The simulation involves several steps, starting with the Hamburg Shelf Ocean Model (HAMSOM) which is used to estimate the tsunami-induced current velocity at 1 km spatial resolution and 1 s time step. This ocean current signal is then superimposed to modelled and measured HF radar backscatter signals using a new modulation technique. After applying conventional HF radar signal processing techniques, the surface current maps contain the rapidly changing tsunami-induced current features, which can be compared to the HAMSOM data. The specific radial tsunami current signatures can clearly be observed in these maps, if appropriate spatial and temporal resolution is used. Based on the entropy of the ocean current maps, a tsunami detection algorithm is described which can be used to issue an automated tsunami warning message.     

Assessing the fidelity of surface currents from a coastal ocean model and HF radar using drifting buoys in the Middle Atlantic Bight

The rapid expansion of urbanization along the world’s coastal areas requires a more comprehensive and accurate understanding of the coastal ocean. Over the past several decades, numerical ocean circulation models have tried to provide such insight, based on our developing understanding of physical ocean processes. The systematic establishment of coastal ocean observation systems adopting cutting-edge technology, such as high frequency (HF) radar, satellite sensing, and gliders, has put such ocean model predictions to the test, by providing comprehensive observational datasets for the validation of numerical model forecasts. The New York Harbor Observing and Prediction System (NYHOPS) is a comprehensive system for understanding coastal ocean processes on the continental shelf waters of New York and New Jersey. To increase confidence in the system’s ocean circulation predictions in that area, a detailed validation exercise was carried out using HF radar and Lagrangian drifter-derived surface currents from three drifters obtained between March and October 2010. During that period, the root mean square (RMS) differences of both the east–west and north–south currents between NYHOPS and HF radar were approximately 15 cm s^?1. Harmonic analysis of NYHOPS and HF radar surface currents shows similar tidal ellipse parameters for the dominant M₂ tide, with a mean difference of 2.4 cm s^?1 in the semi-major axis and 1.4 cm s^?1 in the semi-minor axis and 3° in orientation and 10° in phase. Surface currents derived independently from drifters along their trajectories showed that NYHOPS and HF radar yielded similarly accurate results. RMS errors when compared to currents derived along the trajectory of the three drifters were approximately 10 cm s^?1. Overall, the analysis suggests that NYHOPS and HF radar had similar skill in estimating the currents over the continental shelf waters of the Middle Atlantic Bight during this time period. An ensemble-based set of particle tracking simulations using one drifter which was tracked for 11 days showed that the ensemble mean separation generally increases with time in a linear fashion. The separation distance is not dominated by high frequency or short spatial scale wavelengths suggesting that both the NYHOPS and HF radar currents are representing tidal and inertial time scales correctly and resolving some of the smaller scale eddies. The growing ensemble mean separation distance is dominated by errors in the mean flow causing the drifters to slowly diverge from their observed positions. The separation distance for both HF radar and NYHOPS stays below 30 km after 5 days, and the two technologies have similar tracking skill at the 95 % level. For comparison, the ensemble mean distance of a drifter from its initial release location (persistence assumption) is estimated to be greater than 70 km in 5 days.     

Backtracking drifting objects using surface currents from high-frequency (HF) radar technology

In this work, the benefits of high-frequency (HF) radar ocean observation technology for backtracking drifting objects are analysed. The HF radar performance is evaluated by comparison of trajectories between drifter buoys versus numerical simulations using a Lagrangian trajectory model. High-resolution currents measured by a coastal HF radar network combined with atmospheric fields provided by numerical models are used to backtrack the trajectory of two dataset of surface-drifting buoys: group I (with drogue) and group II (without drogue). A methodology based on optimization methods is applied to estimate the uncertainty in the trajectory simulations and to optimize the search area of the backtracked positions. The results show that, to backtrack the trajectory of the buoys in group II, both currents and wind fields were required. However, wind fields could be practically discarded when simulating the trajectories of group I. In this case, the optimal backtracked trajectories were obtained using only HF radar currents as forcing. Based on the radar availability data, two periods ranging between 8 and 10?h were selected to backtrack the buoy trajectories. The root mean squared error (RMSE) was found to be 1.01?km for group I and 0.82?km for group II. Taking into account these values, a search area was calculated using circles of RMSE radii, obtaining 3.2 and 2.11?km² for groups I and II, respectively. These results show the positive contribution of HF radar currents for backtracking drifting objects and demonstrate that these data combined with atmospheric models are of value to perform backtracking analysis of drifting objects.     

Wind-speed inversion from HF radar first-order backscatter signal

Land-based high-frequency (HF) radars have the unique capability of continuously monitoring ocean surface environments at ranges up to 200 km off the coast. They provide reliable data on ocean surface currents and under slightly stricter conditions can also give information on ocean waves. Although extraction of wind direction is possible, estimation of wind speed poses a challenge. Existing methods estimate wind speed indirectly from the radar derived ocean wave spectrum, which is estimated from the second-order sidebands of the radar Doppler spectrum. The latter is extracted at shorter ranges compared with the first-order signal, thus limiting the method to short distances. Given this limitation, we explore the possibility of deriving wind speed from radar first-order backscatter signal. Two new methods are developed and presented that explore the relationship between wind speed and wave generation at the Bragg frequency matching that of the radar. One of the methods utilizes the absolute energy level of the radar first-order peaks while the second method uses the directional spreading of the wind generated waves at the Bragg frequency. For both methods, artificial neural network analysis is performed to derive the interdependence of the relevant parameters with wind speed. The first method is suitable for application only at single locations where in situ data are available and the network has been trained for while the second method can also be used outside of the training location on any point within the radar coverage area. Both methods require two or more radar sites and information on the radio beam direction. The methods are verified with data collected in Fedje, Norway, and the Ligurian Sea, Italy using beam forming HF WEllen RAdar (WERA) systems operated at 27.68 and 12.5 MHz, respectively. The results show that application of either method requires wind speeds above a minimum value (lower limit). This limit is radar frequency dependent and is 2.5 and 4.0 m/s for 27.68 and 12.5 MHz, respectively. In addition, an upper limit is identified which is caused by wave energy saturation at the Bragg wave frequency. Estimation of this limit took place through an evaluation of a year long database of ocean spectra generated by a numerical model (third generation WAM). It was found to be at 9.0 and 11.0 m/s for 27.68 and 12.5 MHz, respectively. Above this saturation limit, conventional second-order methods have to be applied, which at this range of wind speed no longer suffer from low signal-to-noise ratios. For use in operational systems, a hybrid of first- and second-order methods is recommended.     

A short-term predictive system for surface currents from a rapidly deployed coastal HF radar network

In order to address the need for surface trajectory forecasts following deployment of coastal HF radar systems during emergency-response situations (e.g., search and rescue, oil spill), a short-term predictive system (STPS) based on only a few hours data background is presented. First, open-modal analysis (OMA) coefficients are fitted to 1-D surface currents from all available radar stations at each time interval. OMA has the effect of applying a spatial low-pass filter to the data, fills gaps, and can extend coverage to areas where radial vectors are available from a single radar only. Then, a set of temporal modes is fitted to the time series of OMA coefficients, typically over a short 12-h trailing period. These modes include tidal and inertial harmonics, as well as constant and linear trends. This temporal model is the STPS basis for producing up to a 12-h current vector forecast from which a trajectory forecast can be derived. We show results of this method applied to data gathered during the September 2010 rapid-response demonstration in northern Norway. Forecasted coefficients, currents, and trajectories are compared with the same measured quantities, and statistics of skill are assessed employing 16 24-h data sets. Forecasted and measured kinetic variances of the OMA coefficients typically agreed to within 10–15%. In one case where errors were larger, strong wind changes are suspected and examined as the cause. Sudden wind variability is not included properly within the STPS attack we presently employ and will be a subject for future improvement.     

Wind influence on surface current variability in the Ibiza Channel from HF Radar

Surface current variability is investigated using 2.5 years of continuous velocity measurements from an high frequency radar (HFR) located in the Ibiza Channel (Western Mediterranean Sea). The Ibiza Channel is identified as a key geographical feature for the exchange of water masses but still poorly documented. Operational, quality controlled, HFR derived velocities are provided by the Balearic Islands Coastal Observing and Forecasting System (SOCIB). They are assessed by performing statistical comparisons with current-meter, ADCP, and surface lagrangian drifters. HFR system does not show significant bias, and its accuracy is in accordance with previous studies performed in other areas. The main surface circulation patterns are deduced from an EOF analysis. The first three modes represent almost 70 % of the total variability. A cross-correlation analysis between zonal and meridional wind components and the temporal amplitudes of the first three modes reveal that the first two modes are mainly driven by local winds, with immediate effects of wind forcing and veering following Ekman effect. The first mode (37 % of total variability) is the response of meridional wind while the second mode (24 % of total variability) is linked primarily with zonal winds. The third and higher order modes are related to mesoscale circulation features. HFR derived surface transport presents a markedly seasonal variability being mostly southwards. Its comparison with Ekman-induced transport shows that wind contribution to the total surface transport is on average around 65 %.     

Use of surface currents measured by HF radar in planning coastal discharges

Surface currents measured by HF Radar (OSCR) are analysed to determine near-shore circulation in connection with planning optimum locations for sewage discharges. Dual radar units, each measuring a radial velocity component, were deployed at a total of five sites for three separate periods each of one week's duration (one site being common to two deployments). Observed data from 20 radial ‘bins’ each 1.2 km long were obtained for six beams each 6° wide at all five sites. These data were analysed separately to determine 1. the tidal constituents, 2. the wind driven component and 3. the steady residual component. In deriving the wind-driven response a relationship

was assumed where R(t) is the non-tidal velocity, W_E and W_N east and west components of wind speed, Δt a time lag and α and β empirical coefficients determined by least squares fitting techniques. At positions where beams from any two sites crossed, their separate current components were combined, producing the following spatial distributions of surface current vectors: 1. tidal ellipses, 2. separate responses to unit wind forcing from a. west and b. south and 3. steady residual drifts. The validity of these current distributions is examined by 1. comparison with other observational data, 2. by consideration of their self coherence and 3. in terms of simple theory.     

Range resolution dependence of VHF radar returns from clear-air turbulence and precipitation

With employing 1.5 h of the data observed by the Chung-Li VHF radar, the range resolution dependences of the VHF backscatter from refractivity fluctuation and precipitation are investigated in this article. It indicates that the atmospheric layer structure of refractivity seems to play a role in governing the range resolution dependence of clear-air turbulent echoes. Observations shows that the VHF clear-air echo power ratios for 4 to 2 μs pulse lengths are close to 3 dB in the middle or bottom side of the layer, while the ratios are significantly greater than 3 dB in the top side of the layer. The analysis of the precipitation echo power ratio for 4 to 2 ms pulse lengths shows that basically the ratios above 3.0 km are close to 3 dB, but enormously smaller than 3 dB below 3.0 km. The feature of extraordinarily small echo power ratios below 3.0 km is also observed for the radar returns from refractivity turbulence. The radar recovery effect is thought to be a primary factor responsible for the severe diminution of the echo power ratios at the lower altitudes. In addition, statistical analysis reveals that the range resolution effect on the first and second moments of the Doppler spectra for the radar echoes from clear-air turbulence and precipitation is insignificant and negligible. The dependences of the coefficient A and power B in the power-law approximation V_t=AP^B_r to the terminal velocity V_t and range-corrected echo power P_r are examined theoretically and experimentally. The results show that the coefficient A (powers B) is inversely (positively) proportional to the range resolution, in a good agreement with the observations. Because of the strong dependence of coefficient A and power B on the radar pulse width, it suggests that great caution should be taken in comparing the power-law expressions V_t=AP^B_r established from the radar returns obtained with different range resolutions.     

Transport properties in small-scale coastal flows: relative dispersion from VHF radar measurements in the Gulf of La Spezia

Lagrangian transport characteristics in the Gulf of La Spezia, a 5 × 10-km area along the western coast of Italy, are investigated using data collected from a very high frequency (VHF) radar system with 250 m and 30-min resolution and two clusters of Coastal Dynamics Experiment surface drifters during 2 weeks in the summer of 2007. The surface drifters are found to follow the temporal and spatial evolution of the finite-scale Lyapunov exponents (FSLEs) computed by the VHF radar, indicating the precision of both the radar measurements and the diagnostic FSLE in mapping accurately the transport pathways. In light of this agreement, an analysis of the relative dispersion is conducted. It is found that the average FSLE value varies within a narrow range of 4   day^-1 ￡ l ￡ 7   day^-14 \;\mbox{day}^{-1} \leq \lambda \leq 7 \;\mbox{day}^{-1} and displays an exponential regime over the entire extent of the measurements. The dynamical implication is that relative dispersion is controlled nonlocally, namely by slow, persistent, energetic mesoscale structures as opposed to the rapidly evolving high-gradient small-scale turbulent features. The value of the exponent is about an order of magnitude larger than those found in previous modeling studies and analysis of SCULP data in the Gulf of Mexico but somewhat smaller than that estimated from CLIMODE drifters in the Gulf Stream region. Scaling of the FSLE using a metric of resolved gradients of the Eulerian fields in the form of a positive Okubo–Weiss criterion is useful, but not as precise as in modeling studies. The horizontal flow convergence is found to have a small yet tangible effect on relative dispersion.     

VHF radar observations of surface currents off the northern Opal coast in the eastern English Channel

Two very high-frequency radars (VHFR) operating on the Opal coast of eastern English Channel provided a nearly continuous 35-day long dataset of surface currents over a 500 km² area at 0.6–1.8 km resolution. Argo drifter tracking and CTD soundings complemented the VHFR observations, which extended approximately 25 km offshore. The radar data resolve three basic modes of the surface velocity variation in the area, that are driven by tides, winds and freshwater fluxes associated with seasonal river discharge. The first mode, accounting for 90% of variability, is characterized by an along-shore flow pattern, whereas the second and third modes exhibit cross-shore, and eddy-like structures in the current velocity field. All the three modes show the dominant semi-diurnal variability and low-frequency modulation by the neap-spring tidal cycle. Although tidal forcing provides the major contribution to variability of local currents, baroclinicity plays an important role in shaping the 3D velocity field averaged over the tidal cycle and may strongly affect tracer dynamics on larger time scales. An empirical orthogonal function (EOF) decomposition and a spectral rotary analysis of the VHFR data reveal a discontinuity in the velocity field occurring approximately 10 km offshore which was caused by the reversal in the sign of rotation of the current vector. This feature of local circulation is responsible for surface current convergence on ebb, divergence on flood and strong oscillatory vertical motion. Spectral analysis of the observed currents and the results of the Agro drifter tracking indicate that the line of convergence approximately follows the 30-m isobath. The most pronounced feature of the radar-derived residual circulation is the along-coast intensification of surface currents with velocity magnitude of 0.25 m/s typical for the Regions of Freshwater Influence (ROFI). The analysis has provided a useful, exploratory examination of surface currents, suggesting that the circulation off the Opal coast is governed by ROFI dynamics on the hypertidal background.     

Letter to the editor A comparison of F-region ion velocity observations from the EISCAT Svalbard and VHF radars with irregularity drift velocity measurements from the CUTLASS Finland HF radar

During August 1998, the UK EISCAT special programme SP-UK-CSUB, which combines operation of both the mainland VHF and Svalbard UHF incoherent scatter radars, was run for several hours around magnetic midnight on four consecutive days. The CUTLASS Finland HF coherent scatter radar was, at these times, operating in a discretionary mode, sounding on all 16 beams, one at high-time resolution. This study presents a comparison of the velocities measured by coherent and incoherent techniques during the SP-UK-CSUB experiments. Agreement, particularly between the ion velocities measured by the EISCAT Svalbard radar and irregularity drift measurements by the Finland radar, is remarkable, thereby validating the scientific integrity of both data sets. This work highlights the substantive contribution to our understanding of the solar-terrestrial environment which can be made by use in concert of incoherent and HF coherent scatter radars.     

Rapid variations in echo power maps of VHF radar backscatter from the lower atmosphere

For the first time, echo power maps of aspect-sensitive VHF backscatter are shown, with good time and spatial resolutions, for angles 0°–7° from zenith. Sequences of power maps show large changes in appearance over timescales of a few minutes and height intervals of a few hundred metres. Often, individual power maps are consistent with tilted and distorted specular-type scattering layers, rather than anisotropic turbulence, and the direction of maximum echo power is sometimes several degrees off-vertical. Nevertheless, after time-averaging the variable echo-power patterns, the average pattern can become almost circular and centred on zenith, as has been assumed in the past. Echo power maps measured in strong windshear beneath the jet stream show a skewing of the echo power distribution. However, some power maps in the lower stratosphere, despite stronger wind shear, appear more constrained and their maximum echo power remains closer to zenith.     

Angles of arrival of HF mesosphere radar echoes

The distribution of the angles of arrival of mesosphere radio echoes has been obtained and studied. Special attention has been paid to the region of altitudes in the vicinity of the mesopause. It has been indicated that the angle of arrival distribution in this region has a maximum near the direction toward the zenith and is anisotropic. At the same time, the average intensity of radio echoes substantially exceeds the background level for any angles of arrival.     

VHF radar observations of gravity waves at a low latitude

Wind observations made at Gadanki (13.5°N) by using Indian MST Radar for few days in September, October, December 1995 and January, 1996 have been analyzed to study gravity wave activity in the troposphere and lower stratosphere. Horizontal wind variances have been computed for gravity waves of period (2–6) h from the power spectral density (PSD) spectrum. Exponential curves of the form e^Z/H have been fitted by least squares technique to these variance values to obtain height variations of the irregular winds upto the height of about 15 km, where Z is the height in kilometers. The value of H, the scale height, as determined from curve fitting is found to be less than the theoretical value of scale height of neutral atmosphere in this region, implying that the waves are gaining energy during their passage in the troposphere. In other words, it indicates that the sources of gravity waves are present in the troposphere. The energy densities of gravity wave fluctuations have been computed. Polynomial fits to the observed values show that wave energy density increases in the troposphere, its source region, and then decreases in the lower stratosphere.     

Model of the Long Island Sound outflow: Comparison with year-long HF radar and Doppler current observations

A three-dimensional primitive-equation model is used to simulate the Long Island Sound (LIS) outflow for a 1-year (2001) period. The model domain includes LIS and New York Bight (NYB). Tidal and wind forcing are included, and seasonal salinity and temperature variations are assimilated. The model results are validated with the HF radar, moored acoustic Doppler current profiler (ADCP), and ferry-based ADCP observations. The agreement between simulated and observed flow patterns generally is very good. The difference in seasonal mean currents between the model and moored ADCP is about 0.01 m/s; the correlation of dominant velocity fluctuations between the model and HF radar is 0.83; and the difference in mean LIS transport between the model and shipboard ADCP is about 5%. However, the model predicts a prominent tidally generated headland eddy not supported by the HF radar observation. The model sensitivity study indicates that the tides, winds, and ambient coastal front all have important impact on the buoyant outflow. The tides and winds cause stronger vertical mixing, which reduces the surface plume strength. The ambient coastal front, on the other hand, tends to enhance the plume.     

A statistical study of underestimates of wind speeds by VHF radar

Comparisons are made between horizontal wind measurements carried out using a VHF-radar system at Aberystwyth (52.4°N, 4.1°W) and radiosondes launched from Aberporth, some 50 km to the southwest. The radar wind results are derived from Doppler wind measurements at zenith angles of 6° in two orthogonal planes and in the vertical direction. Measurements on a total of 398 days over a 2-year period are considered, but the major part of the study involves a statistical analysis of data collected during 75 radiosonde flights selected to minimise the spatial separation of the two sets of measurements. Whereas good agreement is found between the two sets of wind direction, radar-derived wind speeds show underestimates of 4–6% compared with radiosonde values over the height range 4–14 km. Studies of the characteristics of this discrepancy in wind speeds have concentrated on its directional dependence, the effects of the spatial separation of the two sets of measurements, and the influence of any uncertainty in the radar measurements of vertical velocities. The aspect sensitivity of radar echoes has previously been suggested as a cause of underestimates of wind speeds by VHF radar. The present statistical treatment and case-studies show that an appropriate correction can be applied using estimates of the effective radar beam angle derived from a comparison of echo powers at zenith angles of 4.2° and 8.5°.