Open access peer-reviewed chapter

Remote Sensing in Land Applications by Using GNSS-Reflectometry

By Yan Jia and Yuekun Pei

Submitted: April 27th 2017Reviewed: December 4th 2017Published: July 25th 2018

DOI: 10.5772/intechopen.72901

Downloaded: 349

Abstract

Global navigation satellite system-reflectometry (GNSS-R) as an efficient tool for remote sensing has gained increasing interests in the last two decades, due to its unique characteristics. It uses GNSS signals as sources of opportunity, providing precise, continuous, all-weather, and 24 hours’ detections, which play a key role in many land applications. The fundamental theoretical part of GNSS-R technique is examined at first. Then, GNSS-R methodologies applied in the soil moisture content, vegetation biomass sensing, and altimetry applications are also detailed. One retrieval method uses only RH (right-hand) reflected data. Another retrieval method for soil moisture content (SMC) aimed to calibrate the measurement by using water reflections, based on the bistatic equations with LH (left-hand) reflected and RH direct signals. The other method for SMC retrieval is related to the polarimetric ratio (PR), the ratio of LH/RH reflected signals can reveal the fluctuations of the SMC. Another vital parameter vegetation biomass was observed by using the variation of reflectivity of the LH and RH reflected components. Finally, the C/A code method was used for exploring the possibility to the altimetry estimation. The features of GNSS-R technique made it a promising remote sensing technique in hydrology, climatology carbon cycles, and other potential applications.

Keywords

  • GNSS-reflectometry (GNSS-R)
  • soil moisture content
  • vegetation biomass
  • altimetry
  • polarimetric measurement

1. Introduction to GNSS

Global navigation satellite systems (GNSSs) are satellite constellations orbiting the Earth to provide continuous positioning and timing information for users on the globe. It mainly includes the United States’ NAVSTAR global positioning system (GPS), European Union’s Galileo, Russia’s Glonass, China’s BeiDou, and Japan’s Quasi-Zenith satellite system (QZSS). For now, such systems are used not only for navigation positioning but also for object detection in various disciplines, such as archeology (a 3-dimension map is generated using the combined signals of GPS, Glonass, and QZSS to detect each artifact), air traffic (the surveillance-broadcast system of air traffic control obtains aircraft positions by using GNSS signals), weather (forecast can be made through the atmospheric bending measurement of GNSS satellite signals), and others. The huge potential of GNSSs has led international governments and researchers made efforts on it [1]. Up to now, the United States’ NAVSTAR GPS and the Russian GLONASS have begun providing global coverage for a long time. In December 2016, the European’s Galileo formally announced to start its initial service to provide its full operational capability. China is in the process of expanding its regional BeiDou Navigation Satellite System into global compass navigation system in the very near 2020. India and Japan are also in the process of developing their regional navigation systems [2].

1.1. GPS system

GPS is the United States’ global positioning system (GPS) developed since 1973. The number of satellites, the positions, and the relative inclinations was optimized to ensure an adequate global coverage. In 1978, the first operational prototype satellite was launched. In 1993, the initial capability and the full GPS system were declared [3]. As a universal positioning system, GPS provides an accurate, continuous, three-dimensional position and velocity information to the user with the appropriate receiver on Earth. Now, it is the most used navigation system in the world and fully operative both for military purpose and civilian aims [4].

The GPS signals are transmitted with right-hand circular polarization (RHCP) so that the ionosphere attenuation does not affect the signal reception. Each GPS satellite transmits microwave carrier signals at three bands: L1 (1575.42 MHz) carries the navigation message and the standard positioning services (SPS), L2 (1227.6 MHz) is used to measure the ionosphere delay by precise positioning services (PPS), and L5 (1176.45 MHz) was introduced to increase the accuracy for safety-of-life transportation applications [5]. The frequency plan and the band of GPS are shown in Figure 1, together with the information of Galileo and Glonass [6].

Figure 1.

GNSS frequency bands.

2. Theory of GNSS-reflectometry

GNSS-reflectometry is usually called for short as GNSS-R. Its main principle is to receive and further extract information from the GNSS signals which were reflected off the Earth surface. It works as a bistatic radar, in which the transmitter and the receiver are separated by a significant distance, comparable to the expected distance to the target. The concept can be expanded that the receiver receives multiple signals simultaneously from different transmitters, in which case, it is multistatic radar. Comparing with the other existing remote sensing tools, GNSS-R has following advantages: 1. The additional transmitter is not needed; 2. Many signal sources can be chosen, including GPS, Galileo, GLONASS, and BeiDou/Compass; 3. Use spread-spectrum communication technology for receivers to receive weak signals; 4. It provides a wide range of applications, for example, soil moisture retrieval, vegetation biomass detection, altimetry estimation, sea wind and wave height, and so on [1, 7, 8, 9, 10, 11, 12, 13].

2.1. GNSS-R geometry

In the case of GNSS-reflectometry, the transmitter and receiver are in the different locations, and the scattered signals can be originated from a variety of geometries, given a target that scatters sufficient power to the direction of the receiver. For the land surface, the typical bistatic geometry consists of a transmitter and receiver above the surface, with scattering taking place mainly from the region of the surface surrounding a specular reflection point. The specific geometry is depicted in Figure 2.

Figure 2.

GNSS-R as a bistatic radar geometry.

Figure 2 depicts signals transmitted from a GPS satellite (about 20,200 km height), scattered from the land surface along the specular direction, and received by a receiver equipped on an aircraft. The direct signal is received by the zenith antenna for processing positioning information or system calibration, and the reflected signal is received by the nadir antenna for remotely sensing the land surface. The incident RHCP GPS signal is predominantly left-hand circular polarization (LHCP) after specular scattering. The specular point is defined as a reflection point characterized that the incident and reflected angles are equal (Snell’s Law). The strong coherent scattering component concentrated about a narrow area called the first Fresnel zone around the specular point. The size of this active region is usually considered to be the first Fresnel zone for which the differential phase change across the surface is constrained to λ/2 radians, where λ is the signal wavelength [14]. Surface roughness causes scattering from a glistening zone around the specular reflection point. The scattering component, also called incoherent or diffuse component, consists of power scattered randomly in all directions, with its magnitude smaller than coherent component.

When the distance from the surface to the transmitter is much larger than the distance from the surface to the receiver (or vice versa), the semi-major axis aand semi-minor axis bof the first Fresnel ellipse are dependent on the incidence angle θ, and minimum height of the transmitter or receiver H was given by [15]:

α=2δHcosϑcosθ2,b=2δHcosθcosθ,E1

where δis the phase change measured in wavelengths (λ/2for the first Fresnel zone) and θis the incident angle.

As the signal scattering, the time delay and the frequency of the received signals change. The narrowly spaced iso-range ellipses (lines of equal delay across the surface) and iso-Doppler hyperbolas (lines of equal Doppler frequency across the surface) can be mapped across the Earth as shown in Figure 3. At each point in the glistening zone, the path delay and reflection angles are different. This results in a range of different path delays (between the transmitter and receiver) and Doppler frequencies at the receiver. Lines of constant delay or iso-range lines can be drawn as ellipses centered at the point of specular reflection. The lines with the same radial velocity are the cuts of the cone surfaces with the earth plane, i.e., parabola. At fixed radio frequency (RF), the radial velocity is proportional to the Doppler frequency. Therefore, the cone can also be called Doppler cone, and the cone cuts with the earth surface called iso-Doppler lines [16, 17].

Figure 3.

Iso range and iso-Doppler lines across the Earth’s surface around the specular reflection point.

In Figure 3, the red ellipses represent the first Fresnel zone. The uneven ring belts indicate equal-code delay lines (yellow lines). The black curves are the Equal-Doppler lines. GPS satellites are far from the Earth and thus the Fresnel belts are small. The surface property can be extracted through analyzing the reflected signal and the direct one.

2.2. Bistatic radar equations

In general, the total scatter power ppqrmeasured by a bistatic radar can be expressed as [18]:

ppqr=ppqc+ppqi,E2

where subscripts qand pdenote the incident and scattered polarizations. ppqcand ppqiare the coherent power and incoherent power, respectively.

2.2.1. Coherent component from smooth surface

For the coherent component in the case of like polarized GPS bistatic radar, the bistatic equation can be written as [18],

Plrc=PrtGt4πRst+Rrs2Grλ24πΓlr.E3

Subscript lrrepresents the scattering when a satellite incident signal (right-hand polarized) is scattered by the surface and inverts the polarization to the left-hand, Prtis the transmitted RHCP signal power, Gtis the transmitter antenna gain, Gris the receiver antenna gain, and λis the wavelength (19.042 cm for GPS L1 signal). Rrsis the distance between the receiver and the specular point, Rstis the distance between the satellite and the specular point, and Γlris the power reflectivity of the reflecting surface.

2.2.2. Coherent component from rough surface

As the smooth surface begins to transition to a rough surface, the coherent component of the reflected power is decreased. Therefore, the term Γlrin Eq. (3) decreased from that of the smooth surface case due to increasing roughness, which was described by [14].

Γlrθ=Rlrθ2χz,E4

where Rlris the Fresnel reflection coefficient and χzis the probability density function of the surface height z.

The Fresnel reflection coefficient Rlrcan be expressed as linearly polarization modes [7]:

Rlr=Rrl=12RvvRhh,E5
Rrr=Rll=12Rvv+Rhh,E6

where Rvvand Rhhare the Fresnel coefficients for horizontal and vertical polarization [19]:

Rhhθ=cosθεrsin2θcosθ+εrsin2θ,E7
Rvvθ=εrcosθεrsin2θεrcosθ+εrsin2θ,E8
εr=ε/ε0j60λσ,E9

where θis the incident angle, in which εis the dielectric constant of the surface, ε0is the dielectric constant of the air, λis the wavelength of the signal, and σis the electric conductivity.

Most of the natural surface can be modeled by a Gaussian height distribution. Thus, the reflectivity can be written as:

Γlrθ=Rlrθ2ehcosϑ2.E10

hwas called roughness parameter, and it is defined as:

h=4k2σ2,E11

where kis the wave number (2π/λ) and σis the standard deviation of the surface. When it is a smooth surface (h=0), the reflectivity Rlris simply the square of the Fresnel reflection coefficient.

2.3. Coordinate reference systems

Latitude, longitude, and height are three coordinates used to characterize the user’s position on Earth. A reference system with respect to the coordinates must be defined correspondingly. In the beginning, the World Geodetic System 1984 ellipsoid (WGS84) in an Earth-centered, Earth-fixed (ECEF) coordinate system is defined for typical GPS measurements. The ECEF coordinates represent positions as (x,y,z). The x-axis points in the direction of 0°longitude (passing through the intersection of the Greenwich meridian with the equatorial plane), and the +y-axis points in the direction of 90°east longitude. The z-axis is chosen to be normal to the equatorial plane in the direction of the geographical north pole, thereby completing the right-handed coordinate system. The ECEF system enables users to know the positions with respect to the center of the Earth. The ellipsoid coordinate origin of WGS84 (latitude, longitude, and attitude) is meant to be located at the Earth’s center of mass. It is a simplest geometric model, which best fits the entire Earth surface [20].

Since the Earth is not spherical and not uniform in density, the Earth Geoid Model 1996 (EGM96) was defined as the locus of all points with the same gravity potential (equipotential surface) best fitting with the mean sea level on the Earth [21]. When the height is indicated relevant to the mean sea level, it refers to the geoid model. The height of the geoid is denoted Nwith respect to the ellipsoid, which is shown in Figure 4. At each point, a geoidal height (N) is defined along the line perpendicular to the ellipsoid. It is measured from the ellipsoid to the geoid. Height his the distance measured from the ellipsoid to the object. Hindicates the height with respect to the geoid called orthometric height.

Figure 4.

Relation of an ellipsoid, geoid, and terrain heights.

Georeferencing specular points on maps were done based on the ECEF coordinate system, by assuming the surface is locally flat and the surface height remains the same with that of the receiver projection position. However, a typical onboard GPS receiver and Google Earth software provide the ellipsoid height hbased on the EGM96 model. Therefore, the EGM96 coordinates information was transferred in WGS84 and further converted to ECEF coordinates, in order to calculate the elevation angle with respect to the satellite position (x,y,z). On the other hand, the orthometric height His useful, such as in the altimetry application. Hence, some computations are needed as Eq. 12 [3].

HhN.E12

3. GNSS-R applications

3.1. Background and objective

After GNSS signals reflected off the Earth surface, the scattered signal received by the GNSS-R system contains the information of the surface. Hence, GNSS-R as an Earth’s surface remote sensing tool has been widely studied for characterizing the reflecting surface. Several retrieval methods for soil moisture content (SMC), altimetry, and vegetation biomass estimation can be used to reveal the geographical and biological parameters of the Earth surface.

For altimetry application, the GPS Bistatic technique for altimetry application was first proposed by Martin-Neira [17]. The aim of the GNSS altimetry is to detect the height between the receiver and the reflecting surface. In principle, the altimetry techniques can be applied to any surface, but their performance will depend on the detecting environment [22, 23, 24].

For soil moisture application, frequencies in the range of 1–3 GHz have been shown to be ideally suited to sense soil moisture due to reduced atmospheric attenuation [25]. Njoku and O'Neill [26] showed that P-band (0.775 GHz) and L-band (1.4 GHz) frequencies are optimal for sensing soil moisture in the top 0–4 and 0–2 cm surface layers, respectively.

For biomass content sensing, a couple of papers mentioned that a retrieval method with very preliminary analysis is available. The one provided by Zavorotny et al. [27] proposed modeling GPS reflected signals having a different polarization state that can be used for sensing a homogeneous soil of rough land surface. The results were obtained by modeling the impact of soil moisture from the ratio between expected power levels in orthogonal polarizations (LR/RR,HR/VR), in function of the elevation angle. By varying the soil moisture level (x-axis) and different incident angles (the angle between the propagation path and the normal to the surface), the ratio of two orthogonal polarizations received power proved to be independent of the surface roughness factor and sensitive to the soil moisture. With this established theoretical hypothesis, the algorithm with respect to the orthogonal polarization was developed and presented as follows.

In the following, the retrieval methods for the SMC, vegetation sensing, and altimetry applications are applied and discussed in the framework of Italian SMAT-F2 (System of Advanced monitoring of the Territory—phase 2) project. First, the SMC retrieval and vegetation biomass sensing methods are detailed. With considering different kinds of measured signals, three kinds of retrieval methods based on the bistatic configuration are presented. Then, the C/A code altimetry method, which provides the possibility for altimetry estimation, will be discussed.

3.2. Methodologies

3.2.1. Overview of soil moisture retrieval methods

Soil dielectric constant and soil moisture retrieval play fundamental aspects in agriculture and water cycle attracting interests from many researchers and have started to produce some results. Three different methodologies were implemented mainly for this kind of applications:

  1. Multipath effect and its relation to soil moisture: it uses a standard ground-based nearly hemispherical RHCP antenna to receive the direct signal (see Figure 5(a)). The reflected signal originated from ground creates multipath effects, since it interferes with the direct signal. The total received signal amplitude has a sinusoid behavior with the change of sine of elevation angle. By measuring the received SNR ratio, it is possible to retrieve soil moisture information [9, 27].

  2. Interference pattern technique (IPT): usually a horizontal-pointing vertically polarized antenna is used (see Figure 5(b)). This technique consists of measuring the power fluctuations related to the signal resulting from the simultaneous reception (and interference) of the direct and the reflected GNSS signals. It is similar to a multipath effect. In this case, the two rays approach is adopted. Soil moisture retrieval is based on finding a specific ‘notch’ point from the interference pattern, versus satellite elevation [22, 28, 29, 30].

  3. Bistatic method: this is based on the separate reception of direct and reflected signals using two different antennas and on the separate measurement of signal powers. Depending on the antenna configuration, three possible observing systems exploiting the bistatic geometry can be further identified:

    1. A down-looking LHCP antenna and an up-looking RHCP antenna: an LHCP antenna receives reflected GPS signal from the surface (see Figure 5(c)). The power reflectivity could be obtained either by using a bistatic radar equation [31] or from the power ratio between the reflected signal and the direct signal [32, 33]. The reflectivity is then a function of the dielectric constant of the soil, the elevation angle, and the surface roughness. By properly choosing the surface roughness parameter (elevation is known from the direct signal) for a certain scattering model (see for example [34]) such as the small perturbation method (SPM) and Kirchhoff approach (KA), the dielectric constant can be retrieved.

    2. One RHCP up-looking antenna and two down-looking antennas (or channels) with one RHCP polarized and the other LHCP polarized [7, 12, 35]. With this configuration (see Figure 5(d)), it is possible to measure both the co-polar component of the terrain reflectivity by using the LHCP signal and the cross-polar component by using the RHCP antenna. The ratio of these two reflectivities was in good correlation with SMC, and it was independent of the surface roughness.

    3. A similar configuration to the one described in (3.a) but with horizontally (H) and vertically (V) polarized antenna for both the up and down-looking directions (see Figure 5(e)): this configuration has been tested through some simulation studies [36]. The ratio between the reflected and the direct power on the horizontal polarization and the same ratio on the vertical polarization depend on the soil reflectivity and surface roughness. If the power ratio between the two channels with orthogonal polarization is considered, the influence of the surface roughness can be canceled. It was verified that the final expression holds under different scattering models, which means that it could be applied to a wide range of surface roughness. This dielectric constant retrieval approach is based on the use of the ratio of power densities scattered at HHand VVpolarizations along the specular direction for different incidence angles. Since the ratio is a function of both the elevation angle and the dielectric constant, a minimum least square technique was applied to better define the dielectric constant, by measuring at least two different elevation angles.

Figure 5.

System schemes for different observing strategies.

The retrieval of the soil moisture from the dielectric constant at microwave band (especially in L band) has been widely investigated, and several well-accepted theoretical and empirical models have been established, such as [25, 28, 37, 38]. The information of soil texture in terms of percentage of clay and sand should be known and provided in input to such models. A further model described by Topp et al. [39] does not depend on any input information, since it models the relative permittivity as a third-order polynomial function of soil moisture.

3.2.2. Bistatic methods for soil moisture retrieval

In this section, the algorithms used for SMC and vegetation sensing are applied based on the (5.a): bistatic method. The retrieval process aims to establish the link between received reflected signals and the dielectric constant of the soil. Moreover, the retrieval of dielectric constant is the key component to the retrieval of soil moisture, since there are several well-established models giving this relationship as introduced in the last section.

3.2.2.1. Only LH reflected GPS signals are measured

Assuming a perfectly smooth surface and regardless the surface roughness and incoherent components (Γlr=Rlrϑ2), the dielectric constant is retrieved from the LHCP-reflected SNR which can be processed by the open loop approach.

The expression of SNR is derived from bistatic radar equation:

SNRpeakreflect=14PtGt4πR1+R22Grλ2GD4πPNRvvRhh2,E13

where Ptis transmitted signal power, Gtand Grare antenna gains of transmitter and receiver, respectively. λis signal wavelength (about 19 cm for GPS L1 signal). R1and R2represent the distances between specular point and transmitter or receiver. The superscript c denotes the coherent component, PNis noise power, and GDis the processing gain due to the de-spread of the GPS C/A code. The value of the correlation gain is determined by the chipping rate of the GPS L1 C/A code (1.023e6Hz) and 1 ms processing interval. It results:

Gpr=10log101.023e61000=30.1dB.E14

Fresnel reflection coefficients of HHand VVare given in Eqs. (7) and (8). The dielectric constant is then easy to be solved by combining the Eqs. (7), (8), and (13). Input parameters such as Pt, Gt, Gr, and GDare usually seen to be constants. Since they are system parameters of GPS, PNshould be estimated or calculated properly. R1, R2, and θare geographic information which can be calculated through receiver positions. SNR is the input parameter obtained from signal postprocessing. Knowing all the above input parameters, the dielectric constant can be solved numerically.

3.2.2.2. Both LH reflected and RH direct GPS signals are measured

In this case, both SNRs are obtained from RH direct signal and LH reflected signal. The dielectric constant can be retrieved from the power ratio of LHCP reflected SNR over direct RHCP SNR, provided that direct RH signals are processed with the same open-loop approach to get SNR.

The SNR expression of the peak power for the direct RHCP signal is:

SNRpeakdirect=PtGt4πR32Grλ2GD4πPN,E15

where R3is the distance between the transmitter and receiver. It must be noted that receiver gain Grand noise power PNgiven in Eq. (15) are not equal to those given in Eq. (13) for the reflected signal path, and this is the reason why calibration is needed. By making the ratio of Eq. (13) to Eq. (15) canceling the same parameters and summarizing the uncertainties of Grand PNinto a single calibration parameter C, we have:

SNRpeakreflectSNRpeakdirect=R324R1+R22RvvRhh2C.E16

Now, input parameters of R1, R2, θand R3are only needed, and they are easy to obtain from positioning information. Calibration parameter Ccan be either ignored if there is no way to calibrate the system or calculated through over-water measurement. The dielectric constant can be obtained from Eq. (16), given all the input parameters defined.

3.2.2.3. Both LH reflected and RH reflected GPS signals are measured

The retrieval process has linked the measured reflected LHCP and RHCP power (or SNR) to the normalized bistatic radar cross section σqpo, which has been widely studied and modeled for different surface roughness and scattering models, such as Kirchhoff Approximation in stationary-phase approximation (Kirchhoff Geometrical Optics, KGO), Kirchhoff approximation in physical optics approximation (KPO), and small perturbation method (SPM) [34].

Applying polarization scattering matrix [40] and combining some basic electromagnetic field theories, the cross-polarization power ratio of LHCP to RHCP has the relation with normalized bistatic radar cross section of linear polarizations:

PlrrPrrr=σohh+σovvσohhσovv.E17

Note here that in Eq. (17), σois approximated by the special case of σqpspecoin the specular direction for simplicity, and therefore, the cross polarization terms of σohvand σovhare computed to be zero for the models of KGO, KPO, and SPM in the specular direction case. Also, σofor hhand vvpolarizations given by these three models are all products of a polarization sensitive parameter, usually the Fresnel reflection coefficient and other polarization independent surface roughness parameters (detailed expressions are given by [Ulaby, 1982]). Then, the expression shown in Eq. (17) can be simplified as:

PlrrPrrr=Rhh+RvvRhhRvv,E18

where Ris the Fresnel reflection coefficient as given in Eqs. (7) and (8). This equation also holds the condition of specular reflection, and the result shows that the cross-polarization power ratio of LHCP to RHCP of reflected signals is independent with surface roughness and is only a function of dielectric constant and incident angle.

The received power for both lrand rris obtained as input parameters after post-processing. The dielectric constant of the soil is solved after taking all the parameters.

3.2.3. Retrieval methods for vegetation biomass sensing

The GNSS signals in bistatic configuration were obtained by the cross-correlation between direct and reflected signals with the pseudo-random noise (PRN) code replica, respectively. The output of the cross-correlation is called Delay-Doppler map (DDM), and Delay Waveforms (DW) is a row from DDM with the Doppler frequency set to a constant value of where the peak is. The phase information of the signal is lost during the integration process, resulting in amplitude <Ypq>or power <Ypq2>waveforms for the incident q-polarized to p-polarized scattered signals. The observable polarimetric parameters Γlrand Γrrrepresent the cross and co-polar reflectivity can be defined as in [12]:

Γpq=Yreflected,p(Δτ,f)peakYdirect,q(0,f)peak2,E19

where Δτis the delay difference between the direct and the reflected paths, and the fis assumed to be the Doppler frequency shift of the signal. With ground-based or low-altitude receivers, it also can be assumed to be the same, fdfr=f. The above equation obtained both the coherent and incoherent component of the received direct and reflected signals. However, the unstable airborne platform, terrain variation, and fading phenomenon may impair the estimation of the coherent scattering component. By adding the time series, an averaging procedure can be used to remove the contribution of the incoherent component on the final observed reflectivity. It was written as:

Γpq=<Yr,ptYd,qt>2=<Yr,ptYd,qt2>δYr,ptYd,qt2.E20

After the reduction of the incoherent component from the received signal, the received signal is mainly composed of coherent components. Therefore, the spatial resolution can be roughly estimated from the first Fresnel zone [31]. For instance, the 600 m flight altitude and 40o incident angle lead to an elliptical active scattering area with a minor axis of 24 m and a major axis of 32 m. Finally, the reflectivities Γlrand Γrrand the ratio of Γlr/Γrr, which was called polarimetric ratio (PR), are the three observable parameters used to investigate the SMC and the vegetation biomass of the terrain.

3.2.4. GPS C/A-code altimetry

The geometry of the reflecting surface with the receiver and transmitter is shown in Figure 6 under the assumption of a flat Earth and ray optics. The path delay (δ) between the direct and the reflected signals can be taken as a delay difference (Nsamples) in the DDM obtained from Eq. 21.

δ=cfsNsamples=cfsNchip10231msfs,E21

where Cis the light speed and fsis the sampling frequency of the GNSS-reflectometer, and Nchipsis the number of chips in one period.

Figure 6.

GNSS altimetry geometry with signals considered in a parallel condition.

Assuming that the signals are parallel as shown in Figure 6, it gives that θ=αfor the incoming signals. Height his possible to be derived from the following equation [41, 42]:

hδ2sinα,E22

where δis the excess path and θis equal to the elevation angle α. Therefore, the estimated altitude of the receiver hdepends on the measured delay between the two signals and the elevation of the satellite.

3.2.5. Data processing

The GNSS-R system used for the acquisition of the received power can measure both the direct GPS signal and the reflected one using two independent acquisition channels: one for the direct signal and the other for the reflected signal. The receiver consists of two commercial front-ends (FEs) connected to two antennas: for instance, a traditional hemispherical GNSS L1 patch antenna pointed upward for the measurement of the direct signal. A high-gain LHCP antenna pointed downward for measuring the reflected signal as shown in Figure 7. Each FE is connected to a PC for data storage. The antennas and the FEs are mounted on a Kevlar bar. A typical tripod is used for more efficient adjustment of the orientation of the antennas.

Figure 7.

GNSS-R system block diagram: Antennas and front-end connections.

The raw data of direct and reflected signals are collected and postprocessed with the software SOPRANO [43] for obtaining the peak power and SNR for each satellite. An open-loop approach is used for processing to obtain delay-Doppler maps (DDMs) and the corresponding delay waveforms. SNR time series were estimated from several noncoherently integrated delay waveforms generated by cross-correlation of the incoming data and local replica C/A codes.

Generally, increasing coherent integration time can help the low SNR to generate the DDM, but the performance may worsen when it exceeds a certain value. In this case, the coherent integration time was 1 ms because of the length of the GPS C/A code (1 ms). On the other hand, since the reflected signal power was attenuated by the surface as well as by the noise effects, it could not detect the real signal power. A noncoherent integration time technique (called averaging) was adopted during the postprocessing in order to mitigate the noise. Generally, it is in the range of 100-1000 ms. The noncoherent integration time, in this case is 500 ms. Moreover, the information of positions and altitude was recalculated, which was used to calculate the specular points on object areas. During the experiment, satellite positions (x,y,z) in an ECEF coordinate were collected by a GPS recorder. In the meanwhile, the receiver position (latitude, longitude, attitude) in geoid coordinate was collected in a file (.kml) which can be opened by the Google Earth software. The flowchart of obtaining the navigation message is shown in Figure 8.

Figure 8.

Flowchart of processing navigation message.

The height coordinate z in the file (.kml) should be transferred into the WGS84 coordinate first, because the height z is in EGM96 coordinate system which already considered the sea surface altitude. The second step is to transfer the WGS84 coordinates (latitude, longitude, height) into the ECEF coordinate that can correlate with the satellite positions (x,y,z). With having receiver and satellite positions both in ECEF coordinates, elevation angle can be calculated with respect to the local tangent plane.

4. Results and discussions

A low-altitude airborne experiment was carried out by a P92 Digisky airplane over the Avigliana lake (45.099°N, 7.369°E) in Italy. An image of the experimental area from Google Map is shown in Figure 9.

Figure 9.

Avigliana area.

In this area, there are two lakes, the size of the northern lake (bigger) is approximately 1 km x 1.3 km, while that of the southern lake (smaller) is 700 m × 1.1 km. The area was selected for the experiment because of the following reasons. First, the presence of two lakes of the experiment can provide the reflections and the known permittivity for calibration. Second, in the North and South directions of the two lakes, there are some grass areas which the flight will go over, and in the West and East directions, the flight coverage areas are full of forest. This kind of geographical environment provides the possibility to investigate the roughness and vegetation biomass at the same time.

The elevation and azimuth angles were obtained, as shown in Figure 8. The sky plot of GPS satellites at the 1570th second of the flight is shown in Figure 10.

Figure 10.

Sky plot of satellites on 18th June 2015.

The reflected LHCP signal was collected through a high gain dual polarization antenna. With the direct RH signal obtained, the reflectivity was computed by the SNR of LH/RH signal component. The SMC was obtained by using the method introduced in Section 3.2.2. The SMC estimation result of PRN 15 with high elevation angle is shown in Figure 11, calibrated by the water reflectivity of 0.67. The color of the specular points is proportional to the estimated soil moisture (the red color denotes totally dry soil, and the green color represents humidity 100%).

Figure 11.

Soil moisture content estimation of PRN 15.

The flight flew covered lake and two forest areas from east to west. The flight route was denoted by the yellow color in Figure 12. Since the reflection point begins from the forest in the east (mid-biomass) to big lake (flat) then stop with the forest in the west (big-biomass) as shown in Figure 12, it can be used to investigate the vegetation biomass sensing.

Figure 12.

Flight route (yellow line) and PR level of reflection point georeferenced in Google Map.

Based on the method introduced in Section 3.2.3, the reflectivities Γlrand Γrrand the ratio of the two polarization (lrand rr) PR can be calculated and shown in Figure 13. From 720 s, the big increment 6 dB of the average lrappears in the lake (L1) compared to the forest. The lrcomponent is quite sensitive with SMC. The average rrin forests (F1 and F2) is 2–4 dB lower than the lake (L1). In addition, considering the average lrand rrin forests (F1 and F2), both average in F2 are higher than F1. It could be considered that the vegetation biomass in F2 is lower than F1. It confirms the fact and further implies that both polarizations are sensitive to vegetation biomass. Moreover, the biggest PR level is the presence of the lake. Therefore, the PR can be taken as a good parameter for soil moisture estimation regardless of the vegetation biomass. The PR level is 5.6 dB in the lake, 3.9 dB in F1, and 2.3 in F2.

Figure 13.

PR level, Γlr and Γrr level.

For the altimetry estimation, the C/A code method (see Section 3.2.4) was applied on the unsmooth surface to obtain the target height. Due to bandwidth and the equipment, the altimetry resolution is limited, according to 1 sample 23.07 m. The altimetry estimation results of satellite PRN 15 with high elevation angle are presented in Figure 14.

Figure 14.

Altimetry estimation of PRN 15 from 1015s (a), 1062s (b).

Figure 14 shows two figures of height estimation of PRN 15 start from 1015s and 1062s. The elevation angle of PRN 15 is 76°. The estimation height is quite same with GPS recorder. The estimation deviation in Figure 14 (a) is bigger than Figure 14 (b). It confirmed that the best condition for altimetry application is when using high elevation angle on the water surface, which can greatly avoid the interference of forest and cities. The accuracy of the altimetry estimation depends on the bandwidth, algorithm, and equipment, and it shows that it is quite affected by the environment as introduced in the beginning of Section 3.

5. Conclusion

In this chapter, the global navigation satellite system-reflectometry (GNSS-R) technique as a method of remote sensing that receives and processes microwave signals reflected from various surfaces to extract useful information is reported. GNSS-R technique has gained widespread attention in a variety of applications. This work assesses the capability of measuring soil moisture and vegetation biomass by means of reflected signals from global navigation satellites from a theoretical and experimental point of view.

The SMC was obtained and georeferenced in Google Maps. The type of terrain ranged from water (on the lakes) to terrain with small bushes to built-up areas. The results showed good correlation with the types of the underlying terrain. It was observed that the LH was rather more sensitive to SMC. The different SMC areas were detected with different PR values. The ratio of the two polarizations was independent of roughness, and it proved to be an optimum parameter for soil moisture estimation. Moreover, it was found that increasing the biomass quantity decreased the Γlrand Γlrpolarization components, since the reflected signal components are attenuated from leaves, trunks, crowns, and so on. The altimetry estimation is based on the calculation of the delay differences between the two signals. In principle, the altimetry techniques can be applied to any surface, but their performance will depend on the scattering behavior. Therefore, GNSS-R altimetry has only been used on strongly reflecting surfaces and geometries, such as waters and smooth ice. The accuracy of the derived height is limited also due to the bandwidth and the sampling frequency of the equipment.

Acknowledgments

Some of the works illustrated in this chapter were developed by the authors and their collaborators at the University of Politecnico di Torino. The authors would like to thank the SMAT-F2 project, ISMB and all the researchers involved. Furthermore, this research is financed by National Science Foundation of China under projects No. 61601076 and sponsored by NUPTSF (Grant No. 217152).

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Yan Jia and Yuekun Pei (July 25th 2018). Remote Sensing in Land Applications by Using GNSS-Reflectometry, Recent Advances and Applications in Remote Sensing, Ming-Chih Hung and Yi-Hwa Wu, IntechOpen, DOI: 10.5772/intechopen.72901. Available from:

chapter statistics

349total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Recent Advances and Applications in Remote Sensing

Edited by Ming Hung

Next chapter

Remote Sensing with Shipborne High-Frequency Surface-Wave Radar

By Junhao Xie, Minglei Sun, Zhenyuan Ji and Guowei Yao

Related Book

First chapter

Application of Spatial Statistics in Transportation Engineering

By Uday R.R. Manepalli and Ghulam H. Bham

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us