The global navigation satellite system-reflectometry (GNSS-R) field has experienced an exponential growth as it is becoming relevant to many applications and has captivated the attention of an elevated number of research scholars, research centers and companies around the world. Primarily based on the contents of two Special Issues dedicated to the applications of GNSS-R to Earth observation, this review article provides an overview of the latest advances in the GNSS-R field. Studies are reviewed from four perspectives: (1) technology advancements, (2) ocean applications, (3) the emergent land applications, and (4) new science investigations. The technology involved in the GNSS-R design has evolved from its initial GPS L1 LHCP topology to include the use of other GNSS bands (L2, L5, Galileo, etc.), as well as consider RHCP/LHCP-receiving polarizations in order to perform polarimetric studies. Ocean applications have included developments towards ocean wind speed retrievals, swell and altimetry. Land applications have evolved considerably in the past few years; studies have used GNSS-R for soil moisture, vegetation opacity, and wetland detection and monitoring. They have also determined flood inundation, snow height, and sea ice concentration and extent. Additionally, other applications have emerged in recent years as we have gained more understanding of the capabilities of GNSS-R.

Global navigation satellite system reflectometry is a multistatic radar using navigation signals as signals of opportunity. It provides wide-swath and improved spatiotemporal sampling over current space-borne missions. The lack of experimental datasets from space covering signals from multiple constellations (GPS, GLONASS, Galileo, and Beidou) at dual-band (L1 and L2) and dual-polarization (right- and left-hand circular polarization), over the ocean, land, and cryosphere remains a bottleneck to further develop these techniques. ³ Cat-2 is a 6-unit (3 × 2 elementary blocks of 10 × 10 × 10 cm ³ ) CubeSat mission designed and implemented at the Universitat Politècnica de Catalunya-BarcelonaTech to explore fundamental issues toward an improvement in the understanding of the bistatic scattering properties of different targets. Since geolocalization of the specific reflection points is determined by the geometry only, a moderate pointing accuracy is only required to correct the antenna pattern in scatterometry measurements. This paper describes the mission analysis and the current status of the assembly, integration, and verification activities of both the engineering model and the flight model performed at Universitat Politècnica de Catalunya NanoSatLab premises. 3Cat-2 launch is foreseen for the second quarter of 2016 into a Sun-Synchronous orbit of 510-km height.

The FSSCat mission was the 2017 ESA Sentinel Small Satellite (S⌃3) Challenge winner and the Copernicus Masters competition overall winner. It was successfully launched on 3 September 2020 onboard the VEGA SSMS PoC (VV16). FSSCat aims to provide coarse and downscaled soil moisture data and over polar regions, sea ice cover, and coarse resolution ice thickness using a combined L-band microwave radiometer and GNSS-Reflectometry payload. As part of the calibration and validation activities of FSSCat, a GNSS-R instrument was deployed as part of the MOSAiC polar expedition. The Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition was an international one-year-long field experiment led by the Alfred Wegener Institute to study the climate system and the impact of climate change in the Arctic Ocean. This paper presents the first results of the PYCARO-2 instrument, focused on the GNSS-R techniques used to measure snow and ice thickness of an ice floe. The Interference Pattern produced by the combination of the GNSS direct and reflected signals over the sea-ice has been modeled using a four-layer model. The different thicknesses of the substrate layers (i.e., snow and ice) are linked to the position of the fringes of the interference pattern. Data collected by MOSAiC GNSS-R instrument between December 2019 and January 2020 for different GNSS constellations and frequencies are presented and analyzed, showing that under general conditions, sea ice and snow thickness can be retrieved using multiangular and multifrequency data.

FSSCAT is an innovative mission concept consisting of two federated 6U Cubesats in support of the Copernicus Land and Marine Environment services. The first CubeSat carries the Flexible Microwave Payload-2 (FMPL-2), a dual microwave payload (a GNSS-Reflectometer and an L-band radiometer with interference detection/mitigation), and the second one a hyper-spectral optical payload. The combined use of these payloads will allow to measure soil moisture, ice extent, and ice thickness, and to detect melting ponds over ice. FSSCAT also includes radio and optical inter-satellite links to test some of the techniques and technologies for the upcoming satellite federations. FSSCAT will be the precursor of a constellation of federated small satellites for Earth observation achieving high temporal resolution and moderate spatial resolution in a cost-effective manner.

The first-ever dual-frequency multi-constellation Global Navigation Satellite Systems Reflectometry (GNSS-R) polarimetric measurements over boreal forests and lakes from the stratosphere are presented. Data were collected during the European Space Agency (ESA) sponsored Balloon Experiments for University Students (BEXUS) 19 stratospheric balloon experiment using the P(Y) and C/A Reflect Ometer (PYCARO) instrument operated in closed-loop mode. Maps of the polarimetric ratio for L1 and L2 Global Positioning System (GPS) and GLObal Navigation Satellite System (GLONASS), and for E1 Galileo signals are derived from the float phase at 27,000 m height, and the specular points are geolocalized on the Earth’s surface. Polarimetric ratio ( ) maps over boreal forests are shown to be in the range 2–16 dB for the different GNSS codes. This result suggests that the scattering is taking place not only over the soil, but over the different forests elements as well. Additionally to the interpretation of the experimental results a theoretical investigation of the different contributions to the total reflectivity over boreal forests is performed using a bistatic scattering model. The simulated cross- (reflected Left Hand Circular Polarization LHCP) and co-polar (reflected Right Hand Circular Polarization RHCP) reflectivities are evaluated for the soil, the canopy, and the canopy–soil interactions for three different biomass densities: 725 trees/ha, 150 trees/ha and 72 trees/ha. For elevation angles larger than the Brewster angle, it is found that the cross-polar signal is dominant when just single reflections over the forests are evaluated, while in the case of multiple reflections the co-polar signal becomes the largest one. The first-ever dual-frequency multi-constellation Global Navigation Satellite Systems Reflectometry (GNSS-R) polarimetric measurements over boreal forests and lakes from the stratosphere are presented. Data were collected during the European Space Agency (ESA) sponsored Balloon Experiments for University Students (BEXUS) 19 stratospheric balloon experiment using the P(Y) and C/A Reflect Ometer (PYCARO) instrument operated in closed-loop mode. Maps of the polarimetric ratio for L1 and L2 Global Positioning System (GPS) and GLObal Navigation Satellite System (GLONASS), and for E1 Galileo signals are derived from the float phase at 27,000 m height, and the specular points are geolocalized on the Earth’s surface. Polarimetric ratio ( ) maps over boreal forests are shown to be in the range 2–16 dB for the different GNSS codes. This result suggests that the scattering is taking place not only over the soil, but over the different forests elements as well. Additionally to the interpretation of the experimental results a theoretical investigation of the different contributions to the total reflectivity over boreal forests is performed using a bistatic scattering model. The simulated cross- (reflected Left Hand Circular Polarization LHCP) and co-polar (reflected Right Hand Circular Polarization RHCP) reflectivities are evaluated for the soil, the canopy, and the canopy–soil interactions for three different biomass densities: 725 trees/ha, 150 trees/ha and 72 trees/ha. For elevation angles larger than the Brewster angle, it is found that the cross-polar signal is dominant when just single reflections over the forests are evaluated, while in the case of multiple reflections the co-polar signal becomes the largest one.

The Flexible Microwave Payload-2 is the GNSS-R and L-band Microwave Radiometer Payload on board ${}^3$Cat-5/A, one of the two 6-unit CubeSats of the FSSCat mission, which were successfully launched on 3 September 2020 on Vega flight VV16. The instrument occupies nearly a single unit of the CubeSat, and its goal is to provide sea-ice extension and thickness over the poles, and soil moisture maps at low-moderate resolution over land, which will be downscaled using data from Cosine Hyperscout-2 on board ${}^3$Cat-5/B. The spacecrafts are in a 97.5${}^{°}$ inclination Sun-synchronous orbit, and both the reflectometer and the radiometer have been successfully executed and validated over both the North and the South poles. This manuscript presents the results and validation of the first data sets collected by the instrument during the first two months of the mission. The results of the validation are showing a radiometric accuracy better than 2 K, and a sensitivity lower than the Kelvin. For the reflectometer, the results are showing that the sea-ice transition can be estimated even at short integration times (40 ms). The presented results shows the potential for Earth Observation missions based on CubeSats, which temporal and spatial resolution can be further increased by means of CubeSat constellations.