The Potential and Challenges of Using Soil Moisture Active Passive ( SMAP ) Sea Surface Salinity to Monitor Arctic Ocean Freshwater Changes

Sea surface salinity (SSS) links various components of the Arctic freshwater system. SSS responds to freshwater inputs from river discharge, sea ice change, precipitation and evaporation, and oceanic transport through the open straits of the Pacific and Atlantic oceans. However, in situ SSS data in the Arctic Ocean are very sparse and insufficient to depict the large-scale variability to address the critical question of how climate variability and change affect the Arctic Ocean freshwater. The L-band microwave radiometer on board the NASA Soil Moisture Active Passive (SMAP) mission has been providing SSS measurements since April 2015, at approximately 60 km resolution with Arctic Ocean coverage in 1–2 days. With improved land/ice correction, the SMAP SSS algorithm that was developed at the Jet Propulsion Laboratory (JPL) is able to retrieve SSS in ice-free regions 35 km of the coast. SMAP observes a large-scale contrast in salinity between the Atlantic and Pacific sides of the Arctic Ocean, while retrievals within the Arctic Circle vary over time, depending on the sea ice coverage and river runoff. We assess the accuracy of SMAP SSS through comparative analysis with in situ salinity data collected by Argo floats, ships, gliders, and in field campaigns. Results derived from nearly 20,000 pairs of SMAP and in situ data North of 50◦N collocated within a 12.5-km radius and daily time window indicate a Root Mean Square Difference (RMSD) less than ~1 psu with a correlation coefficient of 0.82 and a near unity regression slope over the entire range of salinity. In contrast, the Hybrid Coordinate Ocean Model (HYCOM) has a smaller RMSD with Argo. However, there are clear systematic biases in the HYCOM for salinity in the range of 25–30 psu, leading to a regression slope of about 0.5. In the region North of 65◦N, the number of collocated samples drops more than 70%, resulting in an RMSD of about 1.2 psu. SMAP SSS in the Kara Sea shows a consistent response to discharge anomalies from the Ob’ and Yenisei rivers between 2015 and 2016, providing an assessment of runoff impact in a region where no in situ salinity data are available for validation. The Kara Sea SSS anomaly observed by SMAP is missing in the HYCOM SSS, which assimilates climatological runoffs without interannual changes. We explored the feasibility of using SMAP SSS to monitor the sea surface salinity variability at the major Arctic Ocean gateways. Results show that although the SMAP SSS is limited to about 1 psu accuracy, many large salinity changes are observable. This may lead to the potential application of satellite SSS in the Arctic monitoring system as a proxy of the upper ocean layer freshwater exchanges with subarctic oceans.

In situ salinity observations in the Arctic Ocean have been historically very sparse.Spaceborne observation of sea surface salinity (SSS) is anticipated to advance our monitoring capability for the Arctic Ocean.SSS is a good tracer of the freshwater movement in the upper Arctic Ocean.The synoptic large-scale spatial variability of the SSS field that was obtained from the vantage point of space will not only monitor freshwater changes due to regional processes in the Arctic Ocean, but it will also provide a means to improve our understanding of the interaction between the Arctic and subarctic oceans.
Global SSS observations have been available from three satellite missions.The first one is the European Space Agency (ESA) Soil Moisture and Ocean Salinity (SMOS) mission launched on 2 November 2009 [30,31], carrying a passive L-band microwave system.The second one is the Aquarius/SAC-D, a joint venture of National Aeronautics and Space Administration (NASA) and Argentina's Comision Nacional de Actividades Espaciales (CONAE) [32][33][34], which was in operation from 10 June 2011 to 7 June 2015.The NASA Soil Moisture Active Passive (SMAP) observatory, which was launched in 31 January 2015, designed to measure the soil moisture over land [35], is also observing global ocean with an L-band radiometer and is currently being used to estimate SSS [36].
Satellite SSS data have enabled many new scientific results for tropical and subtropical oceans.However, its applications in the polar and sub-polar oceans did not progress as well due to two reasons [37,38].First, satellite SSS retrieval in high latitude seawater is complicated by several factors.L-band instrument sensitivity to SSS is greatly reduced in cold seawater [39,40].Leakage of emissivity from sea ice into the satellite antenna's main lobe or through the antenna side-lobes can contaminate salinity signals if undetected [41].The larger uncertainty in the ocean surface roughness correction in the rough seas frequently encountered under the conditions of high winds and the strong waves of polar oceans influences the accuracy of the SSS retrievals [42,43].The relatively large uncertainties of the ancillary satellite sea surface temperature (SST) data in the Arctic Ocean used in the SSS retrievals also affect the SSS accuracy [44].Second, the paucity of in situ measurements hinders the validation of satellite SSS [45][46][47].The validation of satellite SSS using in situ data can also be complicated by the difference in spatial scale between the satellite measurements over satellite footprints, and the point-wise in situ measurements in regions with strong sub-footprint variability [48][49][50].As a result, quite large gaps near land and ice are masked out (~110 km for Aquarius and ~50 km for previously released SMAP SSS) in order to avoid contaminated retrieval, which could be misinterpreted or interfere with global validation results.
The recent release of SMAP SSS (V4.0) produced at the Jet Propulsion Laboratory (JPL) [51] was based on a more advanced correction strategy, aiming to deliver the valuable information in the critical regions near land and sea ice as much as possible.One of the main improvements is land contamination correction on measured brightness temperature (TB) before SSS retrievals.Land contamination is the leakage of energy from land surface into the radiometer receiver through side lobes or partially through the main lobe.In land-contaminated TB, non-oceanic signatures mixed with oceanic signatures may result in differences between two categories that exceed 100 kelvin (K).If undetected or uncorrected, the non-oceanic signature would cause large bias in the satellite salinity values.The JPL SMAP SSS algorithm has implemented the land contamination correction that is based on two look-up-tables (LUT), both being derived from SMAP data itself.One LUT is the beam-integrated land fraction as a function of footprint position and orientation, while the other is a monthly climatology of land TB extending out to 1000 km from coast, which represents effective contributions of TB from side-lobes over land.Combining the two LUTs, the land contamination contribution for each SMAP footprint is estimated and removed from TB.With improved land correction, the JPL algorithm retrieves SSS within 35 km from the land wherever the sea ice concentration values are less than 3% (based on NCEP SIC analysis [52]).
This study investigates the utility and accuracy of SMAP SSS retrieved in the Arctic Ocean via two complementary approaches: (1) comparing with collocated in situ salinity data, and (2) analyzing the co-variability of SMAP SSS with changes in sea ice and river discharges.We also examined the seasonal and interannual variability of SMAP SSS over major Arctic Ocean gateways to explore the potential of satellite SSS as a proxy to infer the freshwater exchange between the Arctic and sub-oceans.The datasets that were used for the analysis are described in Section 2. The results are presented in Section 3. Section 4 provides further discussion, and Section 5 summarizes the findings.

SMAP SSS
SMAP SSS are produced at Level 2 (L2) and Level 3 (L3).L2 data represents the spatial average of the instantaneous measurements in the swath footprint, delivered for each of the salinity-wind-cell (SWC) along the satellite swath with 1624 × 76 cells along/cross track per satellite revolution.L2 data covers the global ocean in eight days with a spatial resolution of ~60 km, posted on a 25 km grid.There are two Level 3 (L3) products, the first with monthly temporal averaging and the second with eight-day averaging, both on a global 0.25 • × 0.25 • grid.The eight-day product is created daily by averaging eight days of L2 data centered at noon UTC (Coordinated Universal Time) of the day with a search radius of 45 km and a Gaussian weighting half-power distance of 30 km.Also, provided along with L2 and L3 products are the surface salinity output from the Hybrid Coordinate Ocean Model (HYCOM) [53] interpolated to SMAP measurement time and location.
Currently, there are two SMAP SSS products.One is the JPL product [36,51] that is mentioned in Section 1, and the other is based on algorithm independently developed at the Remote Sensing System (RSS) [54].Both datasets are available at NASA Physical Oceanography Distributed Active Archive Center (PO.DAAC) (http://podaac.jpl.nasa.gov).The RSS product leaves large gaps in the Arctic Ocean near the ice edge, and is therefore not suitable for this study.We will focus on analyzing the JPL SMAP SSS V4.0 product (hereafter named SMAP SSS).

In Situ Salinity Data
Individual Argo float profiles are obtained from USGODAE (The United States (US) Global Ocean Data Assimilation Experiment, http://www.usgodae.org/argo/argo.html).The Argo array consists of approximately 3700 floats in the global ocean measuring salinity and temperature [55,56].Argo data were collected and were made freely available by the International Argo Program and the national programs that contribute to it, as part of the Global Ocean Observing System (GOOS).
Salinity data from ships and gliders in the Arctic Ocean are obtained from the Copernicus Marine Environment Monitoring Service (CMEMS), which is an earth observing data center under the European Commission [57].We use the product of Arctic Ocean-In situ Near Real Time Observations [58], which contains data collected from the Arctic ROOS (Regional Ocean Observing System) members, and quality controlled using automated procedures that are available from the In Situ Thematic Assembly Centre (INS TAC).
We also used data sets that were collected in recent field campaigns.During the NASA Ocean Melting Greenland (OMG) mission [10,11] field campaigns that took place in September/October 2016 and October 2017, NASA's G-III research plane, flying along the Greenland coast, deployed Airborne expandable CTD (AXCTD) sensors to measure temperature and salinity profiles from the surface to about 1000 m depth [58].Quality controlled OMG data is made publicly available at (https://omg.jpl.nasa.gov/portal/).University of East Anglia, UK provided us with the salinity data from the seaglider missions MASSMO4, AlterEco and Oman579 (personal communication with Bastien Queste, PI) (http://ueaglider.uea.ac.uk).

Sea Ice Concentration
Sea ice concentration (SIC) data not only identifies ice-free areas where satellite SSS retrieval is possible, but also may influence the SSS values through freshwater exchanges that were associated with sea ice changes [59].We obtain the daily and monthly SIC data from the NASA DAAC at the National Snow and Ice Data Center (NSIDC).The SIC data set [NSIDC-0051] is generated from multi-frequency brightness temperature data and is designed to provide a consistent time series of SIC combining the coverage of several passive microwave instruments [60].The data are provided in the polar stereographic projection at a grid cell size of 25 × 25 km 2 .The SIC data are also available in the Near-Real-Time (NRT) with a daily update [61].The JPL SMAP SSS retrieval algorithm uses the high resolution (1/12 • ) SIC analyses from NCEP [52] to identify ice-free regions (ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/omb/prod/).The SMAP L1B TB footprints falling into bins with SIC < 3% are flagged as ice-free.

River Discharge
Many large rivers in the northern regions produce huge amounts of flow to the Arctic Ocean.River discharge is therefore a dominant component in the Arctic freshwater system.Land hydrology in the polar regions are closely related with climate change and permafrost conditions [62,63].The daily discharge data from six largest Arctic rivers are available from Arctic Great River Observatory (Arctic-GRO) (http://arcticgreatrivers.org), which is a component of the National Science Foundation's (NSF) Arctic Observing Network (AON).Arctic-GRO implemented identical sampling and analytical protocols that were developed by the PARTNERS Project (Pan-Arctic River Transport of Nutrients, Organic Matter, and Suspended Sediments), a coordinated international effort to collect and analyze a time-series of water samples [64,65].We use discharge data from Ob' and Yenisey River during May 2015 to April 2017 for this study [66].Other studies examined the flow regimes and changes for these two rivers [67,68].We focus on the discharge at the downstream stations for basin freshwater input to the Arctic Ocean.

Moorings at Arctic Gateways
The main Arctic gateways are well monitored by in situ mooring arrays [24,69,70].The mooring data over the Bering Strait [69,70] was obtained from the Applied Physics Laboratory (APL) at the University of Washington (http://psc.apl.washington.edu/HLD/Bstrait/bstrait.html).Although the measuring depth of moorings (~50 m) is largely different from the penetration depth (~1-2 cm) of L-band radiometer, and is therefore is not included in the in situ data for validation (Section 3.1), it should provide an indication of the dynamical range of the seasonal variation and interannual trend of SSS.Historical mooring data also exists at the Davis Strait and the Fram Strait.Unfortunately, we were not able to find up-to-date mooring data at either the Davis Strait or the Fram Strait.

Validation with In Situ Salinity
We aggregated in situ salinity data collected by Argo floats, ships, gliders, and in OMG field campaigns to create a gridded in situ data set on 0.  The three-year mean of SSS averaged over collocated grid points from April 2015 to March 2018 are illustrated in Figure 2. In situ data confirms the large-scale salinity feature that was observed in SMAP SSS, i.e., the much saltier seawater on the Atlantic side than on the Pacific side of the Arctic Ocean.The biases with respect to in situ salinity are less than 1 psu in a majority of the areas, which are masked out in white so areas with large biases stand out more clearly (Figure 2, right column).It can be seen that SMAP SSS shows positive biases in a few places near the coast, and negative biases along the east coast of Greenland and near Davis Strait (Figure 2d).On the other hand, HYCOM SSS shows mostly positive biases east of Greenland (Figure 2e).The contrast between SMAP and HYCOM east of Greenland is particularly interesting.It suggests that SMAP SSS retrieval may contain un-detected sea ice effect, while HYCOM is not able to catch the freshening signal due to ice melting where no Argo floats are available for HYCOM's operational data assimilation.
The statistical results of comparing in situ salinity with collocated SMAP SSS and HYCOM SSS are presented in Table 1 and Figure 3.All data are gridded on 0.25° grids daily, covering the period from April 2015 to March 2018.North of 50°N, the bias, standard deviation, and Root Mean Square Difference (RMSD) between 19,738 pairs of SMAP SSS and in situ are 0.442, 2.391 and 2.431 psu with correlation of 0.81.The statistics between HYCOM SSS and in situ is slightly better.The open circles in the scatter plot indicate outliers defined as those pairs where the absolute difference exceeds three times the standard deviation.The number of such defined outliers is less than 1% of the total collocated data set.The entire set of collocated data can be visually divided into two groups: one group of high SSS (>25 psu) and the other of low SSS (<10 psu).It is interesting to note that in the low SSS group, in situ salinity are distinct at three sub-groups that are centered around 2, 5, and 7 psu, with collocated SMAP SSS spread from 0 to 10 psu, while HYCOM SSS almost perfectly aligned with the in situ data.Again, this is expected because most of the in-situ data being used for the evaluation here are from Argo profiling floats, and HYCOM assimilates the Argo float data operationally.The three-year mean of SSS averaged over collocated grid points from April 2015 to March 2018 are illustrated in Figure 2. In situ data confirms the large-scale salinity feature that was observed in SMAP SSS, i.e., the much saltier seawater on the Atlantic side than on the Pacific side of the Arctic Ocean.The biases with respect to in situ salinity are less than 1 psu in a majority of the areas, which are masked out in white so areas with large biases stand out more clearly (Figure 2, right column).It can be seen that SMAP SSS shows positive biases in a few places near the coast, and negative biases along the east coast of Greenland and near Davis Strait (Figure 2d).On the other hand, HYCOM SSS shows mostly positive biases east of Greenland (Figure 2e).The contrast between SMAP and HYCOM east of Greenland is particularly interesting.It suggests that SMAP SSS retrieval may contain un-detected sea ice effect, while HYCOM is not able to catch the freshening signal due to ice melting where no Argo floats are available for HYCOM's operational data assimilation.
The statistical results of comparing in situ salinity with collocated SMAP SSS and HYCOM SSS are presented in Table 1 and Figure      Significant outliers of SMAP SSS are seen as being associated with the high SSS group in two branches, with positive and negative biases, respectively, where positive bias may indicate an over estimation of the surface roughness correction, and the negative bias may suggest un-detected ice contamination.This will be considered for the next version of SMAP geophysical model function to improve the SSS retrieval algorithm.Note that there is no negative bias in HYCOM for the high SSS group.However, there is similar positive bias in the HYCOM versus in situ.After excluding the outliers, the bias, standard deviation, and RMSD between 19,543 pairs of SMAP SSS and in situ are 0.385, 0.987, and 1.060 psu with correlation of 0.82.
North of the Arctic Circle, the number of collocated points drops more than 70%.Table 1 also provides of statistical results North of 65°N.The bias and RMSD between 5712 pairs of SMAP and in situ after removing outliers are 0.339 and 1.227 psu, respectively.
After excluding the outliers, we performed a linear regression of the SMAP (or HYCOM) vs. in situ data.The regression slope for SMAP is close to 0.9 if only the data for in situ greater than 25 psu are considered (Figure 3).If all data are included, then the regression slope is 1.06.Interestingly, Significant outliers of SMAP SSS are seen as being associated with the high SSS group in two branches, with positive and negative biases, respectively, where positive bias may indicate an over estimation of the surface roughness correction, and the negative bias may suggest un-detected ice contamination.This will be considered for the next version of SMAP geophysical model function to improve the SSS retrieval algorithm.Note that there is no negative bias in HYCOM for the high SSS group.However, there is similar positive bias in the HYCOM versus in situ.After excluding the outliers, the bias, standard deviation, and RMSD between 19,543 pairs of SMAP SSS and in situ are 0.385, 0.987, and 1.060 psu with correlation of 0.82.
North of the Arctic Circle, the number of collocated points drops more than 70%.Table 1 also provides of statistical results North of 65 • N. The bias and RMSD between 5712 pairs of SMAP and in situ after removing outliers are 0.339 and 1.227 psu, respectively.
After excluding the outliers, we performed a linear regression of the SMAP (or HYCOM) vs. in situ data.The regression slope for SMAP is close to 0.9 if only the data for in situ greater than 25 psu are considered (Figure 3).If all data are included, then the regression slope is 1.06.Interestingly, HYCOM SSS seems to have a systematic positive bias for in situ SSS in the range of 25 to 30 psu, resulting in a reduced regression slope of 0.467 for the high SSS group of data.The SMAP SSS conditioned on the in situ within this range has a smaller bias, although being more noisy.

The Arctic Ocean SSS and Sea Ice
Figure 4 illustrates the monthly mean SSS from 2015 to 2017 in August, when the Arctic Ocean ice has the minimum extent with maximum ice-free surface areas for satellite salinity retrieval.The large-scale salinity feature, i.e., saltier Atlantic and fresher Pacific is consistently observed in three consecutive years and it agrees with the known contrast between the two ocean basins.North of the Arctic Circle (~66 • N), however, SMAP SSS shows large interannual variations, not only in magnitudes of SSS retrieved, but also in areas with a valid retrieval.For example, North of the Bering Strait, in August 2015, SSS is retrieved in the large areas of the Chukchi Sea and the East Siberian seas; in August 2016 in the Beaufort Sea and the East Siberian Sea; and, in August 2017 areas with valid retrieval extended further north in all three marginal seas.This is due to the differing sea ice extent for the same months among three years, because SMAP SSS retrieval is only possible in open water.For the regions with SSS retrieval in all three years, e.g., in the Hudson Bay or Kara Sea, SMAP SSS shows large interannual differences.In contrast, HYCOM SSS shows a much small variation in the marginal seas within the Arctic Circle.It is known that HYCOM uses climatology river discharge and its SSS is relaxed to a seasonal climatology; both may suppress the magnitude of interannual variations of SSS in HYCOM [71].
SMAP SSS in marginal seas near the ice edge often show very low values.These fresh signatures should be examined and validated very carefully because it could be a mixture of real signal and ice contamination.It is known that in seawater near the ice edge, sea ice formation, and melt have significant contributions to the mixed layer salt budget of the ambient waters, with growing importance toward the ice edge [59].Sea ice concentration (SIC) averaged in the areas North of 65 • N changes dramatically through the seasonal cycle, from 20% in October or November to more than 70% in winter (January-May) (blue curve in Figure 5).It appears that the seasonal cycle of SSS averaged North of 65 • N closely follows that of SIC, with the lowest SSS being observed within one or two months of the minimum SIC.This covariance should be taken with a grain of salt though, since the low SSS in summer also includes the effect of river runoff, and the high SSS in winter largely comes from the Atlantic side when most marginal seas are frozen.However, it is evident that SMAP SSS (black curve in Figure 5) reveals large interannual variations that are consistent with the sea ice conditions.Note that while averaging of SIC was over the entire region including zero or non-zero SIC values, the averaging of SSS is obtained over ice-free areas, which itself changes with time.We found the ice-free area during the ice melting seasonal peak of 2017 is 5% and 10% more than the two previous seasons of 2015 and 2016 (red curves in Figure 5).Correspondingly, the Arctic Ocean averaged SMAP SSS is the lowest in August 2017 despite HYCOM SSS showing essentially no changes.To alienate the possibility that summer SSS minimum was dominated by a few extremely low SSS retrievals or river discharges, Figure 5 also shows the time series of SSS after excluding outliers, which is identified as where the absolute differences from the daily mean exceeding three times of the standard deviation in the averaging domain.After removing the outliers (about 1 to 5% of valid retrievals of the day), the domain averaged SSS increases by 1-2 psu in its summer low for both SMAP and HYCOM, indicating that most outliers are in the low end of SSS.It is important to note that with outliers removed, SMAP SSS still shows clear interannual variability with minimum of 2017 season more than 2 psu lower than that of 2016, in contrast to HYCOM SSS.The consistent seasonal variation of SSS and SIC suggests that SMAP SSS retrieval correctly characterized the effect of sea ice changes.However, this agreement is encouraging but insufficient to completely exclude the effect of un-detected ice-contamination impacting SMAP SSS.The appearance of frazil-pancake ice maybe increasing in the margins of the Arctic during autumn, as identified in a recent field campaign [16].This type of new sea ice forms in the presence of incoming surface waves and it commonly occurs in the wave-dominated marginal ice zone of the Southern Ocean.With the retreat of the Arctic summer ice edge, the increased fetch has led to the increase of surface waves in the summer and fall, particularly in the Beaufort-Chukchi seas.Frazil-pancake ice may present a challenge to identify, say with passive microwave, as the sea ice concentrations may be quite variable over short time scales, which may lead to contamination in the SSS fields.Frazil-pancake ice has long been identifiable in SAR imagery.One recent study [72] describes the validation of a method to estimate the thickness of this young ice type based on wave dispersion.Being able to identify this form of new ice in the Arctic along with other sea ice conditions (e.g., melt, thin ice thickness, eddy formation, presence of waves) and the extent within the marginal ice zone will improve the accurate detection of ice and reduce the uncertainty in the satellite SSS retrieval.

SSS and River Discharge
As indicated in Section 3.1, there are very few in situ salinity data in the marginal seas within the Arctic Circle; for example, in the Hudson Bay and Kara Sea, no in situ data during the three years of SMAP were identified.In this section, we explore an alternative validation by examining the change of SMAP SSS in response to independently measured river discharge data.Massive northern rivers transport huge quantities of water from the continents to the Arctic Ocean.The freshwater inputs that are associated with river runoff should be reflected in the SSS field, particularly near the river mouths and over the shallow shelf areas.Due to its proximity to the two major Arctic rivers (Ob' and Yenisey), the Kara Sea provides an ideal case to examine SSS response to river discharge.
Figure 6 illustrates the evolution of SMAP SSS over the Kara Sea for the warm season (May to October) of 2015 and 2016.No SSS were retrieved in the Kara Sea before June and after October because it was completely covered by sea ice.The first valid SSS retrievals appeared in June, but at However, this agreement is encouraging but insufficient to completely exclude the effect of un-detected ice-contamination impacting SMAP SSS.The appearance of frazil-pancake ice maybe increasing in the margins of the Arctic during autumn, as identified in a recent field campaign [16].This type of new sea ice forms in the presence of incoming surface waves and it commonly occurs in the wave-dominated marginal ice zone of the Southern Ocean.With the retreat of the Arctic summer ice edge, the increased fetch has led to the increase of surface waves in the summer and fall, particularly in the Beaufort-Chukchi seas.Frazil-pancake ice may present a challenge to identify, say with passive microwave, as the sea ice concentrations may be quite variable over short time scales, which may lead to contamination in the SSS fields.Frazil-pancake ice has long been identifiable in SAR imagery.One recent study [72] describes the validation of a method to estimate the thickness of this young ice type based on wave dispersion.Being able to identify this form of new ice in the Arctic along with other sea ice conditions (e.g., melt, thin ice thickness, eddy formation, presence of waves) and the extent within the marginal ice zone will improve the accurate detection of ice and reduce the uncertainty in the satellite SSS retrieval.

SSS and River Discharge
As indicated in Section 3.1, there are very few in situ salinity data in the marginal seas within the Arctic Circle; for example, in the Hudson Bay and Kara Sea, no in situ data during the three years of SMAP were identified.In this section, we explore an alternative validation by examining the change of SMAP SSS in response to independently measured river discharge data.Massive northern rivers transport huge quantities of water from the continents to the Arctic Ocean.The freshwater inputs that are associated with river runoff should be reflected in the SSS field, particularly near the river mouths and over the shallow shelf areas.Due to its proximity to the two major Arctic rivers (Ob' and Yenisey), the Kara Sea provides an ideal case to examine SSS response to river discharge.
Figure 6 illustrates the evolution of SMAP SSS over the Kara Sea for the warm season (May to October) of 2015 and 2016.No SSS were retrieved in the Kara Sea before June and after October because it was completely covered by sea ice.The first valid SSS retrievals appeared in June, but at different locations in 2015 and 2016 due to the different spatial coverage of sea ice.The differences between these two years became more dramatic throughout the season.In the summer of 2015, the freshwater patch first appeared in June, in the area east of the northern tip of the Novaya Zemlya archipelago.It grew and spread to cover almost half of the Kara Sea in August.In contrast, in the season of 2016, the freshwater signature was limited to areas near Ob' and Yenisey Bay.Because the whole region has become almost ice-free after July, the impact of new freshwater inputs from sea ice melt is likely to be minimal.Therefore, the dramatic freshening signature spreading through middle of the Kara Sea from July to September in 2015 and along the Siberia coast in 2016, must have originated from the other freshwater source-river discharge.Below, we show that the contrast in SMAP SSS clearly reflects the differences in river runoff in those two years.As a reference, Figure 7 illustrates the SSS evolving patterns in the Kara Sea from HYCOM SSS.The anomaly in HYCOM SSS is an order of magnitude smaller.As mentioned before, HYCOM is forced by climatological river discharges and its SSS is relaxed to a seasonal climatology [71].These two climatological forcing would suppress the magnitude of interannual variations of SSS in HYCOM for regions without in-situ data or where in-situ data are insufficient to constraint the model SSS, such as the Kara Sea.In the surface layer of Kara Sea ocean currents carry waters from the River Ob' and River Yenisey to the north and northeast.Based on the SMAP SSS observations in the Kara Sea (Figure 6), we find that river discharge adjacent to the Kara Sea must be much larger in 2015 than in 2016.The additional freshwater in 2015 was transported northward and further spread along the path, while the spreading of freshwater discharge in 2016 was limited along the coast.
Indeed, the daily discharge data from the Ob' River and the Yenisey River (Figure 8) show differences that are consistent with SMAP observation.In the first two months of the warm season (May and June), discharge from Ob' shows a similar magnitude for the two years, while Yenisey discharge peaked early and injected 93.4 km 3 more water in 2015 than 2016.In the following four months from July to October, Ob' became the main player, putting 90.9 km 3 more freshwater into the Kara Sea, while Yenisey added another extra 31.5 km 3 (Table 2).Combining the discharges from the Ob' and Yesiney during May to October together, the Kara Sea received more than 210 km 3 extra freshwater in 2015 relative to 2016.It might be useful to roughly estimate the effect of this extra amount of freshwater on the SSS anomaly.About 3 cm of freshwater are needed to dilute 1 m of seawater with a change of salinity by 1 psu.Assuming the extra ~210 km 3 freshwater spread over half of the Kara Sea (total surface area 926,000 km 2 ), it may produce 15 psu salinity changes within top 1 m surface water layer, or 7.5 psu within top 2 m.As seen in Figure 6 (right column), the areas with positive SSS differences that are exceeding 10 psu covering about half of the Kara Sea areas, which is in the same order of magnitude as the freshening effect that was possibly produced by the river discharge difference in the two years.It is understood that large differences in SSS anomaly may depend on how the discharge freshwater transported horizontally and vertically, in terms of depth and spread of the diluted water body.Nevertheless, the agreement in order of magnitude is encouraging.

SSS Variability at Arctic Ocean Gateways
We have demonstrated that SMAP SSS retrieves reasonably good quality data in the Arctic Ocean in terms of comparisons with in situ salinity data (Section 3.1) and in response to river discharge (Section 3.3).Here, we further explore the feasibility of using satellite SSS to monitor the surface salinity variability at Arctic Ocean Gateways.We examine if the SSS variability at the gateways are significantly greater than the retrieval accuracy of 1 psu.
The Arctic Ocean exchanges freshwater with the sub-oceans through four major gateways, as shown schematically in the currents map (Figure 9a from Figure 2a in [14]): the Bering Strait inflows of relative fresh Pacific waters; the Barents Sea Opening (BSO) and part of the Fram Strait inflow of the salty Atlantic water; and, through the Davis Strait and the part of the Fram Strait between Greenland and Svalbard, which comprise the major outflow locations of water modified by the Arctic Ocean in addition to freshwater that is associated with sea ice flux and subsequent melt [12,14].Because the general directions of the ocean currents at gateways are known, the surface salinity variability observed at the gateways carries the information of the freshwater passing through the regions in the upper ocean.
We first describe the procedure to extract SSS information at each of the Arctic gateways.The exact locations of the gateways defined for this study are given in Table 3 and are indicated in Figure 9b.To ensure currents flow through the passage in a roughly unique direction, we divided the Fram Strait into two parts: one between east Greenland coast and 0 • longitude (named EG), the other from 0 • longitude to the west coast of Svalbard (named Fram Strait).We also added a section to extend the narrow Bering Strait from 62 • N to 68 • N so sufficient grid points can be included in the calculation.According to this definition, a maximum of around 250 to 500 grid points can be extracted from SMAP L3 daily maps with 0.25 • spacing, as listed in Table 3.The actual number of grid points with valid SSS retrieval changes with time depending on the ice situation, which varies with the gateway.This is illustrated in terms of the percentage of the ice-free area at each gateway (Figure 10, right column).We note that BSO and the Fram Strait are mostly ice-free year round, while EG and the Davis Strait are only open for about three months in summer, with more than 50% ice cover the rest of the year.The Bering Strait, which is the only gateway connecting the Arctic Ocean and Pacific Ocean, opens about half of the year from June to December.The time series of SSS are obtained by daily averaging over all grid points with valid retrievals at each gateway.In the extremely challenging environment of the Arctic Ocean, we designed two criteria to avoid possible false retrievals.First, we include in the time series calculation only if >30% of the grid points are ice-free.The daily mean and standard deviation were obtained from 30% or more grid points at the gateway with valid SSS data, named SSS Gateway (t) and δ Gateway (t).Examining the time series, we found extremely large values of δ Gateway (t) often occurred at seasonal transitions of ice melting or freezing, suggesting the possibility of undetected ice contamination.Therefore, we defined the second criteria to exclude those outliers.We calculate the standard deviation of SSS Gateway (t) over the whole period (σ Gateway ), and defined outliers as those with δ Gateway (t) exceeding twice of σ Gateway.After filtering with these two criteria, we apply a 30 days moving average on SSS Gateway (t).
Remote Sens. 2018, 10, x FOR PEER REVIEW 15 of 23 The daily mean and standard deviation were obtained from 30% or more grid points at the gateway with valid SSS data, named SSSGateway(t) and Gateway(t).Examining the time series, we found extremely large values of Gateway(t) often occurred at seasonal transitions of ice melting or freezing, suggesting the possibility of undetected ice contamination.Therefore, we defined the second criteria to exclude those outliers.We calculate the standard deviation of SSSGateway(t) over the whole period (Gateway), and defined outliers as those with Gateway(t) exceeding twice of Gateway.After filtering with these two criteria, we apply a 30 days moving average on SSSGateway(t).10).In this study, the Fram Strait is divided into two parts where currents flow in opposite directions: one between east Greenland coast and 0 longitude, the other from 0 longitude to the west coast of Svalbard.Exact coordinates of the four corners of the polygon used for each gateway are given in Table 3.   3.   The time series of SSS at Arctic gateways that were observed by SMAP during the past three summers reveal very rich information (Figure 10, left column).First, it correctly reflects the typical characteristics of the freshwater sources from different ocean basins: relative salty water from the Atlantic and fresh water from the Pacific, with more than 3 psu differences between mean SSSBSO and SSSBering (Table 3).SSSFram is similar to SSSBSO in terms of mean and variability; both bring warm and The time series of SSS at Arctic gateways that were observed by SMAP during the past three summers reveal very rich information (Figure 10, left column).First, it correctly reflects the typical characteristics of the freshwater sources from different ocean basins: relative salty water from the Atlantic and fresh water from the Pacific, with more than 3 psu differences between mean SSS BSO and SSS Bering (Table 3).SSS Fram is similar to SSS BSO in terms of mean and variability; both bring warm and salty Atlantic water into the Arctic Ocean.Also note, SSS Fram is slightly fresher than SSS BSO , possibility due to its vicinity to Greenland.On the other hand, SSS EG and SSS Fram , although being geographically adjacent to each other, bear dramatic contrast in the properties of water transport out from (SSS EG , mean 31.2 psu) and into (SSS Fram , mean 34.8 psu) the Arctic Ocean.The much fresher signature and larger dynamical range of SSS EG suggest the impact of Arctic ice.We note the exceptionally large drop in SSS Fram (~3 psu) and SSS EG (~6 psu) from late 2016 to early 2017.It is unlikely that this large interannual change is caused by the SMAP algorithm or long-term calibration changes because the changes at BSO are much smaller.At this moment, we do not have enough data to pin point the exact mechanism causing such freshening.Also, the magnitude of SSS EG may need further calibration due to less coverage (~40% ice-free area).We speculate that this freshening and recovering in the later part of 2017 is influenced by the Greenland sea ice melting and freezing.This hypothesis may be supported by the observation that the Davis Strait, another passage like EG where outflows of Arctic Ocean modified water mixed with sea ice drift to the northern Atlantic.Different from SSS EG , SSS Davis shows similar seasonal cycles for the three summers that were captured by SMAP.
Consistent seasonal cycles are also observed over the Bering Strait, with a peak-to-peak range of about 4 psu.It is also interesting to note the inter-annual difference.SSS Bering in June-December 2016 (red curve in Figure 10i) is exceptionally low, about 2 psu less than the same period in 2015 (black curve in Figure 10i); while it recovered partially one year later, as seen in Figure 10i, the green curve (June-December 2017) falls between the previous two years.Year 2016 was the strongest El Niño in recent two decades.One well-known characteristic of an El Niño event is the zonal displacement of the western equatorial Pacific warm/fresh pool.The edge of the pool extends eastward during El Niño, retreats westward during La Niña [73][74][75].
Moorings are installed at all major Arctic Ocean gateways, measuring salinity at depth ~50 m year round.At the close of this study, we only obtained mooring data at the Bering Strait up to the summer of 2016 [24,69], with about one year overlapping with SMAP data.Figure 11 illustrates the time series of SMAP SSS Bering (from Figure 10i) along with salinity form June 2015 to June 2016 at three mooring sites: A2 (66 • 19 N, 168 • 57 W) and A3 (65 • 46 N, 168 • 34 W) in the US channel of the strait, and A4 (65 • 44 N, 168 • 15 W) close to the Alaskan coast.All three mooring sites are within the white polygon, defined as the Bering Strait (Figure 9b).They are within 90 km from each other.During August-November, the salinities from the three moorings were mostly within 1 psu from each other, indicating a relatively small spatial variability.However, the A4 mooring data reduced by about 5 psu from November to December, whereas the other two moorings had smaller changes.The SMAP SSS Bering at 1-2 cm depth depicted a change of about 5 psu, well resolved by the 1 psu accuracy (Section 3.1), and appeared to have a consistent trend with the A4 mooring on the amplitude and timing of the freshening peak in December 2015.Nevertheless, the differences with the other two moorings were clearly larger.A2 and A3 moorings revealed less variability relative to SMAP, which is likely due to the vertical and horizontal dilution from the surface to 50 m depth.A longer time series of data is required to bring a more definite conclusion on the consistency or discrepancy, possibly resulting from differences in sensing depth or spatial variability.
accuracy (Section 3.1), and appeared to have a consistent trend with the A4 mooring on the amplitude and timing of the freshening peak in December 2015.Nevertheless, the differences with the other two moorings were clearly larger.A2 and A3 moorings revealed less variability relative to SMAP, which is likely due to the vertical and horizontal dilution from the surface to 50 m depth.A longer time series of data is required to bring a more definite conclusion on the consistency or discrepancy, possibly resulting from differences in sensing depth or spatial variability.

Discussion
This study explores the utility of SMAP SSS in monitoring the Arctic freshwater system.We emphasize the importance to assess the uncertainty and understanding of the limitations of satellite SSS at this explorative stage in the extremely challenging environment of the Arctic Ocean.We discuss several related issues, including the complementarity of satellite SSS with other satellite and in-situ observations.First, we point out the urgent need for in situ measurements of salinity in the polar regions to enhance satellite salinity products.Currently, there are very few in situ salinity data available in the areas North of the Arctic Circle (Section 3.1).A sufficiently large database of in situ salinity data (in conjunction with simultaneous geophysical information, such as SST, wind speed and direction, air temperature, and ocean current) will be invaluable.It will not only provide a robust assessment of satellite SSS uncertainties in the Arctic Ocean, but also support the retrieval algorithm refinement, including the SST effect on the roughness correction in cold water.
Second, improved retrievals of sea ice concentration would be of considerable value to reduce possible ice-contaminated SMAP SSS retrievals.It is known that the quality of sea ice concentration product degrades in thin ice, especially at the ice margins and during summer melt.The method of identifying Frazil-pancake ice based on SAR imagery [72] with the wave dispersion technique may have an advantage over passive microwave, as the sea ice concentrations may be quite variable over short time scales [76].The added information of thickness of this type of young ice [72] will allow quantitative analyses of the SSS response to ice formation and related salt rejection during ice growth on the mixed layer budget.
Third, integrated data analysis is important and necessary to realize the benefit of satellite SSS.One example is using satellite SSS to trace the river discharge.Because satellite observes salinity in the top 1-2 cm of sea surface, SMAP SSS only reveals freshwater variability in the top layer.In stratified regions, SSS alone cannot give the volume of the freshwater distribution, which is traced by satellite from the surface, such as in the Kara Sea (Section 3.3).Quantitative estimation of river discharge impact and distribution over the Arctic Ocean can only be achieved by combining SSS with the salinity stratification structure of the upper ocean layer, based on in situ measurements of salinity profiles and ocean circulation models.Another example is the need of ocean currents in estimating freshwater transport.We demonstrated that the SMAP SSS variability could be used to discern the freshwater transport through Arctic Ocean gateways.However, to quantify the freshwater transport that is associated with those SSS variability, one also needs information of the magnitude and direction of ocean currents in addition to salinity profiles.
Finally, we want to emphasize that the value of SMAP SSS relative to data assimilation products, such as HYCOM, cannot be judged from the comparison with co-located in-situ measurements alone.We have shown (e.g., Section 3.3) that SMAP data reveal much stronger interannual variation of SSS than HYCOM.Given HYCOM's relaxation of SSS to seasonal climatology and the use of seasonal climatology of river discharges, SMAP SSS may provide a more reliable dataset to study interannual freshwater changes of the Arctic Ocean in regions where no or few in-situ data are available to constrain HYCOM.

Conclusions
Here, we demonstrated the potential utility of SMAP SSS in studying Arctic Ocean freshwater changes.The JPL SMAP SSS team adopts an advanced retrieval algorithm near the land and ice edge, aiming to provide the scientific community potentially useful information.We assessed the accuracy of SMAP SSS retrieved North of 50 • N through comparative analysis with in situ salinity data collected by Argo floats, ships, gliders, and in field campaigns.Results that were derived from near 20,000 pairs collocated within a 12.5-km radius and daily time window indicate a RMSD less than ~1 psu and a correlation coefficient ~0.82.We have demonstrated the feasibility of using SMAP SSS to study Arctic Ocean salinity response to changes in sea ice concentration and river discharge.In the Arctic Ocean, SMAP SSS reveals seasonal and interannual variations that are consistent with sea ice concentration.In the Kara Sea, the large SSS contrast between 2015 and 2016 observed by SMAP is corroborated by the anomalous freshwater inputs from the Ob' and Yenisey Rivers.We also demonstrated the feasibility to use SSS as a proxy in monitoring the upper layer freshwater exchanges between the Arctic and sub-Arctic Oceans and through the four major Arctic Ocean gateways.
Results presented in this paper may lead to interesting future studies.Particularly, to understand the surface salinity variability that is observed by SMAP at Arctic Gateways, it will be interesting to examine its relation with other variables, such as SST, wind and current, locally and beyond.In conjunction with in situ data, satellite SSS can provide complementary and critical surface information to advance our knowledge on the linkage between the climate changes in polar regions and sub-oceans.
25 • × 0.25 • grids daily from 1 April 2015 to 31 March 2018.Measurements of different instruments, such as CTD (Conductivity, Temperature and Depth), XCTD (Expandable CTD), and TSG (Thermosalinograph) taken within 5 m of the sea surface are used.Figure 1 shows the distribution of in situ measurements North of 50 • N on 0.25 • grids, as long as there is as least one daily record in three years on the grid point.Areas South of 60 • N are mostly covered by Argo floats, complemented by a few ship tracks.Areas between 60 • N and 80 • N in Greenland area are well covered combining ships, gliders, and OMG data.However, above the Arctic Circle, particularly North of the Bering Strait, there are almost no in situ salinity data.Remote Sens. 2018, 10, x FOR PEER REVIEW 5 of 23 long as there is as least one daily record in three years on the grid point.Areas South of 60°N are mostly covered by Argo floats, complemented by a few ship tracks.Areas between 60°N and 80°N in Greenland area are well covered combining ships, gliders, and OMG data.However, above the Arctic Circle, particularly North of the Bering Strait, there are almost no in situ salinity data.

Figure 1 .
Figure 1.Location where in situ salinity are collected within 5 m from surface by Argo float (green), ship (red), glider (purple), and OMG AXCTD (blue) during the three years from April 2015 to March 2018.

Figure 1 .
Figure 1.Location where in situ salinity are collected within 5 m from surface by Argo float (green), ship (red), glider (purple), and OMG AXCTD (blue) during the three years from April 2015 to March 2018.

23 Figure 2 .
Figure 2. Salinity maps North of 50N averaged from 1 April 2015 to 31 March 2018 from (a) in situ data, (b) Soil Moisture Active Passive Sea Surface Salinity (SMAP SSS), and (c) Hybrid Coordinate Ocean Model (HYCOM) SSS, and the difference of (d) SMAP SSS minus in situ, and (e) HYCOM SSS minus in situ, wherever absolute differences less than 1 psu are masked out in white.

Figure 2 .
Figure 2. Salinity maps North of 50 • N averaged from 1 April 2015 to 31 March 2018 from (a) in situ data, (b) Soil Moisture Active Passive Sea Surface Salinity (SMAP SSS), and (c) Hybrid Coordinate Ocean Model (HYCOM) SSS, and the difference of (d) SMAP SSS minus in situ, and (e) HYCOM SSS minus in situ, wherever absolute differences less than 1 psu are masked out in white.

Figure 3 .
Figure 3. Scatterplot of in situ salinity collected within 5 m from surface and collocated data of (a) SMAP SSS and (b) HYCOM SSS.Collocations are based on gridded data on 0.25 grids averaged daily from April 2015 to March 2018.The three clusters associated with low in situ salinity (2-7 psu) are located in the Gulf of Bothania and the Baltic Sea.Open circles indicate where the absolute difference between the pair exceeds three times of the standard deviation.Linear regression of all data (excluding outliers) bin-averaged in terms of in situ SSS (bin size of 1 psu) is illustrated in red with zero intercepted denote by C0 and slope by C1, while the regression of high SSS (>25 psu, excluding outliers and the low salinity clusters.) is illustrated in black.Red dots indicate bin-average and vertical error bars indicate 1 standard deviation in the bin.Dotted line indicates the diagonal for reference.

Figure 3 .
Figure 3. Scatterplot of in situ salinity collected within 5 m from surface and collocated data of (a) SMAP SSS and (b) HYCOM SSS.Collocations are based on gridded data on 0.25 • grids averaged daily from April 2015 to March 2018.The three clusters associated with low in situ salinity (2-7 psu) are located in the Gulf of Bothania and the Baltic Sea.Open circles indicate where the absolute difference between the pair exceeds three times of the standard deviation.Linear regression of all data (excluding outliers) bin-averaged in terms of in situ SSS (bin size of 1 psu) is illustrated in red with zero intercepted denote by C 0 and slope by C 1 , while the regression of high SSS (>25 psu, excluding outliers and the low salinity clusters.) is illustrated in black.Red dots indicate bin-average and vertical error bars indicate ±1 standard deviation in the bin.Dotted line indicates the diagonal for reference.

Figure 4 .
Figure 4. Salinity maps North of 50N from SMAP (left) and HYCOM (right) for the months of August 2015, 2016, and 2017 (top to bottom).Marginal seas are indicated in the top left panel.Black dots in the top right panel indicate locations of river discharge data are collected.

Figure 4 .
Figure 4. Salinity maps North of 50 • N from SMAP (left) and HYCOM (right) for the months of August 2015, 2016, and 2017 (top to bottom).Marginal seas are indicated in the top left panel.Black dots in the top right panel indicate locations of river discharge data are collected.

Figure 5 .
Figure 5.The time series of SMAP SSS (black) and HYCOM SSS (green) averaged in the domain North of 65N in the ice-free-area (red) over all retrievals (solid line), or with outliers removed (dotted line).Outliers are identified as where the absolute differences from the daily mean exceeding three times of the standard deviation (~1 to 5%).Ice-free-area is normalized (nIFA) by the record maximum (observed in September 2017).Sea ice concentration (SIC) (blue) is averaged daily North of 65N including zero and non-zero values.

Figure 5 .
Figure 5.The time series of SMAP SSS (black) and HYCOM SSS (green) averaged in the domain North of 65 • N in the ice-free-area (red) over all retrievals (solid line), or with outliers removed (dotted line).Outliers are identified as where the absolute differences from the daily mean exceeding three times of the standard deviation (~1 to 5%).Ice-free-area is normalized (nIFA) by the record maximum (observed in September 2017).Sea ice concentration (SIC) (blue) is averaged daily North of 65 • N including zero and non-zero values.

23 Figure 6 .
Figure 6.SSS in Kara Sea for the months from May to October (top to bottom) for the year 2015 (left) and 2016 (middle).The corresponding differences (2016 minus 2015) are shown in the right column.Sampling locations of the Ob and Yenisey River are indicated by black dots.

Figure 6 .
Figure 6.SSS in Kara Sea for the months from May to October (top to bottom) for the year 2015 (left) and 2016 (middle).The corresponding differences (2016 minus 2015) are shown in the right column.Sampling locations of the Ob and Yenisey River are indicated by black dots.

Figure 8 .
Figure 8. Daily river discharge from (a) Ob' and (b) Yenisey Rivers during May 2015 to April 2016 (black) and from May 2016 to April 2017 (red).

Figure 8 .
Figure 8. Daily river discharge from (a) Ob' and (b) Yenisey Rivers during May 2015 to April 2016 (black) and from May 2016 to April 2017 (red).

Figure 8 .
Figure 8. Daily river discharge from (a) Ob' and (b) Yenisey Rivers during May 2015 to April 2016 (black) and from May 2016 to April 2017 (red).

Figure 9 .
Figure 9. (a) From Carmack et al., 2015 [14].Schematic maps of the major ocean currents (long arrows), the four Arctic Ocean gateways in Fram Strait, the Barents Sea Opening, Davis Strait, and Bering Strait (thick bars with red denoting inflow and blue denoting outflow), the gyral circulation patterns (circular arrows), the salt-stratified ocean domains are shown in light blue, and the All Arctic Regions definition of the terrestrial contributing areas shown in white.(b) Locations of the Arctic Gateways indicated on SMAP SSS gridded fields (three years mean), which were used in the calculation of mean SSS and the percentage of ice-free area (Figure10).In this study, the Fram Strait is divided into two parts where currents flow in opposite directions: one between east Greenland coast and 0 longitude, the other from 0 longitude to the west coast of Svalbard.Exact coordinates of the four corners of the polygon used for each gateway are given in Table3.

Figure 9 .
Figure 9. (a) From Carmack et al., 2015 [14].Schematic maps of the major ocean currents (long arrows), the four Arctic Ocean gateways in Fram Strait, the Barents Sea Opening, Davis Strait, and Bering Strait (thick bars with red denoting inflow and blue denoting outflow), the gyral circulation patterns (circular arrows), the salt-stratified ocean domains are shown in light blue, and the All Arctic Regions definition of the terrestrial contributing areas shown in white.(b) Locations of the Arctic Gateways indicated on SMAP SSS gridded fields (three years mean), which were used in the calculation of mean SSS and the percentage of ice-free area (Figure10).In this study, the Fram Strait is divided into two parts where currents flow in opposite directions: one between east Greenland coast and 0 • longitude, the other from 0 • longitude to the west coast of Svalbard.Exact coordinates of the four corners of the polygon used for each gateway are given in Table3.

Figure 10 .
Figure 10.Time series at Arctic Ocean gateways of 1 May 2015-30 April 2016 (black), 1 May 2016-30 April 2017 (red), and 1 May 2017-30 April 2018 (green) of (left) SMAP SSS averaged and (right) percentage of the ice-free area over Arctic gateways defined in Figure 9b.(a,b) The Barents Sea Opening (BSO); (c,d) the Fram Strait east of 0 longitude; (e,f) the Fram Strait opening between east coast of Greenland and Svalbard (EG); (g,h) the Davis Strait; and, (i,j) the Bering Strait.30 days moving average is applied on SSS time series (left column, thick lines).

Figure 10 .
Figure 10.Time series at Arctic Ocean gateways of 1 May 2015-30 April 2016 (black), 1 May 2016-30 April 2017 (red), and 1 May 2017-30 April 2018 (green) of (left) SMAP SSS averaged and (right) percentage of the ice-free area over Arctic gateways defined in Figure 9b.(a,b) The Barents Sea Opening (BSO); (c,d) the Fram Strait east of 0 • longitude; (e,f) the Fram Strait opening between east coast of Greenland and Svalbard (EG); (g,h) the Davis Strait; and, (i,j) the Bering Strait.30 days moving average is applied on SSS time series (left column, thick lines).

Figure 11 .
Figure 11.Daily time series of salinity at the Bering Strait retrieved from SMAP (black), and mooring measurement at mooring sites A2 (red), A3 (greed) and A4 (blue).Daily data are indicated by solid dots and lines illustrate the 8 days moving average.

Figure 11 .
Figure 11.Daily time series of salinity at the Bering Strait retrieved from SMAP (black), and mooring measurement at mooring sites A2 (red), A3 (greed) and A4 (blue).Daily data are indicated by solid dots and lines illustrate the 8 days moving average.

Table 1 .
Results of statistical comparison between in situ salinity and SMAP SSS (JPL V4.0) and HYCOM SSS, for North of 50 • N or North of 65 • N, respectively.All data are gridded on 0.25 • grids averaged daily from April 2015 to March 2018.N is the number of collocated pairs; bias (data minus in situ), standard deviation (Std.) and Root Mean Square Difference (RMSD) are in psu; Corr.stands for the correlation coefficients.

Table 1 .
Results of statistical comparison between in situ salinity and SMAP SSS (JPL V4.0) and HYCOM SSS, for North of 50N or North of 65N, respectively.All data are gridded on 0.25 grids averaged daily from April 2015 to March 2018.N is the number of collocated pairs; bias (data minus in situ), standard deviation (Std.) and Root Mean Square Difference (RMSD) are in psu; Corr.stands for the correlation coefficients.

Table 2 .
Comparison of Ob' and Yenisey River integrated discharge for the warm seasons of 2015 and 2016.

Table 3 .
List of the mean, standard deviation, minimum and maximum values of SMAP SSS retrieved at the Arctic Gateways from 1 May 2015 to 30 April 2018.

Table 3 .
List of the mean, standard deviation, minimum and maximum values of SMAP SSS retrieved at the Arctic Gateways from 1 May 2015 to 30 April 2018.