A Sensitivity Analysis on the Spectral Signatures of Low-Backscattering Sea Areas in Sentinel-1 SAR Images

: Satellite synthetic aperture radar (SAR) is a unique tool to collect measurements over sea surface but the physical interpretation of such data is not always straightforward. Among the different sea targets of interest, low-backscattering areas are often associated to marine oil pollution even if several physical phenomena may also result in low-backscattering patches at sea. In this study, the effects of low-backscattering areas of anthropogenic and natural origin on the azimuth autocorrelation function (AACF) are analyzed using VV-polarized SAR measurements. Two objective metrics are introduced to quantify the deviation of the AACF evaluated over low-backscattering areas with reference to slick-free sea surface. Experiments, undertaken on six Sentinel-1 SAR scenes, collected in Interferometric Wide Swath VV+VH imaging mode over large low-backscattering areas of different origin under low-to-moderate wind conditions (speed ≤ 7 m/s), spanning a wide range of incidence angles (from about 30 ◦ up to 46 ◦ ), demonstrated that the AACF evaluated within low-backscattering sea areas remarkably deviates from the slick-free sea surface one and the largest deviation is observed over oil slicks.


Introduction
Continuous and effective monitoring of the oceans is of paramount importance to improving global marine awareness and the understanding of ocean dynamics, including man-made target surveillance, pollution monitoring, and the impact on climate change [1]. Satellite Earth observation represents a valuable tool that provides extensive data collection over the oceans. An important sensor for ocean observation is the synthetic aperture radar (SAR), an active, band-limited, and coherent microwave imaging sensor that provides day and night imagery in almost all-weather conditions [2]. The exploitation of SAR imagery for marine and maritime applications is now well-established [3][4][5][6] and has been further boosted when the European Space Agency (ESA) started providing Sentinel-1 (S1) SAR satellite measurements free of charge in 2014 [7][8][9][10].
The general physical modeling that rules sea surface scattering in SAR imagery accounts for both the sea surface roughness, i.e., the sea surface spectrum, and the sea permittivity. In the context of marine pollution, low-backscattering areas, whose normalized radar cross section (NRCS) is lower than the surrounding sea one, represent a broad class of targets of interest since they can be often associated to natural or anthropogenic oil slicks. Even though low-backscattering sea areas may have different origins, e.g., organic films, low-wind areas, etc., most analyses on low-backscattering areas have been focused on oil spill detection and classification [11][12][13][14][15][16][17], with important advancements that have been recently accomplished by means of proper physical processing of polarimetric SAR data [18][19][20][21].
Low-backscattering areas are due to the damping of the capillary and small gravity sea waves generated by the local wind and responsible for the measured NRCS, i.e., the centimetric Bragg resonant waves, which result in a reduced sea surface roughness. Accordingly, most of the incident energy is scattered in the specular direction rather than toward the SAR antenna. This is why low-backscattering sea areas appear as dark patches in graytones intensity SAR images [22,23]. However, note that the NRCS values depend on incident wavelength, polarization, and incidence angle [22]. In this paper, a new approach is proposed to exploit VV-polarized, i.e., vertical transmit/vertical receive, S1 SAR data, collected in interferometric wide (IW) swath mode, to characterize different kinds of large low-backscattering sea areas of a known origin. The approach, based on the autocorrelation function evaluated along the azimuth (AACF), relies on the spectral information inherently carried on by SAR measurements and on the SAR peculiar imaging mechanism along the azimuth [24]. The experimental analysis, performed under low-tomoderate wind regime, i.e., wind speed ≤ 7 m/s, and in a broad range of incidence angles (approximately 30 • -46 • ), considers slick-free sea surface as the reference scenario. The deviation of the AACF evaluated over low-backscattering areas of both an anthropogenic and natural origin is quantified using two objective metrics, namely the Euclidean distance and the percentage relative difference. Furthermore, the AACF analysis is compared to the conventional contrast, i.e., the slick-free sea surface to low-backscattering area NRCS ratio, and the effect of incidence angle is also discussed. The experimental results show a pronounced sensitivity of the AACF to low-backscattering sea areas. According to the metrics, the oil AACF differs from the slick-free sea surface AACF and from other natural low-backscattering area AACFs. In addition, the Euclidean distance, if compared to the NRCS contrast, is less dependent on the incidence angle.
The remainder of this paper is organized as follows: The theoretical background is presented in Section 2, while the dataset is described in Section 3. The experiments are presented in Section 4, where the results are discussed, while conclusions are drawn in Section 5.

Physical Rationale
The SAR estimates the sea reflectivity by means of two different scanning mechanisms. The range or across-track direction imaging is done at the speed of light, therefore insensitive to temporal changes of the sea, while the azimuth or along-track direction imaging is done at satellite velocity, therefore sensitive to temporal changes of the sea [24][25][26]. The microwave signal scattered by the sea, under low-to-moderate sea state conditions and in a broad range of incidence angles, can be well described by two-scale scattering models [27]. The class of two-scale scattering models assumes that the full-range sea surface roughness spectrum is artificially split into two parts: The larger-scale roughness, mainly associated to longer surface waves, and the smaller-scale roughness, mainly associated to shorter surface waves, which are responsible for the Kirchhoff and Bragg scattering, respectively [28]. The choice of the K lim , i.e., the wavenumber that splits the sea surface spectrum, has been addressed in several papers and is summarized in [28,29] and references therein.
These two contributions are not independent since, for instance, the non-linear interaction between the capillary (shorter) and gravity (longer) waves is the physical mechanism that underpins the energy transfer from the wind to the waves [30]. However, to achieve a good compromise between accuracy, practical implementation, and interpretation, these two contributions are added incoherently [27]: where σ 0 pq , σ 0 pq,0 , and σ 0 pq,1 are the total NRCS, the zeroth-order Kirchhoff scattering contribution, and the first-order tilted-Bragg scattering contribution, respectively [28]. In (1), the subscripts "q" and "p" stand for transmitted and received polarization, respectively. In the S1 IW SAR imaging mode, the incidence angle is such that the main term contributing to σ 0 pq is the first-order Bragg scattering [28]. The latter, which can be considered as an average of the untilted Bragg scattering over the larger-scale ripple by the long-waves structure, broadens the spectrum of the Bragg resonant waves [28]. Hence, according to the tilted-Bragg model, the microwave signal backscattered to the radar antenna depends directly on the small-scale ripple (through the small-scale sea surface roughness spectrum) and, indirectly, on the longer-wave part of the spectrum through the probability density function of the slopes [28,31,32].
The presence of a surfactant over the sea surface affects both the NRCS and SAR image spectrum [33]. In fact, it modifies the full-range sea surface spectrum [28], i.e., both the short-and the longer-wave part of the sea surface spectrum are affected by the surfactant. The visco-elastic properties of the surfactant have a direct impact on the small-scale part of the spectrum through a damping coefficient that, in the case of monomolecular surface films, can be modeled by the Marangoni damping coefficient [34]. In this case, as expected, the Marangoni damping mainly affects the small-scale part of the sea surface spectrum. The surfactant, reducing the sea surface roughness, also affects the energy transfer from the wind to the sea waves. The latter is typically modeled by the friction velocity, with this phenomenon that is well described by a reduced friction velocity over the surfactant. A reduced friction velocity has a direct implication on the long-wave part of the sea surface spectrum by modifying the peak wavenumber and the significant slope [28].
The above-described physical rationale does not explicitly take into account sea dynamic processes that are of paramount importance when dealing with SAR imaging of sea surface. The effect of the sea dynamics on SAR images is known as the velocity bunching phenomenon [35,36] that is related to the azimuth channel. This imaging is affected by the scene coherence time which, being shorter than the SAR integration time, makes the SAR imaging act as a non ideal filter along the azimuth direction [25,37].
These dynamics processes, which make SAR images appear blurred in the azimuth direction [36], are such that the "optimal" focusing depends on the SAR image patch [36] and it can narrow/broaden according to the sea state [24]. Although all this matter was first seen as a limit in SAR imaging of the oceans [24], it was later considered as a geophysical information to be potentially exploited [36]. This phenomenon is known as azimuth cut-off [38]. In filter theory words, one can say that the azimuth cut-off is a measure of the actual SAR azimuth spatial resolution [24]. In [35] the azimuth cut-off is explained as being due to two main contributions. The radial velocity to the radar of the single moving water particle generates a Doppler shift with respect to a stationary scene and the Doppler shifts of the elementary scatterers in the SAR resolution cell are not all identical, i.e., the single elementary Doppler shift is different from the mean Doppler shift of the SAR resolution cell, producing a "velocity spread" that smears the image of the resolution cell [24,35].
The physical model detailed in [35] shows that the quasi-linear approximation of the SAR image power spectrum P ql (K) is modeled as a perturbation of the linear SAR image power spectrum P l (K) modeled in accordance to the linear imaging theory. Hence, the quasi-linear SAR image power spectrum P ql (K) can be expressed in function of the linear SAR image power spectrum P l (K) as follows [33,35]: where the subscript x refers to the azimuth direction and ξ 2 is the total variance of the azimuthal displacements within the SAR integration time. Equation (2) means that the SAR image spectrum is not able to represent the waves whose wavelength is less than the so-called azimuth cut-off wavelength λ c , which can be interpreted as a measure of the low-pass filtering, in the wavenumber domain, witnessed by the exponential term. In (2), ξ and consequently λ c are related to the sea surface spectrum as follows [38]: where θ is the incidence angle, φ is the wind mean direction relative to the range axis [38], R/V is the ratio between the slant-range distance and the velocity of SAR platform, S(·) is the omni-directional sea spectrum and F(·) is the directional sea spectrum. This reasoning can be exploited over slick-free sea surface, by means of a tailored semi-empirical model and proper estimation procedure, to determine the wind speed [33,39].

Methodology
In this study, the large degree of heterogeneity resulting from the low-backscattering sea class prevents the use of a compact parameter for the spectral analysis as in the case of wind speed estimation, that relies on λ c . Hence, here the analysis is focused on the AACF that is obtained as the inverse Fourier transform of the power spectral density (PSD) following a guideline similar to the one proposed in [39] to estimate λ c . We expect a sensitivity of the AACF to the low-backscattering areas. In fact, the theoretical modeling here outlined shows that the surfactant affects the NRCS and the SAR image spectrum. The relationship between the slick-covered sea surface spectrum and the SAR image spectrum is generally complex and to some extent poorly known. Since the friction velocity reduction affects the dynamic process, it is expected to have an additional impact on the azimuthal SAR spectrum. Hence, we expect that the AACF evaluated over low-backscattering areas exhibits deviation from the reference slick-free sea surface one.
In detail, the methodology proposed to process S1 SAR imagery using the AACF can be summarized as follows, see also the flowchart of Figure 1: • The VV-polarized uncalibrated intensity image of the SAR scene is divided into nonoverlapped square boxes whose size is n = 128 pixels. This value is set according to the S1 pixel spacing, i.e., 10 m, in order to ensure both a satisfactory degree of homogeneity and a consistent number of samples for a reliable AACF estimation; • For each box, the 2-D PSD is evaluated as the square modulus of the Fourier transform of the VV-polarized uncalibrated intensity image; • The 1-D azimuth PSD, PSD x , is evaluated by averaging the PSD along the range direction; • The AACF is obtained by applying the inverse Fourier transform: where PSD = FFT(X), with X is the VV-polarized uncalibrated intensity SAR image; • A smoothing 7 × 1 median filter is applied to the modulus of the azimuth autocorrelation function (AACF) in order to remove the 0-lag contribution.

Dataset
The Copernicus Earth observation S1 mission consists of two polar-orbiting satellites equipped with a C-band (5.405 GHz) SAR. The latter supports dual-polarimetric imaging modes through a single switchable linear polarization transmission (horizontal-H, or vertical-V) while receiving coherently in a linear orthogonal polarization basis H-V. S1 operates in four different acquisition modes: Stripmap, extra-wide swath, wave mode, and IW mode. The IW mode is composed by three sub-swaths obtaing using the terrain observation with progressive scans SAR (TOPSAR) imaging technique. The TOPSAR technique, together with the electronic steering of the beam, result in a high-quality image characterized by no scalloping and homogeneity throughout the swath.
In this study, a SAR dataset consisting of 6 C-band S1 VV+VH ground range detected images collected in IW mode is considered that allows combining a large swath width (250 km) with a moderate geometric resolution (5 m by 20 m). The swath coverage is around 250 km, the pixel spacing is 10 m, and the incidence angle θ spans from about 30 • up to 46 • . Although both co-and cross-polarized channels are available, in this study only co-polarized S1 imagery is used since the cross-polarized channel is significantly affected by noise due to the smaller signal-to-noise ratio (SNR) and, in addition, it exhibits a low sensitivity to low-backscattering areas. All the images are characterized by low-to-moderate wind conditions, i.e., wind speed between 2 m/s and 7 m/s, therefore representing an optimal range for studying biogenic surfactants and oil slicks as suggested in [40]. The wind speeds are provided by space/time co-located ancillary European Centre for Mediumrange Weather Forecasting (ECMWF) information.
The S1 SAR images include low-backscattering areas of a known origin [41][42][43][44][45][46][47] related to both verified oil slicks and other natural phenomena, see Table 1. The corresponding VV-polarized NRCS images are shown in Figure 2, where a graytones dB scale is used.
The regions of interest (ROIs) considered for the experimental analysis are highlighted with dashed boxes that refer to slick-free sea surface (blue and red), oil slick (green), and look-alike (orange).
The image shown in Figure 2a is acquired on 10 August 2017, off the southern coast of Kuwait in the Persian Gulf, where a certified oil spill is present [41]. The oil spill is likely due to an accidental collision between the pipeline laying vessel "DLB 1600" and an old pipeline on the seafloor. A conservative estimation of the oil-covered area is 131 km 2 . Another low-backscattering area, due to very low wind conditions (<3 m/s), is also present. The incidence angle relevant to both low-backscattering areas is equal to 32 • .
The image shown in Figure 2b was acquired on 8 October 2018, over the northern part of the Tyrrhenian Sea between Corsica and Tuscany coasts. A 20-km-long oil spill due to an accidental collision between two cargo ships the day before is present [42,43]. A large low-backscattering area due to very calm sea state is also present. The incidence angles over the two ROIs are 33 • and 41 • , respectively.
The image shown in Figure 2c is acquired on 8 March 2017, off the coast of Fujairah in the western coast of the United Arab Emirates (Persian Gulf). Multiple oil slicks due to a seafloor leakage from the jack-up drill rig "Pasargad 100" are present [44]. The polluted sea area is conservatively estimated to be about 334 km 2 . The incidence angle over the ROIs is 35 • .
The image shown in Figure 2d was acquired on 11 March 2017, over the same area of Figure 2c, during the same oil spill event. Another low-backscattering area due to very low wind conditions (<3 m/s) is also visible [44]. The incidence angle over the ROIs are 35 • and 40 • , respectively.
The image shown in Figure 2e was acquired on 20 July 2019, around the Gotland island in the Baltic Sea between Sweden and Latvia. A diffuse low-backscattering area due to swirling green algae blooms covers most of the observed sea surface [45]. The incidence angle measured at mid-range is equal to 38 • .
The image shown in Figure 2f was acquired on 1 April 2018, over the Balikpapan Bay on the eastern coast of Indonesia. An oil spill due to a 25-m-underwater oil pipeline damaged the day before is present. The oil-affected area was estimated to be about 130 km 2 [46,47]. The incidence angle over the oil slick is 43 • .  Table 1. The regions of interest (ROIs) that refer to slick-free sea surface, oil slick, and look-alike are highlighted with blue-, green-, and orange-dashed boxes, respectively. An extra slick-free sea surface ROI, highlighted with a dashed red box, is also considered to analyze the intrinsic sea surface variability.

Experiments
In this section, the sensitivity of the AACF to the low-backscattering areas highlighted in Figure 2 is analyzed.

Experimental settings
The AACF over each ROI is estimated by averaging at least 10 AACFs evaluated according to the methodology described in Section 2. To allow a fair intercomparison of the estimated AACFs, they are normalized to their maximum value and the pedestal is set to zero, as it will be shown in Section 4.2, where the color coding is in accordance to the ROIs of Figure 2. In all the subsequent experiments, reference is made to slick-free sea surface ROIs (that will be labeled in blue as "Sea ref") at almost the same incidence angle of the low-backscattering areas, see dashed blue boxes in Figure 2. In each SAR image, an extra slick-free sea surface ROI (that will be labeled in red as "Sea"), see the dashed red boxes of Figure 2, is also selected to analyze the intrinsic sea surface variability. In the following, we refer to the slick-free sea surface AACFs as AACF s and to the AACFs estimated within the low-backscattering areas as AACF lb . In all subsequent SAR image analysis, the average contrast ∆, i.e., the difference between the slick-free sea surface and the low-backscattering area NRCS values in dB scale, is evaluated as follows: where the subscripts "s" and "lb" stand for slick-free sea surface and low-backscattering area, respectively. To quantify the deviation of the AACF lb with respect to the AACF s , the Euclidean distance, d E , and the percentage relative difference, D rel , are introduced. The Euclidean distance is defined as: where lb and s refer to AACF lb and AACF s , respectively, while i is the index that refers to the AACF samples of the selected low-backscattering and slick-free sea surface ROIs, each consisting of m = 1000 pixels. The percentage relative difference is defined as: where (d E ) lb is the Euclidean distance between the low-backscattering ROI and the reference sea surface one, while (d E ) s is the Euclidean distance between the two slick-free surface ROIs. This metric, as it is defined, is thought to analyze the deviation of the AACF estimated over low-backscattering areas from the corresponding slick-free sea surface reference one in order to filter out the intrinsic sea surface variability, i.e., induced by NRCS variability with respect to the azimuth angle (wind direction relative to the radar's azimuth look direction). All the above-mentioned quantitative parameters are listed in Table 2. Before proceeding further with the AACF analysis, since we are dealing with low-backscattering areas, a discussion on the effects of additive noise onto the backscattered signal is due. When dealing with S1 SAR data, the worst case nominal noise equivalent sigma zero (NESZ) is −22 dB. Hence, a received signal whose intensity is lower than NESZ may be uninformative. To perform a more accurate analysis, the SNR is used. The latter is not evaluated using the provided (worst case) NESZ but is estimated from the data according to [48,49]: where η and A are the noise calibration parameter and the NRCS calibration factor, respectively, which are pixel-dependent parameters provided in the ESA annotated metadata through a look-up-table. Then, the SNR is evaluated as follows:

Discussion
The SNR images evaluated over the whole SAR dataset are shown, in dB scale, in Figure 3, where land is masked in white. The average SNR values evaluated over the ROIs highlighted in Figure 2 are listed in Table 2.
Considering the reference slick-free sea surface ROIs (dashed blue boxes), the average SNR is always larger than about 7 dB, witnessing that the signal scattered off slick-free sea surface lies well above the NESZ along the whole SAR dataset. Considering the lowbackscattering sea areas, i.e., both oil slicks and look-alikes, the average SNR estimated over the corresponding ROIs significantly varies along the SAR dataset depending on the damping properties and incidence angle. Oil slicks (dashed green boxes), call for an average SNR that lies in the range of 1.5 dB-4.2 dB, while the look-alikes (dashed orange boxes), result in general in lower SNR values falling in the range of 0.5 dB-2.8 dB. Note that these values can be either lower or larger than the oil slick ones. Hence, for the purposes of the spectral analysis, it can be concluded that reference slick-free sea surface samples are noise-free and oil slick samples can be considered to have an average SNR large enough to not affect significantly the spectral analysis, while look-alike samples are partly contaminated by noise and, therefore, particular attention must be paid in their spectral analysis.
The first experiment refers to the SAR image ID 1, see Table 1. The two low-backscattering ROIs are characterized by a great ∆ value, i.e., >6 dB. The corresponding AACFs evaluated over the four ROIs are plotted in Figure 4a. The qualitative analysis clearly shows that the two AACF s are narrow, while the two AACF lb are wider. Although both AACF lb show a distinct behavior with respect to AACF s , the oil-covered one shows the largest broadening due to the oil damping properties and the reduction of the energy transfer from wind to the sea waves. It can be also noted that both ROIs call for a large enough average SNR, i.e., >1.5 dB. The low-wind ROI, whose SNR exhibits large spatial variability, calls for the largest average SNR (2.8 dB). The low-backscattering ROIs are characterized by d E values showing a significant deviation from the reference AACF s , i.e., 3.31 and 1.80 for the oil-covered and the low-wind ROI, respectively. The intrinsic sea variability, measured by computing the Euclidean distance d E between the "Sea ref" and "Sea" AACFs, is equal to 0.30, witnessing a good overlapping of the corresponding AACF s . Note that a similar result is obtained by randomly changing the position of the two slick-free sea surface ROIs. In conclusion, the d E values relevant to the low-wind and oil-covered ROIs are about 6 and 11 times greater than 0.30, respectively. The D rel values are 90.9% and 83.3% over the oil-covered and low-wind ROIs, respectively.       Table 1. Land is masked in white.  Table 1. Note that each AACF is normalized to its maximum value.
The second experiment refers to the SAR image ID 2, see Table 1. The two lowbackscattering ROIs show a great ∆ value, i.e., about 5 dB and 7 dB for the oil-covered and low-wind ROI, respectively. The four AACFs are plotted in Figure 4b. The two AACF lb are wider than the two AACF s , with the oil AACF lb resulting in the largest broadening. d E is equal to 0.55 and 0.14 for the oil-covered and the low-wind ROI, respectively. This occurs even if the oil-covered ROI is characterized by an average SNR significantly larger than the low-wind ROI (4.2 dB versus 0.7 dB, respectively), witnessing that the broadening is mainly due to the different spectral properties of the low-backscattering sea areas rather than to the additive noise. The intrinsic variability of the slick-free sea surface results in d E = 0.01, witnessing a very small AACF s variability. The d E values relevant to AACF lb are at least one order of magnitude larger than the ones relevant to AACF s . The D rel values are 98.2% and 92.9% over the oil-covered and low-wind ROIs, respectively, showing a remarkable sensitivity of the AACF to the low-backscattering ROIs. In fact, AACF lb significantly deviates from AACF s .
The third experiment refers to the SAR image ID 3 and 4, see Table 1. The two scenes refer to the same oil slick observed in two different dates and under different sea state conditions. The SAR image of Figure 2d also includes a low-wind ROI. Again, all the lowbackscattering ROIs are characterized by a great ∆ value, >4 dB. The corresponding AACF evaluated over the considered ROIs are plotted in Figure 4c,d. The AACF s are narrower than the AACF lb , with the oil AACF lb showing the largest broadening. Even in this case, the largest broadening of the oil-covered ROI cannot be attributed to the additive noise since the oil-covered ROI calls for an average SNR which is more than three times larger than the low-wind one (1.7 dB versus 0.5 dB, respectively). For all the low-backscattering ROIs, d E values are in the range of 0.38-2.68, always larger than the intrinsic sea variability that results in d E equal to 0.08 and 0.14 for SAR image ID 3 and 4, respectively. The d E values relevant to the oil and the low-wind AACF lb are more than one order of magnitude and just three times larger than the ones relevant to AACF s , respectively. These results are consistent with the D rel values which are about 94% and 63% for the oil-covered and low-wind ROIs, respectively.
The fourth experiment refers to the SAR image ID 5, see Table 1. The low-backscattering ROI shows a significantly great contrast, i.e., ∆ about 9 dB. The three AACFs are plotted in Figure 4e, where it can be noted that AACF s and AACF lb overlap despite the remarkable contrast. Although it seems that the sensitivity of the AACF to this low-backscattering area is negligible, this is most likely due to the very low wind conditions that apply over the whole sea area. As a result, the longer-wave part of the spectrum does not change significantly from the slick-free sea surface to the low-backscattering ROI. It is interesting to note that AACF s and AACF lb call for almost the same width even if the average SNR of the low-backscattering ROI is more than 10 times lower than reference slick-free sea surface one (0.8 dB versus 8.8 dB, respectively), witnessing that additive noise is not the factor driving the AACF broadening. However, although AACF s and AACF lb appear overlapped, they result in d E = 0.14, i.e., twice the intrinsic sea variability (d E = 0.07). The corresponding D rel value is 50.0% that, even though it is the smallest value among the whole SAR dataset, still represents a remarkable deviation.
The fifth experiment refers to the SAR image ID 6, see Table 1. Even in this case, the oil-covered ROI shows a pronounced contrast, i.e., ∆ = 7 dB. The three AACFs are plotted in Figure 4f. The oil AACF lb clearly deviates from the AACF s . This broadening results in d E = 1.39, 10 times larger than the intrinsic sea variability (d E = 0.14). In this case, the oil-covered ROI calls for D rel = 89.9%.
The experimental results suggest that the NRCS and the AACF call for a different sensitivity to the incidence angle. Hence, a deeper investigation is due. The behavior of ∆ and d E with respect to the incidence angle is analyzed with reference to the scene depicted in Figure 2f. The two metrics are evaluated fixing the oil-covered ROI while moving the slick-free see surface ROI along the range transect to span the available range of incidence angles. As suggested in [48,50,51], for the incidence angles of interest (31 • -43 • ), the slick-free sea surface NRCS variability is expected to be large, i.e., about 10 dB, while the oil NRCS variability is expected to be much smaller. Hence, ∆ is expected to vary significantly with θ. The experimental results are depicted in Figure 5, where a dB scale is used. One can note, as expected, that ∆ significantly varies with the incidence angle. This variability is practically negligible when d E is considered.

Conclusions
In this study, a spectral analysis of low-backscattering sea areas of both anthropogenic and natural origin was addressed. This is of great interest for the marine pollution community. The rationale of the proposed analysis relies on the capability of low-backscattering sea features, including oil slicks, algal blooms, and low-wind regions, to modify the full-range sea surface spectrum. Hence, the sensitivity of the SAR image autocorrelation function evaluated along the azimuth direction, namely AACF, was investigated.
The AACF estimated over low-backscattering areas was analyzed with reference to a slick-free sea surface AACF. The deviation of the low-backscattering AACF from the slickfree sea surface one was quantified according to two objective metrics, i.e., the Euclidean distance d E and the percentage relative difference D rel . A comparison was also made with respect to the intrinsic sea variability, i.e., the AACF deviation between two different slick-free sea surface areas.
Experiments, undertaken on 6 S1 VV-polarized SAR images, collected in IW dualpolarimetric imaging mode, where known oil slicks and low-backscattering sea areas due to natural phenomena are observed under low-to-moderate wind conditions (2 m/s-7 m/s) in a broad range of incidence angles (≈30 • -46 • ), showed that: • The AACF is sensitive to different low-backscattering areas, with d E and D rel values which are at least twice and 50% larger, respectively, than the intrinsic sea surface variability; • Among the low-backscattering sea areas, the oil slicks exhibit the largest AACF deviation with respect to the reference slick-free sea surface, with a maximum of d E = 3.31 and D rel = 98.2%; • The additive noise does not play a key role in broadening the AACF; • The AACF is practically independent on the incidence angle while the backscattering contrast depends on it.

Acknowledgments:
The authors thank the European Space Agency for providing Sentinel-1 SAR data free of charge through the Copernicus Scientific Hub. This paper was written during the COVID-19 pandemic in Italy and is dedicated to all people working in the healthcare system.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: