# Measuring ocean waves from space: objectives and characteristics of the China-France Oceanography SAtellite (CFOSAT)

Danièle Hauser, Céline Tison, Jean-Michel Lefèvre, Juliette Lambin, Thierry Amiot, Lotfi Aouf, Fabrice Collard, Patrick Castillan

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Danièle Hauser, Céline Tison, Jean-Michel Lefèvre, Juliette Lambin, Thierry Amiot, et al.. Measuring ocean waves from space: objectives and characteristics of the China-France Oceanography SAtellite (CFOSAT). 29th International Conference on Ocean, Offshore and Arctic Engineering, OMAE 2010, Jun 2010, Shanghai, China. pp.85-90, 10.1115/OMAE2010-20184. hal-00517959

HAL Id: hal-00517959

https://hal.science/hal-00517959v1

Submitted on 16 Oct 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# MEASURING OCEAN WAVES FROM SPACE: OBJECTIVES AND CHARACTERISTICS OF THE CHINA-FRANCE OCEANOGRAPHY SATELLITE (CFOSAT)

## Hauser Danièle

UVSQ, CNRS, UPMC
LATMOS, 10-12 Avenue de l'Europe
78140- Vélizy, France

## Tison Céline

CNES, "Altimetry and Radar office"
18 avenue Edouard Belin
31 400 Toulouse, France

Lefèvre Jean-Michel
Météo France
Toulouse, France

Aouf Lotfi
Météo France
Toulouse, France

Lambin Juliette
CNES
Toulouse, France

Collard Fabrice
CLS
Brest, France

Thierry Amiot
CNES
Toulouse, France

Castillan Patrick
CNES
Toulouse, France

## ABSTRACT

The Chinese and French Space Agencies are jointly preparing a satellite mission devoted to the monitoring of the ocean surface and related science and applications. This is the so-called "China France Oceanography SATellite" (CFOSAT), to be launched around 2013. This mission will provide simultaneous and collocated observations of wind at the ocean surface and spectral properties of surface ocean waves using two scatterometers, both in Ku-Band: SWIM for measurements of directional wave spectra and SCAT for wind vector measurements. The SWIM instrument will use a real aperture observation technique so as to avoid limitations encountered with SAR systems. This paper describes the main objectives and characteristics of the mission with a focus on the SWIM instrument designed and developed under French responsibility to measure directional spectra of ocean waves.

## THE CFOSAT MISSION

Oceanography greatly benefits from remote sensing satellites for global monitoring and forecast of the sea state. The CFOSAT (China France Oceanography SATellite) mission, whose launch is planned for 2013, is designed to provide at the global scale, observations of wind at the ocean surface and spectral properties of surface ocean waves. It will serve both operational needs for the surface wind and wave forecast (marine meteorology and climatology), and research needs by improving our knowledge on waves hydrodynamics, on interactions between waves and the atmospheric or oceanic boundary layers, and on interactions between electromagnetic signals and the ocean surface.

The CFOSAT mission consists of a polar orbiting system at about 514 km altitude. The platform will embark two

payloads, both are radar at Ku-band (13.2 to 13.6 GHz) scanning around the vertical axis; one is SWIM 'Surface Waves Investigation and Monitoring)', which is a 6-beam radar at small incidence (0 to 10°), with the nadir pointing beam providing like standard altimeter missions wind speed and significant wave height and the 6 to 10° beams devoted to the measurement of directional spectra of ocean waves; it is a wave scatterometer; the other one is SCAT which is a rotating fan-beam wind-scatterometer covering the 20-60° incidence angles. This combination of small to medium incidence angles is appropriate to measure simultaneously the wind vector (from large incidence) and wave properties (from small incidence). The combination of the two instruments will help to separate the different effects governing the evolution of the wave field.

CFOSAT will provide at the global scale, observations of wind vector, directional energy spectra of ocean waves (with wavelength between 70 and 500 m), significant wave height, steepness of long waves, and parameters characterizing the statistics of wave slopes (for waves with wavelength larger than a few tens of centimetres). The mission is designed to get a global full coverage of these parameters. Near-real time transmission of the data to ground station centres, will be used to estimate geophysical parameters which will feed data assimilation into numerical models and forecast procedures.

## THE SWIM INSTRUMENT

SWIM is designed to measure spectral properties of ocean waves based on the Jackson et al.'s concept [1] and used since the 90's on airborne systems [2,3]. The principle is based on the analysis of backscatter modulation due to the tilt of long waves, which is maximum when the look direction is parallel to the propagation direction of the waves. Unlike observations from
SAR, the modulation is integrated over the azimuth dimension of the footprint, and the discrimination in azimuth is obtained through an azimuth scanning. First ideas of designing such a system for a satellite configuration were presented by [4], and first detailed specifications of a satellite mission was proposed by [5].

SWIM is a modified version of the design proposed in [5], now developed as part of the CFOSAT mission. SWIM will operate in Ku-band (13.575 GHz), pointing at incidence angles between 0 and 10° (multi-incidence beams) with a beam scanning and in azimuth (0-360°).

The 6 to 10° beam are designed specifically for measuring directional ocean wave spectra in the wavelength range [70, 500 m] at a scale of about 70 km x 70 km on a 180 km wide swath with a minimum detectable wavelength of about 70m and an angular resolution better than 15°.

The nadir looking beam is designed to measure wind speed and significant wave height with a concept and accuracy similar to standard altimeter missions (e.g. similar to JASON specifications).

All the beams (0 to 10°) will also be used to derive properties of the wave slope statistics from the profiles of radar cross-section.

A scheme of the geometry of observation is shown in Figure 1. For an orbit altitude of 500 km, the footprint will be about 18 km x 18 km. In the look direction, the signal is analyzed with a high range resolution. The footprint will sweep a pseudo-circle with a diameter ranging from 18 to 90 km for incidence angles ranging from 2 to 10°. The surface pattern described by the instrument for all the incidence angles is shown in Figure 2 for successive scans over 360° in azimuth, and for a satellite motion of 7 km/s. The number of elementary measurements (prior to averaging processes) is represented by the colour code. As the feed horns of the antenna will not be aligned on the antenna baseplate, switches of 90° in azimuth view angles occur between some beams, leading to this ground pattern. Nonetheless, at the scale of a resolution cell (70km along-track), all the orientations are covered.

The acquisition durations spent on each incidence angle are called cycles. The global incidence coverage [0°-10°] lasts a “macro cycle” of 218 ms.

The transmitted signal is a chirp, a linearly frequency modulated pulse, with a 320 MHz bandwidth and a 50 µs time length. After reception, the signal is compressed in range onboard, using numerical compression, allowing to reach a range resolution of about 47cm. Besides, each beam's PRF (Pulse Repetition Frequency) is selected in order to maximize the number of independent echoes, under some constraints such as the duty cycle of the high power amplification subsystem and the Tx/Rx chronogram. The pulses acquired during each cycle are averaged in time and in range in order to reduce speckle and thermal noise. For data downlink constraints, some averaging (in time and range) has to be performed on-board.

The algorithm for signal time-integration will compensate in real-time for migrations due to range variations of a surface

cell from one pulse to the other due to satellite advection and antenna rotation [6].

<figure><img src="image_3.png" /><figcaption>Figure 1: Geometry of observation of the SWIM instrument</figcaption></figure>

For the off-nadir beams (6 to 10°), the raw data transmitted to the ground stations will consist in the mean wave form computed in real-time with range migration compensation: signal intensity versus range with a range resolution of about 1.5 m integrated in time over 30 to 35 ms depending of the beam incidence. The geophysical products (directional wave spectra) will be computed at ground.

The validity of the principle to derive the spectra of ocean waves from this concept has been demonstrated several times using airborne systems developed at NASA (ROWS system, [2]) and in France (RESSAC radar, [3]). Simulation approaches have been used first to assess the feasibility of the concept [5] and are presently extended and used to perform detailed analysis of the system performance [7].

<figure><img src="image_4.png" /><figcaption>Figure 2: SWIM sampling at the surface.</figcaption></figure>

The retrieval of the directional spectrum of surface waves is based on the spectral analysis of the radar backscatter modulations in each look direction, providing a modulation spectrum versus wave number of the ocean waves. After taking
into account noise corrections (in particular speckle noise), and normalization factors, the modulation spectrum is related through a linear relationship to the slope spectrum of the waves, $k^2F(k, \phi))$ where $F(k, \phi)$ is the wave height spectrum, $k$ the wavenumber of the waves and $\phi$ the look direction. The azimuth scanning provides the complete information (directional spectrum), with however a 180° ambiguity in the propagation direction. Overlap of the scans of the 6, 8, 10° beams will be used to estimate directional wave spectra within surface boxes of about 70 km x 70 km.

The nadir pointing mode will be used to estimate significant wave height and wind speed as in conventional altimeter missions.

The 0 to 10° radar mean cross-sections will be used to estimate statistical properties of wave slopes (short scales) as inspired by [8,9]

## SIMULATION OF PERFORMANCE OF WAVE SPECTRA ESTIMATES

A simulation tool, called SimuSWIM, has been developed to derive the quality estimation of waves from the SWIM characteristics. SimuSWIM is an end-to-end simulator whose main steps are:

- simulation of the sea surface knowing analytical wave spectrum,

- computation of measured backscattered intensity,

- signal inversion to estimate wave spectrum,

- comparison of the estimation with the input spectrum.

Figure 3 proposes two examples of a theoretical spectrum used as input and the corresponding estimated one, taking into account the up-to-date SWIM definition. The results are accurate as they are very close to the theoretical spectrum used as input.

The modulation spectrum is a random variable. Thus, in [7], the expression of the theoretical standard deviation $\sigma$ has been derived.

$$ \sigma_{\hat{P}_m}(k) = \frac{1}{\sqrt{N_s L_k}} \sqrt{ \frac{1}{N_{imp}^2} \frac{P_{sp}^2(k)}{R^2(k)} + \frac{2}{N_{imp}} \frac{P_{sp}(k)P_m(k)}{R(k)} + \frac{1}{N_{imp}^2} \frac{P_b^2(k)}{R^2(k)} } \quad (1) $$

with $P_{sp}$ the speckle spectrum, $P_m$ the modulation spectrum, $R$ the impulse response, $P_b$ the thermal noise spectrum, $N_{imp}$ the number of impulses averaged in time, $L_{dis}$ the number of samples averaged in space, $N_s$ the number of spectrum averaged, $L_k$ the number of samples averaged in wavenumber.

This enables to define the $3\sigma$ margins around the mean value. The estimation (Fig. 3) is quite entirely within these margins. Nonetheless, the simulation itself induces some artifacts due to digitalization and under-sampling. This may explain why some values are outside the +/-3$\sigma$ margins.

Figure 4 illustrates the simulation on the 2D spectrum in a mixed sea case (non fully-developed wind sea spectrum and

swell) compared to the reference. As automatic comparisons are not easy, some analysis tools are under development to automatically partition the spectrum into main sea states and calculate the associated geophysical parameters such as propagation direction, dominant wavelength, energy, etc.

Until now, theoretical analyses and simulations have confirmed that SWIM design is compliant with scientific requirements.

These results are based on two criteria calculated from the simulated modulation spectrum, $\hat{P}_m$: i) a wavenumber-dependent criteria, $\varepsilon(k)$, and ii) an integrated criterion $\Delta E$:

$$ \begin{gathered} \varepsilon(k) = \frac{\bar{\sigma}_{\hat{P}_m}(k)}{P_m(k)} \\ \Delta E = \left| \frac{E}{E_{ref}} - 1 \right| \\ \text{where } E = \int_{k_{min}}^{k_{max}} \frac{\hat{P}_m(k, \phi)}{k} dk \end{gathered} \quad (2) $$

and $E_{ref}$ is computed with analytical shape of $P_m$. $\hat{P}_m$ is the estimated modulation spectrum. The computation of $\Delta E$ is made for a worthiest case, i.e. assuming that the estimation is equal to the reference spectrum $P_m$ plus the standard deviation:

$$ \hat{P}_m = \bar{\sigma}_{\hat{P}_m} + P_m \quad (3) $$

Figure 5 shows an illustration of these two criteria for the case of a wind sea modeled by a Pierson-Moskowitz function (and a wind speed of 13 m/s). Figure 5a shows that when integrated over the energy containing part of the spectrum (0.04 to 0.12 rad/m), the predicted error on the total energy is less than 15%. Figure 5b shows that for all wavenumbers within $[k_{min}, k_{max}]$, the predicted uncertainty is less than 10%. ($k_{min}, k_{max}$ are wave numbers corresponding to the half-power energy density with respect to the peak).

Table 2 presents a summary of performances obtained for a set of simulations run in different conditions of sea-state and for an optimal set of instrument and processing parameters (parameters affecting the radar signal, $\hat{P}_m$, and its variance). The most important parameters affecting these estimates are the averaging coefficients and the Signal to noise ratio. They are linked to the initial spatial resolution (i.e. the emitted signal bandwidth), the sensor noises and the resolution requirements. According to these simulations, the instrument and processing characteristics chosen for SWIM are appropriate to measure spectral properties of wave in a large variety of conditions.

## ONGOING STUDIES RELATED TO SWIM AND CFOSAT

Several French groups are presently working on the CFOSAT/SWIM preparation.
First assimilation studies have been performed at Meteo-France either using synthetic wave spectra or ENVISAT ASAR level 2 wave data [10]. The main conclusion is that the assimilation of directional information into wave prediction models induces a significant improvement in the estimate of wave parameters, such as the significant wave height and mean wave period. Figure 6 illustrates the impact of assimilation during the forecast period. It shows the index of assimilation, which characterizes the reduction of rms error on significant wave height gained through the assimilation, compared to Jason-1 observations. It shows that the impact of assimilation is larger (higher index of assimilation) and stays longer when directional observations are assimilated. More generally, the assimilation system corrects well the misfit of the wave model in case of severe storm generating high swell travelling through the ocean.

Developments of swell tracking algorithms have been performed and applied on SAR wave mode of ENVISAT [11]. They illustrate the interest of spectral wave observations to improve our knowledge of physical processes impacting swell propagation [12,13]. Such approaches will also be possible for SWIM data analysis.

In order to estimate directional spectra of ocean waves using the specific sampling geometry of SWIM (see Figure 2), work is under progress to prepare level 2 products based on a the multi-beam observations.

Simulations of cloud water and rain impacts on the signal are also under development to prepare rain flags.

## CONCLUSION

SWIM is a new concept for measuring spectral properties of ocean waves from satellite. It is presently under development (phase B realization). It shall provide real breakthroughs in the wave observation field by providing wave spectra under large range of sea-state conditions. SWIM and SCAT combined on the same platform (CFOSAT) will foster new studies on coupled wind/wave phenomena and processes at the air/sea interface affected by wind and/or winds. The mission will be beneficial to operational forecast centers for sea-state forecast.

## REFERENCES

1. Jackson, F. (1981). An analysis of short pulse and dual frequency radar techniques for measuring ocean wave spectra from satellites. *Radio Science*, 16(6):1385-1400.

2. Jackson, F. C., W.T.Walton, and C.Y. Peng, A comparison of in situ and airborne radar observations of ocean wave directionality, *Journal of Geophysical Research*, 90(C1), 1005-1018, 1985

3. Hauser, D., G. Caudal, G.J. Rijckenberg, D. Vidal-Madjar, G. Laurent, and P. Lancelin, RESSAC: A new airborne FM/CW radar ocean wave spectrometer, *IEEE Trans. Geosci. Remote Sensing* 30 (5), 981-995, 1992

4. Jackson F. C., T.W. Walton, and P.L. Baker, Aircraft and satellite measurement of ocean wave directional spectra using scanning-beam

microwave radars, *Journal of geophysical research*, Vol. 90, No. C1, Pages 987-1004, January 20,1985

5. Hauser D., E. Soussi, E.,Thouvenot, L. Rey: SWIMSAT: A real aperture radar to measure directional spectra of ocean waves from space, Main characteristics and performance simulation, *Jour. Atmos. and Oceanic Tech*, vol 18 No3, 421-437, 2001.

6. V. Enjolras V., L. Rey, L. Cros, S. Pouyez, T. Amiot, C. Tison, P. Castillan SWIM: a state-of-the art multi-incidence beams Ku-band waves scatterometer to go beyond current radar systems, *IGARSS'09*, July 2009.

7. C. Tison, T. Amiot, J. Bourbier, D. Hauser, V. Enjolras, L. Rey, P. Castillan, A spaceborne radar for directional wave spectrum estimation: first performance simulations, *IGARSS'09*, July 2009

8. Chapron, B., V. Kerbaol, D. Vandemark, and T. Elfouhaily (2000), Importance of peakedness in sea surface slope measurements and applications, *J. Geophys. Res.*, 105(C7), 17195-17202.

9. Hauser D., G. Caudal, S. Guimbard, et A. Mouche, A study of the slope probability density function of the ocean waves from radar observations, *J. Geophys. Res.*, 113, C02006, doi:10.1029/2007JC004264, 2008

10. Aouf, L., Lefèvre, J-M., Chapron, B. and Hauser, D. Recent improvements of the assimilation of upgraded ASAR Level 2 wave spectra in the wave model. *Proceedings of the SEASAR 08*, Frascati, January 2008.

11. Collard F., F. Ardhuin, B. Chapron (2009), Monitoring and analysis of ocean swell fields from space: New methods for routine observations, *J. Geophys. Res.*, 114, C07023, doi:10.1029/2008JC005215.

12. Collard, F., F. Ardhuin, and B. Chapron, Persistency of ocean swell fields observed from space," J. Geophys. Res., vol. submitted, 2008.

13. Ardhuin F., Chapron B, Collard F., Observation of swell dissipation across oceans, *Geophysical Research Letters*, March 2009; Volume 36 (L06607): P. 1-5

<table><caption>Table 1- Theoretical analysis of SWIM performance for different sea-state conditions.</caption><thead><tr><td>Sea states</td><td>Estimation error (ε) per k on the 3dB wavelength domain</td><td>Energy error (ΔE) on the 3dB wavelength domain</td></tr></thead><tbody><tr><td>Pierson-Moskowitz U=13 m.s-1</td><td>~10-15 %</td><td>~10%</td></tr><tr><td>Pierson-Moskowitz U=10 m.s-1</td><td>~15-20 % (20% for k&lt;0.06)</td><td>~15-20 %</td></tr><tr><td>Jonswap U=13 m.s-1</td><td>~5 %</td><td>~5 %</td></tr><tr><td>Jonswap U=10 m.s-1</td><td>~5 %</td><td>~5 %</td></tr><tr><td>Swell Hs=4m</td><td>~10 %</td><td>~10 %</td></tr><tr><td>Swell Hs=2m</td><td>~15 %</td><td>~15 %</td></tr><tr><td>Swell Hs=1m</td><td>~20 %</td><td>~20 %</td></tr></tbody></table>
<figure><img src="image_5.png" /><figcaption>Figure 3. Simulation of two sea states (wind sea (a) and swell (b)). The black curves are the input spectrum, the blue curves the estimated one taking into account all the properties of the instrument and of the processing chain, the red curves the speckle spectrum. The dotted lines are the +/-3σ margin around the theoretical value (σ begin the standard deviation).</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>Figure 4. Simulated 2D wave spectra for a Jonswap spectrum (U=10 m.s⁻¹fetch = 60 km) after inversion (left) compared to the analytical reference (right). The resulting peak wavelength and peak direction as estimated after partitioning are noted on the figure (left top). The 180° ambiguity in direction is not removed in the example.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>Figure 5. Quality criteria on the inverted modulation spectrum, in a wind-sea case (Pierson Moskowitz spectrum with a 13 m/s wind speed). (a) uncertainty on integrated energy (b) estimation error per wavenumber (the red curve stands for the sum all of sources of errors)</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>Figure 6: Index of assimilation impact versus time during the forecast period. Red line stands for the assimilation of altimeter Ra-2 wave height only, while green line indicates the assimilation of altimeter Ra-2 and ASAR wave spectra.</figcaption></figure>