Comparison of normalized MSE between correlation powers obtained by serial and parallel processing.

## Abstract

In traditional Global Navigation Satellite System (GNSS) application, the reflected GNSS signals from Earth’s surface generally are considered as an interference source to be suppressed or removed. Recently, a new idea which treats the reflected GNSS signal as opportunity source of remote sensing has been proposed to monitor Earth’s physical parameters. This technique is called as GNSS-Reflectometry (GNSS-R) which has the advantages of low-power, -mass and -cost. With the development and modernization of GPS, Galileo, GLONASS, and BeiDou system, spaceborne GNSS could significantly improve the temporal-spatial resolution by receiving and processing the reflected signal from multiple satellites. This chapter mainly describes this new bi-static remote sensing technique. First, basic theories of GNSS-R including spatial geometry, polarization, and scattering model of reflected signal are discussed; second, spaceborne receivers and fast-response processing methods are reviewed and analyzed; finally, the empirical models retrieving wind speed are proposed and demonstrated using the DDM data from the UK-TechDomeSat-1 satellite. Based on the discussion of this chapter, it could be concluded that although GNSS-R still has some key challenges which have to be addressed, it could be an optimal choice of remote sensing in some special conditions, such as the tropical cyclone.

### Keywords

- GNSS-Reflectometry
- remote sensing
- delay-Doppler maps
- wind speed

## 1. Introduction

Sea wind speed is an important parameter which impacts the ocean circulation and global climate. Moreover, tropical cyclones, as one of the most serious natural disasters, interrupt infrastructure and endanger life safety [1]. For these reasons, it is important to monitor sea wind speed to study and forecast some complex weather. Traditional tools, such as buoys, could provide long-term observation, however, have low spatial sampling. The spaceborne microwave remote sensing devices such as radar altimeter [2], scatterometer [3], and radiometer [4], provide all weather and all day observation of global wind speed. However, high-mass, -power, and -cost equipment limit their application in commercial remote sensing based on the single-function microsatellite. GNSS not only provides the positioning, velocity, and timing (PVT), but also could be considered as the opportunity sources of remote sensing to observe Earth physical parameters. This observation is known as GNSS-Reflectometry which uses GNSS signal reflected off Earth’s surface to observe Earth’s physical parameters. It has been demonstrated that GNSS-R could have been utilized to measure sea wind speed [5], sea surface height [6], sea ice [7], and soil moisture [8] from ground-based, airborne, and spaceborne platforms in past 20 years. Compared to traditional remote sensing ways, GNSS-R has the advantages of the low-cost and low-power, because only a receiver is needed. Spaceborne GNSS-R could improve the global temporal-spatial resolution because there are more than 100 GNSS satellites in orbit or planned. In addition, the temporal-spatial resolution could be further improved by microsatellite network.

The idea of using reflected GNSS signal from Earth’s surface for remote sensing was discussed by Hall and Cordey in 1988 [9]. In 1993, Martin-Neria proposed the concept of PAssive Reflectometry and Interferometry System (PARIS) to use GNSS signal reflected off the ocean for the measurement of ocean altimetry [10]. In 1994, Auber et al. detected GNSS signal reflected off the ocean from an aircraft platform [11]. In 1998, Garrison and Katzberg demonstrated that reflected global positioning system (GPS) signal could measure sea wind speed through an aircraft experiment [12]. The first collection of reflected GPS signal from spaceborne platform was reported in 2002 [13]. The use of reflected GPS signal acquired by UK-DMC satellite is recognized as the first spaceborne experiment to measure wind speed [14]. After the success of the UK-DMC, UK TechDemoSat-1 [15], PARIS-IoD [16], GEROS-ISS [17], and CYGNSS [18] were planned and developed to further demonstrate the feasibility of using reflected GNSS signal to observe Earth parameters. Especially, CYGNSS mission aims to improve intensity forecasting by sensing sea wind speed in the inner core of tropical cyclones using a constellation of eight microsatellites.

The basic product of GNSS-R is the so-called Delay Waveform or delay-Doppler map (DDM) of reflected GNSS signal, which describes the distribution of power level in delay or delay-Doppler domain. Based on the bi-static radar equation and the geometric optics limit of the Kirchhoff approximation, Zavorotny and Voronovich developed the theoretical power model of the scattered GNSS signal as a function of the geometrical and environmental parameters [19]. To simulate DDM, the Z-V model can be implemented by defining the scenario in a reference system and evaluating the functions inside the integrand to compute the integral for every delay and Doppler bin. However, such approach is extremely time- and resource-consuming, especially for spaceborne scenario. In [20], a new and efficient method to compute DDM was proposed based on the derivation of explicit expressions of the space coordinates as a function of the delay offset and Doppler shift. The realistic problem that the noisy DDM measurements should be generated motivated the model and the simulator to produce DDM time series based on a multiple-input-multiple-output filter bank [21]. Other approaches of simulating DDM generate a synthetic sea surface with accurate temporal-spatial correlations and model the time evolution of the scattered signal. In [22], Clarizia et al. presented a new facet-based approach to model the scattering of GNSS signal from the synthetic sea surface and used it to compute DDM for the spaceborne scenario. In practice, the DDM is generated by correlating the reflected GNSS signal with the locally generated replicas at different delay offset and Doppler shift. In [23], a delay/Doppler-mapping hardware/software receiver was presented. A real-time instrument called as GPS open-loop differential real-time receiver (GOLD-RTR) was designed and developed to gather the complex correlations of reflected GPS signal in [24]. To conform the development of spaceborne GNSS-R mission, spaceborne GNSS-R receivers have been designed and developed, such as Space GNSS Receiver Remote Sensing Instrument (SGR-ReSI) [25] and GNSS-REflectometry, Radio Occultation, and Scatterometry (GEROS) [26] for UK-TDS-1 satellite and International Space Station.

The purpose of remote sensing is to determine the Earth’s physical parameters, for which it is the key technology to develop the link between GNSS-R observable and the Earth’s physical parameters. Some algorithms have been proposed to retrieve sea wind speed using GNSS-R. The first type of method is to fit the measured Delay Waveform or DDM to theoretical one. In [27], the trailing edge slope and the complete shape of the Delay Waveform were used to fit the analytical models to estimate wind speed. A 2-D least-squares fitting approach was given to retrieve wind field by fitting the simulated DDM to space-based measured one from UK-DMC satellite [28]. These methods based on fitting are time-and resource-consuming. The second type of method is to directly link wind speed and different Delay Waveform or DDM observable by the regression. In 2013, Rodrigues-Alvarez et al. reported their airborne experiment results in which four DDM observable including weighted area, 2-norm Euclidean distance, distance from the geometric center to the maximum value of DDM, and 1-norm Euclidean distance were defined and linked with the wind speed [29]. First results of the retrieved wind speed using the data from UK TechDemoSat-1 satellite were presented in [30], where in-situ wind from Advanced Scatterometer (ASCAT) and the bi-static scattering coefficient derived from peak DDM using bi-static radar equation are utilized to develop a wind speed retrieval algorithm and the accuracy of 2.2 m/s could be obtained in range from 3 to 18 m/s. To improve retrieval accuracy of high wind speed and spatial resolution for CYGNSS mission, Rodriguez-Alvarez defined the generalized linear DDM observable based on maximum signal-to-noise ratio (SNR), minimum variance of wind speed, and principal component analysis (PCA) to develop relationships with wind speed and found that the observable based on PCA had the best performance [31]. The retrieval algorithm proposed by M. P. Clarizia in 2014 was used to produce the Level 2 wind speed data product of CYGNSS mission in [32], in which several additional processing steps were added to correct the influence of incident angle and adaptively select the parameters of the algorithm to improve the retrieving accuracy [32]. The algorithms for CYGNSS mission have limited the spatial resolution in the region of 25 km × 25 km. These retrieval algorithms above all provide a single wind speed descriptor in the overall observation region. The third type of method is based on the mapping between the physical space and the delay-Doppler domain. In [33], an algorithm to retrieve the bi-static scattering coefficient distribution over the observation region from measured DDM using the deconvolution was proposed. D. Schiavulli demonstrated the validity and feasibility of this approach using UK-TechDemoSat-1 data [34].

This chapter will illustrate the GNSS application in retrieving wind speed using spaceborne GNSS-R. Section 2 will discuss the basic theories of GNSS-R including the observation geometry, the polarization, and the model of reflected GNSS signal. Section 3 will present some typical spaceborne GNSS-R receivers and processing methods of reflected GNSS signal. The explicit retrieval algorithms based on the regression will be developed in Section 4. Finally, the conclusion of this chapter will be addressed in Section 5.

## 2. Basic theory of GNSS-R

### 2.1. Geometry

GNSS-R works as a bi-static radar in which the transmitter and the receiver are separated in the spatial distribution. It is noted that when GNSS-R receivers allow to track the scattered signal from different GNSS satellites, the system is called as multi-static radar as shown in Figure 1, in which for each GNSS satellite, bi-static observation occurs at different azimuth angle, and creates a glistening zone. To simplify observation geometry of GNSS-R, a bi-static geometry could be developed as Figure 2, in which (1) the coordinate origin is at Earth’s core; (2) the YOZ plane is in the incident plane of GNSS signal; (3) Z axis has the same direction with the normal of the tangent plane of the specular point; and (4) it is assumed that Earth, the orbits of GNSS and Low Earth Orbit (LEO) satellite are circular. In Figure 2,

where

The coordinate of scattering unit

where

### 2.2. Polarization of reflected GNSS signal

When Right-Hand Circular Polarization (RHCP) GNSS signal occurs reflection phenomenon, the polarization is changed. The Fresnel coefficients of RHCP and LHCP (Left-Hand Circular Polarization) are expressed as the combination of the vertical and horizontal polarization as

The Fresnel coefficients of vertical and horizontal polarization are computed as

Figure 3 gives the simulated Fresnel coefficients of RHCP and LHCP for GPS L1 signal, when sea surface temperature is 25° and sea surface salinity is 35 psu. From the figure, it could be seen that as the incident angle increases, the Fresnel coefficients of RHCP and LHCP show decreasing and increasing trends, respectively. This illustrates that at nadir and small incident angle, the LHCP signal dominates GNSS signal reflected off the sea surface. Therefore, the combination of LHCP and RHCP antenna or only LHCP antenna is used to receive reflected GNSS signal from the sea surface.

### 2.3. Model of reflected GNSS signal

The signal reflected from the sea surface is described as [21]

where

where

where

where

where

where

where

## 3. Processing of GNSS reflected signal

The GNSS-R receiver is a processing unit to cross-correlate reflected GNSS signal with locally generated replicas. Spaceborne GNSS-R receiver not only requires the capacity of processing reflected GNSS signal in real-time, but also with low-mass, -power and -cost. At present, some research groups have developed spaceborne GNSS receivers to meet their spaceborne missions.

### 3.1. Existing receivers

#### 3.1.1. UK-DMC

To promote the development of spaceborne GNSS-R and explore its prospect of commercial application, Surrey Satellite Technology Ltd. (SSTL) designed and developed spaceborne GPS-R receiver which would be carried on UK-DMC. This receiver which was implemented based on spaceborne GPS receiver of SSTL consisted of a zenith antenna to receive direct GPS signal, a LCHP antenna with the gain of 11.8 dBi to downward point to receive reflected GPS signal, a real-time processor unit to perform delay-Doppler mapping of reflected GPS signal as well as log and download raw sampled IF data to the ground for more intensive postprocessing together with a solid state data recorder (SSDR) and high rate downlink (HRD) [14]. Moreover, the processor unit supported to track the direct signal provides the positioning, velocity and timing (PVT). When the receiver was configured as the collection mode, 20 s raw sampled data of the direct and reflected GPS signal could be saved in SSDR to download to the ground.

#### 3.1.2. SGR-ReSI

Following the success of the UK-DMC, SSTL embarked on the new generation spaceborne GNSS-R receiver Space Receiver Remote Sensing Instrument (SGR-ReSI) which had been equipped with UK-TechDemoSat-1 satellite launched in July, 2014 [25]. Compared to the receiver carried on UK-DMC, SGR-ReSI supported to process in real-time and save sampled IF data of multi-frequency GNSS signal. The core component of the receiver is a 24-channel navigation receiver of GPS L1 signal and reprogrammable to a co-processor unit based on Actel ProASIC3 FPGA. The reprogrammable capacity of SGR-ReSI also enables more scientific and remote sensing application. Two types of RF front-ends including Max2769 optimized for GPS L1 signal and Max2112 which is reconfigured to other GNSS bands have been adopted in the receiver. In operation, the receiver could be configured as fixed gain mode (FGM) and unmonitored automatic gain control (uAGC). To allow the storage of both sampled and processed data, a bank of DDR2 memory with a capacity of 1 GB was used. From September 2014 on, the GNSS-R basic observable DDM data formatted as TIFF were downloaded from UK-TechDemoSat-1 to the ground.

#### 3.1.3. GEROS

The main objectives of GEROS mission are to measure the sea surface height and mean square slope, and the secondary ones are to further explore the potential of monitoring vertical profiles of atmospheric using GNSS-RO and to assess the feasibility of GNSS-R for land application [17]. The receiving antenna has

#### 3.1.4. PAU/GNSS-R

Passive Advanced Unit (PAU) is a new instrument concept which aims to monitor ocean parameters. It consists of the PAU/RAD which is an L-band radiometer to measure the brightness temperature of the sea surface, the RAU/GNSS-R which is a reflectometer to measure the roughness of the sea surface and the PAU/IR which is two infrared radiometers used to observe the temperature of the sea surface. PAU/GNSS-R was designed based on FPGA to synchronously process the reflected GPS signal from different satellites in real-time and output corresponding to 2-D DDM [35]. To reduce the hardware resources of the payload, hardware reuse technique was adopted based on two RAM-like registers that change their respective input and output connections to allow the DDM generator to be implemented at a higher clock rate. A software running on MicroBlaze soft processor mainly performs functions including sending and receiving data from the GPS receiver, selecting the available satellites, and computing the DDM parameters transferred to the DDM generator. The computed DDM are transferred to a terminal computer using an USB and an external master based on FPGA.

#### 3.1.5. PYCARO

According to the research and demonstration mission planned by the Remote Sensing Lab and the NanoSat Lab at the Universitat Politècnica de Catalunya-Barcelona Tech, PYCARO will be carried on a 3

### 3.2. General architecture of receiver

According to the description on the spaceborne GNSS-R receivers above, the general architecture of the receiver could be given as Figure 6, which consists of the antennas, the RF front-end, the processor of baseband signal, the computer on-board, and some external components. The antennas include zenith- and nadir-pointing antennas to receive direct and reflected GNSS signal. Usually, the zenith-pointing antenna is an omnidirectional RHCP GNSS antenna, and nadir-pointing antenna is a narrow-beam and high-gain LHCP antenna. However, it is able to be design the different combination of the zenith- and nadir-pointing antennas to perform special spaceborne mission and application. For example, the zenith- and nadir-pointing antennas of PARIS are both the narrow-beam and high-gain antennas to make the incoming signal from the single satellite be received; to significantly measure the soil moisture of the land, the nadir-pointing antenna of GEROS is the combination of the LHCP and RHCP antenna. RF front-ends convert the analog RF signal to the sampled IF one, consisting of Low Noise Amplifier (LNA), bandpass filter, mixer, gain control, and A/D converter. According to the requirement of the mission, the gain control could have different configuration. For altimetry, the Automatic Gain Control (AGC) mode is utilized to dynamically adjust the power level of incoming signal to the most effective input range of A/D converter; however, the scatterometry of GNSS-R generally selects the fixed gain or monitoring AGC to correctly receiving gain to obtain accuracy power level of reflected signal. To improve the generality of the receiver, the gain control could be designed as reconfigurable mode, such as SGR-ReSI of which the gain control supports the fixed gain and unmonitoring AGC mode. The core of the receiver is a baseband processor which consists of controller and processor unit. The controller unit controls the acquirement and tracking of direct signal, estimates the difference of delay and Doppler frequency between direct and reflected signal, and configures the DDM parameters, such as delay and Doppler frequency window range and resolution and others. The estimations of the differences of delay and Doppler frequency between direct and reflected signal are given as

where

### 3.3. Processing method

At present, the processing methods of reflected GNSS signal have serial and parallel correlation. The serial correlation correlates the reflected GNSS signal with the locally generated replicas at the different delay and Doppler frequency one by one in time domain. The parallel correlation first transforms the reflected and local signal to the frequency domain using fast Fourier transform (FFT) and performs a multiplication operation, then transforms multiple results to time domain using Inverse FFT (IFFT). The serial processing is lossless and can be considered as the reference for the assessments of parallel method. The assessment metric is defined as normalized mean squared error (MSE) between serial and parallel correlation as

where

#### 3.3.1. Serial correlation

Serial correlation in time domain is expressed as

where

#### 3.3.2. Parallel correlation in delay domain

Parallel correlation in delay domain parallelly performs correlation operation for all code delay based on Fourier transform in each Doppler channel as

where

#### 3.3.3. Parallel correlation in Doppler domain

Parallel correlation in Doppler domain performs spectrum estimation for the signal removed the modulation of the code at each delay channel as

The block diagram is illustrated as in Figure 9, in which the number of processing channel is

#### 3.3.4. Parallel correlation in delay-Doppler domain

Parallel correlation in delay-Doppler domain [38], as given in Figure 10, is a block processing technique called as double-block zero-padding (DBZP) expressed as

where * i*th blocks of the reflected signal and local replica;

#### 3.3.5. Comparison of computational complexity

Although the computational complexities of the correlation processing have been analyzed above, here, the simulated comparison is conducted. It is assumed that the range and the resolution of Doppler frequency are [−8, 8] KHz and 500 Hz, respectively, and the numbers of delay bin and the DBZP block are same. Figure 11 gives the simulated computational complexities when the numbers of the sampling are 2048, 4096, 8192, and 16,384, from which it is seen that (1) serial correlation needs most multiplier and adder and consumes most hardware resources so that it is unsuitable for spaceborne receiver; (2) parallel correlation methods significantly reduce the number of required multiplier and adder, especially parallel correlation in Doppler domain. It should be noted FFT and IFFT on hardware is a 2-based operation which needs to padding zero to make the number of operating point to be the power of 2. Large number of padding zero causes the increasing of computational complexity; therefore, it has to be considered to choose optimal DDM parameters.

### 3.4. Processing of UK-TDS/UK-TDS-1 data

The raw IF data from UK-DMC and UK-TechDemoSat-1 satellite are used to demonstrate the correlation approaches above. UK-DMC satellite which first carried an experimental GNSS reflectometry aimed to monitor the disaster using optical imaging. In 2003, UK-DMC satellite was launched into the sun-synchronous orbit with 680 km altitude. Some researchers have used GNSS-R data from UK-DMC satellite to demonstrate the feasibility of remotely sensing sea wind speed, sea ice, and others. Following the end of UK-DMC’s operational life in 2011, SSTL launched UK-TDS-1 satellite equipped with new generation spaceborne GNSS-R receiver SGR-ReSI in 2014. The orbit of UK-TDS-1 satellite is also sun-synchronous orbit with the altitude and the inclination of 635 km and 98.4°, respectively. UK-TDS-1 satellite is controlled and operated with a 8-day duty circle shared between 8 payloads. SGR-ReSI can be operated for only 2 days in every 8-day period. Both UK-DMC and UK-TDS-1 satellite support the collection of raw IF data, however, because of the limitation of the storage capacity, only a few raw IF data are available. Here, the data collected by UK-DMC on November 11, 2004 and UK-TDS-1 on January 27, 2015 are used. The corresponding ground track of GNSS, UK-DMC, and UK-TDS-1 satellite, and the positioning of the specular pointing are shown in Figure 12, in which both specular pointings lie to the ocean surface, i.e., the signals collected by UK-DMC and UK-TDS-1 satellite are both reflected from the ocean surface. Figure 13 shows the DDM distribution of reflected GNSS signal collected by UK-DMC and UK-TDS-1 satellite, in which from left to right, the DDMs are computed using serial correlation, parallel correlation in delay domain, parallel correlation in Doppler domain, and parallel correlation in delay-Doppler domain. From Figure 13, it could be seen that (1) the DDMs from UK-DMC and UK-TDS-1 satellites both show the distribution of the horseshoe shape as simulated in Figure 4 of Section 2.3; (2) the DDMs computed by four correlation methods have similar distribution in delay-Doppler domain. The assessment metrics defined by expression (19) are computed for the DDMs in Figure 13. The results are given in Table 1, which indicate that compared to the serial correlation, parallel correlations have negligible MSEs, and in addition, one-dimensional parallel methods have lower MSEs than two-dimensional ones. The results above illustrate that parallel correlations, especially the parallel correlation in Doppler domain, are better chosen to compute the DDM than serial one for spaceborne GNSS-R receiver.

Satellite | Serial | Parallel ( | ||
---|---|---|---|---|

Delay domain | Doppler domain | Delay-Doppler domain | ||

UK-DMC | 0 | 5.81 | 2.62 | 6.94 |

UK-TDS-1 | 0 | 0.09 | 0.32 | 1.17 |

## 4. Retrieving of wind speed

Retrieving wind speed is an inversion problem of remote sensing, of which the key is to find and extract the sensitivity parameters to wind speed from the basic observables of the remote sensing sensors. In fact, L-band signals are less optimal frequency for the measurement of wind speed compared to scatterometers, however, are not significantly attenuated by rainfall so that they are available to observe wind speed in tropical cyclones. As mentioned in Section 1, at present, three types of retrieving approaches for spaceborne GNSS-R have been proposed. Here, the approaches based on the direct mapping from DDM metrics to wind speed are discussed.

### 4.1. Retrieving method

An expression for the distribution of DDM in delay-Doppler domain has been given in (15), in which the parameter related to sea wind speed is bi-static scattering coefficient

The expression computing bi-static scattering coefficient is derived as

In expression (26),

where

As illustrated in Figure 5, the DDM not only gives the distribution of correlation power in delay-Doppler domain, but also presents the spatial distribution of power. The given delay-Doppler interval in expression (28) determines the geometric resolution of spaceborne GNSS-R as [40]

where

### 4.2. Demonstration of UK-TDS-1 and ASCAT data

The gain modes of RF front-end in SGR-ReSI have unmonitoring AGC and fixed gain. After June 2015, SGR-ReSI was configured as the fixed gain of 36 dB. The UK-TDS-1 Level 1 data collected from November 26, 2015 to April 19, 2016 are used to demonstrate the feasibility of retrieving wind speed using the approach above. Wind speed monitored by ASCAT-A/B which have a resolution of 25 km are selected as in-situ data. The matchup between the specular pointing of UK-TDS-1 and the ground track of ASCAT-A/B are set to allow a maximum space and time separation of 0.25° and 1 h. To improve the accuracy of retrieved wind speed, it is necessary to select high-quality L1 data to retrieve wind speed. The first is spatial selection that the latitude range of the specular point is limited from −50°S to −50°N to remove the data reflected off the sea ice, and the distances between the specular pointing and the coastline are controlled to be over 0.5° to remove the impact of the land on the DDM. The second one is the selection based on the antenna gain at the specular point. The nadir-looking antenna gains of SGR-ReSI rapidly attenuate when the incident angles of incoming signal are over 40°. Therefore, the incident angles of GNSS signal are limited within the range of 30°. The reflected GNSS signal arriving at the nadir-pointing antenna is exceedingly weak. The coherent integration of 1 ms and incoherent averaging of 1000 significantly improve the signal-to-noise ratio (SNR); however, some DDMs still are buried in the noise floor. In addition, the SNR is one of the factors influencing on the accuracy of wind speed. Therefore, it is needed to remove DDMs with the low SNR. The SNR of DDM could be defined as

where

After the DDMs are selected, the delay and Doppler interval in expression (28) should be determined. The delay and Doppler interval are limited by the geometric resolution. Figure 14 shows the relationship between the geometric resolution and the incident angle for different delay interval when Doppler interval is [−1, 1] KHz, from which it is seen that when delay interval is [−0.25, 0.25] chips, the geometric resolution is lower than 25 km which is required resolution of the traditional spaceborne remote sensing. Based on the discussion, the delay and Doppler interval in expression (28) are configured as [−0.25, 0.25] chips and [−1, 1] KHz.

The retrieval model is the mapping linking wind speed with the bi-static scattering coefficient computed by expression (26) or (27). Note that because of the unknown transmit power and antenna gain, the influence of

where

Observable | Train set | Test set | |||
---|---|---|---|---|---|

A | B | C | Bias (m/s) | RMSE (m/s) | |

Expression (26) | −0.237 | −0.0115 | 0.024 | 2.11 | |

Expression (27) | −0.2177 | −0.0665 | 0.021 | 2.15 |

The SNR of DDM has an important impact on retrieving result of wind speed as shown in Figure 18. The accuracy of retrieved wind speed is improved by selecting the DDMs with higher SNR. However, it should be noted that the larger the threshold of the SNR is, the lesser the DDMs used to the retrieved wind speed are, i.e., the coverage of the UK-TDS-1 is lower for the higher threshold of the SNR. It has to be noted that no precise calibration and correction for the DDM measurement in SGR-ReSI also is one of the reasons causing the low performance of retrieving wind speed.

## 5. Conclusion

GNSS-R has been a new tool of remote sensing to monitor Earth’s physical parameters, such as sea wind speed, sea surface height, sea ice, soil moisture and others. Following the success of UK-DMC satellite which was carried on the first GNSS-R receiver, the spaceborne GNSS-R which could provide the global and all-days observation has been an hot focus of the research groups. Some spaceborne missions have been proposed and planned, such as CYGNSS of NASA, GEROS-ISS of ESA, ^{3}CAT-2 of the Remote Sensing Lab and the NanoSat Lab at the universitat Politècnica de Catalunya-Barcelona Tech. This chapter was focused on the discussion of spaceborne GNSS-R to monitor sea wind speed. First, the basic theories of spaceborne GNSS-R, including the geometry, polarization, and scattering model were discussed. The geometry of GNSS-R mainly is determined by the incident angle of GNSS signal, and the height of LEO and GNSS satellite. When GNSS signal occurs reflection on the sea surface, the polarization of GNSS signal is changed. At nadir and small incident angle, the LCHP component dominates the reflected GNSS signal. For the diffuse scattering of the GNSS signal, the correlation power level called as DDM spreads to large delay and Doppler frequency and presents the horseshoe shape in delay-Doppler domain. Then, the spaceborne receivers and processing methods were discussed. The results showed that parallel correlation in Doppler domain has lowest computational complexity and optimal performance so that it is optimal to choose to process the reflected GNSS signal in real-time for spaceborne receiver. Finally, the approach based on DDM metric was used to retrieve wind speed. Proposed DDM metrics have inversely proportional relationship with wind speed. By latitude mask and the threshold of SNR to remove low-quality data, proposed approaches could obtain the RMSE of 2.11 m/s and 2.15 m/s. For further improvement of retrieving accuracy, it is necessary to set higher SNR threshold which results in lower spatial coverage. In one word, GNSS-R will become an effective way to observe global sea wind speed, especially low attenuation of GNSS signal for the rainfall will be an optimal choice to monitor tropical cyclone.