This article summarizes some of the activities in which Jordi Font, research professor and head of the Department of Physical and Technological Oceanography, Institut de Ciències del Mar (CSIC, Spanish National Research Council) in Barcelona, has been involved as co-Principal Investigator for Ocean Salinity of the European Space Agency Soil Moisture and Ocean Salinity (SMOS) Earth Explorer Mission from the perspective of the Remote Sensing Lab at the Universitat Politècnica de Catalunya. We have probably left out some of his many contributions to salinity remote sensing, but we hope that this review will give an idea of the importance of his work. We focus on the following issues: 1) the new accurate measurements of the sea water dielectric constant, 2) the WISE and EuroSTARRS field experiments that helped to define the geophysical model function relating brightness temperature to sea state, 3) the FROG 2003 field experiment that helped to understand the emission of sea foam, 4) GNSS-R techniques for improving sea surface salinity retrieval, 5) instrument characterization campaigns, and 6) the operational implementation of the Processing Centre of Levels 3 and 4 at the SMOS Barcelona Expert Centre.
Este artículo resume algunas de las actividades en las que Jordi Font, profesor de investigación y jefe del Departamento de Física y Tecnología Oceanográfica, del Institut de Ciències del Mar (CSIC) en Barcelona, ha estado desarrollando como co-Investigador Principal de la parte de la misión SMOS de la ESA, una misión Earth Explorer, desde la perspectiva del Remote Sensing Lab, de la Universitat Politècnica de Catalunya. Seguramente, estamos olvidando algunas de sus muchas contribuciones a la teledetección de la salinidad, pero esperamos que esta revisión dé una idea de la importancia de su trabajo. Este artículo se focaliza en los siguientes puntos: 1) las medidas de alta calidad de la constante dieléctrica del agua marina, 2) las campañas de medidas WISE y EuroSTARRS que ayudaron a la definición del modelo geofísico relacionando la temperatura de brillo con el estado del mar, 3) la campaña de medidas FROG 2003 que ayudó a entender la emisión de la espuma marina 4) presentación de las técnicas de GNSS-R para la mejora de la recuperación de la salinidad superficial 5) campañas para la caracterización del instrumento y 6) la implantación del centro de procesado operacional de niveles 3 y 4 en el SMOS Barcelona Expert Centre.
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 (
The link between salinity and brightness temperature is through the dielectric constant of sea water. Though the “protected” 1.400-1.427 MHz band is the best choice, these measurements are far from trivial due to the low sensitivity of brightness temperature to SSS, 0.2-0.8 K/psu (psu = practical salinity unit, roughly 1 g of salt per 1 kg of water), the dependence on the physical temperature (
The first SSS airborne map was published in 1996 by Miller and Zaitzeff (
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 (
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
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
The main results are summarized in
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 (
In 1999 we proposed to ESA the possibility to carry out the WInd and Salinity Experiment (WISE) on the REPSOL Casablanca oil rig (40.72°N 1.36°E), ~50 km offshore of the coast of Tarragona (Spain), and ~115 km south of Barcelona (see
The instruments deployed were (
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 (
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 (
(1a) |
(1b) |
(2a) |
(2b) |
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 U10>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.15 K/(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.3 K for low-to-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 ).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 ).![]()
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 (
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
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
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 (
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.
Ultimately, the impact of sea foam can be parameterized as an emissivity increment, multiplied by the foam fraction
(3) |
, during WISE 2000, and , during WISE 2001. |
(4) |
Actually,
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 (
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 (
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 (
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
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 (
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
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 (
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 (
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.
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) (
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 is 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
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 (
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 ).Sample experimental SMOS Mediterranean SSS products produced at the CP34 [ http://cp34-bec.cmima.csic.es/experimental-mediterranean-sss-products-now-available-at-cp34-bec/ ].![]()
– 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 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 ).- 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 ,2015b ).
All these products and the research activities conducted at BEC during the last few years are broadly described in
For his dedication to the pioneer SMOS mission, Professor Font received the Catalan National Art Prize 2011 in the category of Scientific Thought and Culture (see
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 2000 and 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.
This work has been performed under research grants TEC2005-06863-C02-01/TCM, ESP2005-06823-C05 and ESP2007-65667-C04, AYA2008-05906-C02-01/ESP and AYA2010-22062-C05 and ESP2015-70014-C2-1-R, and EURYI 2004 award.