Measurement of ocean wave spectra using a ship-mounted HF radar

The work describes an inversion algorithm for HF radar measurement of nondirectional wave spectra using an omnidirectional receive/transmit antenna. Such a radar would be suitable for deployment on a stationary ship or drill rig. In this approach, wave information is extracted from the radar observations by numerically inverting the integral equation representing the backscatter return from the ocean. Test results of this technique applied to data collected using a 25.4-MHz radar installed on a ship have been very positive. For the two measurements collected, there is a high degree of correlation between the radar wave estimates and those of a WAVE-TRACK buoy     

HF radar ocean current algorithm based on MUSIC and the validation experiments

Wuhan University's ocean state measuring and analyzing radar (OSMAR2000), working at around 7.5 MHz in the low region of the HF band with a 120-m-long linear receiving antenna array, can measure ocean surface current at ranges of up to 200 km. An ocean surface current algorithm based on direction finding (DF) using the multiple signal classification (MUSIC) method is developed for the OSMAR2000 radar. This paper describes the OSMAR2000 ocean surface current algorithm based on MUSIC and the validation experiments in the East China Sea. The results of the ocean surface current measurements demonstrate that the OSMAR2000 ocean surface current algorithm based on MUSIC is feasible for the long range of ocean surface current mapping with a sufficient bearing resolution.     

高频地波雷达生成海洋表面矢量流图

武汉大学研制的双站高频地波雷达系统OSMAR2000利用测得的两幅单站径向海流图生成矢量海流图。经典矢量流图生成方法不能直接应用到OSMAR2000系统中。本文提出一种先在极坐标系下用自然三次样条函数将径向流插值到公共网格上然后直接进行矢量合成的矢量海流图生成方法。OSMAR2000在东海的表面矢量流实测结果与作对比验证的传统海流计测量结果十分吻合。对比数据表明，该方法是可行的，且优于先进行径向流线性插值后矢量合成的矢量流图生成方法。这也是国内首次利用高频地波雷达实现海洋表面矢量流的实时监测。     

Surface current measurements by HF radar in freshwater lakes

HF radar has become an increasingly important tool for mapping surface currents in the coastal ocean. However, the limited range, due to much higher propagation loss and smaller wave heights (relative to the saltwater ocean), has discouraged HF radar use over fresh water, Nevertheless, the potential usefulness of HF radar in measuring circulation patterns in freshwater lakes has stimulated pilot experiments to explore HF radar capabilities over fresh water. The Episodic Events Great Lakes Experiment (EEGLE), which studied the impact of intermittent strong wind events on the resuspension of pollutants from lake-bottom sediments, provided an excellent venue for a pilot experiment. A Multifrequency Coastal HF Radar (MCR) was deployed for 10 days at two sites on the shore of Lake Michigan near St. Joseph, MI. Similarly, a single-frequency CODAR SeaSonde instrument was deployed on the California shore of Lake Tahoe. These two experiments showed that when sufficiently strong surface winds (2 about 7 m/s) exist for an hour or more, a single HE radar can be effective in measuring the radial component of surface currents out to ranges of 10-15 km. We also show the effectiveness of using HF radar in concert with acoustic Doppler current profilers (ADCPs) for measuring a radial component of the current profile to depths as shallow as 50 cm and thus potentially extending the vertical coverage of an ADCP array     

Directional wave information from the SeaSonde

This paper describes methods used for the derivation of wave information from SeaSonde data, and gives examples of their application to measured data. The SeaSonde is a compact high-frequency (HF) radar system operated from the coast or offshore platform to produce current velocity maps and local estimates of the directional wave spectrum. Two methods are described to obtain wave information from the second-order radar spectrum: integral inversion and fitting with a model of the ocean wave spectrum. We describe results from both standard- and long-range systems and include comparisons with simultaneous measurements from an S4 current meter. Due to general properties of the radar spectrum common to all HF radar systems, existing interpretation methods fail when the waveheight exceeds a limiting value defined by the radar frequency. As a result, standard- and long-range SeaSondes provide wave information for different wave height conditions because of their differing radar frequencies. Standard-range SeaSondes are useful for low and moderate waveheights, whereas long-range systems with lower transmit frequencies provide information when the waves are high. We propose a low-cost low-power system, to be used exclusively for local wave measurements, which would be capable of switching transmit frequency when the waveheight exceeds the critical limit, thereby allowing observation of waves throughout the waveheight range.     

A minimum variance unbiased estimator for ocean surface currentsbased on multifrequency microwave radar observations

This paper presents a minimum variance unbiased (MVU) estimator for estimating an ocean surface current using the multifrequency microwave radar technique. In this technique the current information is obtained by finding the difference between the measured phase velocities of some specific surface gravity waves and the phase velocities calculated from the dispersion relation for still water. By defining the problem as a linear estimation problem, we develop an unbiased estimator for the current component along the radar look direction, which has a variance that is inversely proportional to the sum of the squared wavenumbers of the gravity waves used in the measurements. We also study the performance of an MVU vector estimator based on radar observations along two directions. Our analysis shows that the confidence region of this estimator has the shape of an elongated ellipse with semi-axes and orientation which are dependent on the angle between the observation directions, but independent on the true current vector. Furthermore, the theoretical models are thoroughly tested using both simulated and real radar data, and these tests show very good agreement with the model predictions     

Algorithm for HF radar vector current measurements

A new algorithm is proposed, called the stream function method (SFM) for producing vector current maps from radial data measured by dual-site high frequency surface wave radar (HFSWR). In SFM, a scalar stream function is constructed under some oceanographic assumptions. The function describes the two-dimensional (2-D) ocean surface water motion and is used to obtain the distribution of vector currents. The performance of SFM is evaluated using simulated radial data, which demonstrates that SFM has advantages over typical vectorial combination methods (VCM) both in error acceptance and robustness, and excels another method based on least-squares fitting (LSF) in recovering the complicated current models. Furthermore, SFM is capable of providing the total currents based on radials from single-site radar. We also test the assumptions of horizontal non-divergence in the simulation. The new algorithm is applied to the field experiment data of Wuhan University’s ocean state measuring and analyzing radar (OSMAR), collected in the coastal East China Sea during April 11–17, 2004. Quantitative comparisons are given between radar results by three current algorithms and in-situ current meter measurements. Preliminary analysis of the vertical current shear is given based on the current meter measurements.     

Phase Corrections of Small-Loop HF Radar System Receive Arrays With Ships of Opportunity

This paper is an extension of other work that addresses the use of radar echoes from ships of opportunity to determine the proper phase corrections for small-loop phased-array antennas used within high-frequency (HF) ground-wave radar systems. This technique also yields estimates for unknown ship bearings that (for cases where there is adequate signal-to-noise ratio of 20 dB or more) are consistent to within 2deg-3deg among measurements from independent radar frequencies. Within this paper, phase corrections gathered from actual ships of opportunity are compared to phase corrections gathered during a calibrated transponder run, in which the ship bearing is known. The phase corrections derived from the ship of opportunity presented in this paper were consistent with the known phase corrections to within 13.2deg (for the worst case). Furthermore, the estimates of the ship bearings collected from the two usable radar frequencies were consistent to within 1deg of each other     

高频地波雷达的发展与应用现状分析

高频地波雷达作为一种新兴的海洋探测设备,相比传统观测方式,具有大范围、全天候、低成本等优点,因而在世界各地得到了广泛使用,并在近海海洋环境监测中发挥了重要的作用。首先介绍了高频地波雷达探测海洋环境的基本原理,然后概述了阵列式和紧凑便携式两类高频地波雷达的国内外研究现状,接着介绍了高频地波雷达在海洋环境探测中的应用,最后分析了我国与世界发达国家在高频地波雷达海洋环境探测领域的差距并提出了改进建议。     

On the Development of a Second-Order Bistatic Radar Cross Section of the Ocean Surface: A High-Frequency Result for a Finite Scattering Patch

The development of a model for the second-order bistatic high-frequency (HF) radar cross section on an ocean surface patch remote from the transmitter and receiver is addressed. A new approach is taken that allows a direct comparison with existing monostatic cross sections for finite regions of the ocean surface. The derivation starts with a general expression for the bistatically received second-order electric field in which the scattering surface is assumed to be of small height and slope. The source field is taken to be that of a vertically polarized dipole, and it is assumed that the ocean surface can be described, as is usually done, by a Fourier series in which the coefficients are zero-mean Gaussian random variables. Subsequently, a bistatic cross section of the surface, normalized to patch area, is derived. The result is verified by the following two means: 1) the complete form of the bistatic HF radar cross section in backscattering case is shown to contain an earlier monostatic result that has, itself, been used extensively in radio oceanography applications; and 2) the bistatic electromagnetic coupling coefficient is shown to reduce exactly to the monostatic result when backscattering geometry is imposed. The model is also depicted and discussed based on simulated data     

A comparison of multifrequency HF radar and ADCP measurements ofnear-surface currents during COPE-3

A high-frequency multifrequency coastal radar operating at four frequencies between 4.8 and 21.8 MHz was used as part of the third Chesapeake Bay Outflow Plume Experiment (COPE-3) during October and November, 1997. The radar system surveyed the open ocean east of the coast and just south of the mouth of Chesapeake Bay from two sites separated by about 20 km. Measurements were taken once an hour, and the eastward and northward components of ocean currents were estimated at four depths ranging from about 0.5 m to 2.5 m below the surface for each location on a 2 by 2 km grid. Direction of arrival of the signals was estimated using the MUSIC algorithm. The radar measurements were compared to currents measured by several moored acoustic Doppler current profilers (ADCPs) with range bins 2-14 m below the water surface. The vertical structure of the current was examined by utilizing four different radar wavelengths, which respond to ocean currents at different depths, and by using several ADCP range bins separated by 1-m intervals. The radar and ADCP current estimates were highly correlated and showed similar depth behavior, and there was significant correlation between radar current estimates at different wavelengths and wind speed     

Temporal and spatial resolution of HF ocean radars

The spatial and temporal resolutions of the two main types of HF radar are compared, with reference to the phasedarray and the crossed-loop direction-finding systems which make up the Australian Coastal Ocean radar Network. Both genres use a swept frequency “chirp” modulation to define the range of a pixel being observed but the method for determining the azimuth direction of the pixel is a strong point of differentiation. The phased-array systems produce independent maps of surface currents in about 1/7 of the time for the crossed-loop systems because of contrasting noise performance of the antennas. The use of beam-forming analysis in the phased-arrays is shown to give spatial resolutions, for vector currents, of about 10 km close to the shore, and 25 km at ranges of 150 km. The corresponding vector current spatial resolutions for the crossed-loop systems are 40 km and 60 km respectively.     

Characterizing Observed Environmental Variability With HF Doppler Radar Surface Current Mappers and Acoustic Doppler Current Profilers: Environmental Variability in the Coastal Ocean

A network of high-frequency (HF) radars is deployed along the New Jersey coast providing synoptic current maps across the entire shelf. These data serve a variety of user groups from scientific research to Coast Guard search and rescue. In addition, model forecasts have been shown to improve with surface current assimilation. In all applications, there is a need for better definitions and assessment of the measurement uncertainty. During a summer coastal predictive skill experiment in 2001, an array of in situ current profilers was deployed near two HF radar sites, one long-range and one standard-range system. Comparison statistics were calculated between different vertical bins on the same current profiler, between different current profilers, and between the current profilers and the different HF radars. The velocity difference in the vertical and horizontal directions were then characterized using the observed root-mean-square (rms) differences. We further focused on two cases, one with relatively high vertical variability, and the second with relatively low vertical variability. Observed differences between the top bin of the current profiler and the HF radar were influenced by both system accuracy and the environment. Using the in situ current profilers, the environmental variability over scales based on the HF radar sampling was quantified. HF radar comparisons with the current profilers were on the same order as the observed environmental difference over the same scales, indicating that the environment has a significant influence on the observed differences. Velocity variability in the vertical and horizontal directions both contribute to these differences. When the potential effects of the vertical variability could be minimized, the remaining difference between the current profiler and the HF radar was similar to the measured horizontal velocity difference (~2.5 cm/s) and below the resolution of the raw radial data at the time of the deployment     

A Method of Swell-Wave Parameter Extraction From HF Ocean Surface Radar Spectra

A new method for the extraction of swell-wave parameters from high-frequency (HF) radar spectra is presented. The method of extraction of the parameters, period, direction, and height, relies on a frequency-modulation approach that describes the hydrodynamic interaction of the swell waves with the resonant, shorter, Bragg waves. The analysis process minimizes the electromagnetic second-order interaction and a simulation model was used to validate the approach. This simplified method provides a fast means of examining swell conditions over large areas of the ocean surface. Data are acquired using a pair of coastal ocean surface radar (COSRAD) systems deployed at Tweed Heads, Qld., Australia. The radar covers a sweep (approximately 60deg) every 30 min with spatial resolution of the order of 3 km. A sample set of data from this deployment is used in a case study to show the extraction of swell direction and amplitude using these methods. The results support the use of the COSRAD HF radar for mapping swell in the near-shore zone     

Comparison of inversion algorithms for HF radar wave measurements

All ocean wave components contribute to the second-order scattering of a high-frequency (HF) radio wave by the sea surface. It is therefore theoretically possible to estimate the ocean wave spectrum from the radar backscatter. To extract the wave information, it is necessary to solve the nonlinear integral equation that describes the relationship between the backscatter spectrum and the ocean wave directional spectrum. Different inversion techniques have been developed for this problem by different researchers, but there is at present no accepted “best” method. This paper gives an assessment of the current status of two methods for deriving sea-state information from HF radar observations of the sea surface. The methods are applied to simulated data and to an experimental data set with sea-truth being provided by a directional wave buoy     

基于BP神经网络的高频地波雷达海流空间插值

高频地波雷达是海洋环境监测的重要手段,当前已经实现对海流的业务化观测,但是外部因素常引起海流空间探测的不连续性。为解决此问题,尽量保障区域数据的完整性和准确性,本文将BP神经网络技术与空间插值相结合,建立了海流的BP神经网络插值模型,并进行了针对实测数据的缺失插值仿真,通过与反距离权重法和线性插值法插值结果的对比,分析该模型在区域海流大面积缺失、流速整体较大和流速整体较小3个方面的性能。结果表明,BP神经网络插值模型的海流预测效果明显优于其他两种方法,且在流场数据大范围缺失下也取得了良好的效果。     

Measurement of Ocean Surface Currents by the CRL HF Ocean Surface Radar of FMCW Type. Part 2. Current Vector

The Communications Research Laboratory (CRL) has been developing high-frequency ocean surface radars (HFOSRs). The CRL dual-site HFOSR system can clarify the distribution of surface currents with a nominal range of 50 km. This paper presents a theoretical and experimental analysis of the measurement error of the current vector obtained by the CRL HFOSR system, using a comparison of instantaneous current vectors acquired by the HFOSR system and current meters moored at a depth of 2 m, taking account of the vertical current shear. The theoretical analysis shows that the probability distribution of the measurement error of the current vector forms concentric ellipses at a spatial scale that depends on the RMS measurement error of radial current velocity and with an aspect ratio that depends only on the azimuthal difference of the radar beams. When the azimuthal difference is a right angle, the measurement error of the current vector is at a minimum. A comparison between instantaneous current vectors measured by the CRL HFOSR system and moored current meters shows that the distribution of the difference vector between the radar current and the meter current agrees well with the theoretical measurement error of the current vector and that the RMS of difference vector length is about 10 cm s^–1 while the azimuthal difference between two radar beams is between 45 and 135 degrees. The accuracy of current measurement by the dual-site HFOSR system is therefore considered to be less than 10 cm s^–1 in this range of azimuthal difference. The theoretical analysis will be applicable for a wider range of the azimuthal difference of the radar beams.     

Line broadening on HF ocean surface radar backscatter spectra

Experimental results from an array of moored current meters and an HF ocean surface radar support the idea that line broadening on the radar spectra is caused by the velocity distribution within the radar target cell. The experiment was done in the wake of a small island where the velocity variations were severe. An estimate is made of the line broadening which can be expected. In a turbulent flow with dissipation rate of the orderepsilon sim 10^{-10}m^{2}s^{-3}and target cell size 1 3000 m, the line broadening isDeltaf sim 10^{-3}Hz. This would be resolved with a radar time series ofsim 20min and indicates that the HF ocean surface radar technique has potential in the observation of surface velocity distributions.     

HF radar detection of tsunamis

This paper demonstrates that HF radar systems can be used to detect tsunamis well before their arrival at a coastline. We solve the equations of motion and continuity on the ocean surface using models to simulate the signals produced by a tsunami approaching the east U.S. coast. Height and velocity profiles are derived along with expressions for the radar-observed current velocities in terms of bathymetry and tsunami height and period. Simulated tsunami-generated radial current velocities are superimposed on typical maps of radial velocity generated by a Rutgers University HF radar system. A detection parameter is defined and plotted to quantify the progress of the tsunami, which is shown to be detectable well before its arrival at the coast. We describe observations/warnings that would have been provided by HF radar systems at locations in the path of the 2004 Indian Ocean tsunami.     

Ship detection with high-resolution HF skywave radar

This paper presents an overview of ship detection by high-frequency (HF) skywave backscatter over-the-horizon radar (OTHR). Ships have been detected at ranges of 2000 km or more by OTHR that uses sufficient resolution in the radar spatial and Doppler frequency domains. The HF sea-echo Doppler spectrum limits the target signal-to-clutter ratio (SCR), as a function of the ocean wave-height distribution, wind direction, radio frequency, and ship target radial velocity. Maximum sea-clutter spectrum purity, and hence larger SCR, is achieved with the use of stable single-mode ionospheric propagation. Real-time measurement and interpretation of ionospheric propagation features therefore must guide the choice of OTHR operating frequency. Experimental data recorded at the ONR/SR1 Wide Aperture Research Facility (WARF) bistatic OTHR in central California demonstrate reliable ship detection in the Northeast Pacific Ocean. WARF transmits 1-MW average effective radiated power, using a linear frequency-modulated continuous-wave (FMCW) waveform, and receives with a 2.55-km broadside array of vertical monopole element pairs. Swept bandwidths as high as 200 kHz have been used. Sufficient spectral resolution is achieved with a coherent integration time (CIT) of 12.8 s. Longer CIT, and autoregressive (AR) spectral analysis techniques such as Marple's algorithm, have been used to improve Doppler resolution.