A Ka-Band Wind Geophysical Model Function Using Doppler Scatterometer Measurements from the Air-Sea Interaction Tower Experiment

: Physical understanding and modeling of Ka-band ocean surface backscatter is challenging due to a lack of measurements. In the framework of the NASA Earth Ventures Suborbital-3 Submesoscale Ocean Dynamics Experiment (S-MODE) mission, a Ka-Band Ocean continuous wave Doppler Scatterometer (KaBODS) built by the University of Massachusetts, Amherst (UMass) was installed on the Woods Hole Oceanographic Institution (WHOI) Air-Sea Interaction Tower. Together with ASIT anemometers, a new data set of Ka-band ocean surface backscatter measurements along with surface wind/wave and weather parameters was collected. In this work, we present the KaBODS instrument and an empirical Ka-band wind Geophysical Model Function (GMF), the so-called ASIT GMF, based on the KaBODS data collected over a period of three months, from October 2019 to January 2020, for incidence angles ranging between 40 ◦ and 68 ◦ . The ASIT GMF results are compared with an existing Ka-band wind GMF developed from data collected during a tower experiment conducted over the Black Sea. The two GMFs show differences in terms of wind speed and wind direction sensitivity. However, they are consistent in the values of the standard deviation of the model residuals. This suggests an intrinsic geophysical variability characterizing the Ka-band surface backscatter. The observed variability does not signiﬁcantly change when ﬁltering out swell-dominated data, indicat-ing that the long-wave induced backscatter modulation is not the primary source of the KaBODS backscatter variability. We observe evidence of wave breaking events, which increase the skewness of the backscatter distribution in linear space, consistent with previous studies. Interestingly, a better agreement is seen between the GMFs and the actual data at an incidence angle of 60 ◦ for both GMFs, and the statistical analysis of the model residuals shows a reduced backscatter variability at this incidence angle. This study shows that the ASIT data set is a valuable reference for studies of Ka-band backscatter. Further investigations are on-going to fully characterize the observed variability and its implication in the wind GMF development.


Introduction
Microwave wind scatterometers are active instruments designated to estimate ocean surface vector winds based on measurements of ocean surface backscatter [1]. The latter is obtained from the magnitude of the power scattered back to the radar by small-scale capillary waves that respond to the Bragg resonance mechanism (Bragg ocean surface waves) [2], with modifications due to tilting by large waves and potential contributions January 2020, for incidence angles ranging between 40 • and 68 • . A major advantage of the KaBODS system is that it is able to scan in the azimuthal direction, offering a large range of directional modulation. This helps to accurately describe the surface backscatter with respect to the wind direction.
We have compared the ASIT GMF to the Ka-band Dual Co-Polarized Model (KaDP-Mod) developed by Yurovsky et al. [24]. The latter was developed from measurements collected during a tower experiment conducted over the Black Sea. This comparison aims to analyze the KaBODS backscatter sensitivity with respect to the wind speed and wind direction. In addition, the analysis of the difference between the model and the actual data using both GMFs aims to improve understanding of the geophysical mechanisms driving the backscatter at Ka-band. The results of this work will be used as guidance for our future investigations to fully assess the Ka-band wind scatterometry capabilities in order to improve the current DopplerScatt operational GMF and the accuracy of wind retrievals.
An important result that we found is that the observed backscatter is characterized by a geophysical variability and an analysis of the potential source of this variability is discussed in this work. Unfortunately, only single-polarization data could be collected, leading to additional uncertainty in the identification of the sources of back scatter variation mechanisms. It is important to notice that the backscatter variability is not a peculiarity seen only at Ka-band. Several works have attempted to characterize the surface backscatter at Ku-band and C-band frequencies, and it was found that consideration of geophysical components such as, for instance, swell, bound waves, sea spray and breaking waves are needed to reconcile the models with the actual observations [28][29][30].
The KaBODS system and the ASIT anemometer data are described in Sections 2.1 and 2.2, respectively. The methodology used to process the KaBODS backscattered power to convert it into the corresponding ocean surface backscatter is presented in Sections 2.3 and 2.4. The development of the ASIT GMF is described in Section 2.5. The comparison between the ASIT GMF and the KaDPMod, in terms of wind speed and relative wind direction, is shown in Sections 3.1 and 3.2, whereas the analysis of the differences between the models and the data is given in Section 3.3. The discussion of the observed backscatter variability and the conclusions of the results are given in Sections 4 and 5, respectively.

KaBODS Instrument Description
The KaBODS scatterometer is a Ka-band vertically polarized Frequency Modulated Continuous Wave (FMCW) radar. It was located on the main platform of the ASIT at 13.1 m above the mean sea level. This platform pointed toward 245 • (clockwise) relative to the true North. The KaBODS radar has two standard gain pyramidal horn antennas, one of which is used to transmit while the other is used to receive. Each antenna has a boresight angle of 60 • . Through range compression, the radar collects measurements at different incidence angles, ranging between 20 • and 83 • . Originally, both the transmit and receive antennas were configured to be in vertical polarization (VV). On 20 November 2020, a strong wind event rotated the receiving antenna by 45 • , allowing a configuration in vertical-horizontal polarization (VH). Given the small value of the cross-polarized return from the ocean at low to moderate wind speed, we do not expect that these returns will corrupt our co-polarized estimates significantly. Yurovsky et al. [24] showed that the predicted errors associated to a cross-polarized contamination are small for wind speeds lower than 25 m/s, particularly for the VV return. Therefore, in line with these results and considering that the range of wind speeds used in this work does not exceed 18 m/s, we decided to ignore the VH contribution. In addition, the data before and after the polarization switch are consistent.
The system details are presented in Table 1. With a bandwidth of 100 MHz, the radar range resolution, which represents the radar's ability to distinguish between targets, is ∼1.5 m, which is sufficient to resolve surface waves. At the Ka-band radar frequency, the radar is sensitive to Bragg surface wave components of ∼0.8 cm wavelength [2].
The KaBODS radar amplitude measurements were calibrated with calibration targets deployed at various ranges in the far field. The system did not have an internal calibration gain loop, so transmit power variations, due to temperature, for instance, were not corrected, leading to potential uncertainty in the absolute value of the radar cross section, but with a lesser impact on the study of azimuth and angular modulations, since these are relative measurements collected under the same environmental conditions. We note, however, that our absolute calibration is largely consistent with the results reported in [24].
The KaBODS azimuth look direction is controlled by a motor which allows the radar to steer in azimuth with about 160 • range. The radar collects continuous data at each azimuth location for five minutes with a collection cycle of one hour. An example of the nominal KaBODS azimuth positions within a one-hour cycle is shown in Figure 1. The KaBODS dwell time and configuration were set to maximize the KaBODS time within the footprint of a polarimeter, which was installed on the ASIT. This results in a longer dwell time at the radar azimuth position of 203 • , as shown in Figure 1. For this study, we used KaBODS L0 products from October 2019 to January 2020. The L0 products contain the Single Look Complex (SLC) radar data, as a function of time and radar range, as well as the radar azimuth position relative to the normal of the ASIT platform.

Meteorological and Oceanic Data Collection
The wind measurements used in this study were obtained from six sonic anemometers located on the ASIT at different heights. Three anemometers are located on the ASIT Rohn Tower, above the main platform, between 12 and 23 m. The other three anemometers were located on the ASIT mast extending off the end of the main platform. A description of the ASIT instruments with their corresponding location can be found in [27,31]. The wind measurements were converted to a common height of 10 m using the Coupled Ocean-Atmosphere Response Experiment (COARE) bulk algorithm version 3.5 [31]. The measurement heights were adjusted to account for changes in the sea surface height relative to the ASIT, mainly due to tides. To develop the ASIT GMF, we used these 10 m winds, relative to water, adjusted for neutral conditions, known as equivalent neutral wind [31]. To compute the wind speed relative to water, the surface currents available from the Acoustic Doppler Current Profiler (ADCP) located nearby the ASIT were used. The distribution of the equivalent neutral wind speed (U10 N ) and wind direction relative to the radar azimuth (given in meteorological convention) (φ) is shown in Figure 2a,b, respectively. During the experiment, a wide range of wind intensities was encountered, between 1-18 m/s, with the majority of the data ranging between 5-8 m/s. As noted before, thanks to the radar's capability of scanning in the azimuth direction, the data set covers a wide range of relative wind-radar boresight azimuth directions.
For this study, in addition to the wind measurements, we also had measurements of the sea surface temperature (SST), the dominant wave period (T peak ), the significant wave height (H s ) and the precipitation rate (P r ) (Figure 2c-f). T peak represents the period of the most energetic ocean wave component, whereas H s represents the mean height of the highest one-third of the waves [32]. These parameters have been used to filter out the backscatter data that can potentially be compromised by rain, as well as environmental conditions where insufficient data were available for a reliable result. Hereafter, we use only KaBODS data corresponding to T peak < 10 s, H s < 2.5 m and rain rate <1 mm/h. As can be seen from Figure 2, this only excluded a relatively small number of observations.

KaBODS Signal Power Processing
The KaBODS total return power was computed as the magnitude squared of the SLC data, corrected with a calibration constant, as described below. Two major contributions affect the return power, namely, the speckle noise and the thermal noise. Speckle noise is the result of random phase variations in the interaction between the radar signal and the ocean surface [33]. Such phase variations are mainly introduced by the small-scale properties of the surface, the random motion of the point scatters and the variation of the distance between the radar and the point-target. In order to reduce the speckle noise, a multi-looking approach is usually used, averaging independent samples (looks). To this end, we have averaged the KaBODS received power over N L = 50 looks (radar pulses), that results in averaging over 50 ms. The resulting cross section speckle is reduced by a factor of approximately 10 log 10 1 + N −1/2 L ≈ 0.6 dB for each 50 ms measurement. After averaging for 5 min, the estimated speckle variability (not including geophysical variability) would be in the order of ∼0.06 dB.
The thermal noise power comes from the instrument itself and it introduces a bias in the return power [34]. It is dependent on the system temperature and its bandwidth. Figure 3a and b show the KaBODS multi-looked return power as a function of the radar range and incidence angle, respectively. Three different days/times are shown, corresponding to three different wind speed values, 3.7 m/s, 6.4 m/s and 10.6 m/s. As expected, the power increases as the wind speed increases. However, at ranges higher than 86 m, corresponding to incidence angles greater than 82 • , the power does not change with the wind intensity, it appears rather constant, suggesting the presence of a noise component. Therefore, we have estimated the value of the KaBODS thermal noise by averaging the multi-looked power starting at an incidence angle of 82 • . We then computed the power associated only with the signal (echo), by subtracting the thermal noise component (P N ) from the total return power (P r ) [34]: An important aspect shown in Figure (Figure 3b), respectively. These peaks are not significantly affected by the wind speed and they are probably due to a RF-leakage from either the KaBODS system itself or the other instruments on the ASIT, rather than an effect from the ocean surface. Therefore, in this paper, we focus only on the range between 17 m and 35 m, i.e., between 40 • to 68 • , to avoid possible instrument contamination.  Figure 4 shows the KaBODS signal-to-noise ratio (SNR) as a function of the incidence angle for different wind speed ranges. The SNR shows a maximum around 60 • , which is consistent with the radar look angle and main lobe of the antenna pattern. It also shows a second peak between 45 • and 50 • , coincident with an antenna side lobe. Additional details on the antenna pattern are given in Section 2.4.

The KaBODS Surface Backscatter Computation
In order to convert the signal power P s into the corresponding ocean surface radar backscatter (σ • ) we have used the radar equation as shown in (2): where K cal is a radar calibration constant (computed using known corner reflectors) and X is the X−factor. The X−factor is a pre-computed constant that is needed to convert measured return power into ocean surface backscatter. It accounts for both the radar geometry and the antenna pattern. We have modeled the X−factor as shown in (3).
where P t is the transmitted power, G t (θ, α) and G r (θ, α) are the receive and transmit antenna gains as function of elevation θ and azimuth α, respectively. The KaBODS antenna patterns were not calibrated independently, but, since they are commercial standard gain horns used for calibration, their pattern as a function of angle is well known and within the tolerances for this study. The term χ(r − r 0 ) represents the point-target response (PTR), as a function of the range r centered at r 0 . The definition of χ depends on the radar processing of the received signal from a point target. As an FMCW radar, the KaBODS transmits a signal that is frequency-modulated. This allows for pulse compression on the received signal, which is a signal processing technique usually used to improve the radar range resolution as well as the signal-to-noise ratio [1]. However, this technique also enhances the presence of the side lobes in the received compressed pulse in the time domain. These side lobes are in addition to the main pulse from the target and weighting functions are used to reduce their amplitude [35]. For the the KaBODS system, a Hanning weighting function has been used to produce the SLC measurements. Therefore, χ is defined as a Hanning-weighted PTR. Two different configurations of the X−factor were implemented to account for both VV and VH-polarizations, in order to be consistent with the received antenna rotation configuration. The resulting X−factor as a function of the incidence angles for both VV and VH polarization is shown in Figure 5. Here, we show the X−factor for the whole range of incidence angles available from the KaBODS system. The peak corresponding to the antenna pattern main lobe is located at 60 • , whereas two side lobes are located between 45-50 • and 70-75 • . This X−factor model is consistent with the KaBODS power shown in Figure 3b. It also shows that the power spike seen at 23 • does not correspond to a sidelobe of the antenna pattern. The spikes seen at 70 • and 75 • appear to be caused by the energy entering through the sidelobe of the antenna pattern. However, since the power does not significantly change with the wind, we believe that those angles should be excluded from the analysis. The KaBODS calibration constant K cal characterizes the radar system and its estimation is essential in correctly computing σ 0 from the received signal. The KaBODS system was calibrated on the ground by use of the three corner reflectors placed at 12 m, 26 m and 40 m distance from the radar, as described in [36]. These reflectors have a known radar cross section, so that by recording the received power and measuring the transmitted power, the calibration constant can be estimated. The estimated K cal value of the KaBODS system is 182.8 dBW and we have used this value in the data processing.
The measured data were processed to estimate the corresponding σ 0 measurement. We performed a quality check of the obtained σ 0 time series and we found cases with σ 0 spikes exceeding 0 dB. Each spike was visible at different incidence angles and according to the anemometer winds, they were associated with low wind speed conditions. This suggests that such high backscatter values were likely due to particular events happening at the sea surface (e.g., boats), within the radar footprint, rather than wind/wave effects. For this reason, we decided to filter out values of σ 0 higher than −5 dB, to exclude corruption of measurements by bright targets.
An additional mechanism which can also increase the variability of the surface backscatter at the spatial resolution of the KaBODS system is the modulation induced by the long waves. Additional details on the wave-induced backscatter modulation are given in Section 4. The variability caused by large wave modulation is inherently present in large-footprint spaceborne scatterometer measurements. Since we want our results to be applicable to this situation, we replace spatial averaging with temporal averaging, we averaged the surface backscatter over the azimuth dwell. This corresponds to averaging over a time period of about five minutes: this corresponds to spatial averaging over ∼1.6 km to ∼2.6 km for waves ranging between 20 m and 50 m, which is comparable to spaceborne spatial scales. Each averaged σ 0 value (hereafter, <σ 0 >) was then associated to the azimuth mean of the corresponding azimuth subset and it has been collocated with the anemometer measurements. This collocated data set was then used for the GMF development.

Geophysical Model Function Development
We derived the ASIT GMF as a function of the radar incidence angle (θ), the 10 m equivalent neutral surface wind speed (U10 N ) and the wind direction relative to the radar azimuth (φ), by using the functional form as a described in Yurovsky et al. [24], as shown in (4).
where the A i coefficients are presented as functions of θ with a dependence on log U10 N : To build the GMF, we derived the c mik coefficients by mapping the KaBODS < σ 0 > measurements to the collocated anemometer wind speeds and directions. We used a leastsquares fitting method with a Cauchy loss function to weaken the influence of possible outliers. The ASIT anemometer array and ancillary instruments have the large advantage that wind stress can be calculated and converted to an estimate of neutral winds at 10 m [31], and we use these to construct our GMF. We noticed that the functional form given in (4) does not fit well over a large range of wind speeds. In particular, the surface backscatter at very low wind speeds tends to drop faster than the prescribed log parameterization [37]. Therefore, we excluded measurements at wind speeds lower than 3 m/s, in favor of a good fit for moderate and higher winds. With respect to the wind direction, we noticed that the ASIT is a point of distortion for the KaBODS backscatter measurements when the wind is blowing away from the radar. These data correspond to a wind direction of 180 • relative to the radar azimuth (downwind direction). This distortion is due to the fact that the ASIT shadows the water surface, affecting the surface roughness. Therefore, given the few degrees of freedom of the fitting model, we have decided to exclude from the training data set the backscatter measurements with a relative wind direction within 30 • from the downwind direction, in favor of a good fit up to 150 • , which is known to be free from shadowing effects. Moreover, the radar position on the ASIT main platform causes the radar beam, with the azimuth direction of ∼262 • (relative to the true North), to point into the main platform. This results in contaminated backscatter measurements which have also been excluded from the analysis. The values of the obtained ASIT GMF coefficients (c mik ) are shown in Appendix A.

Wind Speed Sensitivity
We first analyzed the sensitivity of the KaBODS backscatter with respect to the wind speed and compared results for the ASIT GMF with the KaDPMod GMF. Figure 6 shows the KaBODS <σ • > binned over 1 m/s range (blue dots with error bars) along with the ASIT GMF (solid orange line) and the KaDPMod GMF (dashed black line) as a function of the wind speed. The results are shown for φ = 0 • (upwind), φ = 90 • (crosswind), φ = 150 • and φ = 180 • (downwind), at three different incidence angles, 40 • , 50 • and 60 • . The KaBODS measurements are modulated by the wind speed and the ASIT GMF is able to describe this wind sensitivity with a residual absolute mean bias (µ ASIT ) ranging between 0.15 dB to 0.95 dB and a residual standard deviation (σ ASIT ) between 1.96 dB and 3.17 dB for relative wind directions up to 150 • . These values were calculated as the mean bias and the standard deviation of the differences between the model and the actual measurements in log space. In downwind direction, the agreement between the ASIT GMF and the measurements is not as good as at 150 • . This could be ascribed to the fact that the ASIT GMF has not been trained at this angle, as explained in the previous section. However, this could also be an indication of the potential shadowing effects. In addition, a disagreement is observed for wind speeds higher than 12 m/s, mostly due to the lack of data used to train the GMF at high wind speed conditions.
Overall, the KaDPMod and the ASIT GMF are consistent, however, the KaDPMod describes a slightly different wind speed sensitivity with respect to the KaBODS backscatter and in turn to the ASIT GMF. The residual absolute mean bias (µ Y16 ) ranges between 0.43 dB to 3.63 dB with a mean standard deviation (σ Y16 ) ranging between 2.05 dB to 3.22 dB. The values of µ Y16 are higher compared to µ ASIT . This is expected since the ASIT GMF is built by fitting the data. Notice that these differences vary with wind speed and incidence angle, and cannot be ascribed to an overall calibration difference between the two system: in general, the absolute calibration seems consistent between the two data sets. The σ Y16 values are mostly consistent with the σ ASIT values for different relative wind directions and incidence angles. An offset is seen between the KaDPMod and the ASIT GMF at 40 • , particularly in the upwind and crosswind directions and this offset decreases at higher incidence angles.  Figure 7 shows the KaBODS surface backscatter with respect to the relative wind direction, at wind speeds of 5 m/s (left), 7 m/s (middle) and 9 m/s (right). The dots and error bars represent the mean with the corresponding standard deviation of the KaBODS backscatter binned over 10 • -direction. To this end, at a given incidence angle, firstly we binned the data over 1 m/s bin centered at 5, 7 and 9 m/s. Bins of 1 m/s are used to ensure enough data points. Each wind speed bin was then sorted in 10 • -direction bins where the corresponding backscatter mean and standard deviation were computed.

Wind Direction Modulation
The KaBODS backscatter shows a clear modulation with the wind direction and as expected, the ASIT GMF is able to reproduce such modulation with a relatively low absolute value of µ ASIT ranging between 0.09 dB to 0.96 dB and σ ASIT between 2.16 dB to 3.71 dB. The KaBODS backscatter shows a stronger directional modulation at an incidence angle of 50 • and 60 • than 40 • . Indeed, as also described by the ASIT GMF, the directional modulation tends to be smoother at 40 • , regardless of the wind intensity.
The comparison between the KaDPMod and the ASIT GMF shows that the KaDPMod GMF does not fully describe the directional modulation of the KaBODS measurements. At 50 • , the KaDPMod GMF shows a slightly smoother modulation compared to the ASIT GMF and the KaBODS data. This effect worsens at 60 • . This uncertainty in the KaDPMod GMF directional modulation can be ascribed to the data used to train this GMF [17]. Due to the nature of the Black Sea tower measurements, the collected data set did not have a wide azimuth variability and this may have affected the model sensitivity with respect to the relative wind direction. In our setup, nearly co-temporaneous azimuth variability was guaranteed by completing azimuth cycles within one hour.
The KaDPMod residuals absolute mean bias µ Y16 ranges between 0.08 dB to 3.25 dB and the standard deviation σ Y16 between 2.48 dB to 3.77 dB. As for the wind speed sensitivity, the models show an offset at 40 • , with an incidence angle dependence, but they show similar values of the residuals mean standard deviations for different wind speeds and incidence angles.

Model Accuracy
In order to further evaluate the accuracy of both models with respect to the KaBODS <σ 0 > data, we computed the statistics of the model residuals for the whole data set along with the root mean square difference and the correlation coefficient. The distribution of the model residuals is shown in Figure 8. The analysis is carried out for a range of incidence angles: 40 • , 45 • , 50 • , 55 • , 60 • and 65 • . The results are shown in Table 2. Overall, both models appear well correlated with the data, with a correlation coefficient ranging between 0.71 and 0.88 for the ASIT GMF and 0.69 to 0.86 for the KaDPMod. As expected, the ASIT GMF better agrees with the KaBODS <σ 0 > measurements then KaDPMod, as shown by lower values of the mean bias. We note that these statistics show a trend with incidence angle. In particular, they tend to improve at 60 • , where we observe the lowest values of mean bias, standard deviation, root mean square errors and the highest correlation. However, they tend to worsen again at 65 • . This particular dependence on the incidence angle could be related to the shape of the antenna pattern. However, because the SNR is high over the whole incidence angle range, these results may also suggest that the surface wave components that better respond to the Ka-band radar are those observed at (and close) an incidence angle of 60 • . It is interesting to notice that the offset between the ASIT GMF and the KaDPMod also shows a similar trend, as it tends to disappear at 60 • (Figure 8). The nature of such incidence angle dependence is currently under investigation.
Both models show a residual distribution that is Gaussian in log-space. It is interesting to see that, consistent with the two previous analyses, the standard deviation values do not significantly change between the two models when using the whole data set to compute the residuals. We believe that such variability is not related to the model accuracy, but rather to a backscatter variability that is not fully captured by the parameterization. Since we have processed the data to reduce the modulation due to the speckle and thermal instrument noise, we believe that such variability has a geophysical source.  Table 2. The mean bias (µ dB) and the mean standard deviation (σ dB) of the residuals of both the ASIT GMF and the KaDPMod from Yurovsky et al. along with the root mean square difference (ε dB) and the correlation coefficient (R).

Discussion
Based on the results of the model residuals, our main focus has been understanding the source of the observed residual backscatter variability.
As mentioned in Section 2.2, we used the winds as water-relative, i.e., corrected by the surface currents available from the ADCP. We also performed the same analysis using the earth-relative winds, but no significant difference was seen as we obtained very similar results. This is due to the fact that the typical currents on the site were lower than 30-50 cm/s, which is much smaller than the wind speeds encountered during the experiment. Figure 9 shows the values of the velocity components obtained by the ADCP. Therefore, this suggests that possible effects due to the surface currents on this site are expected to be small.
An important source of the Ka-band backscatter modulation at moderate and high incidence angles, is hydrodynamic modulation [38]. This phenomenon is caused by the backscatter modulation induced by the long waves. This occurs since the short wave components, locally generated by the wind, are modulated by the long wave orbital velocity [39]. These long waves are usually generated further away, but they propagate into the observed scene and are referred to as swell [32]. According to Yurovsky et al. [38], the backscatter modulation due to long waves tends to decrease as the wind speed increases. As explained in Section 2.3, to account for the effect of the long-wave modulation, we averaged the backscatter over the azimuth dwell time of about five minutes. However, a residual variability could still be associated with longer period waves, where 5 min may not be sufficient to average out their contribution. Therefore, to further assess the impact of the wave-induced modulation, we computed the wave age as c p /U10 N , where c p is the phase velocity of the wave component with period T p [40]. The wave age represents a measure of the wave development status and we used it to identify swell-dominated data [31]. In particular, we used a wave age value of 1.2 to identify the data corresponding to fully developed sea, which is the condition of equilibrium between the wave energy source and dissipation. Wave age values greater than 1.2 identify decaying seas dominated by swell, whereas values of less than 1.2 represent developing seas [31]. As shown in Figure 10, the majority of the ASIT data with wind speed between 5-8 m/s fall under swell conditions. We then computed the model residuals by filtering the data with wave age higher than 2.2, in order to reduce the amount of swell-contaminated data. The results show that the standard deviation of the model residuals does not significantly change. We also used different wave age thresholds, such as, 2, 1.8 and 1.6 to filter the data. However, we observe that even when filtering a larger amount of swell-contaminated data, the standard deviation is not significantly reduced. This suggests that the long-wave modulation is not the primary source of the KaBODS residual backscatter variability. Similar results were also obtained by performing the same analysis over data points within easterly and south-westerly wind directions to provide long fetches and resemble open-ocean conditions.  Past works have shown the importance of including the sea surface temperature (SST) in the Ku-band geophysical model function in order to improve the accuracy of wind retrieval from scatterometers [13]. The SST can modify the roughness of the surface through modification of the surface water viscosity [41]. It has also been shown that SST perturbations affect the wind stress field through the stability of the marine atmospheric boundary layer [42,43]. Vandemark et al. [12], showed that the ocean backscatter at Ka-band for incidence angles near nadir and at low wind speeds varies between cold and warm water. The SST is also expected to have an impact at larger incidence angles. As shown in Figure 2, the SST during the ASIT experiment covered a range between 4 • C to 14 • C. Therefore, in line with these previous studies, we analyzed the averaged KaBODS backscatter with respect to the SST for low-medium wind speed. However, at the KaBODS spatial resolution, we have not found a clear signature of an SST contribution to the observed residual variability. Given the uncertainty in the radar gain as a function of temperature, and the weakness of the observed signal, we are not confident in concluding anything definite at this stage and leave this for future research.
One feature that we observe in the KaBODS backscatter data set is that the KaBODS time series (i.e., the measurements before averaging over the azimuth dwell period) show a spiky structure. As a consequence of this structure, the backscatter distribution is not Gaussian, but rather skewed in linear space. Figure 11a shows the backscatter time series at an incidence angle of 60 • , for two different days/times in similar wind conditions, corresponding to wind speed between 8-9 m/s, in the upwind direction and in the absence of swell. The blue trace shows an example of the typical KaBODS backscatter spiky structure, with the corresponding distribution (Figure 11b). In this case, the distribution skewness is 1.50. We found that in presence of spikes with a magnitude much larger than the background backscatter, as shown from the red trace, the value of the distribution skewness significantly increases. For this case, the skewness increases from 1.50 to 5.46. Similar results are obtained at different wind speed conditions, as shown in Figure 11c,d, for a wind speed between 6-7 m/s. In this case, the skewness increases from 1.45 to 4.3. Yurovsky et al. [44] suggested that the presence of high amplitude spikes, making the distribution strongly skewed in linear space, is likely due to wave breaking events. In particular, they showed that at moderate to high incidence angles, these backscatter spikes represent a non-negligible non-polarized backscatter component associated to quasi specular reflection from breaking facets. Although the contribution of this non-polarized component is shown to be stronger at horizontal polarization than vertical polarization [44], such spikes are also visible in vertical polarization, as shown in [45], and they can also be associated to enhanced roughness of the ocean surface from post wave-breaking events, as described in [46]. Moreover, although the studies of wave breaking effects usually focus on moderate to high wind speed conditions, it was suggested in [38] that the capillary waves observed at Ka-band at low wind speed conditions are generated as parasitic capillary waves from the breaking of short-gravity waves [3], showing a possible breaking contribution at low winds as well. Therefore, in line with these studies, we believe that the large spikes observed in the KaBODS backscatter show a contribution due to wave breaking effects. Further investigations are now ongoing to fully characterize the observed spiky structure as well as the impact that the wave breaking events have on the KaBODS backscatter distribution.

Conclusions
In this work, we present a Ka-band wind geophysical model function based on the data collected during the WHOI ASIT experiment. The ASIT was equipped with a Ka-band Doppler scatterometer along with six sonic anemometers, allowing direct estimation of neutral winds and providing a unique data set with collocated ocean surface backscatter and wind measurements.
The ASIT GMF has been compared to the existing KaDPMod GMF, which was developed from measurements taken during a different tower experiment conducted over the Black Sea. We have analyzed the two GMFs with respect to the wind speed and relative wind direction and we have analyzed the accuracy of both models with respect to the Ka-BODS data set. Both GMFs are well correlated with the data, with a correlation coefficient ranging from 0.7 to 0.9. The correlation coefficient values are slightly larger with the ASIT GMF than the KaDPMod. This is expected by construction since the ASIT GMF was trained using the data. However, the KaDPMod shows a slightly different sensitivity to the wind speed and a smoother directional modulation than the ASIT GMF, particularly at larger incident angles (>50 • ). We believe that the KaDPMod smoother modulation is probably due to the fact that the data used to train the model was sparse in terms of azimuth variability and this may have affected the accuracy of the KaDPMod directional component.
We also observe that the two GMFs have a consistent residual variability that we concluded as being mainly associated with a variability of the surface backscatter. To analyze this residual backscatter variability with respect to the long-wave induced modulation we have taken into account the wave age parameter and we concluded that this modulation is not the primary source of the observed variability. We have also analyzed the possible effect associated to the SST, however, we could not see, at this stage, a clear signature associated with the SST variation at the spatial resolution of the KaBODS system. One feature that we observe is that the surface backscatter has a spiky structure and the presence of spikes with magnitude greater than the background makes the distribution of the backscatter strongly skewed in linear space. Such results are consistent with what has been observed in previous studies, suggesting the presence of a contribution due to wave breaking events.
Moreover, we observe that the statistics of both models as well as the backscatter residual variability show a trend with the incidence angle. They tend to significantly improve at an incidence angle of 60 • . Although the KaBODS antenna has a peak gain at 60 • , it is worth noticing that the SNR is high within the whole incidence angle range, and the variability is a relative measurement over a short time, and not expected to be affected by the absolute gain. Therefore, the observed trend may indicate that the wave components observed at an incidence angle of 60 • (and closer) are those that better respond to the Kaband frequency in terms of wind interaction, according to the Bragg scattering mechanism.
The results of this work show that the ASIT data set is extremely valuable as a reference data set for investigations of new and existing GMFs at Ka-band. The results of the comparison of the GMFs with the actual data provide guidance for future efforts to fully characterize the surface backscatter mechanisms at Ka-band. We are currently investigating the impact of the incidence angle in order to explore the optimal viewing geometry for future mission design. We are also analyzing the potential signature of breaking waves in the surface backscatter distribution and the possible SST effects.