An Analysis of X-Band Calibrated Sea Clutter and Small Boat Reflectivity at Medium-to-Low Grazing Angles

The coherent temporal characteristics of medium-to-low grazing angle sea clutter and small boat reflectivity are considered for different radar waveforms under a range of environmental conditions and geometrical configurations. Accurate empirical modelling of sea clutter enables the inference of the local sea conditions from radar returns, pertinent for port safety and navigation. Understanding the dynamics and associated reflectivity of small boats, in addition to empirical sea clutter models, allows the development of advanced detection and tracking algorithms, which will improve the performance of surveillance and marine navigation radar against small boats. Work presented is based on the empirical analysis of data recorded with two calibrated, coherent, pulsed radar systems at X-band frequencies. Specifically, target echoes from small boats are included in the datasets and subsequent analysis.


INTRODUCTION
There is a growing need for accurate, real-time instrumentation of the sea surface for safe navigation of vessels in and around harbours and shipping lanes.A commercial product, wave and surface current monitoring system (WaMoS) [1], can be connected to a conventional X-band marine radar.WaMoS II processes the unfiltered sea clutter to estimate the wave and surface current parameters in near real-time.According to the manufacturers the instrumented range is 0.1-3 km.At X-band operating frequencies, the main scattering mechanism is that of Bragg scattering, associated with the resonant capillary waves [2].Capillary waves, in turn, are generated by the local near-surface wind and do not propagate beyond the breeze area.Only in a fully developed sea can the wave height be directly related to the present wind speed and can the significant wave height H s be accurately inferred from the average sea clutter reflectivity [3].In transient sea conditions, the best fit for the Georgia Institute of Technology (GIT) mean sea clutter reflectivity model [4] is found if the sea state is related to the local mean wind speed rather than H s [5,6].It can therefore be deduced that it is possible to also infer local wind conditions from X-band sea clutter.This article investigates the temporal characteristics of coherent sea clutter, with specific interest in the Doppler characteristics of sea clutter and the relationship thereof to the local wind and wave conditions (e.g., average wind speed, wind direction, and wind gusts).
The significant increase in sea clutter reflectivity in rough seas with strong winds, together with relatively small radar cross section (RCS) of small boats (e.g., yachts, ski boats, and rigid inflatable boats (RIBs)), has often been blamed for disasters at sea where large ships collided with these small boats [7].In certain cases, the marine radar could not discern the boat signature from the clutter, while in other cases there were too many false tracks established leading to the subsequent disabling of the automatic tracker.For safe navigation it is pertinent that the detection capabilities of marine radar in adverse conditions are improved.With the introduction of cheaper, solid-sate, coherent marine radar [8] a whole new class of coherent detection algorithms has become applicable to marine navigation radar, for example, [9][10][11].Theoretical and first-order empirical analysis suggest subclutter visibility [11].Little work has been presented on the performance of this class of detectors (often referred International Journal of Navigation and Observation to as asymptotically optimal) on measured data of small boats [12].No reference is made to the performance as a function of the type of manoeuvring of the boat and the influence of the boat on the surrounding sea surface.This article investigates the performance of the adaptive linear quadratic (ALQ) detector [11] under different boat manoeuvres for different types of RIBs, including the 4.2 m pencilduck type that is often used for watersport racing.This boat has a very small RCS, but since it has a racing speed of up to 40 kts the local disturbance of the sea surface often exceeds the boat RCS by up to 10 dB.With improved subclutter visibility, the problem arises that first detections are declared not only for small boats, but also for large birds such as seagulls, with a typical RCS of 0.01-0.1 m 2 [3], and angels (flocks of birds flying together).Effective algorithms to discard tracks established on birds have to be developed [13].
Typical scanning surveillance systems (including maritime radar) have to declare a detection using only a limited number of pulses.Due to the long decorrelation time of sea clutter [14] and the more often than not spiky amplitude statistics [15], detection is quite difficult due to the short dwell time.Persistent, ubiquitous surveillance has become a top priority internationally.Typical entities of interest range from small recreational watercraft to large tanker ships.One of the characteristics of such systems, for instance, AwareNet [16], is the ability to employ long dwell times at specific areas of interest through the utilization of multiple, electronically steered receiver channels.Improved discernibility of small boats with long dwell times is therefore investigated in this article.
Two sea clutter and boat reflectivity measurement trials were conducted in 2006 and 2007 on the south western coast of South Africa.The aim of these trials was firstly to record datasets of sea clutter returns at different frequencies, range resolutions, grazing angles, look angles, and environmental conditions to validate current sea clutter models.Secondly, the aim was to record boat reflectivity datasets for a number of small boats to investigate its detectability with open literature detectors that will hopefully lead to the development of improved detection algorithms for radar systems employing adaptive dwell times.
The layout of the article is as follows.Section 2 presents an overview of the two measurement trials.The results presented in this article were obtained from the analysis of the data recorded during these trials.A description is given of the different radar systems, the experimental setup, as well as the system and data integrity verification procedures.Section 3 investigates a subset of the sea clutter measurements, focusing on the amplitude statistics and temporal characteristics for fixed frequency and frequency agile waveforms.The RCS and temporal fluctuation of a variety of small boats are investigated in Section 4 for different manoeuvres.Of particular interest is the effect of the boat manoeuvring on the local sea surface and its subsequent reflectivity.Section 5 presents an analysis of the detectability of these small boats, seagulls, as well as angels.

OVERVIEW OF MEASUREMENT TRIALS
The first measurement trial was conducted with the Fynmeet dynamic RCS measurement facility (Figure 1) at the Overberg Test Range (OTB).The site provided azimuth coverage of 135 • predominantly up-swell with well-developed waves and with a significant variation in wind direction.Sea clutter at grazing angles 0.3-3 • were recorded.The second measurement trial was conducted with an experimental, monopulse, X-band radar (Figure 5) deployed on top of Signal Hill in Cape Town.This site provided azimuth coverage of 140 • from up-to cross-swell, but with only two predominant wind directions for the duration of the trial.The sea was more representative of open sea conditions.Sea clutter and littoral clutter at grazing angles 0.3-10 • were recorded.The experimental radar uses pulse compression to increase the system gain and subsequently yields extended range capabilities compared to Fynmeet.

Radar and experimental set-up
The radar was deployed at OTB at location 34 • 36 56.52 S, 20 • 17 17.46 E, 67 m above mean sea level (AMSL).The shortest distance to the coastline was 1.2 km due south.A plan overview of the deployment site is depicted in Figure 2. The important specifications of Fynmeet are listed in Table 1.
Local wind speed and direction (Figure 4(a)) were logged with two weather stations separated by 1 km.The local wave direction φ wave , significant wave height H s , maximum wave height H max and wave period T wave (Figure 4  with a directional recording wave buoy.The local wave structure is influenced greatly by significant weather patterns from the south west and further perturbed by the diffraction patterns due to the cape, southwest of the deployment site which is located in a small bay area with a sea bed depth varying between 10-30 m at ranges of 3-10 km.The ground truth tracks of the boats were estimated using a differential processing global positioning system (GPS) receiver.

System and data integrity verification
For absolute RCS calibration the response from a sphere suspended below a helicopter, tracked in range with a typical α − β tracker and in angle with a video centroid tracker, was measured and the calibration coefficient was calculated as where N is the number of pulses transmitted, M(n) the range gate with maximum return for the nth pulse, A the receiver attenuation, and σ cs the sphere RCS.A standard deviation of 1 dB was achieved.Daily stability verification measurements were done with a corner reflector, exhibiting variations in the order of about 1 dB across the measurement period.Linearity of the quadrature receiver channels were ascertained by the analysis of calibrated noise source, receiver noise, and blue sky measurements taken throughout the trial.This analysis included the estimation of channel skewness, kurtosis as well as the 2nd to 4th normalized intensity moments, I 2 − I 4 .The amplitude and phase imbalance of the quadrature channels were estimated as 0.03 dB and 1 • resulting in a negative Doppler image of ≤−41dBc.Additionally, there were also harmonically related spurious responses at a level of ≤ −50dBc.A 5 MHz leak-through signal was identified and removed from the data.The signal phase was nondeterministic and therefore the amplitude and phase were estimated from the dataset itself.Sea clutter with a strong steady state component biases this estimate and the best results were obtained by using a censored mean level technique in the estimation process.The percentage of the dataset censored was chosen such that the resultant estimate yielded the lowest variance.Applying the discrete fourier transform (DFT) to the corrected data, a 0 Hz frequency bin with a comparable power density to adjacent frequency bins was obtained without suppressing the steady state clutter response.

Trial summary
Sea clutter datasets were recorded over a period of 11 days, with 112 fixed frequency and 38 stepped frequency datasets centred at four transmit (Tx) frequencies over an azimuth angle range from 90 • N to 225 • N and a grazing angle range 0.3-3 • .The local weather pattern may be described as roughly following a 6-day cycle [17] as cold front systems pass by from the west to the east.Over the trial period, the average wind speed ranged between 1-20 kts, with a maximum gust of 40 kts.Wind direction spanned 360 • , but the high wind speeds were mainly from the south west.H s varied between 1-3.8 m, with a maximum recorded wave height of 7.31 m.The wave direction was roughly 180 • N, and slowly changed direction toward the end of the trial to 135 • N.
Even with a low wind speed of 1 kt, the significant wave height was ≥1 m.This is due to the strong incoming south westerly swell combined with the diffraction patterns of the close-by cape and the reduction in sea depth.The illuminated sea area can therefore not be defined as a fully developed sea [3] and the standard tables relating wind speed and wave height to sea state do not apply.In terms of the average wind speed, sea states 1 to a low 5 were observed.In terms of the wave height, sea states from a high 2 to a high 5 were observed.
A 5.7 m RIB, a glass fibre ski boat and a wooden chokka fishing vessel (Figure 4) were deployed on 4 days with conditions ranging from calm to rough seas.The boats sailed a number of manoeuvres at different ranges and azimuth angles.A total of 55 fixed frequency and 43 stepped frequency datasets were recorded centred at four Tx frequencies.
The recorded datasets have been made available to the international research community.For information on the available datasets and how to access these datasets, refer to http://www.csir.co.za/small boat detection/.

Radar and experimental set-up
The experimental, X-band, monopulse radar was deployed on Signal Hill at location 33  Bay and around Robben Island.The tracks of the instrumented boats (Figure 8) were estimated using a differential-processing GPS receiver.

System and data integrity verification
Due to the similarity of the two radar systems, similar system and data integrity verification process were followed as for the Fynmeet radar described in Section 2.1.2.The experimental radar employs matched filter pulse compression, where pulse compression codes c pc are designed to yield a specific pulse compression gain, sidelobe levels and blind range.In the calibration procedure the height and range of the helicopter carrying the calibration sphere over the sea were restricted.The above-mentioned restrictions result in not all codes being calibrated.It is possible however, to estimate C cal for the uncalibrated codes from the calibrated codes by adding the relative pulse compression gain for the uncalibrated code whereN is the uncalibrated and K the calibrated code lengths.Equation ( 2) is valid for a matched filter with unit noise gain.Similarly, the Doppler processing gain can be estimated for the uncalibrated Doppler processing coefficients c dp , where M is the uncalibrated and L the calibrated Doppler processing coefficient lengths.For the experimental radar, the image rejection and spurious response is sufficiently below the noise floor.Pulse compression codes utilized during the measurement trial yielded range sidelobe levels in the order of −35 dB.

Trial summary
Sea clutter datasets were recorded on eight different days over a period of thirteen days.The predominant wind direction was northwestern, but with southeastern intervals.The average wind speed varied between 0 kts and 40 kts, with a maximum gust of 60 kts.The significant wave height ranged in 1-4.5 m, whilst the swell direction varied between 230 • N and 270 • N.
Datasets of the instrumented boats depicted in Figure 8 were recorded on five different days.
The Nadine Gordimer is a 10 m Class A deep sea rescue vessel with two MTU 1000 turbo diesel inboard motors and a range of communication antennas.The Rotary Endeavour is a Class 3 5.5 m RIB with two 60 hp Yamaha outboard motors and a single VHF antenna.The South African National Parks (SANParks) RIB is a 4.8 m RIB with a 60 hp Yamaha outboard motor.The 4.2 m pencilduck has a single 50 hp outboard motor with no antennas.In addition, datasets were recorded for a large variety of noncooperative boats of opportunity.Recordings were made using a range of fixed frequency and stepped frequency waveforms.

SEA CLUTTER ANALYSIS
Various statistical properties are evaluated in this section for seven sea clutter datasets recorded during the OTB 2006 measurement trial.These datasets represent low and high sea states at grazing angles 1 • and 0.5 • for a single Tx frequency and pulse widths of 100 nanoseconds and 300 nanoseconds.The range-time intensity plot for the high sea state, 1 • , 100 nanoseconds dataset CFC16-001 is presented in Figure 9.
The strong underlying modulation caused by the well developed waves is clearly visible for this up-swell configuration dataset.For cross-swell configurations and further ranges coupling between the waves and the underlying modulation becomes less pronounced as multiple waves are contained within a resolution cell, which is defined by the azimuth beamwidth and radar range resolution at low grazing angles.The underlying modulation also becomes less pronounced in weaker swell conditions.

Mean reflectivity and amplitude statistics
The mean reflectivity σ 0 and clutter-to-noise ratio (CNR) for the different OTB 2006 datasets are tabulated in Table 2 and compared to the GIT and the hybrid (HYB) models [4].In both models the sea state S was derived from the mean wind speed v wind using the empirical relation From Table 2 it can be concluded that there is good agreement between empirical σ 0 and the HYB model, with values generally between that of the HYB and the GIT.The GIT typically underestimates σ 0 at low grazing angles for low sea states by a significant margin, as discussed in detail in [4].Of particular interest is the good fit found by matching    The 2nd to 4th normalised intensity moments I 2 − I 4 and the estimated shape parameter ν are tabulated in Table 3, assuming a k-distributed envelope process and compared to the shape parameter model ν model [15].In this estimation, the theoretical relationship between the actual shape parameter ν and the effective shape parameter ν eff in the presence of noise is used [18]: ( The non-Rayleigh envelope statistics is evident.Theoretically spikiness should increase with a decrease in grazing angle.However, this is contradicted in the empirical analysis where a decrease in spikiness is observed with a decrease in grazing angle.This may be due to well-developed waves located closer to the radar (3 km) at the high grazing angles yielding increased spikiness, with less developed waves at the far-out ranges (8 km) where the low-grazing-angle sea clutter data was recorded.
These observations are indicative of the highly complex scattering environment and illustrate a still incomplete understanding of sea clutter.

Average doppler characteristics
As the capillary waves are the main scattering mechanisms at X-band and they are directly influenced by the nearsurface local wind [2], it can be expected that the Doppler and autocorrelation properties of sea clutter are directly influenced by the local wind.The sea clutter speckle autocorrelation r(τ) [14] is plotted in Figure 10 for four different datasets.From the magnitude response it is evident that the speckle decorrelation time is 10-20 milliseconds, which is consistent with literature [14].It also indicates that the decorrelation time is affected by sea state, where the decorrelation time decreases as the sea state (roughness of the sea) increases.An explanation for this may be the more rapid deformation of the capillary waves in rough seas.The complex autocorrelation is strongly coupled to the Doppler characteristics of the sea clutter.Evaluation of the real and imaginary components of r(τ) reveals that the second zerocrossing of I{r(τ)} approximates 1/2 of the mean projected Doppler period, This together with the empirical model [3] f d φ wind ≈ 2v wind cos φ wind f 0 4c (7) enables the estimation of the local projected wind speed from an analysis of the estimated autocorrelation.With a complete azimuth scan of the radar it would be possible to infer both wind speed and direction.This lies outside the scope of this paper and will be the subject of future research.of Doppler frequency, yielding an indication of the spikiness per Doppler resolution cell.The spectrogram in Figure 11 is representative of an up-swell sea (H s = 3.4 m) with strongly developed waves and up-wind (v wind = 15 kts) configuration at a range of 3.8 km.From the spectrogram the different individual waves can be distinguished as they are propagating through the given range cell.The different individual waves have very different Doppler spectra, which results in a significantly raised I 2 at the Doppler velocities associated with localised wind gusts.Thus in the Doppler domain these echoes will compete with those of real targets and hence may have an adverse effect on the false alarm rate.Since the individual waves are resolved, I 2 is also higher than 2 at the mean Doppler frequency.The spectrogram in Figure 12 is representative of a 70 • cross-swell sea (H s = 2.8 m) and down-wind (v wind = 15 kts) configuration at a range of 5.3 km.The sea was more representative of open sea conditions.From the spectrogram it is clear that the shorttime Doppler spectrum is much more homogeneous and it is impossible to distinguish individual waves or events.At the mean Doppler frequency I 2 tends to the theoretical value of 2 and is only slightly raised at the average Doppler spectrum edges.The spectrogram in Figure 13 is representative of a 20 • up-swell sea (H s = 2.5 m) and up-wind (v wind = 7 kts) configuration at a range of 5.6 km.Once again the sea was representative of open sea conditions.For this dataset the short-time Doppler spectrum is inhomogeneous compared to the previous dataset, but not as severe as the first dataset analysed in this subsection.For this up-swell configuration there is evidence of individual waves propagating through the range cell, but it is clear that more than one wave are contained within the resolution cell.The events associated with the broadened Doppler response may be associated with whitecaps blown off the top of the waves by the higher wind and/or gusts.The spikiness at the Doppler frequencies associated with the local maximum wind speed is confirmed by a significant raise in I 2 .

Spectrally inhomogeneous sea clutter
This brief analysis of the sea clutter Doppler spectrum and I 2 ( f d ) suggests that it is possible to also infer the existence and severity of whitecaps from the sea clutter.

Frequency agility
It is generally accepted that sea clutter speckle decorrelates with frequency agility when the frequency step size exceeds the pulse bandwidth, Δ f c ≥ B [3].The correlation coefficient ρ( f 0 , f n ) is plotted for a coherent processing interval (CPI) of 100 milliseconds at a fixed range cell over a period of 60 seconds-depicting the correlation between the base Tx frequency f 0 and an offset frequency of up to f 0 + 130 MHz for a pulse bandwidth of 10 MHz.From Figure 14 it can be concluded that in general the sea clutter speckle decorrelates whenever Δ f c ≥ B. However, there are a number of CPIs where the speckle only decorrelates after a step size of 40 MHz.Discrete spike events can also be identified where there is strong correlation for a step size of up to 130 MHz.Most radar detection mechanisms will declare these spikes as targets.It is important to note that these discrete spike events have a typical lifetime of 0.5-2 seconds.
Thus, overall these observations show broad agreement with those reported elsewhere and with the GIT and HYB models.They also re-emphasise the complexity of the scattering environment in which a target is required to be detected.This aspect is examined further in the following two sections.

SMALL BOAT REFLECTIVITY ANALYSIS
This section presents the results of the analysis of a range of small instrumented boats deployed during the two trials detailed in Section 2.

Small boat RCS and amplitude statistics
The mean RCS of the small boats, averaged over aspect angle, have been estimated from the measured data and tabulated in Table 4. Isolation of the boat signature from the sea clutter was obtained by Doppler filtering, using the GPS-estimated Doppler frequency as input.The responses of the threerange cells closest to the GPS range were coherently added to counter range-gate straddling losses.The effects of multipath fading from smooth sea surfaces and of shadowing of the  boat in higher sea states on the RCS values were not corrected for in the results presented in Table 4.
From Table 4 it is clear that the RCS of an RIB is related to its physical size and ranges between 1-5 m 2 .The RCS of solid boats are larger in general, with the mean RCS of the chokka fishing vessel up to 16 m 2 .In addition to the mean RCS of the small boat it is also critical for detection performance calculation to have a good model for its RCS fluctuation.The isolated boat RCS signature for the chokka fishing vessel is plotted in Figure 15, with the empirical probability density function (PDF) plotted for both the chokka fishing vessel and the WaveRider RIB in Figure 16.
Slight differences between the PDFs of the two boats are visible as well as the inability of the Swerling models to accurately describe their amplitude distributions.A characteristic of small boats in heavy sea is the fading of the RCS as the boat steers into the troughs of the waves, as indicated by the close-up view about 25-29 seconds in Figure 15.This suggests a strong correlation of boat RCS with the local sea waves and explains the poor fit of the empirical PDFs to the Swerling models at low RCS values.The noncoherent autocorrelation of a boat steering directly into the waves (Figure 17) shows periodicity with the same period as the mean sea wave period.This correlation of boat RCS and sea clutter echo strength further complicates an already challenging detection problem.

Small boat doppler bandwidth
The Doppler bandwidth of a target is an important design parameter for optimal coherent detection.Coherent processing gain is only achieved by an increase in the CPI whilst it still remains less than or equal to the inverse of the target Doppler bandwidth.High Doppler resolution spectrograms have been computed for the different small boats, with the results for CPI's of 10 milliseconds and 200 milliseconds plotted in Figure 18 for the WaveRider RIB.
The Doppler bandwidth for the WaveRider RIB can be observed from the high-resolution spectrogram.As for most of the other small boats, this is approximately 10 Hz.This yields an optimal CPI ranging between 100-200 milliseconds.As the CPI increases beyond this value, the boat energy will start to spread over multiple Doppler resolution cells, yielding no additional coherent processing gain.This is a key consideration in the design of an optimal detector.

Highly manoeuverable small boats
Due to its small size, light weight and powerful engines, the RIB class of small boats is highly manoeuvrable with the ability to reach speeds of up to 40 kts for even the small 4.2 m pencilduck.Especially in the design of coherent detection and tracking algorithms, it is important to have a good understanding of the anticipated manoeuvrability of the target.In addition, the disturbance of the manoeuvring boat on the local sea surface may also greatly influence its detectability either adversely or positively.The high Doppler resolution spectrograms of two RIB's are plotted in Figure 19-21 for three different manoeuvres.
The narrow Doppler response of the drifting pencilduck (Figure 19) is evident, as well as the slight movement of the pencilduck due to the local waves.The drifting pencilduck caused little disturbance on the local sea surface.The WaveRider RIB steering radially outbound at a speed of about 10 kts (Figure 20) had a narrow Doppler response, with a local disturbance of the sea surface visible when the RIB was crashing through the crests of the waves.This local disturbance is observed as quite broad Doppler bandwidth noise with Doppler velocities ranging from slightly higher than the speed of the boat down to the Doppler velocity of the local sea clutter speckle.The spectral density of the disturbance is 20 dB lower than the boat signature.The pencilduck racing at 40 kts radially outbound (Figure 21) still had a narrow Doppler response for the body of the boat, but caused a significant local disturbance of the sea water (e.g., splashing waves and water spray by the propeller), decreasing the localised signal-to-interference ratio (SIR) to less than −10 dB.With such a low SIR, it becomes increasingly difficult to detect the boat with clutter suppression algorithms, even though the local disturbance of the sea surface may be detected by a basic envelope thresholding detector.However, there is still ample Doppler separation between the boat and interference, and in principal a long dwell time range-Doppler detector could be constructed that will consistently declare detections for this fast moving boat.Also of interest is the case where the boat is racing crossrange.This still yields strong self-induced interference, but the Doppler response of the body of the boat will be buried within the interference and it will become extremely difficult to detect.From this subsection it can be concluded that the exact manoeuvre of the small boat has a great influence on its detectability, especially due to its potential disturbance to the local sea surface.It is also clear that not only the speed, but also the heading of small boats has to be modelled for accurate performance prediction.

Frequency agility correlation for small boats
Theoretically the correlation coefficient ρ( f 0 , f n ) for the RCS of a point scatterer for different frequencies should be unity.For a target that can be approximated as a point scatterer it is assumed that |ρ( f 0 , f n )| → 1.In the presence of clutter and multipath fading, the correlation will be adversely effected.ρ( f 0 , f n ) is plotted in Figure 22 for a CPI of 100 milliseconds for the two range cells containing most of the energy for the chokka fishing vessel.
The frequency agility correlation for the boat exceeded 0.5 with Δ f c = 100 MHz for 84% of the total time period, compared to only 10% for sea clutter only as represented in Figure 14.There were time periods when the correlation coefficient dropped to similar levels as the sea clutter, which coincided with low levels of SCR and/or SNR.For a 10 MHz pulse bandwidth empirical evidence suggest that small boats exhibit significant frequency agility correlation for frequency step sizes up to and beyond 130 MHz.It is possible to design a detection algorithm that uses ρ( f 0 , f n ) as the basis for its test statistics.This detector may still detect discrete spikes, since they have similar characteristics as the small boats (e.g., at t ∈ (6, 51) s in Figure 14).Once again, these discrete spikes only have a limited lifetime of less than 2 seconds.

DETECTABILITY OF SMALL BOATS
This section presents the detectability of small boats undergoing different manoeuvres using the ALQ detector as an example of the asymptotically optimal class of detectors.

Overview of the ALQ detector
The ALQ detector is designed by extending the generalized likelihood ratio test approach, as suggested by Kelly [9] for Gaussian interference, to the spherically invariant random process model for non-Gaussian interference [19].Assume that the radar transmits a coherent train of m pulses.The associated m received complex samples can be constructed as a vector z = [z(1), . . ., z(m)] T .Under the assumption that M is known exactly, the ALQ detector can be expressed mathematically as where p is the steering vector typically constructed with elements p i = e j2πi fdT [20], T the radar PRI and f d the target Doppler frequency [11].It is generally accepted that M is highly dependent on the radar configuration, geometry, and the environmental conditions and has to be estimated from adjacent range gates that are not contaminated by the boat itself.Various estimation techniques have been proposed [21].Gini and Greco [11] describe one such technique that makes a good compromise between detection losses and hardware processing requirements: for i = 0, 1, 2, . . ., N it .During each iteration the approximately maximum likelihood (AML) estimation is normalized such that its trace is equal to m.Since the ALQ detector involves inversion of M, care has to be taken to ensure that the matrix does not become singular.This can be ensured by setting the number of independent sea clutter time vectors at different range gates k equal to at least the length of the test vector m, k ≥ m [9].Detectability can be improved by increasing this ratio, but at the expense of increased hardware processing requirements.

ALQ performance for different RIB manoeuvres
In the first dataset evaluated, the 4.2 m pencilduck was floating close the southwestern shore of Robben Island, with H s = 3 m and v wind = 6 kts NE.The radar look angle was 343 • N at range R = 11 km with grazing angle θ = 1.5 • .SCR and CNR were 6 dB and 24 dB, respectively.Figure 23( of 10 kts at range R = 3.3 km.The radar look angle was 166 • N, whilst H s = 3.2 m and v wind = 16 kts SSE.SCR and CNR were 4 dB and 17 dB, respectively.The performance of the ALQ detector for P FA = 10 −4 is plotted in Figure 24, with E{P d } = 62%.Even though the SCR is lower than in the previous case, a significant increase in P d is observed.This is most probably due to the increased Doppler separation. Figure 24(c) indicates that low P d is in general associated with low SCR and/or low SNR.
With the observed increase in sensitivity with an increase in Doppler separation, it can be expected that detectability will increase even further for the high-speed pencilduck.The performance of the ALQ detector is plotted in Figure 25 for the pencilduck racing at 40 kts radially outbound at a range of 21.5 km.Even though the boat can be distinguished in the high Doppler resolution spectrogram (Figure 21) the ALQ detector only manages very intermittent detections.The sea clutter and localised disturbance are whitened over all Doppler, effectively masking the boat.In this case the ALQ cannot be classified as an asymptotically optimal, since a range-Doppler cell-averaging CFAR detector can be configured to steadily detect the boat due to the separation in Doppler of the interference and the boat signature and its narrow Doppler spectrum.

Detection of seagulls and angels
With improved subclutter visibility, the problem arises that first detections are declared not only for small boats, but also for large birds such as seagulls.ALQ detections for the last dataset in Section 5.2 are overlaid on the rangetime intensity plot in Figure 26.Even though only intermittent detections were declared for the racing pencilduck, a large number of detections were declared for the entire dataset.Examination of these detections reveals that they coincide with scatterers yielding an RCS of approximately 0.01-0.1 m 2 and a narrow Doppler spectrum with very fast acceleration.Comparing this to the Doppler signatures of birds in [13] and observation by the copilot of the pencilduck led to the conclusion that these scatterers are indeed seagulls.
Of particular interest in Figure 26 is the consistent detections declared from a range of 21.9 km at t = 0 second closing in to a range of 21.4 km.Closer inspection revealed that this was a flock of about 6 seagulls flying in formation.The combined RCS of this "angel" was 10 dB lower than the pencilduck, but still yielded significantly higher P d .The main reason for this is the Doppler separation and that the seagulls caused no local disturbance of the sea surface.The resultant SIR for the angel was indeed higher than for the pencilduck.As radar sensitivity is increased, this will become a more and more significant problem.

CONCLUSIONS
Current commercial products provide near real-time estimation of basic wave and surface current parameters using the video output of standard X-band marine radar.This paper investigated sea clutter and small boat reflectivity in the littoral and proved that sea clutter reflectivity is related to v wind rather than H s .Temporal characteristics of sea clutter were investigated, with empirical results suggesting that v wind and φ wind can be estimated from the sea clutter speckle autocorrelation.The spectral inhomogeneity of sea clutter was investigated for different sea conditions.The brief analysis of the sea clutter Doppler spectrum and I 2 ( f d ) suggested that it is possible to also infer the existence and severity of whitecaps from the sea clutter.Discrete spikes in sea clutter were clearly visible when the frequency agility decorrelation was estimated.
For safe navigation, it is pertinent that the detection capabilities of marine radar in adverse conditions are improved, especially for small boats.This requires an in-depth understanding of the dynamics and associated reflectivity of these boats.The absolute RCS, amplitude statistics, and temporal characteristics of a range of small boats have been analysed using a comprehensive set of recorded datasets.Of particular interest were the dependency of the boat reflectivity on the local sea, deviation from the Swerling RCS models, the perceived persistence of reflectivity for short periods of time and the distinguishable pulse-to-pulse frequency agility correlation properties of small boats.It was shown using real data that the ALQ detector can, under certain conditions, be subject to selfmasking.A definite contribution to the knowledgebase is the importance of not only modelling the sea clutter and boat reflectivity accurately, but also to model the local disturbance caused by small boats, especially during fast manoeuvring.

Figure 2 :
Figure 2: Radar deployed on Signal Hill with open view of sea.

Figure 6 :
Figure 6: Plan overview of radar deployment site.

Table 1 :
Fynmeet system and performance specifications.

Table 3 :
Empirical versus model amplitude statistics.
sea state to local wind speed rather than H s in transient sea conditions (not fully developed).
Section 3.2 provided empirical evidence that the average wind speed can be inferred from the sea clutter speckle autocorrelation, which is strongly correlated to the aver- age Doppler response thereof.Experimental data suggests however that the sea clutter spectrum is inhomogeneous in both range and time in general.High Doppler resolution spectrograms of three different geometrical configurations and environmental conditions are presented in Figures11-13.In addition to the spectrograms I 2 is plotted as a function