Institute of Oceanography

University of Hamburg

 | search | contact | imprint
Deutsche version
Institute of Oceanography
University of Hamburg
Bundesstraße 53
D-20146 Hamburg
Tel.: +49 40 42838-2605 / -5449
Fax: +49 40 42838-7488
E-Mail:  waltraut.domke-sommer(at)zmaw.de

SAR-signature interpretation of Oceanic Features

Spaceborne synthetic aperture radar (SAR) images are high resolution microwave radar images which are usually acquired from satellites at altitudes on the order of 800 km which orbit around the Earth within about 100 minutes. SAR images have widths between some tens of kilometers and 500 km and resolutions between meters and tens of meters. Well known satellites with SAR systems have been the American Seasat (1978), the European satellites ERS-1 (launched 1991), ERS-2 (launched 1995), and Envisat (launched 2002), and the Canadian Radarsat-1 (launched 1995).

Fig. 1: Envisat in space (artist's view) and in a thermal test facility at ESA-ESTEC
(Noordwijk, Netherlands). © ESA.

Since a radar is an active instrument and microwaves penetrate clouds, SAR images of the ocean can be acquired at any time of the day and under any weather conditions, which is a major advantage compared to passive instruments working at smaller wavelengths such as, for example, optical instruments or infrared radiometers. Since the SAR image intensity depends mainly on the sea surface roughness, one can derive information on surface winds and waves as well as on current features that modulate the surface roughness by hydrodynamic wave-current interaction.


a

b

c

d
Fig. 2: Examples of ERS-1 / ERS-2 SAR images (100 km × 100 km) showing signatures of (a) underwater bathymetric features off the coast of China, (b) oceanic internal waves at the Strait of Gibraltar, (c) the plume of the Kutai river, Borneo, (d) atmospheric convection cells and other features north of the Strait of Messina (Mediterranean Sea). © ESA.

The Remote Sensing group of the Institute of Oceanography has a long tradition in the analysis and interpretation of SAR signatures of oceanic features such as fronts, internal waves, and the underwater bathymetry in coastal waters, as well as wind features near the water surface such as atmospheric convection cells. We have participated in the development of theoretical models and algorithms as well as in field experiments and laboratory experiments within the framework of numerous national and international projects, and we have published many articles on these subjects.

Within the framework of MARSAIS we have implemented an algorithm for the retrieval of information on current fronts from their SAR signatures. The following series of 12 figures with explanations illustrates how such algorithms work.

 1  2  3  4  5  6  7  8  9 10 11 12

This is a 48 km × 48 km section of an ERS-1 SAR image of the mouth of the Rhine river at the Dutch coast, reduced in size by averaging over 8 × 8 pixels and contrast optimized. The image was acquired on 13 October 1993, 21:46 UTC (orbit 11740, frame 1035, © ESA). It exhibits clear signatures of a freshwater plume released by the Rhine river, which is floating on the seawater of the North Sea. The arrows indicate the location of the Rhine plume front, whose radar signatures are most pronounced in the northeast. The image intensity variations reflect variations in the surface roughness due to hydrodynamic modulation of waves by spatially varying currents.

 1  2  3  4  5  6  7  8  9 10 11 12

For a quantitative interpretation of the radar signatures, we need to know some radar parameters and environmental parameters. The frequency of the ERS-1 SAR is 5.3 GHz (C band), the polarization is VV, and the incidence angle is about 23°. The Rhine plume image was acquired during an ascending overpass, looking to the right, as indicated by the red arrows. The backscattered power detected by the SAR is basically proportional to the intensity of ocean waves with wavelengths of about 7 cm which propagate in and against this direction. The satellite's heading was 344° i.e., North is 16° clockwise from the vertical axis of the image.

 1  2  3  4  5  6  7  8  9 10 11 12

According to weather observations at Hoek van Holland, the wind speed was about 11 m/s from 210°. In principle, this information could also be retrieved from the SAR image itself, using a wind retrieval tool.

 1  2  3  4  5  6  7  8  9 10 11 12

The time of the SAR image acquisition was 2:07 hours before high water in Hoek van Holland. According to a current atlas, the tidal currents at this phase are quite small and converge in the test area, making the water level rise. At the Rhine plume front itself we expect some current shear and convergence, since certain circulation patterns, including a slow rotation, must be established to maintain a clear separation between a freshwater lens and the ambient seawater for some time.

 1  2  3  4  5  6  7  8  9 10 11 12

We will now have a closer look at the region marked by the white square.

 1  2  3  4  5  6  7  8  9 10 11 12

This is the main region of interest at full resolution (pixel size = 12.5 m × 12.5 m, area size = 6000 m × 6000 m). The front is a pronounced bright line. Unfortunately, SAR images at full resolution are always noisy. This so-called speckle noise results from coherent processing of the radar signal, which is required to obtain the high spatial resolution of SAR images. One can reduce the noise by a factor of N1/2 by averaging over N pixels (assuming that the intensity values represent independent samples of the same expectation-value image intensity).

 1  2  3  4  5  6  7  8  9 10 11 12

 

In this area, the front forms an angle of about 36° with the vertical axis of the SAR image. We cut a transect across the front to determine the magnitude and width of the image intensity variations. Considering just one single line of pixels, we obtain a very noisy result.

 1  2  3  4  5  6  7  8  9 10 11 12

This animation shows how we obtain a better estimate of the mean intensity profile across the front by averaging over up to 80 transects (1000 m along the front). The image intensity at the front is about 2.2 times the mean image intensity. The width of the SAR signature is on the order of 100 to 150 m.

 1  2 3  4  5  6  7  8  9 10 11 12

In order to determine the current variations at the front, we will now try to reproduce this observed radar signature with a numerical SAR imaging model, which can simulate theoretical SAR signatures of a given current field. For example, we can try a simple convergent current pattern (modeled by a hyperbolic tangent function) as given by the blue line, and we obtain the simulated SAR signature given by the red line. The radar imaging model M4S simulates such radar signatures by computing the spatial variations of the surface wave spectrum, which result from hydrodynamic wave - current interaction, converting the resulting wave spectra into Doppler spectra of the backscattered radar signal, and generating expectation-value SAR signatures from the Doppler spectra according to SAR imaging theory.

 1  2  3  4  5  6  7  8  9 10 11 12

This animation shows results of test simulations with pure current shears in the vicinity of the front. Assuming a cyclonal rotation of the Rhine plume, we expect a (small) positive current in y direction on the left and a (small) negative current on the right of the front. This would result in a dark signature! However, even if we reverse the shear current pattern, the simulated SAR signatures are much smaller than the observed one, even for unrealistically strong current shears. We conclude that the current shear at the front cannot be the dominant source of the observed SAR signature.

 1  2  3  4  5  6  7  8  9 10 11 12

Better agreement between model results and observations is obtained with a pure convergence at the front, as shown in this animation, which illustrates the iterative current field optimization procedure. Hydrodynamic modulation of the surface waves by convergent currents appears to be the main source of the observed bright SAR signature of the Rhine plume front. According to M4S, the current component normal to the front changes by about 0.4 m/s within about 20 m. This would be the main information extracted from the given SAR image by using MARSAIS tools. However, quantitative results obtained in this context are somewhat questionable, since the model validation is very difficult.

 1  2  3  4  5  6  7  8  9 10 11 12

Finally, we can have another look at the effect of shear currents in combination with the best-fit solution for the convergent cross-front current component. The dependence of the simulated SAR signature on the current shear is still quite weak. Accordingly, we will usually not be able to determine the strength of shear currents at oceanic fronts from their SAR signatures.

For more information see our other webpages, particularly the one on the modeling of SAR signatures and the one on MARSAIS, as well as the corresponding articles listed in the publication list of Roland Romeiser.