Assessment and Improvement of Sea Surface Microwave Emission Models for Salinity Retrieval in the East China Sea

: Accurate prediction of sea surface emission is the key for sea surface salinity retrieval from satellite microwave radiometer. In order to retrieve salinity from satellite observation, several sea surface microwave emission models have been developed based on theoretical or empirical methods and validated by in-situ measurements in di ﬀ erent regions. However, their performances are still unclear in the Chinese coastal waters. In this study, based on two cruises measurements in the East China Sea (ECS), including the brightness temperature measured by a shipborne microwave radiometer and other auxiliary data (sea surface salinity, sea surface temperature and wind speed), the performances of di ﬀ erent sea surface emission models are tested. The results showed that the developed models provide fairly good accuracy in predicting brightness temperature; for example, the accuracy of small perturbance / small scale approximation model (SPM / SSA), two-scale model (TSM) and empirical model is in the range from 0.6 K to 3 K. Moreover, the TSM and empirical models are further improved by optimizing the model parameters in the ECS. Finally, the sea surface salinity were retrieved from shipborne measured data based on the improved models, and the results show that the root mean square (rms) di ﬀ erences between retrieved and in-situ sea surface salinity is about 0.4 psu, indicating the signiﬁcant improvement by the regional model parameters.


Introduction
Salinity is one of the most important marine parameters, which controls many processes such as physical circulations, biogeochemistry dynamics from regional to global ocean [1,2]. Although drifters and buoys, together with cruises, have accumulated large amount in-situ water salinity data in different regions, it is still difficult to monitor global ocean salinity with high temporal-spatial resolution by the in situ measurements. In the past decade, with the successful launch of the soil moisture and ocean salinity (SMOS) satellite by the European Space Agency (ESA) [3], the Aquarius/SAC-D satellite [4], and the soil moisture active passive (SMAP) satellite [5,6] by the National Aeronautics and Space Administration (NASA), global sea surface salinity (SSS) observations from space have become possible and a significant improvement has been made to understand the ocean dynamics and climate change.
The early concept of remote sensing of SSS has been demonstrated in the late 1970s with observations by Skylab [7] and two airborne L-band radiometry experiments [8,9]. At the end of 1990, two airborne microwave interferometers, the electronically scanned thinned array radiometer The radiometer was mounted on the stabilized platform, which was about 6 m high above the sea surface and provided dual polarization measurements of brightness temperature at 5 fixed incidence angles (16°, 26°, 36°, 45° and 60°), and the circular-shaped field of view on the sea surface has the radius in the range of 3 m to 6 m according to different incidence angles (for example, the radius of field of view is about 3 m with incidence angle at 16°). During the cruises, some stations were measured in daytime, which were hampered by sunlight. In order to avoid radiation coming from the Sun, the antenna of the radiometer was pointed away (backwards) from the sun and the radiometer's field of view was in the shadow of the ship. At few stations, the measurements were affected by the Sun, which was removed easily because the intensity of radiation emitted by the Sun was high (~10 K or higher) and the reflectivity of the sea surface was large (in the order of 0.7). Although the effect of radio-frequency interference (RFI) is another possible source of error for the measured brightness temperature, the shipborne radiometer was deployed about 6 m above the sea surface (unlike the airborne radiometer), which could avoid potential sources of RFI. Moreover, the radiometer was pointed towards the east direction at most of the time to avoid the RFI effect from the land. In addition, it was observed that, even under the wind speeds up to 11 m/s, the foam coverage in the radiometer's field of view was less than 5%, and thus, the foam effect was neglected in this study. During the cruises, the following instruments were deployed: a dually polarimetric L-band radiometer designed and produced by Radiometrics Corporation (Figure 2a), and a stabilized platform designed and produced by Watson Corporation to keep radiometer stable during the whole cruise ( Figure 2b). Moreover, to obtain the auxiliary data matched with brightness temperature measurements, two CTD (Conductivity Temperature Depth) sensors from SeaBird Electronics and RBR Corporation and a meteorological station from RM Young anemometer were used to measure SSS, SST and wind speed. The details of these instruments/experiments are briefly given below.
(1) L-band Radiometer PR-1475: The PR-1475 is a dually polarimetric L-band radiometer designed and produced by Radiometrics Corporation. The antenna is a potter horn, which is a narrow band low sidelobe antenna, with a half-power beam width of 30 • and measured side lobe levels at −23 dB. The potter horn design characteristically has sidelobes down more than 30 dB. The radiometric sensitivity is 0.1 K to 1 K depending on the integration time from 10 to 2500 msec and bandwidth of the radiometer is 25 MHz. (2) Stabilized platform: The stabilized platform is a particular model of the Watson product intended to hold the PR-1475 payload aligned to downward direction on a moving vehicle. This system includes a mounting bracket suited to hold the payload for rotation about its center of gravity, a mounting post made for 35 degrees of physical motion clearance and a custom counterweight designed to compensate for the PR-1475 payload.
(3) Salinity and temperature sensors: Two different CTD sensors (SBE 37-SI MicroCAT from Sea-Bird Instruments and XR-620 from RBR Corporation) were placed on the stern and the prow of the ship. The SBE 37-SI was placed at about 1 m below the sea surface on the stern of the ship. The XR-620 was deployed to obtain the profiles of salinity and temperature at sampling stations, and it was programmed for a sampling rate of 6 Hz. Since the sampling rate was higher compared to the deployment speed of XR-620, the measured salinity in the first effective record can be recognized as sea surface salinity.

Radiometric Calibration
The calibration equations (horizontal and vertical channels) are determined from two points: a hot and a cold targets. In this study, the sky was used as "cold target", and an ambient temperature blackbody target as the "hot target". The "hot target" consists of pyramidal open-cell blackbody foam in which two precision temperature sensors are embedded ( Figure 3). Since the "hot target" is seen as a near-blackbody target, the value of the emitted brightness temperature is equal to its physical temperature. The brightness temperature for "cold target" is nominally 3.5 K, provided that the field of view of the radiometer is away from the galactic center and the sky is clear. To accomplish this, according to the galactic noise map [34,35], the antenna was pointed to the north, and at an elevation angle approximately 65° above the horizon. The radiometer was mounted on the stabilized platform, which was about 6 m high above the sea surface and provided dual polarization measurements of brightness temperature at 5 fixed incidence angles (16 • , 26 • , 36 • , 45 • and 60 • ), and the circular-shaped field of view on the sea surface has the radius in the range of 3 m to 6 m according to different incidence angles (for example, the radius of field of view is about 3 m with incidence angle at 16 • ). During the cruises, some stations were measured in daytime, which were hampered by sunlight. In order to avoid radiation coming from the Sun, the antenna of the radiometer was pointed away (backwards) from the sun and the radiometer's field of view was in the shadow of the ship. At few stations, the measurements were affected by the Sun, which was removed easily because the intensity of radiation emitted by the Sun was high (~10 K or higher) and the reflectivity of the sea surface was large (in the order of 0.7). Although the effect of radio-frequency interference (RFI) is another possible source of error for the measured brightness temperature, the shipborne radiometer was deployed about 6 m above the sea surface (unlike the airborne radiometer), which could avoid potential sources of RFI. Moreover, the radiometer was pointed towards the east direction at most of the time to avoid the RFI effect from the land. In addition, it was observed that, even under the wind speeds up to 11 m/s, the foam coverage in the radiometer's field of view was less than 5%, and thus, the foam effect was neglected in this study.

Radiometric Calibration
The calibration equations (horizontal and vertical channels) are determined from two points: a hot and a cold targets. In this study, the sky was used as "cold target", and an ambient temperature blackbody target as the "hot target". The "hot target" consists of pyramidal open-cell blackbody foam in which two precision temperature sensors are embedded ( Figure 3). Since the "hot target" is seen as a near-blackbody target, the value of the emitted brightness temperature is equal to its physical temperature. The brightness temperature for "cold target" is nominally 3.5 K, provided that the field of view of the radiometer is away from the galactic center and the sky is clear. To accomplish this, according to the galactic noise map [34,35], the antenna was pointed to the north, and at an elevation angle approximately 65 • above the horizon.

Radiometric Calibration
The calibration equations (horizontal and vertical channels) are determined from two points: a hot and a cold targets. In this study, the sky was used as "cold target", and an ambient temperature blackbody target as the "hot target". The "hot target" consists of pyramidal open-cell blackbody foam in which two precision temperature sensors are embedded ( Figure 3). Since the "hot target" is seen as a near-blackbody target, the value of the emitted brightness temperature is equal to its physical temperature. The brightness temperature for "cold target" is nominally 3.5 K, provided that the field of view of the radiometer is away from the galactic center and the sky is clear. To accomplish this, according to the galactic noise map [34,35], the antenna was pointed to the north, and at an elevation angle approximately 65° above the horizon. Then the radiometer measurements were performed at different observation angles from 16° to 60° ( 16 , 26 , 36 , 45 and 60 i θ =      ). The total number of the measurements was 101 for each polarization with five different incidence angles. The data possibly influenced by the RFI effect were first sorted out, and the measurements corresponding to the incidence angles outside the range of 1° larger or smaller than the incidence angles were also removed. For example, the measured brightness Then the radiometer measurements were performed at different observation angles from 16 • to 60 • (θ i = 16 • , 26 • , 36 • , 45 • and 60 • ). The total number of the measurements was 101 for each polarization with five different incidence angles. The data possibly influenced by the RFI effect were first sorted out, and the measurements corresponding to the incidence angles outside the range of 1 • larger or smaller than the incidence angles were also removed. For example, the measured brightness temperature exceeding the range of 35 • to 37 • for a 36 • incidence angle would be rejected. After data quality control by the RFI effect and incidence angle abnormal, 18 measurements were rejected for each polarization locating at six stations, and the remain 83 high quality measurements were used for further analysis (Table 1).

Correction for Sea Surface Reflection of Sky Radiation
The reflection of downward radiation (including atmospheric emission and extraterrestrial radiation) is a significant contribution to the radiometer measured brightness temperature. The extraterrestrial radiation is about 2.725 K caused by the cosmic microwave background radiation Remote Sens. 2019, 11, 2486 6 of 18 related to the Big Bang nucleosynthesis and continuum radiation [36]. Although the impact of galactic radiation (hydrogen emission and continuum radiation) could be strong when these radiation sources are limited to the specular reflection into the main antenna beam, these contaminated observations could easily be discarded because reflectivity of the sea surface is large (in the order of 0.7) and the galactic signal scattered from the sea surface is near specular reflection [10,17,37]. Thus, the steady and continuous downward radiation could be due to the cosmic background radiation plus the emission from the atmosphere. The sea surface reflected sky radiation could be estimated by the L-band radiative transfer model developed in our previous work [38]. In this model (Figure 4.), the two-scale model was used to approximate the sea surface reflected radiation towards the radiometer at different observation angles, and the downwelling atmospheric emission and radiation attenuation was obtained by using the millimeter propagation model developed by Liebe [39,40]. In addition, the profiles of the atmosphere properties were obtained from the National Center for Environmental Prediction (NCEP) FNL data, in which the whole atmosphere was divided into 26 layers, and for every layer, the gas extinction and emission were calculated with several input data including the air temperature, pressure, and relative humidity.

Correction for Sea Surface Reflection of Sky Radiation
The reflection of downward radiation (including atmospheric emission and extraterrestrial radiation) is a significant contribution to the radiometer measured brightness temperature. The extraterrestrial radiation is about 2.725 K caused by the cosmic microwave background radiation related to the Big Bang nucleosynthesis and continuum radiation [36]. Although the impact of galactic radiation (hydrogen emission and continuum radiation) could be strong when these radiation sources are limited to the specular reflection into the main antenna beam, these contaminated observations could easily be discarded because reflectivity of the sea surface is large (in the order of 0.7) and the galactic signal scattered from the sea surface is near specular reflection [10,17,37]. Thus, the steady and continuous downward radiation could be due to the cosmic background radiation plus the emission from the atmosphere. The sea surface reflected sky radiation could be estimated by the L-band radiative transfer model developed in our previous work [38]. In this model (Figure 4.), the two-scale model was used to approximate the sea surface reflected radiation towards the radiometer at different observation angles, and the downwelling atmospheric emission and radiation attenuation was obtained by using the millimeter propagation model developed by Liebe [39,40]. In addition, the profiles of the atmosphere properties were obtained from the National Center for Environmental Prediction (NCEP) FNL data, in which the whole atmosphere was divided into 26 layers, and for every layer, the gas extinction and emission were calculated with several input data including the air temperature, pressure, and relative humidity.

Auxiliary Data Preparation
The high accuracy of the input auxiliary data such as sea surface temperature, wind speed and significant wave height is required as these input data can significantly affect the accuracy of SSS retrieval. It is also necessary to obtain the brightness temperature data and auxiliary data simultaneously. In this study, together with the radiometer measurements, most of the auxiliary data were obtained by in-situ measurements with a high accuracy. The descriptions of the data processing for the auxiliary data are given below: (1) The SSS and SST were measured by the SeaBird CTD (SBE 37-SI MicroCAT) and RBR XR-620, respectively. These two measurements were compared to ensure that the surface water properties were measured instead of the subsurface water properties. When the two measurements were not the consistent (mostly the differences less than 0.2 • C and 0.1 psu), an average value was then used.
(2) The wind speed and direction were measured by a wind monitor (RM Young anemometer, Model 05106) with the accuracy of 0.3 m/s and 3 • . The anemometer was installed at a height of 8 m above the sea surface, thus the wind speed at 10 m above the sea surface was used as an input for the Remote Sens. 2019, 11, 2486 7 of 18 sea surface emission model (i.e., U 10 ), as approximated by Equation (1). The wind speed given at any elevation can be calculated from the friction velocity U * by: where U * is related to Z 0 by: and the azimuth angle is defined as φ w = 0 with the wind blowing towards the radiometer.
(3) Significant wave height (SWH) was obtained from the ERA5 hourly reanalysis data. Although the remote sensed SWH could achieve better accuracy than the numerical model results, only few stations could be matched to the remote sensed SWH data due to the lack of satellite SWH data in coastal regions. Thus, the SWH data from ERA5 hourly reanalysis data were used instead of the remote sensed SWH products. In fact, the ERA5 is the fifth generation ECMWF atmospheric and oceanic reanalysis production with the improved temporal-spatial resolution of 0.5 degrees in space and 1 h in time.

Sea Surface Emissivity Models
(1) Small perturbance/small scale approximation model (SPM/SSA): In the SPM/SSA model, the Stokes vector of the emission from the sea surface is expressed as an integral over the surface directional spectrum [21] as: where R vv are the horizontally and vertically polarized Fresnel reflection coefficients, respectively; W(k, φ) is the surface directional spectrum; g h and g v are the weighting functions in the formulation for correcting the reflectivity of a flat sea surface induced by small-scale waves. The weighting functions in Equation (3) are based on some theoretical understanding [18] and may require modifications based on measurements. For this, the incoherent Bragg scatter and second order correction terms were combined; the developed coordinates in [18] were modified by [21] and induced the modification of scattering coefficients. In this study, the surface directional spectrum used in the SPM/SSA model was developed by Durden and Vesecky [41] (hereafter DV spectrum). For well-developed ocean waves, the sea surface wave spectrum can be obtained as: where S(k) is the omnidirectional part of the surface wave spectrum and Φ(k, φ) is the angular portion of the spectrum. Following Durden and Vesecky, the omnidirectional part of the surface spectrum takes the form: where g * = 9.81 + 7.25 × 10 −5 k 2 ; k is the wavenumber and U * is the friction velocity, which is defined in Equation (2) from 0.004 to 0.008 and b 0 is selected in such a way that S(k) is continuous at k = 2 rad·m −1 for a given a 0 . The angular part is represented by: where s is a constant set as 1.5 × 10 −4 , and c is given by: with R = 0.003 + 0.00192U 12.5 0.00316U 12.5 .
Following [23], some misprints of coefficients in the DV spectrum described in [41] is corrected.
(2) Two-scale model: The two-scale model assumes that small-scale waves are rippled and superimposed on to the large-scale waves. In this model, a cutoff wavenumber k c to determine the large-scale waves and small-scale waves is given as: As described in the SPM/SSA model, the DV spectrum was used in this study. For small-scale waves, the Stokes vector of the emission from a local surface patch is exactly the same as Equation (3), and only the cutoff wavenumber is applied in the second terms of Equation (3). For large-scale waves, the integration is used to obtain the Stokes vector as: where R is the rotation matrix, which rotates the coordinate from every local patch to the global coordinate, and the explicit expression of R is given in [23]; h is the hydrodynamic modulation (HM), which induces upwind/downwind asymmetry, and is also given in [23] and P is the probability density function of large-scale waves, which can be expressed as the Gaussian function as [42]: where σ u and σ c are variances of wave slope for upwind and crosswind directions, respectively. The variances can be calculated by a double integral over the wave and azimuth angle as: (3) Empirical model: Unlike the TSM and SPM/SSA models, the empirical model decouples the effects of SST, SSS and sea surface roughness by expanding the sea surface emission with the Taylor series function and only the first order term of the series is used to describe the sea surface roughness effect. In that case, Remote Sens. 2019, 11, 2486 9 of 18 the observed brightness temperatures at the sea surface were firstly divided into two parts due to the perfectly flat sea surface and roughness sea surface conditions as: where where R p is the Fresnel reflection coefficient, subscript "p" stands for the polarization. It is widely accepted that ∆T B,rough is a function of WS. Camps [13] also showed that the ∆T B,rough is determined by both WS and SWH. In addition, the incidence angle is another parameter that is to be considered. Thus, the emission vector is expressed as: where ∆T h (θ i , U 10 , SWH) and ∆T v (θ i , U 10 , SWH) represent the emitted brightness temperature variations due to the sea surface roughness, and subscripts "h" and "v" denote horizontal and vertical polarizations, respectively. There are two similar methods to estimate sea surface roughness terms [25,43], and the only difference between these two methods is that the significant wave height (SWH) is introduced rather than using wind speeds (WS) directly to obtain ∆T(θ i , U 10 ) [44]. The sea surface roughness terms developed by Camps (hereafter Emp1) are given by: where θ i is the incidence angle and U 10 is the wind speed at 10 m above the sea surface. The other one (hereafter Emp2) is: These two empirical models are both based on the wind and salinity experiment (WISE), which were performed around the Casablanca oil rig platform, and the validity of these models in East China Sea waters remains to be explored.

Validation of the Sea Surface Emission Models in the ECS
It is known that the prediction discrepancies of the different sea surface emission models (especially, TSM model and empirical model of [25]) are quite significant for SSS retrieval. For the TSM model, the choice of the cutoff wavenumber and the factor a0 can affect the prediction of brightness temperature significantly. The a0 value is always assumed to vary from 0.004 to 0.008, and in some studies, the DV spectrum is multiplied with a factor of 2 [16,32,44]. In this paper, the TSM model with four typical combinations of factors was used to predict the brightness temperature and compare with the in-situ measured data, and the choices of the factor and a0 values are shown in Table 2. The influence of different factor and a0 values on the predicted brightness temperature is shown in Figure 5, in which the SST and SSS were fixed as 15 • C and 35 psu, and the cutoff wavenumber was set as k c = k 0 /3. Clearly, the predicted brightness temperature obtained from the model M1 (DV spectrum multiplied with the factor 2, hereafter DV2) is higher than those from the models M2, M3 and M4, especially for the horizontal polarization component. Comparing to vertical polarization, the horizontal polarization is sensitive to the wave spectrum, and the differences among the models M1, M2, M3 and M4 can reach up to 2 K at nadir, and increase to 10 K with an incidence angle at 60 • , which is too large for SSS retrieval. The differences between the models M1, M2, M3 and M4 are less than 2 K for vertical polarization, which is reasonable because the vertical polarized brightness temperature is less influenced by the sea surface roughness effect. Although the differences are small for vertical polarization, it is still too large for SSS retrieval. Thus, these factors need to be determined by fitting the model predictions to the in-situ measurements. On the other hand, due to the limited brightness temperature measurements in coastal waters of China, the regional bias of the empirical model in the ECS is still unclear and needs to be examined. wavenumber was set as 0 /3 c k k = . Clearly, the predicted brightness temperature obtained from the model M1 (DV spectrum multiplied with the factor 2, hereafter DV2) is higher than those from the models M2, M3 and M4, especially for the horizontal polarization component. Comparing to vertical polarization, the horizontal polarization is sensitive to the wave spectrum, and the differences among the models M1, M2, M3 and M4 can reach up to 2 K at nadir, and increase to 10 K with an incidence angle at 60°, which is too large for SSS retrieval. The differences between the models M1, M2, M3 and M4 are less than 2 K for vertical polarization, which is reasonable because the vertical polarized brightness temperature is less influenced by the sea surface roughness effect. Although the differences are small for vertical polarization, it is still too large for SSS retrieval. Thus, these factors need to be determined by fitting the model predictions to the in-situ measurements. On the other hand, due to the limited brightness temperature measurements in coastal waters of China, the regional bias of the empirical model in the ECS is still unclear and needs to be examined.    Table 3 presents the results for different sea surface emission models obtained by comparing the model-predicted and in situ brightness temperatures. The M4 model had the largest root-mean-square difference (rms; up to 3 K) for horizontally polarized brightness temperature, which could lead to 3 psu or a higher bias in SSS retrieval. Although the wave spectrum was multiplied with a factor of 2 to fit with the measured brightness temperature at higher frequencies [18], as done previously [16,32], the DV2 spectrum would overestimate the brightness temperature values. The rms values were the lowest for horizontally polarized brightness temperature with the M2 model and for vertically polarized brightness temperature with the M3 model. According to the mean biases, the M2 model slightly underestimated the brightness temperature, whereas the M3 model overestimated the brightness temperature. The underestimation or overestimation of the brightness temperature by these models might be attributed to the choice of cutoff wavenumber. Thus, the cutoff wavenumber needs to be considered in the model (as discussed in a later section). Furthermore, the empirical models developed based on the WISE campaigns show the regional biases in the ECS, and the coefficients in these models also need to be optimized. The choice of k c is rather arbitrary in previous studies; for example, k c = k 0 /4, k c = k 0 /5 and k c = k 0 /10 according to [32], [45] and [24], respectively. The kc values were given in the range of k 0 /10 to k 0 /3. The reason is that, as shown in [32], for the range of k 0 /10 to k 0 /3, the effect of k c on the predicted brightness temperature is small or even negligible. In addition, when using the DV2 spectrum, the difference caused by the cutoff wavenumber is less than 0.4 K with the incidence angle at 30 • and 40 • . Thus, k c = k 0 /3, k c = k 0 /5 and k c = k 0 /10 were tested in this study.
As shown in Figure 6, the horizontally polarized brightness temperatures had higher correlations with a0 and cutoff wavenumbers than the vertically polarized brightness temperatures. Moreover, the brightness temperature was underestimated when a0 was set as 0.004, especially for horizontal polarization, and when a0 was 0.006, or a0 = 0.008 together with k c = k 0 /3, the errors were relatively small compared to other combinations of a0 and kc values. The performances of the TSM model with different a0 and kc values are shown in Table 4, by comparison between the predicted and in situ brightness temperatures. The results show that a0 = 0.006 with k c = k 0 /10 fit the in situ values best with the rms values of 0.39 K and 0.37 K for horizontally polarized and vertically polarized brightness temperatures, respectively. In addition, the biases and standard deviations were also small when a0 = 0.006 and k c = k 0 /10. Although the rms values were similar for different a0 and k c values, the biases for horizontally and vertically polarized brightness temperatures had a significant influence on the SSS retrieval using these models.

Adjustment of the Empirical Model for the ECS
Following Equations (17)- (19), for every incidence angle, the sensitivity of ∆T B,rough to WS (or SWH) is represented by the slope of linear regression of ∆T B,rough to WS (or SWH). Generally, ∆T B,rough can be expressed as: where S i (θ i ) is the sensitivity of ∆T B,rough to WS (or SWH), i = 1, 2 and 3 and c is the offset. In this study, although all the in situ data were used, the number of measured brightness temperature was still small for different incidence angles. Thus, the wind speed sensitivity was obtained based on the previous studies. It should be noted that the fitted sensitivity of vertically polarized brightness temperature to wind speed from different campaigns was quite different. For example, at 20 • and 30 • incidence angles, the vertical polarization sensitivity to wind speed, as derived by the WISE, was half as it was obtained from the PALS and Hollinger's data. It decreased to −0.03 with the incidence angle at 50 • , whereas it was 0.22 from the PALS dataset. This could be attributed to the discrepancies in the dielectric constant model used to compute the flat sea surface brightness temperature (Equation (16)). The dielectric constant model used in this study was from Klein and swift [46], which is the same used in [13]. Thus, the wind speed sensitivity at different incidence angles was exactly the same as [13]. A linear regression method was used to obtain the coefficients in Equations (20) and (21). The relationships of ∆T B,rough with wind speeds and incidence angles are shown in Figure 7. Overall, the ∆T h,rough shows a good correlation with sensitivity and wind speed, and the vertical polarization variations show a lower sensitivity to wind speed compared to the horizontal polarization, which were consistent with the previous studies [10]. This might be attributed to the fact that the vertically polarized brightness temperature was less influenced by the WS and roughness. It should be noted that, for both polarizations, the brightness temperature variations did not fit well with incidence angle at 60 • (green points), this was particularly more evident for horizontal polarization. When the data at a 60 • incidence angle were excluded, the correlation coefficients between ∆T h,rough and wind speed (multiplying sensitivity) increased from 0.78 to 0.88. Although the brightness temperature measured at high incidence angles were less used for SSS retrieval, further study is still needed to examine the variation of sea surface roughness induced brightness temperature at low incidence angles and high incidence angles separately. From this analysis, the optimized ∆T h,rough and ∆T v,rough models for the ECS are given by: where S is the sensitivity of different parameters (i.e., WS and SWH), subscripts "h" and "v" represent horizontal and vertical polarizations, respectively. It should be noted that the sensitivity of horizontally polarized brightness temperature to SWH obtained from [13] is notably different from that of [25] (θ/142 comparing to 1 − θ/50). For example, when the incidence angle at 60 • , the sensitivity was 0.4225 and −0.2 from these two models. In order to fit well with the measurements, these two sensitivities were both tested in this study, and the results showed that the sensitivity from [13] fitted well with the measurements, and finally θ/142 was used in Equation (22).

Comparison of the Retrieved Sea Surface Salinity using Different Models
For SSS retrieval, the iterative Levenberg and Marquard method was used, which is fully described in the literature [47,48]. The L-M method was performed by minimizing the cost function defined as:  23)). It is known that using the first Stokes parameter can achieve a higher sensitivity to SSS, however, when the forward model overestimates or underestimates the sea surface emission, the error of retrieved SSS would be larger. In this paper, three models were used to retrieve SSS: 1) TSM with a0 = 0.006 and 0 /10 c k k = , 2) modified EMP1 and 3) modified EMP2 (the results are reported in Table 5). It can be seen that the RMS of the retrieved SSS was about 0.4 psu for all three models, but the modified TSM model had the highest accuracy. The lower performance of the EMP models compared to the TSM model might be due to the fewer in-situ measurements used to improve the empirical model. As mentioned in Section 4.2, the , B rough T Δ needs to be calculated separately at low and high incidence angles in order to achieve the higher accuracy. The spatial distributions of retrieval errors are shown in Figure 8. Clearly, the accuracy of retrieved SSS using the modified TSM model was higher than the two empirical models, especially

Comparison of the Retrieved Sea Surface Salinity using Different Models
For SSS retrieval, the iterative Levenberg and Marquard method was used, which is fully described in the literature [47,48]. The L-M method was performed by minimizing the cost function defined as: where T mea Bi and T mod  23)). It is known that using the first Stokes parameter can achieve a higher sensitivity to SSS, however, when the forward model overestimates or underestimates the sea surface emission, the error of retrieved SSS would be larger. In this paper, three models were used to retrieve SSS: (1) TSM with a0 = 0.006 and k c = k 0 /10, (2) modified EMP1 and (3) modified EMP2 (the results are reported in Table 5). It can be seen that the RMS of the retrieved SSS was about 0.4 psu for all three models, but the modified TSM model had the highest accuracy. The lower performance of the EMP models compared to the TSM model might be due to the fewer in-situ measurements used to improve the empirical model. As mentioned in Section 4.2, the ∆T B,rough needs to be calculated separately at low and high incidence angles in order to achieve the higher accuracy. The spatial distributions of retrieval errors are shown in Figure 8. Clearly, the accuracy of retrieved SSS using the modified TSM model was higher than the two empirical models, especially at stations in the vicinity of the coast. For these stations, the shallow water was one of the factors that affected the sea surface roughness contributions, since shallow water waves were significantly affected by the seabed topography in these coastal areas. Moreover, the error of the SWH was another source, which hampered the accuracy of model prediction and SSS retrieval. In this study, the spatial resolution of the SWH data obtained from the ECMWF ERA-5 reanalysis dataset was 0.5 • × 0.5 • , which was too coarse to match the field of view of radiometer. Thus, in order to derive an empirical model with high accuracy, the SWH should be measured simultaneously with the radiometer measurement. at stations in the vicinity of the coast. For these stations, the shallow water was one of the factors that affected the sea surface roughness contributions, since shallow water waves were significantly affected by the seabed topography in these coastal areas. Moreover, the error of the SWH was another source, which hampered the accuracy of model prediction and SSS retrieval. In this study, the spatial resolution of the SWH data obtained from the ECMWF ERA-5 reanalysis dataset was 0.5° × 0.5°, which was too coarse to match the field of view of radiometer. Thus, in order to derive an empirical model with high accuracy, the SWH should be measured simultaneously with the radiometer measurement.

Conclusions
Satellite remote sensing of SSS is based on the variation of sea surface emission caused by the change of salinity, but this method generally assumes that the sea surface is perfectly flat. However, in the real world, the sea surface is always rough, and sea surface emission is significantly affected by the sea surface roughness. Thus, one of the most important issues in remote sensing of SSS is to remove the contribution of sea surface roughness on sea surface emission. In this study, two cruise experiments were performed in the ECS to assess and improve the sea surface microwave emission models in the ECS. To obtain the brightness temperature of the sea surface, the contributions of surface reflected downwelling sky radiations including cosmic microwave background radiation and atmospheric emission were first removed by using a developed radiative transfer model. Based on

Conclusions
Satellite remote sensing of SSS is based on the variation of sea surface emission caused by the change of salinity, but this method generally assumes that the sea surface is perfectly flat. However, in the real world, the sea surface is always rough, and sea surface emission is significantly affected by the sea surface roughness. Thus, one of the most important issues in remote sensing of SSS is to remove the contribution of sea surface roughness on sea surface emission. In this study, two cruise experiments were performed in the ECS to assess and improve the sea surface microwave emission models in the ECS. To obtain the brightness temperature of the sea surface, the contributions of surface reflected downwelling sky radiations including cosmic microwave background radiation and atmospheric emission were first removed by using a developed radiative transfer model. Based on the in situ brightness temperature and auxiliary data (including SSS, SST, wind speed and direction and SWH), the performances of different emission models including TSM, SPM and SSA and empirical models were examined. The results showed that the accuracy of these models varied from 0.6 K to 3 K. The accuracy of TSM was influenced by the choices of parameters in the surface wave spectrum. The error of TSM could reach up to 0.6 K for a0 = 0.006 and decrease to 3 K with the wave spectrum multiplied with a factor 2, and the accuracy of SSA/SPM and empirical models were quite similar, which were 1.1 K and 0.8 K for horizontally and vertically polarized brightness temperatures, respectively.
The validation in this study showed that the existing models performed poorly in the ECS and required to be optimized for regional applications. For the TSM, it is well known that the accuracy is affected by the choice of the parameters a0 and cutoff wavenumber k c . It was found that with a0 = 0.006 and k c = k 0 /10, the brightness temperature could be fairly accurately predicted with an error less than 0.4 K for both polarizations. Furthermore, the empirical models were improved to be applicable for the ECS by optimizing the model parameters, and the brightness temperature variation as a function of both wind speed and SWH was also established. The results showed that the sea surface roughness induced brightness temperature variations were well fitted by wind speed at low incidence angles, especially for horizontal polarization, with the correlation coefficient of 0.88. At large incidence angles, the brightness temperature variations also showed a good correlation with sensitivity and wind speed. However, when all data were used, the accuracy of model prediction became lower, which needs further study to improve the empirical model performance. Based on SSS retrievals with the improved models using shipborne brightness temperature data, it was found that that the modified TSM retrieved SSS with an RMS and standard deviation around 0.4 psu. Compared to the modified TSM model, the accuracy of retrieved SSS using the EMP model was relatively low, which might be caused by the influence of seabed topography on the surface wave spectrum and the coarse resolution of the input SWH data.