From field experiments to salinity products : a tribute to the contributions of Jordi Font to the SMOS mission

Adriano Camps 1,2,3, Carolina Gabarró 3,4, Mercè Vall-llossera 1,2,3, Sebastià Blanch 1, Albert Aguasca 1, Francesc Torres 1,3, Ignasi Corbella 1,3, Nuria Duffo 1,3, Antonio Turiel 3,5, Marcos Portabella 3,4, Joaquim Ballabrera-Poy 3,4, Verónica González-Gambau 3,4, Justino Martínez 3,4, Ramón Villarino 5, Luís Enrique 6, Alessandra Monerris 7, Xavier Bosch 8, Roberto Sabia 9, Marco Talone 10, Maria Piles 1,2,3, Míriam Pablos 1,2,3, Enric Valencia 1


INTRODUCTION
The Remote Sensing Laboratory of the Department of Signal Theory and Communications at the Universitat Politècnica de Catalunya -BarcelonaTech, has been involved since 1993 in the technological aspects of the Microwave Imaging Radiometer with Aperture Synthesis (MIRAS), the single payload of the European Space Agency Soil Moisture and Ocean Salinity (SMOS) Earth Explorer Opportunity mission.We first met Professor Jordi Font in 1998, during the preparation of what later became the SMOS Mission.Since then, we have enjoyed a fruitful scientific collaboration, which has also led to a longstanding friendship.This paper is a tribute to Professor Font's achievements in remote sensing of the ocean over the last 17 years.
Sea surface salinity (SSS) can be remotely measured using microwave radiometry at L-band (Swift et al. 1983).During the late 1940s and early 1950s, laboratory experiments were conducted to quantify the relationship between salinity and brightness temperature.Until SMOS, only two (experimental) space-borne radiometers at L-band had been launched, in Cosmos 243 (1968) and in Skylab S-194 (1973), but the first evidence from space came only from the Skylab observations (Fig. 1).
The link between salinity and brightness temperature is through the dielectric constant of sea water.Though the "protected" 1.400-1.427MHz band is the best choice, these measurements are far from trivial Editor: J.L. Pelegrí.
Copyright: © 2016 CSIC.This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-by) Spain 3.0 License.Fig. 3. -A, first sea surface salinity map derived from L-band microwave radiometry in 1996 (Miller et al. 1996); B, first sea surface salinity map derived from a synthetic aperture radiometer in 2000 (LeVine et al. 2000).
due to the low sensitivity of brightness temperature to SSS, 0.2-0.8K/psu (psu = practical salinity unit, roughly 1 g of salt per 1 kg of water), the dependence on the physical temperature (Fig. 2), incidence angle, polarization, surface roughness (induced by local wind or by far waves), the technical difficulties of achieving a precise calibration and stability, and especially the presence of radio-frequency interference (RFI).
The first SSS airborne map was published in 1996 by Miller and Zaitzeff (Miller et al. 1996) using the scanning low frequency microwave radiometer (SLFMR), a six-beam real aperture radiometer (Fig. 3A).In 2000, LeVine et al. obtained an SSS map using the electronically steered thinned array radiometer (ESTAR), the first in-orbit 1D synthetic aperture radiometer (LeVine et al. 2000) (Fig. 3B).
In 1995, the ESA organized the Soil Moisture and Ocean Salinity Workshop to define the roadmap to remotely measure these two geophysical variables, and it was concluded that the most promising technique was aperture synthesis radiometry, which had been successfully demonstrated a few years earlier with the ESTAR instrument (Ruf et al. 1998).In 1998, an international team of scientists proposed the SMOS mission to ESA.It was selected as the second Earth Explorer Opportunity Mission (Silvestrin et al. 2001) within the ESA Living Planet Programme and was launched on 2 November 2009 from Plesetsk, Russia (Figs 4 and 5).

TOWARDS A BETTER UNDERSTANDING OF SEA SURFACE EMISSION
In 1998, when the SMOS proposal was submitted, as in any mission, there were many uncertainties and many lessons to be learnt.For ocean salinity retrieval, three key goals related to the geophysical model function had to be met: -to more accurately characterize the impact of salinity and temperature on the dielectric constant, -to account for the impact of the sea state (wind speed, significant wave height, etc.) on the brightness temperatures, and -to account for the impact of sea foam on the brightness temperatures.

Towards a more accurate determination of the sea water dielectric constant
Two of the most widely used dielectric constant models are those of Klein and Swift (1977) and Ellison  (1998).Although they show similar results (see Fig. 6), their differences are too great to meet the accuracy requirements, i.e. for an absolute accuracy of 0.2 psu, the brightness temperature must be in the order of 0.1 K for cold waters (Fig. 2).
Most of the methods used to determine the dielectric constant are based on reflectivity measurements in a transmission line, such as the open-ended coaxial sensor.However, these methods require well-known calibrators and they are not accurate enough for this application.Other techniques based on the change in the resonant properties of a cavity are very accurate but provide more limited information (single frequency, temperature and salinity).Finally, transmission methods provide accurate and more comprehensive information from the measurement of the propagation changes of the electromagnetic signal (S-parameters measured with a vector network analyser).This was the technique selected in Blanch and Aguasca (2004) to derive the dielectric constant of sea water.The test device is shown in Figure 7.It consists of a metallic waveguide through which temperature controlled salty water can circulate by means of a pump.The dimensions were calculated so that a single mode propagates even when the device is filled with salt water.
The main results are summarized in Figure 6.The Klein and Swift (KS) and Blanch and Aguasca (BA) models show similar dielectric constant results for medium temperatures, whereas the Ellison model is slightly different, especially in the real part of the permittivity.Also, the KS and BA models give similar results on brightness temperature, while the Ellison model tends to overestimate it.However, it is difficult to ensure absolute accuracy because of the lack of absolute references.

Towards a better understanding of the impact of sea state
Since 1971, no systematic measurements had been performed to assess the impact of the sea state (parameterized in terms of the wind speed, the significant wave height or other descriptors).Additionally, only a few airborne experiments had been carried out, mostly with the SLFMR real aperture radiometer (Miller et al.The instruments deployed were (Fig. 9): LAURA (L-band AUtomatic RAdiometer), the first ever fully polarimetric (4 Stokes elements) L-band radiometer (UPC), a Ka-band polarimetric radiometer (UMass), a stereo-camera to determine surface topography and root mean square slopes of the sea surface (CETP), four oceanographic and climatological buoys for nearsurface salinity, temperature, wind speed and direction, wave height and period, etc. (ICM-CSIC and LODYC), a portable meteorological station with atmospheric pressure, air temperature, relative humidity and rain rate (UPC), a video camera mounted on the L-band radiometer pedestal to determine sea surface foam coverage (UPC), and an IR radiometer to determine sea surface temperature estimates (UV).Additionally, during WISE 2000, temperature and salinity were recorded from the platform at 5 m below sea level, and simultaneously ocean colour, wind vector and sea surface temperature were acquired by different satellites.
Professor Font and Agustí Julià designed an ad hoc buoy for WISE 2000, to collect conductivity and temperature data near the sea surface close to the radiometer field of view and send them to a data logging station installed on the platform using a real-time link (Figs 9E and 9F).The buoy was built by EMS Environmental Monitoring Systems S.L. and modified for WISE 2001 mainly to host extra power batteries.It was a toroidal body with a stainless steel structure to house the signalization elements (flash, radar reflector and satellite ARGOS beacon) and the measuring and trans- mitting instruments.The main instrument was a Sea-Bird MicroCAT system (model SBE37-SM) used to record water temperature and conductivity for further salinity determination.An RS-232 interface allowed real-time data transmission by an external UHF link.An additional submersible pump was added to ensure a constant water flow through the conductivity cell.The water inlet was situated at 20 cm below sea level in the central part of the toroid, to minimize the effect of waves that might cause air bubbles to be introduced into the measuring cell.Salinity was computed according to established standards with 0.003 psu/month stability, and 0.0002 psu resolution.The conductivity cell was equipped with a chemical poison device to avoid biofouling and degradation of the conductivity measurement.
The main accomplishment of the WISE campaigns was probably the determination of the brightness temperature dependence with respect to the 10 m height wind speed and the significant wave height required to perform the necessary corrections to retrieve SSS (Camps et al. 2002, Camps et al. 2004, Gabarró et al. 2004), including a weak azimuthal signature, which was more clearly seen during a very intense storm that occurred in November 2001: The processing of the data revealed the following: -A sensitivity to wind speed extrapolated at nadir of ~0.23 K/(m s -1 ), or a little bit higher, ~0.25 K/(m s -1 ), when the atmospheric instability or only the measurements corresponding to U 10 >2 m s -1 are accounted for.This sensitivity increases at H-polarization up to ~0.5 K/(m s -1 ) at 65°, and decreases at V-polarization down to ~-0.2 K/(m s -1 ) at 65°, with a zero-crossing around 50°-55°, in agreement with the small slope approximation method using Durden-Vesecky times 2 and Elfouhaily et al. sea spectra.It is very likely that the computed wind speed sensitivities below 2 m s -1 are erroneous.
-A modulation of the instantaneous brightness temperatures due to wave slopes (and also foam), which makes the standard deviation of this modulation increase with wind speed at a rate of ~0.1-0.15K/(m s -1 ), depending on polarization, and very weakly on incidence angle.
-A sensitivity to significant wave height extrapolated to nadir of ~1 K/m, increasing at H-polarization up to ~1.5 K/m at 65°, and decreasing at V-polarization down to -0.5 K/m at 65°. -A small azimuthal modulation of ~0.2-0.3K for lowto-moderate wind speeds, in reasonable agreement with the numerical models, but very large (peak-to-peak modulations of 4-5 K) during a strong storm recorded on 10 November 2001 (Fig. 10).
The brightness temperature sensitivity to wind speed and significant wave height compared satisfactorily with numerical models; these sensitivities were used until calibrated SMOS data were available, from which new roughness models were formulated (Gabarró et al. 2004, Guimbard et al. 2012).
EuroSTARRS was an ESA campaign in support of the SMOS mission that took place from 17-24 November 2001, while the WISE 2001 campaign was taking place.The US STARRS (Salinity Temperature and Roughness Remote Scanner) L-band radiometer was mounted on a DLR plane and overflew several regions of Europe (see Fig. 11A).The main objective of Eu-roSTARRS 2001 was to acquire SMOS-like observations to address a range of critical issues relevant to the SMOS mission.Several flights were dedicated to ocean salinity experiments and were coordinated by Jordi Font.The experiment was defined, among other objectives, to allow the data acquired to estimate multiangular L-band microwave signatures from low-altitude flights; to compare data acquired at low altitude with data acquired at high altitude; to investigate emission from urban areas and the level of RFI, to validate ocean salinity retrievals along salinity gradients from the coast out to the sea; and to investigate azimuthal effects on the emission over sea.One of the sites was around the Casablanca oil rig, near the mouth of the Ebro River (Tarragona), mainly to study the effect of wind speed on the salinity measurements and to investigate azimuth dependence.In the Casablanca area, the flight acquisition plan consisted of three phases (see Fig. 11B, C).The first one was a straight transect from the coast towards the platform, which was selected to avoid interferences from the platform.The second consisted in performing ten circles at a constant bank angle of 22° to obtain a large range of incidence angles.The last was the reverse course along the first transect flown at low altitude.The measurements were made after sunshine, in order to avoid interference from sun glint.
The Institut de Ciències del Mar (ICM) participated actively during the flight over the Casablanca oil rig on 21 November.Simultaneously to the airborne flight, two types of measurement were made over the area.The CSIC's R/V Garcia del Cid carried out a survey in a rectangular area around the Casablanca rig coincident with the plane overflight, collecting oceanographic data, such as underway near-surface temperature and salinity records, CTD profiles spread across the site, acoustic Doppler current profiles and meteorological measurements on-board in a continuous mode.The results were published (Font et al. 2003, Gabarró et al. 2003, Emelianov 2003).

Understanding the impact of sea foam
One of the open questions from WISE 2000 and 2001 was the evident modulation of the brightness temperatures when a patch of foam appeared in the antenna beam (e.g.Fig. 10b).This was confirmed by the infrared images, in which the foam patches appeared as hot spots.To quantify this effect, as well as the impact of rain and oil spills, and to validate the existing numerical models, the FROG 2003 (Foam, Rain, Oil spills and GPS-reflections) experiment was performed at the IRTA facilities at Poble Nou del Delta (Tarragona, Spain) to study specifically the impact of foam, rain, and oil spills on sea surface emissivity (Camps et al. 2005).Figure 12A shows the experimental set-up in which the LAURA radiometer was installed on the side of a 3×7 m pool filled with a mixture of fresh and sea water.A total of 104 air diffusers distributed over the pool produced a fairly uniform foam coverage.Foam structure was analysed using video images (e.g.Fig. 12B) that were processed to infer the bubble size and distribution, and an array of conductivity sensors was used to infer the air/ water fraction as a function of depth (Camps et al. 2005).
Results indicated that foam impact is not negligible and should be included in the emission model at L-band.Rain (Fig. 12C) and a thin layer of oil on the surface produce a modification of the water surface roughness that is sufficient to be detectable at L-band.
Ultimately, the impact of sea foam can be parameterized as an emissivity increment, multiplied by the foam fraction F(U 10 ): Actually, F(U 10 ) is also dependent on the air-sea stability and salinity, among others, and can be quite variable (see Fig. 13).
The WISE and FROG campaigns were very useful before the SMOS launch to confirm the general behaviour of the models as a function of the wind speed and foam, but once the satellite had been launched the models needed to be recomputed using the new SMOS data, and the BEC group worked hard to improve the emissivity model (Guimbard et al. 2012), which is still being used as one of the three parameterizations in the official processing chain.

GNSS-R TECHNIQUES TO IMPROVE SEA SURFACE SALINITY RETRIEVALS FOR SMOS FOLLOW-ON MISSIONS
As explained, the capability to properly correct for the sea state impact on the brightness temperature is limited; that is, the remaining uncertainty is too large to meet the scientific requirements.In this Section the main outcome of a project conducted to assess the reliance of GNSS-R (Global Navigation Satellite Systems Reflectometry) techniques to perform a better sea state correction is presented.We take advantage of the closeness in frequency of the protected passive microwave band (1400-1427 MHz) to the radio-navigation bands (L1, L2, and L5, mainly), since both emission and scattering should be affected by the same range of mechanical wavelengths (Camps et al. 2009).
The underlying concept in GNSS-R is quite simple: when the electromagnetic wave is scattered over the sea surface, the scattered signal changes its polarization (from right-hand circular to mostly left-hand circular) and arrives at the receiver mainly from the specular reflection point, determined by the shortest distance between the transmitting GPS satellite (Glonass, Beidou or Galileo) and the receiver.However, as the sea gets roughened, the scattered signals come from a wider region (known as the "glistening zone") that increases in size with increasing roughness, in a similar manner as the Sun reflecting over the sea.
In classical GNSS-R, observables are derived by cross-correlating the received scattered signal with a local replica of the pseudo-random noise signals (i.e. the open C/A code in the case of GPS) transmitted by the navigation satellites, generated at different delays and Doppler frequency shifts.Actually, the values of this two-dimensional cross-correlation are called the Delay-Doppler Map (DDM).When the sea gets rougher, the volume below the normalized DDM (peak equal to 1) down to a predetermined threshold larger than the noise level increases.
The field experiments ALBATROSS 2008 and 2009 performed in the Gran Canaria island already demonstrated that the changes in the brightness temperatures at different polarizations and incidence angles were correlated to the changes in the volume of the DDM, and that this information could be used synergetically to correct for the sea state.Actually, it was also found that the length of the tail of the waveform (the cut of the DDM for Doppler shift equal to 0 Hz) at 1/e was also highly correlated to the sea state (Camps et al. 2009, Marchan et al. 2009, Talone et al. 2009).This last result was also used in the COSMOS-OS experimental data to improve the retrieval of the salinity.The experiment conditions were a number of flights perpendicular to the coast in the Gulf of Finland passing from fresh water in a lake close to the coast to (a slightly) salty water.The results are shown in Figure 14.In Figure 14A the peaks of the brightness temperature correlate well with the peaks in the normalized waveform area, suggesting a sea state correction without other auxiliary data.Figure 14B shows the improvement produced by these corrections (green dots), as compared with the retrievals using raw data.

THE SMALL AIRBORNE MIRAS CAMPAIGN
The Small airborne MIRAS (AMIRAS) was a reduced version of the MIRAS: Four elements per arm with a unit spacing of 0.875 λ and a noise injection radiometer in the centre.Most of the hardware units, including front-ends, correlators and antennas, were first-version prototypes in the framework of the SMOS project.AMIRAS had been specifically designed mechanically and thermally to fly on board the Short SC-7 Sky-van of the Helsinki University of Technology (TKK).The demonstrator was first installed in an anechoic chamber at UPC, where it underwent a number of characterization tests, including measurement of the antenna patterns.
A ground-based experiment took place with the instrument placed in the open air pointing to the cold sky (IRTA facilities, Spain).This experiment was used to measure the instrument's response to a constant target and also to make an image of the Galaxy over-pass (Martin-Neira et al. 2008).Figure 15 shows the instrument placed at the IRTA facilities and the first image of the Milky Way obtained by AMIRAS.Finally, the demonstrator was installed on board the TKK sky-van and was flown three times in order to acquire dual-pol images of coastlines and islands.A sequence of snapshots produced during a flight over an island was presented at IGARSS'06 and is shown in Figure 16.It clearly reproduced the coastal geometry, qualitatively demonstrating the imaging capability of the system.
The limited spatial resolution of the AMIRAS instrument and the small number of snapshots over a single point produced a large amount of dispersion of the results, but they were extremely useful for understanding the specificities of the SMOS field of view (Corbella et al. 2008).

ON-GROUND CHARACTERIZATION OF THE SMOS PAYLOAD
The instrument was successfully assembled by EADS-CASA Espacio (Spain) and it was subjected to an extensive testing process aimed at checking the hardware and software operation as well as validating the calibration procedures and the imaging capability.MIRAS was fully deployed inside the Maxwell anechoic chamber at ESA-ESTEC (Fig. 17) and some tests were performed on it (Fig. 18).They included the following: calibration sequence evaluation, polarization check, stability, inter-element phase retrieval, correlation offset, electromagnetic compatibility and image validation (Corbella et al. 2009).
MIRAS was successfully characterized on-ground in terms of basic operation, thermal cycling, image validation internal phase and RFI.The results showed that the instrument's performance was excellent, having all parameters according to specifications, except for an anomaly detected in the temperature dependence of the phase of correlation and a minor interference of the nominal X-band transmitter.The internal calibration concept and operation was validated, providing highly repeatable and accurate calibration parameters with low temperature sensitivity.Absolute amplitude calibration was demonstrated in the anechoic chamber using the average of all Payload Modules, PMS values, showing that the instrument could measure absolute brightness temperature with an accuracy of 0.1 K. Images of the brightness temperature of the empty chamber showed radiometric accuracies of less than half a kelvin and radiometric sensitivities according to the theoretical expectations as a function of integration time and bandwidth.

BARCELONA EXPERT CENTRE AND PROCESSING CENTRE FOR LEVELS 3 AND 4
Some of the above experimental activities and results were the embryo of the so-called empirical model (Model 3) in the SMOS L2 OS Algorithm Theoretical Baseline Document (ATBD) (SMOS ESL, 2014), which describes all the steps required to perform the retrieval of the ocean surface salinity from the calibrated and geolocated brightness temperatures (Level 1C).
Professor Font's group was one of the three Expert Support Laboratories giving advice to ESA on how to retrieve SSS (the "holy grail" of physical oceanography) from the instrument measurements from 2005.
The SMOS Barcelona Expert Centre on Radiometric Calibration and Ocean Salinity (SMOS-BEC) was created in 2007 as a joint initiative of the Spanish National Research Council (CSIC) and the Universitat Politècnica de Catalunya (UPC) to contribute to the ground segment of ESA's SMOS mission, gathering in a single entity the full spectrum of expertise required within the SMOS data processing, ranging from instrument engineers to oceanographers.The UPC's participation is made through the Passive Remote Sensing Group/Remote Sensing Lab of the Department of Signal Theory and Communications.The main CSIC actor in SMOS is the Physical and Technological Oceanography Department of the ICM, together with the Earth Observation Group of the Institute of Space Sciences (ICE).The SMOS-BEC is installed at the Institute of Marine Sciences (ICM) building, Barcelona, and Professor Font was the head of the group.
BEC's missions were to provide assessment to ESA as a Level 2 Ocean Salinity Expert Support Laboratory, to contribute to SMOS radiometric calibration and validation activities, and to develop and validate new algorithms for generating added-value products at Levels 3 and 4, helping the Spanish Processing Centre, devoted to Levels 3 and 4 (CP34), to maintain the operational generation and distribution of SSS and soil moisture (SM) products.
ESA's mandate is to process instrument data up to Level 2, (geophysical products), while Levels 3 and 4 (spatio-temporal averages and synergistic products) are produced by national agencies.Spain pushed for the CP34 (Processing Centre for Levels 3 and 4), which was initially installed at ESAC (the European Space Astronomy Centre of ESA, at Villafranca del Castillo, Spain).BEC was nominated to define and coordinate the SMOS high level products and distribution system of the CP34.INDRA was the company in charge of developing the operational software to generate the products and the dissemination tools, while the company INSA was in charge of the CP34 maintenance at ESAC.All these activities where supported by several National Plan projects led by Professor Font, called MIDAS.
In 2013, the operations of CP34 at ESAC were discontinued and BEC took over.As BEC started to generate and distribute SSS and SM products, some deep, useful changes with respect to the structure of the former CP34 were introduced: a new, easy to navigate webpage was put in place; the structure of the products was simplified and unified under the netCDF format in MF convention to better meet the needs of the scientific community; and new products and new applications were made available: the new CP34-BEC service was created.BEC's visibility increased enormously thanks to those changes, and now BEC is a reference in the processing and distribution of marine and soil moisture remote sensing products in Spain.
The data produced at BEC are served in netCDF format and are freely distributed to the scientific community by means of OPENDAP, HTTPServer and NetcdfSubset services through a THREDDS server.Maps produced by BEC can also be visualized online by means of a Web Map Service.All products are served at http://cp34-bec.cmima.csic.es/data/available-products.
BEC creates new salinity global maps as new data reprocessing campaigns are performed by ESA and exhaustive quality control reports are made.The product description and quality reports are published at BEC-CP34 (http://cp34-bec.cmima.csic.es/doc/BEC-SMOS-0001-PD.pdf).
Currently, the BEC-CP34 web portal distributes the following products: -Salinity products according to the direction of the SMOS orbit passes: ascending, descending and both products.These products are created in a variety of averaging periods: 3 days and 9 days generated every 3 days, monthly, quarterly and annually.The spatial averaging is computed by default in a regular lat-lon grid of 0.25×0.25°.Several processes have been applied: -Binned products are constructed by simple weighted averaging of the filtered L2 SSS values.The weighted average of SSS follows (Boutin et al. 2012).
-Optimal interpolation (OI) products: Binned products are optimally interpolated (Objective Analysis) to produce maps of higher consistency and fewer gaps as compared with the L3 binned products.The OI is performed using monthly WOA 2009 data as background field.
-L4 Fused products using singularity analysis techniques are obtained with a singularity-analysis-based fusion technique.A template variable of good quality (sea surface temperature) is used as a template to restore the multifractal structure of singularity fronts in a noisy variable (SSS in our case).Further information on the multifractal structure of ocean scalars can be found in Turiel et al. (2009).
-SMOS Mediterranean SSS products: Objectively Analysed 9-day 0.25×0.25°L3 SSS product and Fused 1-day 0.25×0.25°L4 SSS product (Fig. 19).-Land products according to ascending and descending passes -Global Soil Moisture binned products are constructed by simple spatial averaging taking into account only the last orbit measurements in each cell.They are operationally generated by averaging daily, 3-day, 9-day, monthly and annual data on an EASE 25 km grid.
-Soil Moisture Level 4 products: By combining SMOS Level 1C TBs, with Terra/Aqua MODIS land surface temperature and normalized difference vegetation index data, high resolution soil moisture maps (1 km) can be produced over the Iberian Peninsula with a delay of less than 12 hours.(Piles et al. 2014).
-Forest fires risk assessment and senescence: As shown in (Chaparro et al. 2015), fire ignition is more likely to occur in the northeast of the Iberian Peninsula under dry soil conditions (as seen by SMOS).Since 2012, BEC SMOS high resolution soil moisture maps in near real time are being routinely used by the Diputacio de Barcelona (DiBa) for fire risk assessment and prevention during the forest fire season.-Other studies and products in development at BEC are: -Sea ice concentration maps with SMOS data: At 1.4GHz, the brightness temperature values are very different for sea and sea ice, allowing an index of ice concentration to be established.
A new statistical algorithm to retrieve this index with SMOS data is being used to produce daily maps (Gabarró et al. 2015).
-SSS maps at high latitudes (>50ºN) have been produced.It is well known that the sensitivity of the TB to SSS decreases for cold waters (Fig. 2), but by using empirical bias correction technics, it is possible to retrieve salinities from cold waters with reasonable accuracy.
-A sea surface density product has been implemented based on the new equation of state for sea water, the thermodynamic equation of seawater, TEOS10.The product is ready to go into operation, although validation is still needed.-The ice sheet layer in Antarctica with Aquarius and SMOS data: A good agreement was observed between the brightness temperatures and the ice thickness variations over selected transects of eastern Antarctica.The presence of subglacial lakes, which are generally located at ~3 km depth, may influence the physical temperature profile and/or the dielectric properties of the ice layers above, affecting the observed brightness temperatures (Pablos et al. 2015a(Pablos et al. , 2015b)).All these products and the research activities conducted at BEC during the last few years are broadly described in Turiel et al. (2016).The paper reviews the main scientific activities and achievements of the Barcelona Expert Centre under the direction of Jordi Font, as well as the future directions, highlighting the importance of the oceanographic applications of the mission For his dedication to the pioneer SMOS mission, Professor Font received the Catalan National Arts Prize 2011 in the category Scientific Thought and Culture (see Fig. 20).

CONCLUSIONS
This article has summarized some of the contributions of Jordi Font, co-Principal Investigator for Ocean Salinity of the ESA SMOS Earth Explorer Mission.
Among his many direct contributions, and all those triggered by his work, we have emphasized the ICM/CSIC contributions to the field experiments WISE 2000and 2001, and FROG 2003, which helped define the geophysical model function relating the brightness temperature to the sea state.We have also stressed his efforts to establish a new ocean salinity dielectric constant.Finally, we have highlighted his drive towards the creation of the Barcelona Expert Centre, the implementation of the CP34, and the development of new products.We will miss his leadership, but his approach to oceanography will always remain.

Fig. 4 .
Fig. 4. -A, SMOS artist's view in orbit; B, Jordi Font with the ROCKOT launcher behind with SMOS inside, ready to be launched at Plesetsk, Russia.

Fig. 6 .
Fig. 6. -Comparison between the Klein and Swift, Ellison and Blanch and Aguasca dielectric constant models.A, real and imaginary part of the dielectric constant; B, brightness temperature at nadir for 37.5 psu as a function of the temperature (Gabarró 2004).

Fig. 8 .
Fig. 8. -A, location of the REPSOL Casablanca oil rig in the Mediterranean coast of Catalonia; B, image of the north side showing the location of the instruments; C, heavy equipment was transported to Casablanca oil rig by ship, while delicate goods and personnel were flown by helicopter.

Fig. 9 .
Fig. 9. -Some instrumentation deployed during WISE 2000 and 2001: A, LAURA installed and acquiring data from the 32 m height deck; B, LAURA (L-band AUtomatic polarimetric Radiometer) during installation (UPC); C, Carolina Gabarró and Agustí Julià checking the radio link with the Aanderaa buoy; D, Aanderaa CMB3280 (meteorological buoy, ICM/CSIC); E, Prof. Font preparing the buoy designed at ICM with the CT recorder to sample near-surface salinity; F, underwater view of the buoy.

Fig
Fig. 10.-A, picture taken from 32-m-high deck of the Casablanca oil rig, waves destroyed the 7-m deck and seriously harmed the 13-m deck; B, vertical polarization (upper row), horizontal polarization (central row), and azimuth angle (bottom row) for two azimuthal scans during the storm (Camps et al. 2004).

Fig. 11
Fig. 11.-A, picture of the STARRS instrument mounted on the plane; B, flight acquisition plan; C, image of the plane flying over the Casablanca oil rig.

Fig. 13 .
Fig. 13.-Foam emissivity increase per foam thickness unit in mm at 100% coverage as a function of SSS and incidence angle: A, fit at horizontal polarization; B, fit at vertical polarization (modified from Camps et al. 2005).

Fig. 14
Fig.14.-A, instantaneous variations of the first Stokes parameter (nadir-looking antenna) and area under the normalized waveform (same philosophy as explained for the DDM volume, but in the experiment GOLD-RTR was not computing DDMs); B, retrieved sea surface salinity values using GNSS-R data from GOLD-RTR from IEEC -Institute for Space Sciences/CSIC to correct for sea state(Valencia et al. 2011).

Fig. 15
Fig. 15.-A, AMIRAS instrument on IRTA facilities for obtaining the first image of the Galaxy; B, first image of the Milky Way by AMIRAS (April 27-28, 2006).

Fig. 16
Fig. 16.-First snapshot images obtained with AMIRAS during an over-pass of an island in a lake in Finland.
The maps are based in the combination of OI SSS L3 and OSTIA SST.-TB Intercomparison between SMOS and Aquarius: Several studies on SMOS and Aquarius intercomparison have been performed for several Earth regions: Antarctica, the open ocean, dry and wet lands, the Amazon, etc. (Pablos et al. 2014).