Internal Wave Dark-Band Signatures in ALOS-PALSAR Imagery Revealed by the Standard Deviation of the Co-Polarized Phase Difference

: Analysis of synthetic aperture radar (SAR) images in L-band of short-period internal waves (IWs), and classification of their radar signatures is presented by means of a polarimetric data set from ALOS-PALSAR mission. We choose the polarimetric feature named standard deviation of the co-polarized (std) phase difference (CPD) to identify fundamental differences in SAR signatures of internal waves, and divided them into three different classes, according to their backscattered modulation depths and morphology as well as the std CPD, namely: double-signed, single-negative, and single-positive signatures, for IW normalized image transects that display, respectively, signatures in the form of bright/dark, dark, and bright bands that correspond to positive/negative, negative, or positive variations of radar backscatter. These radar power types of signatures have a counterpart in the std CPD normalized transects, and in this paper we discuss those correlations and decorrelations. We focus in the single-negative type of signature, that is dark bands on gray background, and show that the std CPD is greatly enhanced over the troughs and rear slopes of those IWs. It is suggested that such behavior is consistent with the presence of surface slicks owing to enhanced surfactant concentration. Furthermore, those single-negative SAR signatures appear at locations where and when biological productivity is enhanced. It is found that the modulation depths associated to the std CPD is higher than the one associated to the HH-polarized radar backscatter for single-negative signatures propagating in the range direction, while the reverse occurs for the other types of signatures.


Introduction
Internal waves (IWs) are widespread, distributed around the world's oceans as the result of the variety of ways they can be generated (including from nonlinear disintegration of internal tides, lee waves, resonance, tidal beams, buoyant near-surface plumes, etc; see e.g., [1]). The diversity of the generation mechanisms, environmental conditions, and locations produce waves with wavelengths ranging from few tens of meters to kilometers; however, all reasonably described by some solitary wave theory. The IWs propagate in the ocean along the pycnocline, i.e., a water column portion characterized by a sharp change in density resulting, typically, from temperature and/or salinity vertical gradients. They connect large-scale tides to smaller scale turbulence and retain their form for a long time period after their generation [1].
The IWs are important for heat exchange in the vertical, as well as nutrients vertical exchange and other properties between shelf and the open ocean, playing a key role in biological primary radar modulation depths as well as its form (bright and dark bands that we denominate double-sign signatures in this article), do capture most of the IW signatures that have been observed in L-band SARs. However, there are many instances when the L-band signatures of IWs appear as dark bands on a gray background, as in the case of Figure 1, the usually companion bright bands being either absent or strongly attenuated. This fact is consistent with film slick effects proposed in [6] for C-band SARs, but it is also applicable to L-band (see [6] and their Figure 13). Here we discuss in some detail dark-band signatures such as those in Figure 1, which we denominate single-negative signatures, using a method known as the standard deviation of the co-polarized phase difference [17,18]. In this framework, it is well known that the SAR polarimetry provides extra benefits for a range of marine applications [17,18] and in the literature several polarimetric features measuring the departure from the almost deterministic sea surface Bragg scattering have been proposed [18]. The polarimetric feature named standard deviation (std) of the co-polarized phase difference (CPD) is related to the scattering properties of the observed scene and it is considered a measure of the complex correlation between the co-polarized channels. The latter polarimetric feature has been recognized to be very effective for surface slick observation purposes discriminating an oil slick from a weak-damping surfactant film [19]. Furthermore, the standard deviation of the CPD has been shown almost independent of the estimation window size and unbiased when low correlation of co-polarized channels applies [19,20]. It is important to highlight that the std CPD behavior is affected by incidence angle, noise equivalent sigma zero (NESZ), and meteo-marine conditions as discussed by [21][22][23].
In this paper, further insights on the effects of the surface films on the IW patterns in SAR imagery are analyzed by means of the std CPD for the first time in literature. The theoretical background is presented in Section 2. In Section 3, the study area is presented as well as the SAR data set and the ancillary information are described. The experimental results and discussion are presented in Section 4. Finally, the conclusion is shown in Section 5.

Theoretical Background
The complete Scattering matrix S is measured by a fully-polarimetric SAR. Taking into account the horizontal-vertical linearly polarized electric fields, S can be expressed as in [24] The scattering amplitude, i.e., each complex element of S, is expressed as where {sq} = {HV}. Upon the assumption of reciprocity (i.e., the two cross-polarized terms are equal, S HV = S V H [24]), we can write Equation (1) as S = e jφ VV S HH e jφ CPD S HV e jφ XPD S HV e jφ XPD S VV where and are, respectively, the CPD and cross-polarized phase difference, where and * stand for the phase part and complex conjugate operator, respectively. For most natural targets, the XPD probability density function (pdf) follows approximately a uniform distribution over [−π, π] [25], therefore, it contains less specific information related to the target. The analytical expression of the CPD pdf is given by [26,27] where and N is the number of looks, ρ is the correlation of the co-polarized channels, L (·) and Γ (·) are the first kind-Legendre and Gamma functions, respectively. Considering 0 < ρ < 1, the CPD pdf resembles a Gaussian bell with mean represented by µ φ CPD and standard deviation henceforth denominated by σ φ CPD . The increasing N or ρ implies that the CPD pdf becomes narrower. Thus, when ρ tends to 0, the CPD pdf is uniformly distributed between [−π, π]; otherwise, when ρ tends to 1, the CPD pdf tends to a Dirac delta function. The σ φ CPD is inversely proportional to ρ [19,26,27].
Considering a homogeneous sea surface where the Bragg scattering applies, the cross-polarized scattering amplitudes are almost negligible, a low depolarizing effect is expected as well as a high correlation between the co-polarized channels. Thus, a low σ φ CPD value is expected, i.e., σ φ CPD tends to 0 • . A surface slick damps the small scale ocean wave structures causing a departure from the conventional sea surface Bragg scattering. In this scenario higher depolarization effects occur as well as a lower correlation between the co-polarized channels are expected. Thus, higher σ φ CPD values apply [19,28]. Please note that according to a different theory [23] the behavior of a wide range of polarimetric features (including the σ φ CPD ) over both the slicked and unperturbed sea surface can be adequately explained accounting only to Bragg scattering and the noise effect, i.e., the "non-Bragg" theory is not needed to be evoked.
As reported in literature [21,23], the decreasing signal to noise ratio, SNR, (both signal-to-additive and -multiplicative) plays an important role in broadening the p (φ CPD ), i.e., the σ φ CPD increases with decreasing SNR. However, it was pointed out by [21,22] that the σ φ CPD is less influenced by environmental parameters such as wind (under low to moderate) and sea state conditions.

Study Area and Data Set
The study areas are located in the Andaman Sea of the Indian Ocean (offshore the coast of Myanmar and Malaysia) and the Galapagos Marine Reserve (GMR) of the eastern Pacific Ocean. In the Andaman Sea, IW signatures have been observed on both SAR and optical imagery. The region shows complex solitons as well as soliton interactions as a result of widely IW potential sources that are operating simultaneously in close proximity (e.g., internal tides, sills between islands and underwater volcanic seamounts) [29,30]. According to [31], the IWs interact strongly with surface waves, resulting in surface rips (short, choppy, and breaking waves). The phytoplankton biomass shows a variable spatial distribution, where the higher concentration is found near the shore and river mouths (for example, Malacca Strait and Myanmar estuaries), as well as a seasonal pattern, i.e., the highest concentration occurs during the winter (November-February) and the lowest during the summer (March-May) [32,33].
The GMR of the Pacific Ocean is known for being a complex transition zone that is under the influence of tropical, subtropical and upwelled water, where slicks, riplles and IWs are widely observed [34,35]. Based on satellite imagery, a complex distribution of IWs both between the Galapagos Islands and away from the archipelago is reported in [29]. The IWs are pointed in [35] as an important source of upwelling in the archipelago. In [34], upwelling nutrient-rich water zones are correlated to productive habitats, with higher phytoplancton biomass concentration than the surrounding areas. According to these authors, from 2002 to 2007, a seasonal pattern in the phytoplankton biomass peak is observed in the region, where, in south of 0.75 o N, a peak is observed around September and another peak in March. Furthermore, the phytoplankton biomass is also affected by the El Niño-Southern Oscillation events, i.e., the GMR experiences a reduction in the phytoplankton biomass during the El Niño events (reduction in topographic upwelling) and an increase during cooling events (sea surface temperature anomaly at 0.8 • C).
The SAR data set is composed by 9 Single Look Complex (SLC) L-band ALOS-PALSAR (APL) polarimetric scenes acquired between April 2007 and March 2011 (i.e., 2 images acquired in Andaman Sea and the others in GMR), with incidence angles ranging from about 22 • to 25 • . All images were acquired in ascending passing. The APL system was launched in January 2006 by the Japan Aerospace and Exploration Agency (JAXA) and operated at center frequency of 1.27 GHz with near range NESZ of −29 dB. The APL was operational until May 2011. The APL data scenes are provided by the Alaska Satellite Facility Distributed Active Archive Center (ASF DAAC) supported by NASA (https://search.asf.alaska.edu). Figure 2 shows 5 APL excerpt images where IW signatures are present. Furthermore, general information about the APL system and an overview of the SAR data set are presented in Table 1.
The wind information is acquired from the NCEP-DOE AMIP-II Reanalysis (R-2) data, i.e., National Centers for Environmental Prediction-Department of Energy Atmospheric Model Intercomparison Project (https://psl.noaa.gov/data/gridded/data.ncep.reanalysis2.html). The NCEP-DOE AMIP-II Reanalysis (R-2) is an updated 6-hourly global analysis series from 1979 to present and uses an improvement forecast model and data assimilation system. The 10 m height wind reanalysis u and v components are provided by NOAA Earth System Research Laboratory's (ESRL) Physical Science Division (PSD) with spatial resolution of 2.00 • latitude × 1.75 • longitude [36]. The SAR images were acquired under moderate wind conditions (i.e., wind speed ranging from about 4.5 to 8.5 m/s). The wind information is organized on Table 1.   Table 1.
The chlorophyll-a (CHL-a) concentration is acquired from the NASA Goddard Space Flight Center, Ocean Ecology Laboratory, Ocean Biology Processing Group (https://oceandata.sci.gsfc.nasa.gov/ MODIS-Aqua/Mapped/Monthly/4km/chlor_a/). The CHL-a concentration is calculated using an empirical relationship derived from in situ measurements of CHL-a and remote sensing reflectances (Rrs) in the blue-to-green region of the visible spectrum from the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument aboard the AQUA satellite. We used the MODIS Aqua Level 3 monthly products (binned, with 4 km in spatial resolution) [37]. The CHL-a concentration in the study areas associated to the SAR images are within the range from 0.11 to 2.12 mg/m 3 .

Experimental Results and Discussion
The APL scenes are radiometrically calibrated and the φ CPD is calculated according to Equation (4). The σ φCPD is estimated using a 5 × 5 sliding window. The full-polarized intensity channels are pre-processed to reduce speckle using a 5 × 5 average moving window.
For each scene of our data set, one IW profile is selected. However, since the IW patterns of two particular scenes differ significantly in type within a single scene in different locations, two profiles are selected for those scenes (i.e., Image ID 5 and Image ID 9), in total 11 IW signature profiles.
In order to compare the IW profiles and to study the IWs signatures quantitatively, the IW contrast, δ, is defined as follows where δ is calculated for each selected IW profile, considering both HH-and VV-polarized intensity and σ φCPD (respectively, referred in the text as δ σ 0

VV
, and δ σ φCPD ). Thus, l assumes the HH-and VV-intensity and σ φCPD information within the IW profile, respectively. l 0 assumes the mean value associated to the background, which was calculated selecting a region of interest (ROI) representing the unperturbed sea surface for each SAR image composing the data set. Inside each ROI, we randomly selected 1000 pixels for calculating the mean HH-and VV-intensity and σ φCPD , respectively.
In addition, we propose a simple figure of merit defined as where δ max and δ min are the maximum and minimum δ values in the IW profile. This feature is proposed in order to measure and compare the IW modulation associated to the HH-and VV-intensity and σ φCPD , respectively, ∆ σ 0 HH , ∆ σ 0

VV
, and ∆ σ φCPD . The flowchart of the methodological approach is shown in Figure 3. As represented in the flowchart, we conduced three different experiments, i.e., thread 1, 2, and 3 in the flowchart. In [6], the authors have identified three different kinds of IW contrast profiles, i.e., double and single-negative and -positive signatures. A scheme representing the three IW signature types is shown in Figure 4. Thus, since we identify the same type of signatures in our data set, the first experiment (i.e., thread 1) consists of splitting the IW contrast profiles into those three different classes and discuss the results following this rationality. The second experiment (i.e., thread 2) consists of calculating the SNR for each IW profile of our data set. In order to accomplish the analysis, for each IW profile, we split the SNR values into two different classes, i.e., IW crests and troughs. The third experiment (i.e., thread 3) comprises the analysis of chlorophyll-a concentration maps associated to each SAR image of our data set.

IW Signature Types Analysis
In this subsection, we analyze the results associated to the IW contrast profiles, which were classified into three different IW signature types, i.e., double sign, single-negative and positive. Please, see the thread 1 in the flowchart (Figure 3). The calculated ∆ values for each IW contrast profiles composing our data set is shown in Table 2. Please, note that the ∆ σ 0 HH is higher than the ∆ σ 0 VV considering all profiles except the profile ID 3, which have ∆ σ 0 . Considering the tilted Bragg scattering (two-scale) model, higher IW HH-intensity modulations are expected [8,38]. Consequently, taking into account this rationality, the follow results are discussed considering only the HH-intensity channels. Furthermore, considering each IW profile, the correlations between δ σ 0 HH and δ σ φCPD are calculated and the results are also presented in Table 2.   Table 2. General view of the IW signature patterns and associated parameters. The double sign and single-negative and -positive signatures are presented, respectively as + \ −, −, and +. profiles and their respective δ σ φCPD are shown in Figure 5. Further details are listed in Table 2. The δ σ 0 HH profiles have positive signatures more than 33% stronger than the negative ones. Please note that the δ σ φCPD peaks are associated to the δ σ 0 HH pits. As explained by the hydrodynamic theory related to IWs [3], the front and rear slope of the solitary wave are associated, respectively, with increasing and decreasing sea surface roughness (i.e., respectively, bright and dark band in SAR imagery). Considering the IW rear slope, the decreasing in the signal backscattered of the sea surface causes a consequent decreasing in the SNR values. Thus, since, as pointed by [21][22][23]39], the σ φCPD over the sea surface tends to increase with decreasing SNR, the IW front and rear slopes are associated, respectively, with lower and higher δ σ φCPD values (i.e., dark and bright bands in the δ σ 0 •

single-negative IW signatures
Considering the single-negative signatures, 4 HH-polarized contrast profiles of our data set correspond to this class (or kind) of signature. The δ σ 0 HH profiles and their corresponding δ σ φCPD are depicted in Figure 6. The δ σ 0 HH are mostly negative, as expected; while, the δ σ φCPD are mostly positive. The ∆ σ 0 HH varies between 1.02 (Image ID 5) and 1.47 (Image ID 8). Please, see Table 2. The ∆ σ φCPD varies between 0.73 and 2.59. Note that, the ∆ σ φCPD are more than 22% higher than the corresponding ∆ σ 0 HH considering the range-propagating IW signatures, (i.e., Image ID 3, 8 and 9). The possible explanation for the higher ∆ σ φCPD associated to those single-negative IW signatures is the coupled effect of the hydrodynamic modulation (as discussed previously for the IW double sign signatures) and the role of surfactant films that enhance the co-polarized channel decorrelation, increasing the σ φCPD values associated to the IW rear slope signatures. As demonstrated in [6], the IW negative contrast in radar power is stronger as surfactant film concentration (or film elasticity) increases. In another way, Image ID 5 has ∆ σ φCPD 30% lower than its corresponding ∆ σ 0 HH . It is important to point out that the IW signature in Image ID 5 is propagating in azimuth direction. As discussed by [6], the azimuth-propagating IW signatures are many times dominated by single-negative signatures. We recall that we use the same model in [6] to explain the class of signatures reported in this paper. In the particular case of azimuthally propagating IWs, positive contrast variations are absent in the model, and the backscatter contrast is expected to be negative, in agreement with the observation in Figure 6b (top panel). This is a consequence of assuming Bragg scattering and hence the IW hydrodynamic modulation does not change significantly range propagating Bragg waves for azimuth propagating IWs. Thus, the lower values of ∆ σ φCPD in this case can be explained by the fact that the increasing decorrelation between the co-polarized channels associated to the IW rear slopes are due only to the surface film modulation effect.
Please note that considering all single-negative IW signatures, the maximum and minimum δ σ φCPD values are associated to the IW troughs and crests, respectively, and the correlation between the δ σ 0 HH and δ σ φCPD values are always higher than 0.81, for single-negative IW contrast.
• single-positive IW signatures We identify 2 HH-polarized contrast profiles that follow the single-positive IW signature. The δ σ 0 HH profiles and their corresponding δ σ φCPD are depicted in Figure 7. As expected, the δ σ 0 HH values are mostly positive; while, considering image ID 2, the δ σ φCPD signatures have a double sign pattern and, for image ID 5, the δ σ φCPD values are mostly negative. The ∆ σ 0 HH are 0.80 (Image ID 2) and 0.47 (Image ID 5). Please, see Table 2. The ∆ σ φCPD for Image ID 2 and Image ID 5 are, respectively, 0.61 and 0.20, i.e., more than 23% lower than the ∆ σ 0 HH . According to [6], when the wind speed is very low (<2 m/s), the IW can be imaged as bright bands only since the expected dark bands are lost in the dark SAR image background. Another complementary theory presented in [40] associates the IW single-positive signatures to either or both the following mechanisms: (1) generation of bound centimeter-decimeter waves with Bragg wavelengths; and (2) wave breaking. The latter mechanisms are the responsible by a positive contrast even on the IW rear slope due to the indirect contribution of meter-scale waves to the backscattered signal [41]. We note in passing that some advanced radar imaging models use a composite surface expansion, which account for long-wave-short-wave interaction terms resulting in upwind-downwind differences of the backscattered signal, and hence could explain single-positive signatures [42]. Thus, no decreasing in the sea surface backscattering related to the IW rear slope is expected as well as no increasing decorrelation between the co-polarized channels (consequently, no decreasing in the SNR values are expected), resulting in no clear modulation of the σ φCPD . Please note that, since Image ID 2 and Image ID 5 are not acquired under low wind condition (see Table 1), the second theory (i.e., the one accounting for wave breaking) is more reliable. Considering both images there is no clear correlation between δ σ 0 HH and δ σ φCPD . This is confirmed by the lowest correlation between the δ σ 0 HH and δ σ φCPD (lower than 0.20).
In order to better understand the behavior of the std CPD associated to each type of IW contrast profiles, normalized histograms were computed and the results are shown in Figure 8. It is important to highlight that, for computing the histogram related to the IW single-negative sign signature, the azimuth-propagating IW image was not considered (because of its anomalous behavior as pointed out previously). It can be noted that the empirical statistical distribution of std CPD is different for each kind of IW contrast profile. For the double sign signatures (please, see Figure 8a), more than 50% of the std CPD values are concentrated between the range of 20 • and 35 • , being the mean and std values, respectively, 30 • and 10 • . Considering the single-negative sign signatures (see Figure 8b), the histogram is shifted to higher std CPD values, i.e., more than 50% of the std CPD values are concentrated within the range of 22 • to 52 • . Furthermore, the histogram becomes wider, with mean and std values of, respectively, 47 • and 20 • . The single-positive sign signatures (see Figure 8c) are associated to the most left-shifted histogram (i.e., shifted to lower std CPD values, with more than 50% of the std CPD values being within the range of 22 • to 26 • ) as well as it has the narrowest histogram (mean and std values of, respectively, 24 • and 4 • ). Since higher std CPD values are expected over the sea surface under the influence of slicks [19][20][21]43], the results related to the histograms analysis confirm the previous outcomes in this paper, principally, regarding the possible influence of surfactant films on the IW single-negative sign signatures.

Signal to Noise Ratio Analysis
In this subsection, we present the analysis of the SNR for each IW contrast profile of our data set. Please see the thread 2 in the flowchart presented in Figure 3. It was discussed in previous studies the manifestations of film slick formation related to IWs as well as the impact of those films on the IW signatures in SAR imagery [6,9]. In order to interpret the SAR imagery of IWs and to understand the role of film slicks in their signatures, the SNR must be accounted for. Hence, the SNR is calculated using the nominal NESZ associated to the APL sensor at off nadir angle 21.5 • [44]. For each IW profile, we split the σ 0 values into two different classes, i.e., crests and troughs. Figure 9 shows the SNR associated to the IW crests and troughs for each profile. The results are presented in the graphical format known as box plot, where the median SNR is represented by the horizontal line inside the box, while the first (Q 1 ) and third quartiles (Q 3 ) are, respectively, the lower and upper lines of the box. The minimum and maximum values excluding outliers are, respectively, the lower and higher short horizontal bars. The minimum and maximum values are calculated, respectively, using Q 1 − 1.5IQR and Q 3 + 1.5IQR, where IQR is the interquartile range (i.e., IQR = Q 3 − Q 1 ). Since we have identified three different types of HH-polarized contrast profiles, δ σ 0 HH , in our data set, i.e., double (Figure 9a) and single-negative and -positive signatures (respectively, Figure 9b,c), the results are organized according to the type of IW signature. Please note that the SNR associated to the VV channel is also shown, as well as the dashed red line representing the 3 dB threshold (below this threshold we consider the measurements as very low signal samples).
As expected, the median HH-channel SNR associated to the IW troughs are lower than the crests, consistently for all profiles. Considering the IW double sign signatures (Figure 9a), the SNR values are well above the 3 dB threshold for both IW crests and troughs. The SNR values associated to the IW single-negative signature profiles (Figure 9b) are well above the threshold for IW crests as well as for the IW troughs associated to Image 5 and 9; while, the SNR values are near the threshold for Image ID 8 and fall below the threshold for Image ID 3. The SNR values associated to single-positive signatures (Figure 9c) are well above the threshold considering both crests and troughs.
The highest separation in SNR statistics between IW crests and troughs are found for single-negative sign signatures; while, for single-positive signatures, the SNRs are very near for both crests and troughs. Note that, as expected, higher SNR values separability between the IW crest and troughs classes are associated to higher δ σ φCPD modulation and higher correlation between the δ σ 0 HH and δ σ φCPD values. The SNR variability (i.e., the SNR standard deviation) associated to the IW troughs is higher than to the crests considering all profiles. The difference between the SNR variability associated to the IW crests and troughs is more evident for the single-negative signatures (Figure 9b), excluding Image ID 5, where IWs travel in azimuth direction. The difference is less evident for the single-positive signatures. As discussed previously maybe the large SNR variability associated to the single-negative signatures is mainly due to the damping properties of the surfactant layer that varies in space along the IW trains, and consequently along IW propagation direction (i.e., the surfactant concentration is different for the different solitons in the same IW train [9]) causing the large SNR variability related to the single-negative IW trough's signatures, see Figure 9. This variability associated to the different damping properties of the surfactant layer is also discussed by [45]. The same analysis applies for VV channel. Please note that, as pointed out in other studies [23,45], the VV channel is less noise corrupted considering all profiles.

Biomass Production Validation
We present in this section the analysis of the chlorophyll-a maps associated to each SAR image of our data set as presented in the thread 3 in the flowchart, Figure 3. The biomass production by marine organisms in the oceanic environment, principally phytoplankton, is pointed out by marine biologists as the largest source of natural surfactant films at the sea surface [12,46]. Several studies discuss the capability of phytoplankton for producing surface-active materials as metabolic by-products [12,13]. The in situ production of organic substances by plankton and bacteria are the most important source of surface-active material in the marine environment and only about 10% of surfactants own their existence to anthropogenic input [47]. Laboratorial experiments conducted by [12] showed that surfactants are produced in culture media of species of marine phytoplankton depending on the age of culture and species. These authors conducted field experiments showing a positive correlation between surfactant activity in seawater column and phytoplankton production at several northern Adriatic Sea stations. Furthermore, space borne radar scatterometers have shown a strong reduction of the Normalized Radar Cross Section (NRCS) (6 to 20 dB) in areas where ocean color data show a high chlorophyll-a concentration [48]. In [49] it has been suggested that, natural slick coverage can be used as a proxy for ocean productivity. Here, we consider chlorophyll-a concentration as a proxy to the enhancement of surfactant film concentration (in more technical language, enhanced surface film pressure). Figure 10 shows the CHL-a concentration associated to each SAR image of our data set. The results are presented in a graphical format known as box plot. The outliers are represented by the red plus signs. The outliers are considered when the elements are >1.5IQR below the first quartile or above the third quartile. The sections with gray background are the ones identified as single-negative signatures of IWs propagating in range direction. The yellow background sections mean that the chlorophyll-a images are not useful owing to cloud coverage. Considering the double sign and single-positive IW signature types, the maximum CHL-a concentration excluding outliers is associated to Image ID 1 (0.87 mg/m 3 ). It is important to highlight that the CHL-a product related to Image ID 1 is very cloud-covered, with contaminated pixels being about 50%. Therefore the value should be considered with caution. In fact we discard this measurement. For all the other images following these signature types, excluding outliers, the maximum CHL-a concentration is always equal or lower than 0.57 mg/m 3 . The values of CHL-a concentrations associated to the single-negative signatures with IWs traveling in the range direction (i.e., Image ID 3, 8, and 9, gray sections in Figure 10) have maximum CHL-a concentrations excluding outliers always equal or higher than 0.63 mg/m 3 (Image ID 8). Excluding outliers, the maximum CHL-a concentration is found for Image ID 3, i.e., 1.1 mg/m 3 . Furthermore, the statistical distribution of the CHL-a concentration associated to these images are more evidently right-skewed than the others. Image ID 5, where the IW is traveling in the azimuth direction, is associated to the lower values of CHL-a concentration for the single-negative signatures. The latter image has maximum CHL-a concentration of 0.41 mg/m 3 . However, it is important to point out the CHL-a product for this image is jeopardized due to intense cloud-coverage, since the contaminated pixels are about 60%.
We note that, in general, the single-negative IW signatures (gray sections in Figure 10) present higher chlorophyll-a values, which also have a much larger number of outliers. We attribute the outlier values to possible phytoplankton blooms, which are translated into higher chlorophyll-a values. Note that the images are Level 3 monthly, indicating, in our view, that the outliers are true values, but not bad data. Hence, it is suggested that the elevated values of chlorophyll-a concentration in Image IDs 3, 8, and 9 are an indication of enhanced concentration of surfactant films. We stress that, this is consistent with our hypothesis for explaining the single-negative "dark-banded" IW signatures. Figure 11 shows the frames of the SAR images associated to the IW single-negative sign signature overlaid onto the productive habitats of the GMR mapped by [34]. Please, note that the SAR scenes were acquired near the productive habitats. The IW profiles that have higher ∆ σ φCPD values (i.e., Image IDs 3, 8, and 9) are associated with SAR images acquired in March. To the south of 0.75 • N, a biomass peak is observed around March [34], as presented previously in Section 3. These SAR images were also acquired during a weak and strong La Niña event, respectively, 2008/2009 (associated to Image ID 3) and 2010/2011 (associated to Image IDs 8 and 9).

Conclusions
A study aimed at investigating the IW signatures in SAR imagery by means of the polarimetric feature named std CPD was performed for the first time in literature. Experiments on actual SAR measurements were accomplished considering 9 L-band ALOS-PALSAR polarimetric scenes. The IW patterns were classified according to 3 different types of IW contrast profiles, i.e., double sign and single-negative and -positive signatures. The main outcomes of this study are summarized as follows: • Considering the IW double sign signatures, the modulation associated to the std CPD is lower than the one associated to the HH-and VV-polarized intensity channels. The decreasing correlation between the co-polarized channels on the IW rear slope (higher std CPD values) is presumably due to the lower sea surface roughness (caused by hydrodynamic modulation) and, consequently, lower SNR.

•
Taking into account the IW single-negative sign signatures, the modulation associated to the std CPD is higher than the one associated to the HH-and VV-polarized intensity channels for the range-propagating IW signatures. Probably, the reason is the coupled effect of the hydrodynamic modulation and the surfactant films associated to the IW rear slopes and over the IW troughs, which admittedly decreases the correlation between the co-polarized channels (causing the raise of the std CPD values). In other way, for the azimuth-propagating IW signature, the std CPD is lower than the one associated to the HH-polarized intensity channel. The likely explanation is that the decreasing correlation between the co-polarized channels (and, consequently, the increase of the std CPD values) are mostly due to only one effect, i.e., the surfactant films present in the IW rear slope and over the IW trough.

•
For the IW single-positive sign signatures, the modulation associated to the std CPD is lower than the one associated to the HH-and VV-polarized intensity channels. Since no decreasing in the IW rear slope sea surface backscattering is expected, no clear modulation of the std CPD is observed.
In this study we apply the same classification scheme used in [6] for C-band SAR images to IW signatures observed with the L-band ALOS-PALSAR. In [6], it was demonstrated that IWs appear as dark bands on a gray background in the presence of surface films, which were classified as single-negative sign signatures. In this study too, we observe dark signatures on a gray background for L-Band, and use the same classification scheme. We note in passing that the single-negative sign signatures were found in highly productive hot-spots off the Galapagos Islands, i.e., in the island-wakes prone to upwelling and in certain seasons. Nevertheless, this possible connection between surfactants from biologically productive areas and IW signatures needs further investigations.
The outcomes described in this paper may support the design of algorithms aimed at IWs automatic extraction and classification based on co-polarized coherent SAR measurement and the design of polarimetric SAR architectures that aim at observing a wide range of ocean features (e.g., ripples, oceanic fronts, river plumes, and IWs). In future work we may include the organization of field campaigns coincident with satellite SAR acquisitions, aimed at sampling surfactant films on the sea surface associated to IW features, and the extension of such analysis to a larger data set. The latter should include SAR images acquired with higher incidence angles and under different wind regimes such as very low and/or very high wind speed, as well as the use of SAR scenes acquired by different polarimetric sensors working in different frequencies (e.g., X-and C-Bands).

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations
The following abbreviations are used in this manuscript: