2-D and 3-D processing and interpretation of multi-fold ground penetrating radar data:a case history from an archaeological site

A 2.5-D and 3-D multi-fold GPR survey was carried out in the Archaeological Park of Aquileia (northern Italy). The primary objective of the study was the identification of targets of potential archaeological interest in an area designated by local archaeological authorities. The second geophysical objective was to test 2-D and 3-D multi-fold methods and to study localised targets of unknown shape and dimensions in hostile soil conditions. Several portions of the acquisition grid were processed in common offset (CO), common shot (CSG) and common mid point (CMP) geometry. An 8×8 m area was studied with orthogonal CMPs thus achieving a 3-D subsurface coverage with azimuthal range limited to two normal components. Coherent noise components were identified in the pre-stack domain and removed by means of FK filtering of CMP records. Stack velocities were obtained from conventional velocity analysis and azimuthal velocity analysis of 3-D pre-stack gathers. Two major discontinuities were identified in the area of study. The deeper one most probably coincides with the paleosol at the base of the layer associated with activities of man in the area in the last 2500 years. This interpretation is in agreement with the results obtained from nearby cores and excavations. The shallow discontinuity is observed in a part of the investigated area and it shows local interruptions with a linear distribution on the grid. Such interruptions may correspond to buried targets of archaeological interest. The prominent enhancement of the subsurface images obtained by means of multi-fold techniques, compared with the relatively poor quality of the conventional single-fold georadar sections, indicates that multi-fold methods are well suited for the application to high resolution studies in archaeology.

Radar interferometer measurements of space debris using the Evpatoria RT-70 transmitter

The ability of the Evpatoria RT-70 radar complex to perform research on space debris was investigated in four trial experiments during 2001-2003. The echo-signals of 25 objects at geostationary, highly elliptical and medium-altitude orbits were recorded on magnetic tapes at radio telescopes in Russia, Italy, China and Poland. The multi-antenna system configuration gives potential to supplement the classic radar data with precise angular observations using the technique of Very Long Baseline Interferometry. The first stage of such processing was fulfilled by the correlator in N. Novgorod, Russia. The cross-correlation of transmitted and received signals was obtained for the 11 objects on the Evpatoria-Bear Lakes, Evpatoria-Urumqi and Evpatoria-Noto baselines. This activity also promoted developing the optical observations of geostationary objects, conducted for the improvement of the radar target ephemerides. (C) 2004 COSPAR. Published by Elsevier Ltd. All rights reserved.

Measurements of ocean wave and current field using dual polarized X-band radar

A new ocean wave and sea surface current monitoring system with horizontally-(HH) and vertically-(VV) polarized X-band radar was developed. Two experiments into the use of the radar system were carried out at two sites, respectively, for calibration process in Zhangzi Island of the Yellow Sea, and for validation in the Yellow Sea and South China Sea. Ocean wave parameters and sea surface current velocities were retrieved from the dual polarized radar image sequences based on an inverse method. The results obtained from dual-polarized radar data sets acquired in Zhangzi Island are compared with those from an ocean directional buoy. The results show that ocean wave parameters and sea surface current velocities retrieved from radar image sets are in a good agreement with those observed by the buoy. In particular, it has been found that the vertically-polarized radar is better than the horizontally-polarized radar in retrieving ocean wave parameters, especially in detecting the significant wave height below 1.0 m.

Remote sensing of ocean waves by along-track interferometric synthetic aperture radar

A new wave retrieval method for the Along-Track Interferometric Synthetic Aperture Radar (AT-InSAR) phase image is presented. The new algorithm, named parametric retrieval algorithm (PRA), uses the full nonlinear mapping relations. It differs from previous retrieval algorithms in that it does not require a priori information about the sea state or the wind vector from scatterometer data. Instead, it combines the observed AT-InSAR phase spectrum and assumed wind vector to estimate the wind sea spectrum. The method has been validated using several C-band and X-band HH-polarized AT-InSAR observations collocated with spectral buoy measurements. In this paper, X-band and C-band HH-polarized AT-InSAR phase images of ocean waves are first used to study AT-InSAR wave imaging fidelity. The resulting phase spectra are quantitatively compared with forward-mapped in situ directional wave spectra collocated with the AT-InSAR observations. Subsequently, we combine the parametric retrieval algorithm (PRA) with X-band and C-band HH-polarized AT-InSAR phase images to retrieve ocean wave spectra. The results show that the ocean wavelengths, wave directions, and significant wave heights estimated from the retrieved ocean wave spectra are in agreement with the buoy measurements.

A new model to estimate significant wave heights with ERS-1/2 scatterometer data

A new model is proposed to estimate the significant wave heights with ERS-1/2 scatterometer data. The results show that the relationship between wave parameters and radar backscattering cross section is similar to that between wind and the radar backscattering cross section. Therefore, the relationship between significant wave height and the radar backscattering cross section is established with a neural network algorithm, which is, if the average wave period is <= 7s, the root mean square of significant wave height retrieved from ERS-1/2 data is 0.51 m, or 0.72 m if it is >7s otherwise.

Wave parameters retrieved from QuikSCAT data

A new algorithm is proposed to estimate significant wave height from QuikSCAT scatterometer data. The results show that the relationship between wave parameters and the radar backscattering cross section is similar to that between wind and the radar backscattering cross section. Therefore, the relationship between significant wave height and the radar backscattering cross section is established with a neural network algorithm. If the average wave period is less than or equal to 7 s, the root mean square errors of the significant wave height retrieved from QuikSCAT data are 0.58 m for HH polarization (HH-pol) and 0.60 m for VV polarization (VV-pol). If the average wave period is greater than 7 s, the root mean square errors of the significant wave height retrieved from QuikSCAT data are 0.83 m (HH-pol) and 1.10 m (VV-pol), respectively.

A parametric method of retrieving ocean wave spectra from synthetic aperture radar images

A parametric method that extracts the ocean wave directional spectra from synthetic aperture radar (SAR) image is presented. The 180 degrees ambiguity of SAR image and the loss of information beyond the azimuthal cutoff can be overcome with this method. The ocean wave spectra can be obtained from SAR image directly by using iteration inversion mapping method with forward nonlinear mapping. Some numerical experiments have been made by using ERS-1 satellite SAR imagette data. The ocean wave direction retrieved from SAR imagette data is in agreement with the wind direction from the scatterometer data.

A new method for the estimation of oceanic mixed-layer depth using shipboard X-band radar images

Shipboard X-band radar images acquired on 24 June 2009 are used to study nonlinear internal wave characteristics in the northeastern South China Sea. The studied images show three nonlinear internal waves in a packet. A method based on the Radon Transform technique is introduced to calculate internal wave parameters such as the direction of propagation and internal wave velocity from backscatter images. Assuming that the ocean is a two-layer finite depth system, we can derive the mixed-layer depth by applying the internal wave velocity to the mixed-layer depth formula. Results show reasonably good agreement with in-situ thermistor chain and conductivity-temperature-depth data sets.