Article Analytical Computation of the Spatial Resolution in GNSS-R and Experimental Validation at L1 and L5
Adriano Camps 1,2,* and Joan Francesc Munoz-Martin 1,2
1 CommSensLab-UPC, Department of Signal Theory and Communications, UPC BarcelonaTech, c/Jordi Girona 1-3, 08034 Barcelona, Spain; [email protected] 2 Institut d’Estudis Espacials de Catalunya-IEEC/CTE-UPC, Gran Capità, 2-4, Edifici Nexus, despatx 201, 08034 Barcelona, Spain * Correspondence: [email protected]
Received: 3 November 2020; Accepted: 24 November 2020; Published: 28 November 2020
Abstract: Global navigation satellite systems reflectometry (GNSS-R) is a relatively novel remote sensing technique, but it can be understood as a multi-static radar using satellite navigation signals as signals of opportunity. The scattered signals over sea ice, flooded areas, and even under dense vegetation show a detectable coherent component that can be separated from the incoherent component and processed accordingly. This work derives an analytical formulation of the response of a GNSS-R instrument to a step function in the reflectivity using well-known principles of electromagnetic theory. The evaluation of the spatial resolution then requires a numerical evaluation of the proposed equations, as the width of the transition depends on the reflectivity values of two regions. However, it is found that results are fairly constant over a wide range of reflectivities, and they only vary faster for very high or very low reflectivity gradients. The predicted step response is then satisfactorily compared to airborne experimental results at L1 (1575.42 MHz) and L5 (1176.45 MHz) bands, acquired over a water reservoir south of Melbourne, in terms of width and ringing, and several examples are provided when the transition occurs from land to a rough ocean surface, where the coherent scattering component is no longer dominant.
Keywords: GNSS-R; spatial resolution; diffraction; experiment; airborne; L1; L5
1. Introduction Global navigation satellite systems reflectometry (GNSS-R) is a relatively new remote sensing technique that uses navigation signals as signals of opportunity in a multi-static radar configuration [1]. The spatial resolution of GNSS-R instruments was first analyzed in [2], and it was estimated from the dimensions of the area determined by the intersection of the first iso-delay and iso-Doppler lines. In this case, for a Low Earth Orbit (LEO) satellite at 700 km, the predicted spatial resolution for an instrument using the GPS L1 C/A (Coarse/Acquisition) code would be in the range of 36 to 54 km, for elevation angles from 90◦ to 40◦. Later on, the definition of spatial resolution was revised and “an effective spatial resolution derived as a function of measurement geometry and delay–Doppler (DD) interval, and as a more appropriate representation of resolution than the geometric resolution previously used in the literature” was proposed in [3]. This effective spatial resolution was defined following an approach similar to the definition of effective bandwidth used in filter theory. An increase in the effective spatial resolution is observed for increasing incidence angle, receiver altitude, and delay–Doppler window. With this definition, the effective spatial resolution is estimated to range from 25 to 37 km for a 700 km height receiver (Figure 4, [3]).
Remote Sens. 2020, 12, 3910; doi:10.3390/rs12233910 www.mdpi.com/journal/remotesensing Remote Sens. 2020, 12, 3910 2 of 15
The above definitions implicitly assume that the scattering is incoherent, as the radar equation used is the one that corresponds to surface scattering. However, the development of advanced GNSS-R receivers,Remote Sens. capable 2020, 12 of, x FOR processing PEER REVIEW the coherent component and separating it from the incoherently2 of 15 scattered component, showed that much finer features can be distinguished [4,5]. In fact, when coherent GNSS-R receivers, capable of processing the coherent component and separating it from the scattering occurs, the GNSS reflectivity values have a very high spatial resolution, as most of the power incoherently scattered component, showed that much finer features can be distinguished [4,5]. In fact, is mainly coming from the first Fresnel zone [6], an ellipse of semi-major axis a (Equation (1a)) and when coherent scattering occurs, the GNSS reflectivity values have a very high spatial resolution, as semi-minor axis b (Equation (1b)): most of the power is mainly coming from ther first Fresnel zone [6], an ellipse of semi-major axis a RT RR (Equation 1a) and semi-minor axis b (Equationa = 1b):λ · (1a) RT + RR 𝑅a ·𝑅 𝑎=b = 𝜆 (1a) 𝑅 +𝑅, (1b) cos(θi) where λ is the electromagnetic wavelength, R𝑏=T and RR, are the distances from the transmitter(1b) and receiver to the specular reflection point, and θi is the incidence angle. This means that even from an where λ is the electromagnetic wavelength, RT and RR are the distances from the transmitter and LEO satellite, the spatial resolution can be as small as ~400–500 m. receiver to the specular reflection point, and θi is the incidence angle. This means that even from an In [7] the spatial resolution under coherent scattering was theoretically analyzed in detail, including LEO satellite, the spatial resolution can be as small as ~400–500 m. the effectIn of [7] the the di ffspatialerent Fresnelresolution zones. under It coherent was shown scattering that up was to ~45theoretically◦ most of analyzed the power in isdetail, coming fromincluding an area determined the effect of bythe thedifferent first Fresnel Fresnel zone,zones. although It was shown a power that up equivalent to ~45° most to the of the one power received is in free-spacecoming conditions from an area is actually determined coming by the from first a regionFresnel ~0.58–0.62 zone, although times a thepower first equivalent Fresnel zone. to the However, one significantreceived contributions in free-space fromconditions regions is actually farther coming away than from the a region first Fresnel~0.58–0.62 zone times also the contribute first Fresnel to the receivedzone. power,However, making significant its interpretation contributions from more regions difficult. farther This away region than wasthe first ~800 Fresnel m away zone atalso a 60◦ elevationcontribute angle to for the an received LEO satellite power, atmaking 500 km its height. interpretation more difficult. This region was ~800 m awayIn this at work,a 60° elevation we have angle gone for one an step LEO further satellite and at 500 analyzed km height. the spatial resolution from a theoretical point of view,In this using work, known we have principles gone one of electromagnetism. step further and analyzed The use ofthe these spatial results resolution makes from it possible a to computetheoretical in apoint very of compact view, using form known the expression principles of of electromagnetism. the ringings and The oscillations use of these observed results makes when the it possible to compute in a very compact form the expression of the ringings and oscillations observed specular reflection point transitions from the land into the ocean or vice versa. Analytical results are when the specular reflection point transitions from the land into the ocean or vice versa. Analytical validated with airborne experimental data gathered at L1 and L5 with the Microwave Interferometric results are validated with airborne experimental data gathered at L1 and L5 with the Microwave ReflectometerInterferometric (MIR) Reflectometer instrument (MIR) when instru crossingment the when coastline crossing [8 ].the coastline [8]. 2. Methodology 2. Methodology TheThe scattering scattering geometry geometry is shownis shown in in Figure Figure1, where1, where the the receiver receiver (R) (R) and and transmitter transmitter (T) (T) lielie inin the YZ plane,the YZ and plane, are and at heights are at heights hR and hR h andT over hT over a flat a flat surface. surface. RT RandT and R RRare are the the distancesdistances to the the origin origin of coordinatesof coordinates where where the specular the specular reflection reflection point is,point and is,θ iandis the θi is incidence the incidence angle, angle, which which for the for specular the reflectionspecular point reflection (origin point of coordinates) (origin of coordinates) is equal to is the equal reflected to the anglereflectedθr .angle The diθr.ff Theerent different disks represent disks the Fresnelrepresent zones, the Fresnel with thezones, odd with ones the (1, odd 3, 5,ones etc.) (1, colored3, 5, etc.) dark colored gray, dark and gray, the and even the ones even (2, ones 4, 6,(2, etc.) colored4, 6, white.etc.) colored white.
FigureFigure 1. Global 1. Global navigation navigation satellite satellite systems systems reflectometryreflectometry (GNSS-R) (GNSS-R) scattering scattering geometry geometry showing showing the the transmitter,transmitter, receiver, receiver, specular specular reflection reflection point point (origin(origin of of coordinates), coordinates), and and the the different different Fresnel Fresnel zones. zones.
Remote Sens. 2020, 12, 3910 3 of 15 Remote Sens. 2020, 12, x FOR PEER REVIEW 3 of 15
InIn what follows,follows, a a few few principles principles of electromagnetismof electromagnetism are reviewed. are reviewed. They willThey be will used be to computeused to computea closed-form a closed-form solution solution of the GNSS-R of the GNSS-R response response to a step to function, a step function, such as insuch the as transition in the transition between betweenwater and water the ocean.and the ocean.
2.1. Review of Some Basic Principles of Electromagnetism 2.1. Review of Some Basic Principles of Electromagnetism 2.1.1. Principle of Equivalence and Image Theory 2.1.1. Principle of Equivalence and Image Theory The principle of equivalence [9] allows us to compute the solutions in the region of interest The principle of equivalence [9] allows us to compute the solutions in the region of interest (half- (half-space in front of the conducting plane) by replacing the plane conductor with the images of the space in front of the conducting plane) by replacing the plane conductor with the images of the dipoles. The image sources must satisfy the boundary condition that the tangential electric field at the dipoles. The image sources must satisfy the boundary condition that the tangential electric field at conducting surface should be zero. The total electric field above the conducting plane (Figure2) is then the conducting surface should be zero. The total electric field above the conducting plane (Figure 2) computed as the sum of the one created by the original currents and the image ones, while the total is then computed as the sum of the one created by the original currents and the image ones, while the electric field below the conducting plane is zero. Applying this principle to our problem, the transmitter total electric field below the conducting plane is zero. Applying this principle to our problem, the can be replaced by its “image” underneath the plane, as illustrated in Figure3. transmitter can be replaced by its “image” underneath the plane, as illustrated in Figure 3.
Figure 2. Left: linear and loop currents over a perfect conductor. Right: the conductor has been Figure 2. Left: linear and loop currents over a perfect conductor. Right: the conductor has been replaced by the image currents. Note that the sign reversal of the currents parallel to the conducting replaced by the image currents. Note that the sign reversal of the currents parallel to the conducting plane are responsible for the polarization change of the “reflected” wave. plane are responsible for the polarization change of the “reflected” wave.
FigureFigure 3. 3. GNSS-RGNSS-R scattering scattering geometry geometry in in Figure Figure 11 modifiedmodified according according to to the the principle principle of of equivalence: equivalence: thethe transmitter transmitter is is replaced replaced by by its its image. image. Computed Computed electric electric fields fields ar aree only only valid valid in in the the region region z z >> 0.0. 2.1.2. The Huygens Principle and Knife-Edge Diffraction 2.1.2. The Huygens Principle and Knife-Edge Diffraction The Huygens principle states that every point (red dots in Figure4) on a primary wave front The Huygens principle states that every point (red dots in Figure 4) on a primary wave front (blue curve in Figure4) may be regarded as a new source of spherical secondary waves, and that at (blue curve in Figure 4) may be regarded as a new source of spherical secondary waves, and that at each point in space these waves propagate in every direction at a speed and frequency equal to that of each point in space these waves propagate in every direction at a speed and frequency equal to that of the primary wave, and in such a way that the primary wave front at some later time is the envelope of these secondary waves (red curve in Figure 4).
Remote Sens. 2020, 12, x FOR PEER REVIEW 4 of 15
The Huygens principle can be applied to compute the electric fields when part of the trajectory is partially obstructed by a conducting plane as illustrated in Figure 5, leaving a clearance h with respectRemote Sens. to the2020 line, 12, 3910of sight. This academic case is known as the “knife-edge diffraction”, and it is 4used of 15 in radio communications to assess the propagation losses with respect to the free space caused by the Remoteobstruction Sens. 2020 by, 12 a, xmountain, FOR PEER REVIEW building, etc. In the latter this will represent the transition from4 ofone 15 the primary wave, and in such a way that the primary wave front at some later time is the envelope of medium to another [6]. theseThe secondary Huygens waves principle (red can curve be inapplied Figure to4). compute the electric fields when part of the trajectory is partially obstructed by a conducting plane as illustrated in Figure 5, leaving a clearance h with respect to the line of sight. This academic case is known as the “knife-edge diffraction”, and it is used in radio communications to assess the propagation losses with respect to the free space caused by the obstruction by a mountain, building, etc. In the latter this will represent the transition from one medium to another [6].
Figure 4. Graphical representation of the Huygens diffraction principle. Figure 4. Graphical representation of the Huygens diffraction principle. The Huygens principle can be applied to compute the electric fields when part of the trajectory is partially obstructed by a conducting plane as illustrated in Figure5, leaving a clearance h with respect to the line of sight. This academic case is known as the “knife-edge diffraction”, and it is used in radio communications to assess the propagation losses with respect to the free space caused by the obstruction by a mountain, building, etc. In the latter this will represent the transition from one medium to anotherFigure [6]. 4. Graphical representation of the Huygens diffraction principle.
Figure 5. Huygens principle applied to the “knife-edge diffraction” case. The light-gray area partially blocks the propagation of the signal, hiding parts of some Fresnel zones. h is the clearance of the obstruction.
The development of an analytical solution can be found in many text books, and only the final
result is provided below: FigureFigure 5. 5. HuygensHuygens principle principle applied applied to to the the “knife-edge “knife-edge di diffraction”ffraction” case. case. The The light-gray light-gray area area partially partially blocksblocks the the propagation propagation of of the the signal, signal, 𝐸hiding hiding=𝐸 parts parts of ofsome some·𝐹 Fresnel 𝑣 Fresnel zones. zones. h is htheis clearance the clearance of the of (2) the obstruction. where:obstruction. The development of an analytical solution can be found in many text books, and only the final The development of an analytical solution can be found in many text books, and only the final result is provided𝐹 𝑣 = below: 𝑐𝑜𝑠 𝑡 𝑑𝑡 − 𝑗 𝑠𝑖𝑛 𝑡 𝑑𝑡 = −𝐶 𝑣 −𝑗 −𝑆 𝑣 , (3) result is provided below: E = E F(v) (2) 𝐶 𝑣 and 𝑆 𝑣 are the so-called Fresneldi f fintegrals,f ree space which· are given be: 𝐸 =𝐸 ·𝐹 𝑣 (2) where: where: "Z Z # 1 + j ∞ π ∞ π 1 + j 1 1 F(v)=ˆ cos t2 dt j sin t2 dt = C(v) j S(v) , (3) 𝐹 𝑣 = 2 v 𝑐𝑜𝑠 2𝑡 𝑑𝑡− − 𝑗 v 𝑠𝑖𝑛 2𝑡 𝑑𝑡 = 2 −𝐶2 − 𝑣 −−𝑗 2−𝑆− 𝑣 , (3) 𝐶 𝑣 and 𝑆 𝑣 are the so-called Fresnel integrals, which are given be:
Remote Sens. 2020, 12, 3910 5 of 15
C(v) and S(v) are the so-called Fresnel integrals, which are given be: Remote Sens. 2020, 12, x FOR PEER REVIEW 5 of 15 Z ν π C(v) = cos t2 dt, (4a) 0 2 𝐶 𝑣 = 𝑐𝑜𝑠 𝑡 𝑑𝑡, (4a) Z ν π 2 S𝑆( v𝑣) == 𝑠𝑖𝑛sin 𝑡t dt𝑑𝑡,, (4b) (4b) 0 2 andand v𝑣is is the the Fresnel–Kirchho Fresnel–Kirchhoffff parameter: parameter: s 22·(R T𝑅++𝑅RR) v𝑣=ℎ·= h · (5(5)) · λ𝜆·𝑅RT R·𝑅R · · FigureFigure6 6 shows shows thethe absoluteabsolute valuevalue of 𝐹F (𝑣v ) inin Equation Equation (3) (3) in in linear linear units. units. When When h =h 0,= half0, half of the of theelectric electric field field is blocked, is blocked, which which translates translates into into a apower power loss loss of of 6 6 dB dB with with respect to thethe free-spacefree-space conditions.conditions. Note the ringing for h << 0 due to the didifferentfferent FresnelFresnel zoneszones beingbeing blockedblocked oror notnot andand contributingcontributing toto thethe totaltotal electricelectric fieldfield withwith didifferentfferent phases. phases.
Figure 6. Absolute value of F(v). Negative v indicates that the line-of-sight (LoS) trajectory is Figure 6. Absolute value of 𝐹 𝑣 . Negative 𝑣 indicates that the line-of-sight (LoS) trajectory is unobstructed, while positive v indicates an LoS obstruction (see Figure5). unobstructed, while positive 𝑣indicates an LoS obstruction (see Figure 5). 2.1.3. Babinet’s Principle 2.1.3. Babinet’s Principle Babinet’s principle [10] states that the diffraction pattern from an opaque body is identical to that fromBabinet’s a hole of theprinciple same size [10] and states shape, that exceptthe diffraction for the overall pattern forward from an beam opaque intensity. body The is identical concept isto expressedthat from graphicallya hole of the in Figuresame 7size[11 ].and shape, except for the overall forward beam intensity. The conceptBabinet’s is expressed principle graphically can be used in Figure to compute 7 [11]. the electric field diffracted when the signal emitted by the transmitter’s image is diffracted by the light-gray opaque area in the z = 0 plane (the white area being “transparent”), and when it is diffracted by the white opaque area in the z = 0 plane (the light-gray area now being “transparent”). Remote Sens. 2020, 12, x FOR PEER REVIEW 5 of 15