Abstract

The main objective of this paper is to demonstrate the capability of compact polarized (CP) synthetic aperture radar (SAR) to retrieve ocean wave field parameters. Souyris’ and Nord’s algorithms are used to carry out the reconstruction of CP SAR pseudo quad-polarized data for the ocean surface under both the circular transmit linear receive (CTLR) and π/4 mode. The results show that, for the CP reconstruction, Nord’s algorithm has a better convergence ability than Souyris’. In addition, the investigation of the reconstruction accuracy shows that the CTLR mode is superior to the π/4 mode, in terms of ocean surface reconstruction. It is, therefore, concluded that the reconstructed parameters of CP CTLR mode data by Nord’s algorithm adapt to retrieve ocean wave information. The ocean wave slope spectrum and other main wave parameters are also calculated from reconstructed CP data and compared with measurements from in situ National Data Buoy Center (NDBC) matched-up buoys. Comparison of CP SAR-based wave field information with buoy outputs also shows good agreement in the case of dominate wave height, wave direction, and wave period, with biases of 0.36 m, 17.96°, and 0.88 s, respectively.

1. Introduction

Ocean waves are the most important ocean surface features in the fields of ocean physics and coastal engineering. Ocean surface parameters, such as wave height, wave direction, and wave length, are important for sea-ice monitoring, ship detection, marine pollution forecasting, and in other areas [1, 2]. The wave slope spectrum, which provides information about the 2-D directional distribution of the wave energy density, is important to ocean wave information retrieval. In recent years, remote sensing as an efficient tool for earth surface monitoring gradually plays a key role in ocean information retrieval. Of all remote sensing sensors, SAR systems are the active satellite sensors with a resolution high enough to provide information about the ocean wave spectrum and so these systems play a key role in the measurement of ocean surface waves [35]. In contrast to in situ buoys, SAR can operate in all weather conditions or in areas far from the shore and can measure ocean waves across the globe continuously [6]. In addition to improving the understanding of the properties of the ocean, SAR observations can be used to develop the forecast accuracy of the operational wave model. It is, therefore, important to investigate the capability of SAR to carry out ocean surface monitoring, especially in terms of ocean wave retrieval.

In recent decades, satellites such as ERS-1, ERS-2, Envisat, Radarsat-1/2, Sentinel-1 A/B, and Chinese GF-3 have carried on many SAR systems for use in ocean surface monitoring. Several wave retrieval algorithms have also been developed for ocean information retrieval. For example, K. Hasselmann and S. Hasselmann [7, 8] proposed the Max Planck Institute (MPI) algorithm based on the closed nonlinear transformation relation between the ocean spectrum and the SAR image spectrum; the output from the wave model (WAM) is input as a first-guess spectrum in this method. Mastenbroek and de Valk [9] developed a semiparametric retrieval algorithm (SPRA) based on Hasselmann’s nonlinear theory [7] by separating the wind-driven wave and swell—the wind-wave information is provided by synchronous scatterometer measurements and the swell information is derived from the rest of the SAR signal. Wang [10] proposed a semiempirical algorithm which has the objective to estimate the wave height from Envisat ASAR wave mode imagery without any prior knowledge. Shao et al. [11] develop a semiempirical function for significant wave height (Hs) and mean wave period (Tmw) retrieval from C-band VV polarization Sentinel-1 SAR. With the development of polarization technology, quad-polarized SAR has also been applied in ocean surface monitoring. Zhang et al. [12] present a synergistic method to retrieve both ocean surface wave and wind fields from spaceborne quad-polarization SAR, with consideration to the nonlinear mapping relationship between ocean wave spectra and SAR image spectra, in order to synergistically retrieve wind fields and wave directional spectra. Schuler and Lee [13, 14] proposed an algorithm for retrieving ocean wave spectra using quad-polarized SAR imagery and also investigated the results using L- and P-band SAR data. However, this method is only suitable for the L and P bands. By modifying the polarization orientation angle, He et al. [15, 16] developed a new retrieval algorithm so that it could be applied to C-band quad-polarized SAR data. He et al.’s method was also validated using Radarsat-2 quad-polarized SAR data by Zhang et al. [17, 18].

In addition, although the quad-polarized SAR can provide abundant polarization information to develop the process of ocean information retrieval, the limited swath width has also become the disadvantage of quad-polarized SAR for large-scale monitoring. For overcoming the limitation of CP SAR, a new type of SAR sensor called compact p(CP) SAR has been invented, which provides a new way to measure ocean waves and its potential for use in ocean surface monitoring has gradually become clear [19]. Compared with quad-polarized SAR, CP SAR operates at lower power but still has the capability to reconstruct quad-polarized data, thus overcoming the disadvantages of the narrow swath and high power of quad-polarized SAR and the insufficient polarized information provided by single- or dual-polarized SAR [20, 21]. CP SAR can, in fact, be interpreted as a compromise between dual-polarized and quad-polarized SAR so that swath width is increased and the operating power is reduced at the same time. To date, India’s RISAT-1, launched in April 2012, and the Japanese ALOS-2 SAR satellite, launched in May 2014, have carried CP SAR. In the future, the planned Radarsat-2 constellation of satellites will also support CP SAR. Although CP SAR is carried on existing SAR platforms, there is still no well-calibrated data that can be used in research, especially in the case of ocean wave measurements. It is, therefore, important to study the potential of CP SAR for the retrieval of ocean wave information. In a previous study, Nunziata et al. [22] investigated sea-surface scattering with a polarized model and exploited the hybrid polarized (HP) mode to interpret ocean surface polarization parameters. Yin et al. [23] proposed a new extended X-Bragg model and also used C-band quad-polarized SAR data acquired by both SIR-C/X-SAR and Radarsat-2 to emulate the HP and π/4 CP modes. Geldsetzer et al. [24] proposed that the dependence between incidence angle, wind speed, and wind direction of CP SAR parameters are presented for open water. Selected 20 CP parameters are related to C-band geophysical model function (CMOD) output. CP parameters are simulated from polarized Radarsat-2 data and emulate data available on the pending Radarsat Constellation Mission (RCM). Li et al. [25] present a new empirical reconstruction algorithm for the modified N based on a data set of more than 2000 Radarsat-2 quad-polarization images and collocated buoy observations. It is proven that cross-polarization facilitates wind speed retrieval and improves the accuracy of the results, especially with respect to high wind speeds. In summary, the capability of CP SAR to reconstruct pseudo quad-pol data and ocean wind retrieval has been investigated in detailed; however, its potential for ocean surface reconstruction has not been discussed yet, let alone wave information retrieval.

In this paper, two quad-polarized SAR images are used to reconstruct the CP pseudo quad-polarized data and demonstrate the potential of CP SAR for ocean wave information retrieval. The data set used in this paper is illustrated in Section 2. In Section 3, we briefly introduced the CP theory and the ocean wave retrieval method. In Section 4, we investigate two CP algorithms used for ocean surface reconstruction and evaluate the optimum CP mode for ocean surface reconstruction. The FP and CP ocean wave retrieval results are compared with NDBC buoy data. Conclusions are given in Section 5.

2. Data Sets

2.1. SAR Data

Two scenes of Radarsat-2 quad-polarized image are used in this study. Both of two image acquisition parameters are listed in Table 1, and the study areas are shown in Figures 1(a) and 1(b). HH-, HV-, and VV-polarized single look complex image (SLC) SAR images are also shown in Figures 1(a) and 1(b).

2.2. Buoy Data

In order to validate the CP SAR wave retrieval results, NDBC buoy data available from the National Oceanic and Atmospheric Administration (NOAA) were collected along with the SAR imagery. The main wave parameters were extracted from the original measurements and are listed in Table 2. Among all the wave measurement data, the dominant wave height (DWH), the dominant wave period (DWP), and dominant wave direction (DWD) were used for validating the retrieval results.

Because of limited collection of NDBC buoy, only one buoy was collocated with SAR imagery of scene 2. So the imagery of scene 1 was performed to CP reconstruction and ocean slope spectrum retrieval methods, and scene 2 was added to validate the method.

3. Theory and Method

3.1. CP Reconstruction

CP theory is a newly developed method that can produce pseudo quad-polarized data from dual-polarized SAR data. The polarized information can be reconstructed in a satisfied precision with the decrease of sensor operation power [26, 27]. The CP has three main modes: the CTLR, dual circular polarized (DCP), and π/4 modes.

3.1.1. CTLR Mode

In the CTLR mode, the microwaves are transmitted in left-handed or right-handed circular polarization and received in HH and VV polarization. The scattering vector is expressed as

The relevant Hermitian covariance matrix is where stands for matrix transpose operations, stands for matrix conjugate operations, 〈〉 stands for average operations, is the real part of complex number, is polarization scattering matrix and p, q are horizontal or vertical polarized channel, is the imaginary unit, and is the imaginary part of complex number. It is also noticeable that this matrix can be simplified based on the assumption that the reflection is symmetrical.

In order to reconstruct a pseudo quad-polarized covariance matrix for the CP CTLR mode, also needs to be iteratively solved. This covariance matrix can be expressed as

It is noticed that in (4) is element from CP 2 × 2 covariance matrix.

3.1.2. DCP Mode

In DCP mode, the microwaves are transmitted in right-handed circular polarization and received in left-handed or right-handed circular polarization. The scattering vector is expressed as

The relevant Hermitian covariance matrix is

It should be noted that the last term can be set to zero based on the assumption of reflectional symmetry and that the covariance matrix can be simplified to

In order to reconstruct a pseudo quad-polarized covariance matrix for the CP DCP mode, needs to be iteratively solved. This covariance matrix can be expressed as [28]

It is noticed that in (8) is element from CP 2 × 2 covariance matrix.

3.1.3. π/4 Mode

In π/4 mode, the microwaves are transmitted in linear polarization with a 45° inclination and received in HH and VV polarization. The scattering vector is expressed as

The relevant Hermitian covariance matrix is

It should be noted that the last term can be set to zero based on the assumption of reflectional symmetry and that the covariance matrix can be simplified to

In order to reconstruct the pseudo quad-polarized covariance matrix for the CP π/4 mode, needs to be iteratively solved. The pseudo quad-polarized covariance matrix for the π/4 mode can be expressed as [28]

It is noticed that in (12) is element from CP 2 × 2 covariance matrix.

It should also be noted that the scattering vectors and are related by

The polarization information for the CTLR and DCP modes will be similar because of the above linear relationship, which means that the CP reconstruction result for these two modes will be the same [29, 30].

3.2. Wave Spectrum Retrieval Algorithm

For the linearly polarized backscattering cross section, the quad-polarized covariance matrix is used for the ocean surface wave polarization signatures, as follows [31]:

It is noticed that in (14) is element from FP 3 × 3 covariance matrix. Because the maximum contribution to the wave slope is obtained using only linear polarization, here, we set the ellipticity to zero and so the linearly polarized backscattering cross section at any polarization orientation angle can be derived as [14, 15]: where ,, and are the linearly polarized backscattering cross sections for the horizontal, cross, and vertical polarizations, respectively. is the correlation between the horizontal and vertical polarizations, which can be directly derived from the covariance matrix.

In the linear modulation theory of ocean waves, the ocean surface height can be regarded as a superposition of sine waves and the radar backscattering cross section expressed as [16, 32, 33] where is the ocean surface height, is the wavenumber, is the wave space vector, is the wave angular frequency, andis the complex conjugate of the series. is the modulation transfer function for SAR, which includes the tilt modulation, hydrodynamic modulation, polarization angle modulation, range and azimuthal direction modulation, and velocity bunching modulation. In all modulation functions, only the tilt and polarization orientation angle modulation depend on the radar polarization [16]. Moreover, the wave slope induced by the polarization orientation angle modulation can be neglected for horizontal or vertical polarization [16]. By (18), the terms that are not sensitive to the polarization vanish and the equation can be simplified to [16] where and are the complex conjugates of series, respectively; and are given by [16] where is the radar incidence wavenumber and is the incidence angle. and are given by (16) in [16]. By consolidating the equations, the straightforward algebraic calculation can be summarized as [16] where and are the range and azimuthal wave slope, respectively. Parameters and are also given by (8) [16]. With an appropriate choice of , the range wave slope and azimuthal wave slope can be directly obtained. The wave slope spectrum and the other wave parameters can be estimated from wave slope.

4. Results and Discussion

4.1. Comparision of Reconstruction Algorithms

Besides the reflectional symmetry hypothesis, the relation between the copolarized and cross-polarized responses is the main condition limiting the solution of the CP reconstruction equations [21, 34].

Although both Souyris et al.’s [21] and Nord et al.’s [34] algorithms are used to reconstruct pseudo quad-polarized data, they differ in their definitions of . Souyris et al. use 4 as the value of , while Nord et al. calculate the value as

We, therefore, used the two reconstruction algorithms to calculate the slope rate of curve between and (Figure 2). It is important to note that in Figure 2(a) is equal to 4 while in Figure 2(b) is given by , the value of which was calculated from FP data. We found that the convergence of Figure 2(a) was much inferior to Figure 2(b); the correlation coefficients, , were equal to 0.53 and 0.93, respectively. To determine the optimum algorithm, we also calculated statistics for the nonconvergence of pixels during the reconstruction process. The results for this are shown in Table 3. The number of nonconvergent pixels found when using the two subimages of the ocean surface was much lower for Nord’s algorithm than for Souyris’ algorithm. Nord’s algorithm is, therefore, more suitable for reconstructing pseudo quad-polarized data of the ocean surface.

4.2. Evaluation of CP Reconstruction Accuracy

CP SAR has three main modes: CTLR mode, DCP mode, and π/4 mode. Each CP mode has different orientation and ellipticity angles, which describe the transmitted and received wave status. Because of the linear relationship between the DCP and CTLR modes, these two modes have similar polarization information; therefore, we discuss only the π/4 mode and the CTLR mode in this paper. Figure 3 shows scatter plots detailing how well the derived pseudo quad-pol results fit the original quad-pol ones. The figure shows that the values of and between FP and CP are concentrated along the one-to-one line, which means that for and , the reconstruction of the pseudo quad-polarized data was good. However, the reconstruction of shows a slight underestimation in the CTLR mode and an overestimation in the π/4 mode, which means that the reconstruction of was not as good as for and . The reason for the low reconstruction accuracy for may be related to the weak volume scattering at the ocean surface. It has been previously documented that the CP reconstruction algorithms are generally applied to fields and forests, where the total scattering contains a strong volume scattering component. However, the cross-polarization information is not directly used in our paper, Therefore, the inferior reconstruction performance of will not make any influence on retrieval results.

To examine the reconstruction performance in more detail, we also selected 7 types of polarization parameters, including copolarization backscattering coefficient,, cross-polarization backscattering coefficient, copolarization correlation coefficient, and coherence matrix eigenvalues,, and, for each parameter, calculating the reconstruction accuracy for the CTLR and π/4 modes. The reconstruction accuracy was used to describe the amount of deviation in the CP reconstruction. The reconstruction accuracies for , , , , , , and are listed in Table 4. We used the descriptions “high,” “medium,” and “low” to describe the deviations, where “high” was defined as less than 10%, “medium” corresponded to 10%−30%, and “low” was defined as larger than 30%. Table 4 shows that, for both the CTLR and π/4 modes, the reconstruction accuracy for the parameters , , , and was high, with the values for the CTLR mode being superior to those for the π/4 mode. We, therefore, decided to use the CRTL mode for the subsequent ocean surface parameter reconstruction and ocean wave slope spectrum retrieval.

4.3. Retrieval of Ocean Wave Slope Spectra

In this section, we describe how Nord’s algorithm was used to reconstruct the pseudo quad-pol data in the CTLR mode. The retrieved ocean wave slope spectra based on the reconstructed pseudo quad-polarized data are shown in Figure 4. Compared with the ocean wave slope spectrum derived from the FP data, the results for the CP CTLR mode obtained in this study are satisfied. It is also clear that the wave slope spectrum derived from FP presents more details compared to that derived from the CP data. However, the two wave slope spectra are similar on the whole.

In order to make a quantitative assessment of the wave retrieval results, the main wave parameters of dominant wave height, wave direction, and wave period were also extracted from the wave slope spectra and compared with the NDBC buoy data. The expressions of main wave parameters are where represents the root mean square of wave slope, represents the dominated wave length, and , respectively, represent the components of in azimuth and range direction, is the dominated wavenumber, is acceleration of gravity, and is the water depth.

Figure 5 shows four subimages’ position which is performed to retrieve result statistics, and the 180° ambiguity of wave slope spectrum has been removed according to the measured wind direction of buoy. The results shown in Table 5 show that the average errors in the dominant wave height, wave direction, and wave period derived from the FP data are 0.14, 17.12, and 0.61, while those derived from the CP data are 0.36, 17.96, and 0.88, respectively. It should be noted that, although the average errors in the dominant wave height, wave direction, and wave period derived from the FP data were much larger than for the CP, the retrieved CP SAR results are also satisfactory for ocean wave measurements. It is also true that the retrieval accuracy obtained using CP SAR in this study was limited by the reconstruction of the CP pseudo quad-pol data—this was especially due to the loss in precision for . In future studies, the retrieval accuracy obtained using CP SAR will need to be improved.

5. Conclusion

In this paper, we investigated the capacity and potential of CP for ocean wave measurement. A Radarsat-2 quad-polarized image was used to reconstruct CP pseudo quad-polarized data and validate the potential of this for ocean slope spectrum retrieval. The main conclusions of this study are as follows.

For ocean surface waves, Nord’s reconstruction algorithm is much more effective than Souyris’. Because of its convergence capacity and iteration time, we recommend adopting Nord’s algorithm for the reconstruction of CP pseudo quad-polarized data.

For both the CTLR and π/4 modes, the reconstruction accuracy for the parameters , , , and was high whereas for , , and it was low. In terms of the retrieval of ocean wave slope spectra using the CP mode, the CTLR mode is superior to the π/4 mode and, therefore, for the retrieval of ocean wave slope spectra, we recommend using the CTLR mode.

The potential of CP for ocean wave slope spectrum retrieval has been demonstrated in this paper. The values for the main wave parameters were also validated using in situ NDBC buoy data. A comparison of the CP SAR-based wave field information with the buoy outputs showed good agreement in terms of the dominant wave height, wave direction, and wave period, with biases of 0.36 m, 17.96°, and 0.88 s, respectively. Overall, the results for the retrieval of ocean wave slope spectra and ocean wave parameters using CP SAR were satisfactory. Moreover, considering only single SAR imagery is performed to wave retrieval algorithm study, the universality validation of ocean wave information retrieval from CP SAR will still need more results to support. In our future study, we will collect more quad-pol SAR to match up with buoy for plotting the wave parameter difference between FP and CP SAR precisely.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this article.

Acknowledgments

The authors gratefully acknowledge the financial support provided by the National Key Research and Development Program of China (2016YFB0502504, 2016YFB0502500) and National Science Foundation of China under Grant nos. 41301500, 41431174, 61471358, and 41401427. This article is partly sponsored by the scholarship funding from the Chinese Academy of Sciences. The authors would also like to acknowledge Professor Zhang from the First Institute of Oceanography, State Oceanic Administration, for supporting research on CP theory. The authors would also like to acknowledge Dr. Bian and Dr. Liu from Institute of Remote Sensing and Digital Earth for polishing the original manuscript.