Next Article in Journal
Weighting of Multi-GNSS Observations in Real-Time Precise Point Positioning
Previous Article in Journal
The Benefits of the Ka-Band as Evidenced from the SARAL/AltiKa Altimetric Mission: Quality Assessment and Unique Characteristics of AltiKa Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Performance and Requirements of GEO SAR Systems in the Presence of Radio Frequency Interferences

1
School of Information and Electronics, Beijing Institute of Technology, Beijing 100081, China
2
Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milano, Italy
3
Key Laboratory of Electronic and Information Technology in Satellite Navigation (Beijing Institute of Technology), Ministry of Education, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2018, 10(1), 82; https://doi.org/10.3390/rs10010082
Submission received: 27 September 2017 / Revised: 1 January 2018 / Accepted: 3 January 2018 / Published: 9 January 2018

Abstract

:
Geosynchronous Synthetic Aperture Radar (GEO SAR) is a possible next generation SAR system, which has the excellent performance of less than one-day revisit and hundreds of kilometres coverage. However, Radio Frequency Interference (RFI) is a serious problem, because the specified primary allocation frequencies are shared by the increasing number of microwave devices. More seriously, as the high orbit of GEO SAR makes the system have a very large imaging swath, the RFI signals all over the illuminated continent will interfere and deteriorate the GEO SAR signal. Aimed at the RFI impact in GEO SAR case, this paper focuses on the performance evaluation and the system design requirement of GEO SAR in the presence of RFI impact. Under the RFI impact, Signal-to-Interference-plus-Noise Ratio (SINR) and the required power are theoretically deduced both for the ground RFI and the bistatic scattering RFI cases. Based on the theoretical analysis, performance evaluations of the GEO SAR design examples in the presence of RFI are conducted. The results show that higher RFI intensity and lower working frequency will make the GEO SAR have a higher power requirement for compensating the RFI impact. Moreover, specular RFI bistatic scattering will give rise to the extremely serious impact on GEO SAR, which needs incredible power requirements for compensations. At last, real RFI signal behaviours and statistical analyses based on the SMOS satellite, Beidou-2 navigation satellite and Sentinel-1 A data have been given in the appendix.

Graphical Abstract

1. Introduction

L band, C band and X band spaceborne microwave systems play important roles in lots of remote sensing applications, such as the surface deformation measurements of landslides and earthquakes and vegetation mapping of forests. Nevertheless, the available microwave frequencies allocated for the spaceborne systems, called the active primary allocation frequencies, are given by the International Telecommunications Union’s (ITU’s) Earth Exploration-Satellite Service (EESS) [1]. Active primary allocations within the microwave frequency extension from L band to X band are provided in Table 1 [1]. Beyond the primary allocation frequencies in different radio frequency bands, there are protected frequency bands or void frequency bands. Taking a classical protected frequency band in L band as an example, the frequency interval 1400–1427 MHz is protected, within which unwanted emissions of active service operations are forbidden and even the emissions from the nearby bands should be guaranteed to be lower than an expected value [2].
However, as more and more microwave devices are implemented, the primary allocation frequencies are becoming increasingly crowded. Different systems, which utilize the same or adjacent microwave frequency bands, will interfere with each other. Since only the active primary allocation frequencies can be used by the spaceborne microwave systems, the systems cannot avoid being impacted by the Radio Frequency Interference (RFI). Based on L band Aquarius scatterometers and the SMOS radiometer, strong RFI emissions, which might come from the terrestrial sources (e.g., military radars), have been detected [2,3,4]. Besides those detected by the radiometers, in 2016, Pascual et al. discussed the crosstalk RFI in the interferometric Global Navigation Satellite System Reflectometry (iGNSS-R), when the Delay-Doppler Map of the desired satellite is overlapped by other satellites [5]. Other detailed discussions about the classification of RFI sources, impacts analysis and mitigation methods in the global navigation satellite system (GNSS) system case were given in [6].
RFI is a common and serious problem in Synthetic Aperture Radar (SAR) as well. RFI decreases the SAR image quality by generating the noise-like artefacts, which is shown in Figure 1. RFI can also cause the saturation of the SAR echo and even produce spurious targets in the images [7]. Lots of RFI has been detected by SAR systems of different working bands. With respect to the L band SAR, the impact of RFI has been observed in ALOS-1 PALSAR data, which showed that RFI could not only cause haze-like image artefacts, linear pattern or blurs in the focused SAR image but reduce the capability of obtaining polarimetric signatures and the phase accuracies of SAR images [8]. Likewise, for C band systems, the RFI deriving from the Radio Local Area Networks (RLAN) could degrade the signals of the C band meteorological radars when the RFI power was comparable to that of the signal obtaining the nominal Signal-to-Noise Ratio (SNR) [9,10]. RFI signals could be observed affecting a few range lines in the C band Sentinel-1A SAR image [11]. Many low power RFI sources with the whole bandwidth and a few very high-power temporal RFIs are detected by Sentinel-1A SAR image as well [12]. Significantly, Radarsat-2 and Sentinel-1A have mutual interferences because they share the same working band and have the almost same equatorial crossing time and multiple repeat periods of each other. Although RFI signals in X band are not as obvious as those in C band, increased RFI mitigation capability is needed due to the wide designed bandwidth, as based on the analysis of the C- and X- band Windsat microwaves radiometer data [13,14]. Moreover, a small RFI modulated with a pulse wave has been observed in the X band TerraSAR-X system, which if notched out, will give rise to sidelobes in the SAR image [15]. Therefore, RFI is an important factor needed to be considered, including the future launched Geosynchronous orbit SAR (GEO SAR) system.
To improve the revisit and the coverage performance of the spaceborne SAR system, GEO SAR has been proposed [16]. Unlike a totally geostationary communication satellite, GEO SAR satellite runs in a geosynchronous orbit with non-zero inclination and eccentricity. Thus, a GEO SAR can obtain the relative movement with the targets on the ground, which ensures sufficient Doppler bandwidth after long time integration for realizing two-dimensional SAR imaging. It can achieve a revisit of the target region of one day and a wide beam footprint of more than hundreds of kilometres [17]. According to the inclination, the designed GEO SAR systems can be classified into the inclined orbit GEO SAR system (large inclination) and the quasi-geostationary system (small or zero inclination, small eccentricity but not zero) [18]. To satisfy the SNR requirement, inclined orbit GEO SAR system needs a higher power and a larger antenna size [19]. It achieves the designed azimuth resolution by a relative shorter integration time and it is sensitive to ionosphere and cluster disturbances [20,21,22,23,24,25,26,27,28]. In contrast, the geostationary system has lower requirements of the power and the antenna size [20]. Moreover, it has better revisit and observation performance for the target [29]. Nevertheless, it needs an integration time of hours. Thus, atmosphere distortion is a serious problem in this design [30]. GEO SAR has a longer integration time and a larger coverage compared with a Low Earth Orbit SAR (LEO SAR). Thus, it has an increased probability to be impacted by the RFI and the related impact will be more severe as well. Shown in Figure 2, the high orbit of GEO SAR gives rise to a large swath of about 400 km, while the swath of a LEO SAR is only around tens of kilometres in across track. Resultantly, more RFI signals will come into the beam of GEO SAR, compared with that in a LEO SAR case. Although the RFI impact depends on its time-frequency characteristics, such as the bandwidth and the transmitting duration, it is general noise-like, which can finally decrease the imaging quality of GEO SAR images by reducing the Signal-to-Interference-plus-Noise Ratio (SINR). Therefore, we need to consider and evaluate the RFI impact in the GEO SAR system design.
The paper is organized as follows. The RFI impact on the SINR and the required average transmitted power to compensate RFI impacts in the GEO SAR case are discussed in Section 2. In Section 3, focusing on the bistatic scattering RFI case, the RFI impact on SINR and the required transmitted powers in the specular scattering and non-specular scattering cases are studied. Moreover, under the assumption of the point-like target case, the Doppler filtering effect in the LEO SAR bistatic scattering RFI case is specially analysed. In Section 4, in the presence of the RFI, simulations are conducted to analyse the GEO SAR system design requirements in detail. Finally, Section 5 concludes this paper. In addition, some important real RFI measurement results in different frequencies are given in the Appendix.

2. Ground RFI Impacts on GEO SAR Signal

The ground RFI is the most common case which often derives from ground implemented radio electronic devices and ground-based radar systems. Thus, we analyse its impact on GEO SAR in this section. A sketch of a GEO SAR system in the presence of RFI is shown in Figure 3. The received RFI signal will raise the noise floor and degrade the Signal-to-Interference Ratio (SIR). Generally, the definition of the SIR is [32].
S I R = P s P R F I ,
where P R F I is the power of the RFI and P s is the power of the received effective GEO SAR signal. Based on the radar equation [29], P s can be written as
P s ( θ d , ψ ) = P a A t 2 F 2 ( θ d , ψ ) L l o s s 2 4 π λ 2 R 4 σ 0 ρ a ρ g ,
where P a is the average transmitted power of the GEO SAR signal, A t is its antenna aperture area, F represents the gain, θ d and ψ represent the pointing direction of the antenna, L l o s s is the one-way total loss, σ 0 is the backscatter coefficient of the target, λ is the wavelength, R is the slant range and ρ a and ρ g are azimuth resolution and the ground range resolution, respectively.
Considering the thermal noise, we can define the SINR as
S I N R = P s P R F I + P t h ,
where P t h is the power of the thermal noise, which can be written as [32]
P t h = k B T t h B g ,
where k B = 1.38 × 10 23 J / K is the Boltzmann constant, T t h is the thermal noise temperature and B g is the bandwidth of the GEO SAR system.

2.1. Point-Like RFI Source

If the ground RFIs are contributed by the point–like RFI sources (strong, localized and continuous emitting) [8,33], the powers of all the received point-source RFIs have to be summed together. Referring to Figure 3, the total power of RFI at the receiver of GEO SAR can be expressed as
P R F I = A t 4 π R 2 i = 1 N P e , i p e , i F ( θ e , i , ψ e , i ) ,
where P e is defined as the average Effective Isotropic Radiated Power (EIRP) of each single point-like RFI source, p is the probability of each point-like RFI source, i is the index, N is the number of the point-like RFI sources, θ e and ψ e indicates the position of each point-like RFI source. Because (5) relates to the antenna pattern, the position of the point-like RFI source will impact the power of P R F I . A high-power point-like RFI source located far away from the centre of the beam of GEO SAR may have a smaller impact than a lower power point-like RFI source which is near the centre of the beam of GEO SAR.

2.2. Distributed RFI Source

Practically, many RFI sources on the ground are similar and independent distributed within a fixed area (e.g., personal electronic devices). Therefore, in this case, the emission of the RFI is noise-like, which is very similar to a black-body.
According to the basic theory of the black-body in microwave spectrum, under the assumption that the antenna receives the full EM field illuminated by the black-body, the antenna directivity within the solid angle is uniform and the bandwidth is small compared with the carrier frequency, we have [32]
P R F I = k B B g T R F I ,
where T R F I is the brightness temperature contributed by the distributed RFI sources. As the antenna area and the area of the beam foot-print are in the inverse proportion, P I will be independent on the antenna area. It only depends on the brightness temperature contributed by the RFI and the bandwidth of the receiver.

2.3. Average Transmitted Power Requirements in the Presence of RFI

RFI will degrade the SINR performance of a fixed GEO SAR system (e.g., specified average transmitted power, antenna size). To meet the requirements for achieving the good performances of remote sensing applications, we should determine the proper average transmitted power.
Considering the incoherent RFI and combining (2), the Noise-Equivalent Sigma Zero (NESZ) N E S Z T can be expressed as
N E S Z T ( θ d , ψ ) = 4 π R 4 λ 2 k B ( T R F I + T t h ) P a A t 2 F 2 ( θ d , ψ ) L l o s s 2 T a ρ a ρ g ,
where, especially, in the multi point-like RFI source case,
T R F I = h I A t 4 π k B R 2 ,
is the RFI equivalent brightness temperature, h I = i = 1 N P e , i p e , i F ( θ e , i , ψ e , i ) B I , i is the power spectrum density, B I is the bandwidth of the point-like RFI source and T a is the integration time of GEO SAR. In the distributed RFI source case, T R F I is the brightness temperature caused by the noise-like RFI.
Based on (7), in both the point-like RFI source case and the distributed RFI source case, the SINR will have the following relationship in dB
S I N R = σ 0 N E S Z T ,
According to (9), with respect to the fixed target (specified σ 0 ), to satisfy the SINR requirement for SAR applications with a good performance (e.g., to keep the SINR above 10 dB in SAR interferometry applications for making the coherence of the interferogram above 0.9 [34]), we should design the average transmitted power of the GEO SAR system as
P a r = 4 π λ 2 R 4 k B ( T R F I + T t h ) S I N R r A t 2 L l o s s 2 T a ρ a ρ g σ 0 ,
where P a r is the required average transmitted power of GEO SAR and S I N R r is the required SINR.

3. Bistatic Scattering RFI Impact on GEO SAR

3.1. Impact of Bistatic Scattering RFI on GEO SAR

Besides the RFI impacts coming from the ground RFI directly, some other received RFI signals may derive from the signals emitted from spaceborne radio systems (e.g., LEO SAR and GNSS satellite system) with the same or adjacent radio frequencies and bistatic scattering from the ground [6,35,36,37]. The sketch of the impact of the bistatic scattering RFI case on GEO SAR is shown in Figure 4. P i a is the average transmitted power of the spaceborne radio system (RFI source), O represents the position of the target, L is the position of the spaceborne system, R s is the slant range of the spaceborne systems, R ^ s is the unit vector from O to L , S is the position of GEO SAR, R is the range of the GEO SAR, R ^ is the unit vector from O to S , θ 1 is the incidence angle of the spaceborne system and θ 2 is the corresponding RFI bistatic scattering angle to GEO SAR, z ^ is the normal unit vector of the scene, x ^ is the unit vector of the projection of R ^ s on the scene plane, y ^ is determined by the cross product of orthonormal basis z ^ and x ^ , ψ s is the out-of-plane angle of the bistatic scattering direction, β is the bistatic angle which relates to θ 1 , θ 2 and ψ s , σ B 0 is the bistatic scattering coefficient, θ R F I is the beam width of the spaceborne system, θ G E O is the beam width of the GEO SAR.

3.1.1. Non-Specular Bistatic Scattering Case

Generally, if the bistatic scattering signal does not form a specified angle to produce the specular scattering case, non-specular scattering occurs (out-of-plane case ( θ 1 = θ 2 , ψ s 0 ) and in-plane case ( θ 1 θ 2 , ψ s = 0 )). For the out-of-plane case, σ B 0 turns minimum when ψ s approaches 90° and σ B 0 is almost same with that in the in-plane case when ψ s is near 0° or 180° [38]. Thus, if we consider the in-plane case for the worst condition, based on the constant- γ bistatic-scatter-region model [38], in-plane land clutter bistatic scattering coefficient can be approximately expressed as
σ B 0 = γ ( sin θ 1 sin θ 2 ) 1 2 ,
where γ is the normalized reflectivity parameter, which can be estimated by
γ = σ M 0 sin θ 1
where σ M 0 is the monostatic scattering coefficient.
Regardless of the large-bistatic-angle case (e.g., forward-scatter region), we simply focus on the bistatic scattering case with the relative small bistatic angle for the above analysis. Such as for ALOS-2 PALSAR, the bistatic angle formed by the LEO SAR and the GEO SAR is general less than 30° when the LEO SAR signal is received by the GEO SAR within its antenna mainlobe. For the LEO SAR signal received outside the GEO SAR antenna mainlobe, as there will be a large attenuation for the RFI signal, we ignore the large-bistatic-angle case. Moreover, although the GNSS system and the GEO SAR can generate the bistatic system with large bistatic angle, because σ B 0 is generally smaller than σ M 0 [38], a small bistatic angle assumption is better for analysing the worst cases.
In the non-specular case, the received RFI signal is weak as the gain in the scattering direction of the GEO SAR is low.

3.1.2. Specular Bistatic Scattering Case

At times, if θ 1 = θ 2 and ψ s = 0 , a specular reflection happens. In this case, the power of the RFI will concentrate in the scattering direction of GEO SAR and the received RFI signal from the spaceborne system is very strong. In the specular scattering case, like a mirror reflection [32,39], σ B 0 can be expressed as
σ B 0 = 4 π A g s λ 2 ,
where A g s is the specular scattering area on the ground. Here and in the following part, we also use λ as the wavelength of the spaceborne radio system for simplicity. As the targets within the same range resolution cannot be distinguished by the range delay difference, ideally, the maximum area of A g s can be approximately determined within the range resolution bin ρ r by the cross-section area between the ellipsoid formed by GEO SAR and the spaceborne system and the Earth ellipsoid, which has been shown in Figure 5.
Based on the above analysis, the considerably different behaviour of the two cases derives from the different bistatic scattering coefficients which are shown in (11) and (13). In the specular scattering case, the bistatic scattering coefficient is large, while it will be very small in the non-specular scattering case.

3.1.3. SINR and Required Power Analysis

In the bistatic scattering RFI case, the total brightness temperature T R F I _ s a can be shown as
T R F I _ s a = 1 k B k = 1 N b 1 T o T w ( P i a , k ( θ R F I R s ( t ) ) 2 σ B , k 0 ( t ) A s A t 4 π R ( t ) 2 ) p b , k L l o s s b , k F ( θ d i , k , ψ i , k , t ) d t B I , k ,
where P i a , k is the average transmitted power of the kth spaceborne system, A s is the coverage area of the spaceborne system, p b , k is the corresponding illuminating probability for the bistatic scattering case, T o is the time period of one orbit for the spaceborne system, T w is the operation time of the spaceborne system per orbit, L l o s s b is the transmitting loss of the RFI signal, t is the slow time, k is the index and N b is the number of the spaceborne systems. The integration in (14) refers to the variation of the bistatic scattering coefficient and the gain of GEO SAR towards the spaceborne system within the integration time of GEO SAR.
A s changes in different cases, when θ R F I > θ G E O (e.g., GNSS system) or θ R F I < θ G E O (e.g., LEO SAR system), which can be shown as
A s = { ( θ G E O R ) 2 , θ R F I > θ G E O ( θ R F I R s ) 2 , θ R F I < θ G E O .
Based on (6) and (14), the SINR in the bistatic scattering RFI case can be written as
S I N R ( θ d , ψ ) = 4 π λ 2 R 4 k B ( T t h + T R F I _ s a ) P a A t 2 F 2 ( θ d , ψ ) L l o s s 2 T a ρ a ρ g σ 0 .
Likewise, being the same with (10), we can use (16) to obtain the average transmitted power requirement of GEO SAR P a r for meeting the required SINR in the presence of the spaceborne system RFI.

3.2. Doppler Filtering

Specially, if our interesting target is point-like and isolated, the received pulses will have deterministic phases. In this case, Doppler filtering on RFI signal will have a good performance to weaken the RFI impact. Figure 6a shows the geometry difference between GEO SAR and LEO SAR. v s and v g are the velocity vectors of the LEO SAR system and the GEO SAR system, respectively. Vectors R s and R represent the incidence direction (line of sight (LOS)) and the bistatic scattering direction, respectively. ϑ is the LOS angle between v s and R s . GEO SAR has a larger beam footprint (e.g., typically 400 km) than that of LEO SAR (e.g., typically 10 km in along-track). Moreover, the velocity of GEO SAR is small compared with that of a LEO SAR. Resultantly, the Doppler rate of GEO SAR is very small compared with that of a LEO SAR. To achieve the same azimuth resolution, the Doppler bandwidth of GEO SAR is much smaller (only dozens of Hz at most) than the Doppler bandwidth of a LEO SAR (hundreds to thousands Hz), which is shown in Figure 6b. Moreover, as the pulse repetition frequency (PRF) of GEO SAR is still smaller than the Doppler bandwidth of a LEO SAR, the Doppler frequency of a LEO SAR will be aliased. When doing the Doppler filtering in the imaging processing, a GEO system achieves the azimuth focusing by exploiting the small Doppler bandwidth. Therefore, the components of the LEO SAR RFI Doppler spectrum out of the Doppler bandwidth of GEO SAR could be filtered, as shown as the yellow regions in Figure 6b. It should be noticed that the aliased Doppler frequencies of the LEO SAR within the GEO SAR Doppler bandwidth cannot be filtered out.
As the integration time of GEO SAR (hundreds of second or even hours) is far longer than that in a LEO SAR (typically 1 s), if we assume the LEO SAR observes the target within the integration time of the GEO SAR in a worst case (like a staring mode), the LEO SAR signal can be simplified as weighted by the envelope of GEO SAR. Moreover, the worst geometry relationship between the GEO SAR and LEO SAR is assumed as well, which makes the RFI signal of LEO SAR be weighted by the GEO SAR antenna pattern. Under the assumption, the received RFI signal from the LEO SAR is expressed as
s l ( t ) = A sin c 2 ( v g L g λ R 0 t ) exp ( j 4 π R s ( t ) λ ) = A sin c 2 ( v g L g λ R 0 t ) exp [ j π ( 4 R 0 , l λ + f d c , l t + f d r , l t 2 + ) ] ,
where A (relates to SINR) is the amplitude of the RFI, L g is the antenna size of GEO SAR, v g is the velocity of GEO SAR, R 0 , l is the shortest slant range of R s , f d r , l is the Doppler rate of the LEO SAR, f d c , l is the Doppler centroid frequency of the LEO SAR. Sinc function represents the envelope of the GEO SAR signal.
The received GEO SAR signal is
s l ( t ) = sin c 4 ( v g L g λ R 0 t ) exp ( j 4 π R ( t ) λ ) = sin c 4 ( v g L g λ R 0 t ) exp [ j π ( 4 R 0 λ + f d c , g t + f d r , g t 2 + ) ] ,
where R 0 is the shortest slant range of GEO SAR, f d r , g is the Doppler rate of the GEO SAR, f d c , g is the Doppler centroid frequency of the GEO SAR.
The energy of the LEO SAR RFI signal E after Doppler filtering is
E = A 2 T a 2 T a 2 | T a 2 T a 2 sin c 2 ( v g L g λ R 0 u ) exp [ j π ( 4 λ R 0 , l + f d c , l u + f d r , l u 2 ) ] r e c t ( t u T a ) exp [ j π f d r , g ( t u ) 2 ] d u | 2 d t .
Expanding the exponent terms and combining them, ignoring the impact of the antenna pattern for simplicity, we have
E A 2 | f d r , l f d r , g | T a 2 T a 2 | exp [ j π ( f d r , l f d r , g ) ( f d c , l f d c , g 2 + f d r , g t ) 2 ] | 2 d t A 2 T a | f d r , l f d r , g | .
In the meantime, the energy of the effective GEO SAR signal E m after Doppler filtering is expressed as
E m = 1 f d c , g T a 1 f d c , g T a | T a 2 T a 2 sin c 4 ( v g L g λ R 0 u ) exp [ j π ( 4 λ R 0 + f d c , g u + f d r , g u 2 ) ] r e c t ( t u T a ) exp [ j π f d r , g ( t u ) 2 ] d u | 2 d t .
Simplifying the equation, we have
E m 0.9 T a | f d r , g | .
The ratio between the energy of the GEO SAR signal and the energy of the RFI after Doppler filtering is shown as
= E m E 0.9 A 2 | f d r , l f d r , g | | f d r , g | P R F g f d r , l T a ,
where P R F g is the pulse repetition frequency of GEO SAR. If GEO SAR uses a small PRF, the Doppler spectrum of the RFI will be aliased, leading to a small value of , that, in turn results in a small SINR. The difference of the Doppler rates of GEO SAR and LEO SAR is the parameter to describe the attenuation of LEO SAR RFI signal after Doppler filtering. A higher difference between the Doppler rates of GEO SAR and LEO SAR will be helpful to raise . A higher indicates a better imaging performance of GEO SAR in the presence of LEO SAR RFI.
In the case of the distributed target, the LEO SAR illuminates it with a PRF that causes the backscattered returns to be independent, white. The LEO-GEO bistatic scattering will then be nearly white and performance of the Doppler filtering will lose.

4. Performance Evaluation and Discussion: Example of GEO SAR Design

In this section, we evaluate the performance of GEO SAR in the presence of RFI by using System Tool Kit (STK) software to generate the GEO SAR orbit, LEO SAR orbits and the GNSS orbits for the simulation. The High-Precision Orbit Propagator (HPOP) in STK is adopted to obtain more accurate orbit simulation. The input parameters of our inclined orbit GEO SAR system are given in Table 2 [19,21]. We consider the classical “Figure 8” inclined GEO SAR orbit with the small eccentricity and the large inclination. The designed GEO SAR system can work in three frequency bands (L, C and X). The large antenna size is utilized to obtain the higher gain of the GEO SAR signal.
The performance of a GEO SAR defined in Table 2, is shown in Table 3. We assume that both the antenna size and the swath are kept constant with different frequencies (from L band to X band), by using scan SAR mode. From Table 2, when the working frequency of GEO SAR increases from the L band to X band, the number of beams will accordingly change from 1 to 64 to keep the coverage requirement. We design the basic SINR requirement as 10 dB to make the system reach the imaging requirement for interferometric applications. Proper signal bandwidth and integration times make the GEO SAR system have range and azimuth resolutions better than 20m.

4.1. Impact of Ground RFI

In the simulation of the ground RFI case, the location of the beam centre of the GEO SAR is set in Beijing (40°N, 116°E), when the GEO SAR is at the orbit position of 162° true anomaly with a 2.6° look angle. Moreover, we assume the one-way signal loss of GEO SAR is 3 dB. According to the analysis of the real RFI data obtained by the spaceborne systems in the appendix, as the RFI brightness temperature is general thousands of Kelvin, the scope of the RFI brightness temperature in the analysis is from 0 to 10,000 K.
According to the SINR requirement in Table 3, the average transmitted power requirement of the GEO SAR system with the variation of RFI intensity and the distance from the position of the point-like RFI source to the beam centre of GEO SAR is shown in Figure 7. The average transmitted power requirements in all three working frequencies will increase with the increasing of the brightness temperature. In the meantime, the average transmitted power requirements are also very sensitive to the relative distance between the beam centre of GEO SAR and the position of the point-like RFI. When the distance increases, the average transmitted power requirement will decrease, because the power of the RFI after weighting by the antenna pattern of GEO SAR becomes small. Taking the L band system and the brightness temperature of 5000 K as an example, when the distance changes from 0 to 60 km, the average transmitted power requirements varies from more than 2.4 kW to less than 2.3 kW. At last, because of the shorter wavelength, the smaller beam width and the higher backscatter coefficient, comparing the three sub-figures in Figure 7, a higher frequency is helpful to reduce the required average transmitted power. Under the case of a 5000 K brightness temperature RFI and locating at 60 km away from the centre of the GEO SAR beam, the average transmitted power requirement changes from about 2.3 kW to only 0.2 W from L band to X band. Thus, an L band GEO SAR system needs a higher transmit power to compensate the RFI impact.
Likewise, shown in Figure 8, in the case of distributed RFI, the average transmitted power requirement will increase when the RFI brightness temperature turns higher. In the meantime, similarly, the GEO SAR system with a higher working frequency also has a much lower average transmitted power requirement in the distributed RFI case, which suggests the system will be easier to be implemented. Nevertheless, as the interference is uniform, the impact does not depend on the positions of the RFI and the beam centre of GEO SAR.

4.2. Impact of Bistatic Scattering RFI

4.2.1. LEO SAR Bistatic Scattering RFI Case Analysis

In this part, we analyse the bistatic scattering RFI impacts of LEO SAR on GEO SAR. In Table 4 [40,41,42,43,44,45], the utilized system parameters of the main L band, C band and X band LEO SAR systems in the simulation are given, including most important parameters about the power and the orbit. It can be noted that the number of the LEO SAR systems is large, because some of the systems are formed by not only a single satellite but many satellites as a constellation.
The orbits of the ALOS-2 PALSAR, Sentinel-1 A and TerraSAR-X have been used as the examples of L-, C- and X band LEO SAR systems to calculate the occurrence probabilities of non-specular scattering case, specular scattering case and mainlobe to mainlobe specular scattering case by STK software. The analysis result of the occurrence probabilities is provided in Table 5. The occurrence probabilities are obtained based on the orbit simulation of the GEO SAR and the LEO SAR with different periods, which are the repeat orbit cycles of the three systems (14-day, 12-day and 11-day, respectively). For all three considered working frequencies, because the GEO SAR can often ‘see’ the LEO SAR due to its high running orbit, the non-specular scattering case has a high occurrence probability of about 40%. In contrast, due to the strict condition of the occurrence of the specular scattering case, the occurrence probabilities in all working frequencies are very low, less than 0.1%. In fact, the calculation for the specular scattering case gives out the worst case, because the practical occurrence of the specular scattering is also determined by many factors, such as the incidence angle and the type of the scene. As the beam width of the LEO SAR turns smaller with the increasing of the frequency, the occurrence probability decreases further for the X band system compared with the L band or C band systems. Owing to the very little occurrence probability, mainlobe to mainlobe specular scattering case has never happened during the simulation.
The average transmitted power requirements in GEO SAR in the presence of the LEO SAR bistatic scattering RFI are shown in Table 5. Some geometry parameters used in the simulation are given in the caption. To analyse the worst case, we consider the received LEO SAR signal is weighted by the envelope of the GEO SAR one dimensional antenna pattern within the integration time. According to their working frequencies, the LEO SAR satellites in Table 5 are classified and used as the input RFI sources. To simplify the analysis, different satellites are assumed having the same illuminating geometry for GEO SAR.
Under the non-specular reflection case, the average transmitted power requirements in all three frequencies of GEO SAR are very low. Even for an L band GEO SAR, the required average transmitted power is only 0.7 kW. However, although the probability of the specular scattering case is very small, very large average transmitted power requirements are needed for GEO SAR to reach the designed SINR regardless of the working frequencies when the specular scattering case happened. Especially for the L band system, average transmitted power requirement even turns to 1.1 × 109 kW, which is impossible to be achieved. Deriving from the higher gain, the stronger backscatter and lower occurrence probability of the specular scattering case, a high frequency GEO SAR system will be more resistant to LEO SAR bistatic scattering RFI impact. Moreover, the results also depend on the distance between the LEO SAR beam centre and the GEO SAR beam centre. When the distance increases, the average transmitted power requirements will decrease. When the distance increases to about 5000 km, the average transmitted power requirement turns to about 2.4 × 104 kW rather than 1.1 × 109 kW for an L band GEO SAR. In the non-specular scattering case, as the RFI signal is very weak, the distance has the negligible impact on the design requirement.

4.2.2. GNSS Bistatic Scattering RFI Case Analysis

GNSS has been designed to work only in L band. Nevertheless, although Global Positioning System (GPS), Global Navigation Satellite System (GLONASS) and Galileo satellites are Medium Earth orbit (MEO) satellites, Beidou-2 satellites have different orbit types, including the MEO, geostationary orbit (GEOS) and inclined geosynchronous orbit (IGSO). Some important utilized system parameters of the main GNSS are given in Table 6 [46,47,48,49,50]. Orbiton software is used to calculate the number of GNSS satellites simultaneously in the sky above Beijing.
GPS satellite, Beidou-2 IGSO satellite and GEOS satellite are used as the examples of MEO satellite, IGSO satellite and GEOS satellite to calculate the occurrence probabilities of non-specular scattering case, the specular scattering case and the mainlobe to mainlobe specular scattering case in the GNSS bistatic scattering RFI case. The occurrence probabilities are given in Table 7. We obtain the occurrence probabilities by virtue of 10-day orbit simulation of GEO SAR and GNSS satellites. Because of the higher orbit and the wider beam width of GNSS, GNSS bistatic scattering RFI cases have very high occurrence probabilities. Probabilities for non-specular and specular scattering cases in GEOS and IGSO cases are 100%. The MEO satellite even has nearly 100% for non-specular and the specular scattering cases. Because of the very large beam width of the GNSS antenna, mainlobe to mainlobe specular scattering case can even occurrence.
The average transmitted power requirements in GEO SAR in the presence of GNSS bistatic scattering RFI is shown in Table 7. The geometry parameters used in the simulation are marked in the caption. As the mainlobe to mainlobe specular scattering case is the worst case, we discuss this case. The GNSS satellites given in Table 6 are used as the input RFI sources. To simplify the analysis, different satellites are regarded to have the same illuminating geometry for GEO SAR.
From Table 7, it can be noted that the required average transmitted power is very small in the non-specular scattering case and the GNSS bistatic scattering RFI has little impact on GEO SAR. However, because of the large number of the GNSS satellites and the high gain, much worse in the specular scattering case, the required average transmitted power is very high in the specular case, if we consider the maximum Ags (about 30 km wide). In fact, due to the roughness, the specular area can be very small (e.g., 100 m × 100 m), the required average transmitted power of GEO SAR is less than 1 kW, which can be achieved to mitigate the RFI impact.

4.2.3. Doppler Filtering Analysis

As we discussed in Section 3.2, we can conduct the Doppler filtering for the LEO SAR signal. Shown in Figure 9 (assuming the power of the LEO SAR signal is the same with that of the GEO SAR signal), a large Doppler difference bias will be helpful for reducing the LEO SAR RFI energy after imaging. When the Doppler rate difference is large, shown in Figure 9a, the LEO SAR signal after filtering is noise-like with only about −40 dB level. Otherwise, the peak of the LEO SAR signal after filtering will be comparable with the amplitude of the GEO SAR signal shown in Figure 9b. The relationship between the various LOS angles, the PRF of the GEO SAR, the working frequency and the ratio are given in Figure 10. From the figure, we can find that a larger LOS angle, a higher PRF of the GEO SAR and the lower working frequency are the effective methods to improve the performance of Doppler filtering in the presence of the LEO SAR RFI. The improvement derives from that they are helpful to filter the RFI energy outside the GEO SAR Doppler bandwidth by increasing the difference of the Doppler spectra between the LEO SAR and GEO SAR, as well as decreasing the alias. However, as the PRF and the working frequency also depend on the imaging swath and applications, they should be traded-off during the GEO SAR system design. In the meantime, LOS angle is determined by the geometry of the LEO SAR system. Because the LEO SAR generally works in broadside imaging, LOS angle is often large enough for Doppler filtering.

5. Conclusions

In this paper, the analysis of SINR and the average transmitted power requirement of GEO SAR system in the presence of RFI impact are discussed. Two RFI cases, the ground RFI and the bistatic scattering RFIs from the LEO SAR and GNSS systems, are considered.
In the ground RFI case, the average transmitted power requirement depends on the RFI intensity and system frequency. A lower RFI intensity and a higher frequency can help to reduce the required average transmitted power, which makes the system more practical to be implemented. Moreover, for the point-like RFI sources, the relative distance between the beam centre of GEO SAR and the position of the point-like RFI also affects the required average transmitted power of GEO SAR. The average transmitted power requirement will decrease when the position of the point-like RFI is far away from the position of the GEO SAR beam centre.
With respect to the bistatic scattering RFI, in LEO SAR RFI case, the non-specular scattering case has a high occurrence probability of about 40%, while the occurrence of the specular scattering case is very low, less than 0.1%. Under the non-specular scattering case, no special requirements for average transmitted power are needed. However, when the specular scattering case happens, we can obtain the achievable average transmitted power to reach the required SINR only for C and X band systems when the LEO SAR beam centre is thousands of kilometres away from GEO SAR beam centre. For the point-like target, Doppler filtering will be an effective method to resist the LEO SAR RFI when there is a large Doppler rate difference between GEO SAR and LEO SAR.
In the GNSS RFI case, it has very high probabilities (even 100%) of the non-specular and the specular scattering cases due to its high orbit and the wide beam width. The required average transmitted power is still very small in the non-specular scattering case, while an impossible average transmitted power requirement is needed in the specular scattering case.
In the future, some RFI mitigation methods based on post signal processing are needed to address the RFI impact, especially for the ground RFI with a very high intensity and the specular scattering cases for the bistatic scattering RFI.

Acknowledgments

This work is supported by National Natural Science Foundation of China (Grant no. 61471038, 61501032, 61427802), Beijing Natural Science Foundation (Grant no. 4162052) and Chang Jiang Scholars Program (T2012122) and 111 project of China under Grant B14010.

Author Contributions

Y.H. Li and A. Monti Guarnieri conceived and designed the methods and wrote the paper; Y.H. Li performed the simulation and analysed the data. F. Rocca commented about the paper and helped editing. C. Hu gave many suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Measuring and Modelling RFI

Appendix A.1. Radio Frequency Interferences in L Band

Appendix A.1.1. SMOS Satellite RFI Data

SMOS satellite is an L band (1413.5 MHz central frequency) instrument to monitor the soil moisture and ocean salinity [4]. Importantly, SMOS can provide the land brightness temperature file which can be used to estimate the radiation of the RFI on the ground.
We can utilize the SMOS data to evaluate the impacts of RFI in China. The SMOS L3 products, the RFI probability product and the Level 1C Browse product over Land (BWLD1C), are utilized [51]. In Figure A1a, RFI probability in parts of China during 15 February to 16 March 2011 shows that China suffers high probability RFI impact. Most parts of China are interfered by the RFI (with the probability 1) and only a small part in the north of China has a relative lower RFI probability. L1C browse products contain brightness temperature values per pixel at a fixed incidence angle of 42.5, which are identical to the full L1C products. The SMOS BWLD1C on 28 February 2011 is utilized in the following analysis.
Figure A1. (a) RFI probability in parts of China (15 February–16 March 2011); (b) Detected distribution of the RFI sources based on the cluster selection.
Figure A1. (a) RFI probability in parts of China (15 February–16 March 2011); (b) Detected distribution of the RFI sources based on the cluster selection.
Remotesensing 10 00082 g0a1
On the basis of the above data, we study the SINR. The detected distribution of the RFI sources based on the cluster selection RFI sources positioning [2] is given in Figure A1b, in which we have detected 161 RFI sources. We select the RFI sources located within the mainlobe of GEO SAR and obtain the average brightness temperature, which is 5203 K. Based on the brightness temperature and a 2250 w average transmitted power of GEO SAR, the calculated SINR is given in Figure A2. The maximum SINR will be 9.5 dB, which suggests the RFI make the system performance deteriorate slightly (compared with a 10 dB SINR requirement).
Figure A2. GEO SAR SINR simulation based on the SMOS data (28 February 2011).
Figure A2. GEO SAR SINR simulation based on the SMOS data (28 February 2011).
Remotesensing 10 00082 g0a2

Appendix A.1.2. Beidou-2 IGSO Navigation Satellite RFI Data

Beidou-2 IGSO navigation satellite runs in the geosynchronous orbit and works in L band. The L band receiver is set at the top of the building in Beijing Institute of Technology, China. The building is marked as the red point in Figure A3. Considering the sheltering effect, the possible coming area of RFIs obtained by the receiver is within the area marked as the blue lines in Figure A3.
Figure A3. Possible coming area of RFI.
Figure A3. Possible coming area of RFI.
Remotesensing 10 00082 g0a3
We obtain the Beidou-2 IGSO navigation satellite from the No.2 IGSO satellite. When there is no RFI, the received signal is shown in Figure A4a. The corresponding spectrum is shown in Figure A4b. In the absence of RFI, the signal is random, deriving from thermal noise. In Figure A4c,d, it shows the RFI signal in time and frequency domains, respectively. An obvious periodic RFI emerges in the figure. The detected RFI is a kind of pulse interference. Its pulse width is about 100 us and the pulse repeated frequency is nearly 900 Hz. The bandwidth of the RFI is about 4 MHz (wide spectrum part in Figure A4d). The equivalent antenna temperature is about 1573 K. There are some other high amplitude RFI signals existing in the spectrum (e.g., 100 Hz bandwidth). However, as their bandwidths are very small, their energies are small compared with the noise floor.
Figure A4. Beidou-2 IGSO navigation satellite signal analysis. (a) signal without RFI in time domain; (b) signal without RFI in frequency domain; (c) signal with RFI in time domain; (d) signal with RFI in frequency domain.
Figure A4. Beidou-2 IGSO navigation satellite signal analysis. (a) signal without RFI in time domain; (b) signal without RFI in frequency domain; (c) signal with RFI in time domain; (d) signal with RFI in frequency domain.
Remotesensing 10 00082 g0a4
If we add the detected RFI of Beidou-2 IGSO navigation satellite data into the GEO SAR case with a 2250 w average transmitted power, the SINR in the presence of the RFI is shown in Figure A5. The maximum SINR is about 15.7 dB, which shows little impact on GEO SAR imaging (compared with a 10 dB SINR requirement).
Figure A5. GEO SAR SINR based on the Beidou-2 IGSO navigation satellite data.
Figure A5. GEO SAR SINR based on the Beidou-2 IGSO navigation satellite data.
Remotesensing 10 00082 g0a5

Appendix A.2. Radio Frequency Interferences in C Band

The analyses of RFI in C band have been carried out by exploiting Sentinel-1 data according to the procedure described in [12]. The procedure is general and can be applied worldwide, thanks to the systematic acquisition of S1. The evaluation of RFI EIRP transmitted in the direction of the satellite has been achieved by averaging RFI over a time interval of 10 pulses (6 ms) each 20 km. S1 receiver bandwidth is 40–70 MHz, depending on the swath and the thermal noise is about 800 °K.
As example of detection of RFI, achieved by processing 400 datasets, each of them spanning an area 250 × 250 km wide and acquired in 2016–2017 is shown in Figure A6. In the figure, all the detected RFI, affecting half of the datasets, are stacked. The vertical axis is the received power spectrum, normalized to S1 noise floor. The spectra show many distributed sources and a few, strong RFI. An evaluation of the RFI probability of power distribution with frequency is in Figure A7. In the figure on the left, each column is an estimate of the cumulated density of probability (CDF) that RFI exceeds a certain power threshold. Notice that in most of the cases, this probability of exceeding a value of around 3, that roughly corresponds to 2400 °K, is less than 1/100. In the same figure, on the right, the mean and the standard deviations are plotted for the RFI stacked in Figure A6 and accounting for the cases in which no RFI was detected. Here again, one can assume that the RFI level does not exceeds a threshold of 2400 °K. However, this holds only by assuming random distributed RFI and not—for example, for the very strong values observed in some specific geographic locations.
Figure A6. Stacking of all the RFI power spectra evaluated in Europe by Sentinel-1. Vertical axis is power, normalized to Sentinel-1 noise power spectrum.
Figure A6. Stacking of all the RFI power spectra evaluated in Europe by Sentinel-1. Vertical axis is power, normalized to Sentinel-1 noise power spectrum.
Remotesensing 10 00082 g0a6
Figure A7. Left: cumulated histograms (right) of RFI power distributions with frequency, normalized to S1 noise, estimated by processing data from all-Europe. Vertical scale is probability in log10. Right: Plots of average (above) and standard deviation (below) of the RFI mean normalized power for each frequency bin.
Figure A7. Left: cumulated histograms (right) of RFI power distributions with frequency, normalized to S1 noise, estimated by processing data from all-Europe. Vertical scale is probability in log10. Right: Plots of average (above) and standard deviation (below) of the RFI mean normalized power for each frequency bin.
Remotesensing 10 00082 g0a7
The case of point-like RFI needs different analyses. Several repeated acquisitions centred in location in Middle East where the strongest RFI are measured, have been exploited leading to the PDF and CDF shown in Figure A8. The plots show that strong RFI repeats systematically in the same bandwidth, that is however only a portion of the entire spectrum available. This suggests mitigation with an agile frequency allocation.
Figure A8. Histograms (left) and cumulated histograms (right) of RFI power distributions with frequency, normalized to S1 noise estimated by processing data coming from the highest RFI in Figure A6 averaged in different times. Vertical scale is probability in log10 (left) and linear (right).
Figure A8. Histograms (left) and cumulated histograms (right) of RFI power distributions with frequency, normalized to S1 noise estimated by processing data coming from the highest RFI in Figure A6 averaged in different times. Vertical scale is probability in log10 (left) and linear (right).
Remotesensing 10 00082 g0a8
To evaluate the RFI impact in C band GEO SAR with a 2250 w average transmitted power, we utilize the Sentinel-1 RFI data. Based on the Sentinel-1 RFI data in Beijing area, we integrate the RFI energy within the 400 km swath in the GEO SAR case. According to the previous analysis, the equivalent antenna temperature is about 2400 K (total energy received by the antenna). The GEO SAR SINR map by using the Sentinel-1 data is given in Figure A9. The maximum SINR is more than 24 dB. Under the proposed case, the impact of the C band RFI has no impact on GEO SAR imaging (compared with a 10 dB SINR requirement).
Figure A9. GEO SAR SINR based on Sentinel-1 data.
Figure A9. GEO SAR SINR based on Sentinel-1 data.
Remotesensing 10 00082 g0a9

References

  1. Ulaby, F.; Avery, S. A Strategy for Active Remote Sensing Amid Increased Demand for Radio Spectrum; The National Academies of Sciences, Engineering and Medicines, Ed.; The National Academies Press: Washington, DC, USA, 2015. [Google Scholar]
  2. Oliva, R.; Nieto, S.; Félix-Redondo, F. RFI Detection Algorithm: Accurate Geolocation of the Interfering Sources in SMOS Images. IEEE Trans. Geosci. Remote Sens. 2013, 51, 4993–4998. [Google Scholar] [CrossRef]
  3. Freeman, A.P.; Fischman, M.A.; McWatters, D.A.; Spencer, M.W.; Piepmeier, J.R. The Detection and Mitigation of RFI with the Aquarius L-Band Scatterometer. In Proceedings of the IEEE IGARSS, Boston, MA, USA, 7–11 July 2008. [Google Scholar]
  4. Oliva, R.; Daganzo, E.; Kerr, Y.H.; Mecklenburg, S.; Nieto, S.; Richaume, P.; Gruhier, C. SMOS radio frequency interference scenario: Status and actions taken to improve the RFI environment in the 1400–1427-MHz passive band. IEEE Trans. Geosci. Remote Sens. 2012, 50, 1427–1439. [Google Scholar] [CrossRef] [Green Version]
  5. Pascual, D.; Park, H.; Onrubia, R.; Arroyo, A.A.; Querol, J.; Camps, A. Crosstalk Statistics and Impact in Interferometric GNSS-R. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 4621–4630. [Google Scholar] [CrossRef]
  6. Musumeci, L. Advanced Signal Processing Techniques for Interference Removal in Satellite Navigation Systems. Ph.D. Thesis, Politecnico di Torino, Turin, Italy, 2014. [Google Scholar]
  7. Yang, L.; Zheng, H.F.; Feng, J.; Li, M.; Chen, J.Q. Detection and suppression of narrow band RFI for synthetic aperture radar imaging. Chin. J. Aeronaut. 2015, 28, 1189–1198. [Google Scholar] [CrossRef]
  8. Meyer, F.J.; Nicoll, J.B.; Doulgeris, A.P. Correction and Characterization of Radio Frequency Interference Signatures in L-Band Synthetic Aperture Radar Data. IEEE Trans. Geosci. Remote Sens. 2013, 51, 4961–4972. [Google Scholar] [CrossRef] [Green Version]
  9. Brandao, A.L.; Sydor, J. 5 GHz RLAN Interference on active meteorological radars. In Proceedings of the 2005 IEEE Vehicular Technology Conference, Stockholm, Sweden, 30 May–1 June 2005. [Google Scholar]
  10. Tercero, M.; Sung, K.W.; Zander, J. Impact of aggregate interference on meteorological radar from secondary users. In Proceedings of the 2011 IEEE Wireless Communications and Networking Conference, Cancun, Mexico, 28–31 March 2011; pp. 2167–2172. [Google Scholar]
  11. Meadows, P. Sentinel-1A and -1B Annual Performance Report 2016; MPC-0366; Collecte Localisation Satellite: Ramonville Saint-Agne, France, 2016. [Google Scholar]
  12. Monti-Guarnieri, A.; Giudici, D.; Recchia, A. Identification of C-Band Radio Frequency Interferences from Sentinel-1 Data. Remote Sens. 2017, 9, 1183. [Google Scholar] [CrossRef]
  13. John, J.T.; Chen, C.C.; O’Brien, A.; Smith, G.E.; McKelvey, C.; Ball, M.A.C.; Misra, S.; Brown, S.; Kocz, J.; Jarnot, R. The CubeSat Radiometer Radio Frequency Interference Technology Validation (CubeRRT) mission. In Proceedings of the IEEE IGARSS, Beijing China, 10–15 July 2016. [Google Scholar]
  14. Ellingson, S.W.; Johnson, J.T. A polarimetric survey of radio-frequency interference in C- and X-bands in the continental United States using WindSat radiometry. IEEE Trans. Geosci. Remote Sens. 2006, 44, 540–548. [Google Scholar] [CrossRef]
  15. Hamdadou, N.; Chen, J.; Kamel, H.; Hui, K. Bidirectional notch filter for suppressing pulse modulated radio-frequency-interference in SAR data. In Proceedings of the IEEE IGARSS, Quebec, QC, Canada, 13–18 July 2014. [Google Scholar]
  16. Tomiyasu, K. Synthetic aperture radar in geosynchronous orbit. In Proceedings of the Digest International IEEE Antennas Propagation Symposium, College Park, MD, USA, 15–19 March 1978. [Google Scholar]
  17. Madsen, S.N.; Edelstein, W.; Didomenico, L.D.; LaBrecque, J. A geosynchronous synthetic aperture radar; for tectonic mapping, disaster management and measurement of vegetation and soil moisture. In Proceedings of the IEEE IGARSS, Sydney, Australia, 9–13 July 2001. [Google Scholar]
  18. Hobbs, S.; Mitchell, C.; Forte, B.; Holley, R.; Snapir, B.; Whittaker, P. System design for geosynchronous synthetic aperture radar missions. IEEE Trans. Geosci. Remote Sens. 2014, 52, 7750–7763. [Google Scholar] [CrossRef]
  19. Bruno, D.; Hobbs, S.E.; Ottavianelli, G. Geosynchronous synthetic aperture radar: Concept design, properties and possible applications. Acta Aston. 2006, 59, 149–156. [Google Scholar] [CrossRef]
  20. Monti Guarnieri, A.; Hu, C. Geosynchronous and geostationary SAR: Face to face comparison. In Proceedings of the EUSAR, Hamburg, Germany, 6–9 June 2016. [Google Scholar]
  21. Hu, C.; Li, Y.H.; Dong, X.C.; Long, T. Optimal data acquisition and height retrieval in repeat-track geosynchronous SAR interferometry. Remote Sens. 2015, 7, 13367–13389. [Google Scholar] [CrossRef]
  22. Hu, C.; Li, Y.H.; Dong, X.C.; Cui, C.; Long, T. Impacts of temporal-spatial variant background ionosphere on repeat-track GEO D-InSAR system. Remote Sens. 2016, 8, 916. [Google Scholar] [CrossRef]
  23. Hu, C.; Li, Y.H.; Dong, X.C.; Wang, R.; Ao, D.Y. Performance Analysis of L-band Geosynchronous SAR Imaging in the Presence of Ionospheric Scintillation. IEEE Trans. Geosci. Remote Sens. 2017, 55, 159–172. [Google Scholar] [CrossRef]
  24. Hu, C.; Li, Y.H.; Dong, X.C.; Wang, R.; Cui, C.; Zhang, B. Three-Dimensional Deformation Retrieval in Geosynchronous SAR by Multiple-Aperture Interferometry Processing: Theory and Performance Analysis. IEEE Trans. Geosci. Remote Sens. 2017, 55, 6150–6169. [Google Scholar] [CrossRef]
  25. Recchia, A.; Monti Guarnieri, A.; Broquetas, A.; Leanza, A. Impact of scene decorrelation on geosynchronous SAR data focusing. IEEE Trans. Geosci. Remote Sens. 2016, 54, 1635–1646. [Google Scholar] [CrossRef]
  26. Hu, B.; Jiang, Y.; Zhang, S.; Zhang, Y.; Yeo, T. Generalized omega-K algorithm for geosynchronous SAR image formation. IEEE Geosci. Remote Sens. Lett. 2015, 12, 2286–2290. [Google Scholar] [CrossRef]
  27. Wang, C.; Chen, L.; Liu, L. A New Analytical Model to Study the Ionospheric Effects on VHF/UHF Wideband SAR Imaging. IEEE Trans. Geosci. Remote Sens. 2017, 55, 4545–4557. [Google Scholar] [CrossRef]
  28. Wang, C.; Zhang, M.; Xu, Z.W.; Chen, C.; Guo, L.X. Cubic phase distortion and irregular degradation on SAR imaging due to the ionosphere. IEEE Trans. Geosci. Remote Sens. 2015, 53, 3442–3451. [Google Scholar] [CrossRef]
  29. Monti-Guarnieri, A.; Broquetas, A.; Recchia, A.; Rocca, F.; Ruiz-Rodon, J. Advanced radar geosynchronous observation system: ARGOS. IEEE Geosci. Remote Sens. Lett. 2015, 12, 1406–1410. [Google Scholar] [CrossRef]
  30. Ruiz-Rodon, J.; Broquetas, A.; Makhoul, E.; Monti-Guarnieri, A.; Rocca, F. Geosynchronous SAR focusing with atmospheric phase screen retrieval and compensation. IEEE Trans. Geosci. Remote Sens. 2013, 51, 4397–4404. [Google Scholar] [CrossRef]
  31. SMOS Blog. Available online: http://www.cesbio.ups-tlse.fr/SMOS_blog/smos_rfi/?q=image/12894-proba-rfi-d-20170826-20170909-bb (accessed on 15 September 2017).
  32. Curlander, J.C.; McDonough, R.N. Synthetic Aperture Radar: Systems and Signal Processing; Wiley: New York, NY, USA, 1991. [Google Scholar]
  33. Meyer, F.J.; Nicoll, J.; Doulgeris, A.P. Characterization and extent of randomly-changing radio frequency interference in ALOS PALSAR data. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Vancouver, BC, Canada, 24–29 July 2011. [Google Scholar]
  34. Ferretti, A.; MontiGuarnieri, A.; Prati, C.; Rocca, F.; Massonnet, D. SAR Principles: Guideline for SAR Interferometry Processing and Interpretation; ESA: Noordwijk, The Netherlands, 2007; pp. 57–71. [Google Scholar]
  35. Rosen, P.A.; Hensley, S.; Le, C. Observations and mitigation of RFI in ALOS PALSAR SAR data: Implications for the DESDynI mission. In Proceedings of the IEEE Radar Conference, Rome, Italy, 26–30 May 2008. [Google Scholar]
  36. Kankaku, Y.; Osawa, Y.; Suzuki, S.; Watanabe, T. The Overview of the L-band SAR Onboard ALOS-2. In Proceedings of the Electromagnetics Research Symposium, Moscow, Russia, 18–21 August 2009. [Google Scholar]
  37. Huang, S.W.; Tsang, L.G. Bistatic scattering, backscattering and emissivities of randomly rough soil surfaces at L band based on numerical solutions of Maxwell equations of 3 Dimensional simulations. In Proceedings of the IEEE IGARSS, Honolulu, HI, USA, 25–30 July 2010. [Google Scholar]
  38. Willis, N.J. Bistatic radar, chapter 25. In Radar Handbook, 2nd ed.; Skolnik, M.I., Ed.; McGraw-Hill Education: New York, NY, USA, 1990. [Google Scholar]
  39. Skolnik, M. Introduction to Radar Systems; McGraw-Hill Education: New York, NY, USA, 2001. [Google Scholar]
  40. ESA Sentinel-1. Available online: https://sentinel.esa.int/documents/247904/349449/S1_SP-1322_1.pdf (accessed on 8 September 2017).
  41. EoPortal Directory. Available online: https://directory.eoportal.org/web/eoportal/satellite-missions/a/alos-2 (accessed on 8 September 2017).
  42. Microwave Radar and Radiometric Remote Sensing. Available online: http://mrs.eecs.umich.edu/sensors.html (accessed on 8 September 2017).
  43. Bordoni, F.; Younis, M.; Rodriguez-Cassola, M.; Prats-Iraola, P.; López-Dekker, P.; Krieger, G. SAOCOM-CS SAR imaging performance evaluation in large baseline bistatic configuration. In Proceedings of the IEEE IGARSS, Milan, Italy, 26–31 July 2015. [Google Scholar]
  44. EoPortal Directory. Available online: https://directory.eoportal.org/web/eoportal/satellite-missions/r/radarsat-2 (accessed on 8 September 2017).
  45. DLR. Available online: http://www.dlr.de/hr/Portaldata/32/Resources/dokumente/tdmx/Tandem-X-RADAR-2004.pdf (accessed on 8 September 2017).
  46. GPS Information. Available online: http://gpsinformation.net/main/gpspower.htm (accessed on 8 September 2017).
  47. Wikipedia. Available online: https://en.wikipedia.org/wiki/GLONASS (accessed on 8 September 2017).
  48. European GNSS Agency. Available online: https://www.gsa.europa.eu/galileo/programme (accessed on 8 September 2017).
  49. WEEBAU. Available online: http://weebau.com/satellite/B/beidou2.htm (accessed on 8 September 2017).
  50. Steigenberger, P.; Thoelert, S.; Montenbruck, O. GNSS satellite transmit power and its impact on orbit Determination. J. Geod. 2017. [Google Scholar] [CrossRef]
  51. SMOS Blog. Available online: http://www.cesbio.ups-tlse.fr/SMOS_blog/?page_id=4087 (accessed on 8 September 2017).
Figure 1. RFI in ENVISAT image before removal (left) and after removal (right) (data are courtesy of European Space Agency (ESA)).
Figure 1. RFI in ENVISAT image before removal (left) and after removal (right) (data are courtesy of European Space Agency (ESA)).
Remotesensing 10 00082 g001
Figure 2. Comparisons between the GEO SAR swath and the LEO SAR swath and their covered RFI signals (SMOS RFI occurrence probability data in East Asia for 20170902 ± 7 days [31], red regions stand for the high probability of RFI occurrence).
Figure 2. Comparisons between the GEO SAR swath and the LEO SAR swath and their covered RFI signals (SMOS RFI occurrence probability data in East Asia for 20170902 ± 7 days [31], red regions stand for the high probability of RFI occurrence).
Remotesensing 10 00082 g002
Figure 3. Sketch of GEO SAR system in the presence of RFI.
Figure 3. Sketch of GEO SAR system in the presence of RFI.
Remotesensing 10 00082 g003
Figure 4. Sketch of the bistatic scattering RFI geometry.
Figure 4. Sketch of the bistatic scattering RFI geometry.
Remotesensing 10 00082 g004
Figure 5. Formation of the specular scattering area.
Figure 5. Formation of the specular scattering area.
Remotesensing 10 00082 g005
Figure 6. (a) Geometry relationship between LEO SAR and GEO SAR; (b) Doppler bandwidths of GEO SAR and LEO SAR.
Figure 6. (a) Geometry relationship between LEO SAR and GEO SAR; (b) Doppler bandwidths of GEO SAR and LEO SAR.
Remotesensing 10 00082 g006
Figure 7. Average transmitted power requirements in GEO SAR systems of different working frequencies. (a) L band; (b) C band; (c) X band, to ensure a SINR better than 10 dB.
Figure 7. Average transmitted power requirements in GEO SAR systems of different working frequencies. (a) L band; (b) C band; (c) X band, to ensure a SINR better than 10 dB.
Remotesensing 10 00082 g007
Figure 8. Average transmitted power requirement in the distributed RFI case.
Figure 8. Average transmitted power requirement in the distributed RFI case.
Remotesensing 10 00082 g008
Figure 9. GEO SAR and LEO SAR RFI signals after Doppler filtering: (a) about 300 Hz/s2 Doppler rate difference; (b) about 0.024 Hz/s2 Doppler rate difference.
Figure 9. GEO SAR and LEO SAR RFI signals after Doppler filtering: (a) about 300 Hz/s2 Doppler rate difference; (b) about 0.024 Hz/s2 Doppler rate difference.
Remotesensing 10 00082 g009
Figure 10. (a) ratio and the LOS angle ϑ (L band); (b) average ratio (averaged by ϑ within 0°–90°) and the PRF (L band); (c) ratio and the working frequency ( ϑ = 90°).
Figure 10. (a) ratio and the LOS angle ϑ (L band); (b) average ratio (averaged by ϑ within 0°–90°) and the PRF (L band); (c) ratio and the working frequency ( ϑ = 90°).
Remotesensing 10 00082 g010
Table 1. Active EESS primary allocations in L, C and X bands.
Table 1. Active EESS primary allocations in L, C and X bands.
BandFrequency Range (MHz)Services in Band
L1215–1300Radio Navigation Satellite System (RNSS) and radiolocation
C5250–5570Aeronautical RNSS and radiolocation
X9300–9900Radiolocation and radio navigation
Table 2. Classical parameters of the inclined orbit GEO SAR system.
Table 2. Classical parameters of the inclined orbit GEO SAR system.
ParameterValueParameterValue
Carrier central frequency (GHz)1.25/5.4/9.6Inclination (°)10–60
Look angle (°)2–7Argument of perigee (°)270
Eccentricity0.07Antenna area (m2)531
Duty cycle (%)15Incidence angle (°)10–60
Bandwidth (MHz)18Integration time (s)30–250
PolarizationHHEquivalent thermal noise (K)879
Table 3. Performance of GEO SAR system.
Table 3. Performance of GEO SAR system.
ParameterValueParameterValue
Number of beams1 (L band)
16 (C band)
64 (X band)
SINR (dB)10
Resolution (m)<20 × <20 (Single look)Coverage (km)400 × 400
Table 4. Utilized system parameters of the main LEO SAR (L, C and X bands) in the simulation.
Table 4. Utilized system parameters of the main LEO SAR (L, C and X bands) in the simulation.
FrequencyL BandC BandX Band
SARALOS-2Tandem-LSAOCOM-1A/BSentinel-1A/BRadarsat-2RCMTanDEM-X/PAZCOSMO-SkyMedSAR-LupeKOMPSAT-5
Peak power (W)51003622310043682280160022607600250 (total)1700
Duty cycle of power (%)845121213.75181135
Orbit period (min)98.599.797.298.5100.796.495.020% (max)20% (assumed)95.8
Working time per orbit (min)49305–2025, 74 (wave mode)281532
Satellite number1222133451
Average transmitted power (W)20344323947634131675013
Bandwidth (MHz)42 (middle)8450100100100150/300400300 (assumed)120
Table 5. Occurrence probability of the LEO SAR bistatic scattering RFI and the average transmitted power requirement in GEO SAR (HH polarization, Shrub scene, in the specular scattering case: θ 1 = θ 2 = 30 , in the non-specular case: θ 1 = 30 , θ 2 = 45 ).
Table 5. Occurrence probability of the LEO SAR bistatic scattering RFI and the average transmitted power requirement in GEO SAR (HH polarization, Shrub scene, in the specular scattering case: θ 1 = θ 2 = 30 , in the non-specular case: θ 1 = 30 , θ 2 = 45 ).
ParametersLCX
Occurrence probability (%)Non-specular scattering44.838.139.0
Specular scattering5.0 × 10−21.5 × 10−22.0 × 10−3
Mainlobe to mainlobe specular scattering000
σ0 (dB)−14.8−9.7−7.6
Required average transmitted power (kW)Non-specular scattering00.76.0 × 10−22.0 × 10−2
5000 km0.75.2 × 10−21.8 × 10−2
Specular scattering01.1 × 1093.2 × 1074.7 × 105
5000 km2.4 × 1047.83.7 × 10−2
Mainlobe to mainlobe specular scattering---
Table 6. Utilized system parameters of the main GNSS in the simulation.
Table 6. Utilized system parameters of the main GNSS in the simulation.
Navigation SystemGPSGLONASSGalileoBeidou-2
Transmit power (W)240135 265130 (MEO)
185 (IGSO/GEO)
Satellite number above the scene131186 (GEO)
8 (IGSO)
5 (MEO)
Bandwidth (MHz)210420
Table 7. Occurrence probability of the GNSS bistatic scattering RFI case and the average transmitted power requirement in GEO SAR (HH polarization, Shrub scene, in the specular scattering case: θ 1 = θ 2 = 30 , in the non-specular case: θ 1 = 30 , θ 2 = 45 ).
Table 7. Occurrence probability of the GNSS bistatic scattering RFI case and the average transmitted power requirement in GEO SAR (HH polarization, Shrub scene, in the specular scattering case: θ 1 = θ 2 = 30 , in the non-specular case: θ 1 = 30 , θ 2 = 45 ).
ParametersValue
Probability (%)Non-specular scatteringMEO96.9
IGSO100.0
GEO stationary100.0
Specular scatteringMEO96.9
IGSO100.0
GEO stationary100.0
Mainlobe to mainlobe specular scatteringMEO9.0
IGSO22.5
GEO stationary48.4
Required average transmitted power (kW)Non-specular scattering0.7
Specular scattering-
Mainlobe to mainlobe specular scattering2.2 × 108 ( A s ≈ 30 km × 30 km (maximum)) 0.8 ( A s = 100 m × 100 m)

Share and Cite

MDPI and ACS Style

Li, Y.; Monti Guarnieri, A.; Hu, C.; Rocca, F. Performance and Requirements of GEO SAR Systems in the Presence of Radio Frequency Interferences. Remote Sens. 2018, 10, 82. https://doi.org/10.3390/rs10010082

AMA Style

Li Y, Monti Guarnieri A, Hu C, Rocca F. Performance and Requirements of GEO SAR Systems in the Presence of Radio Frequency Interferences. Remote Sensing. 2018; 10(1):82. https://doi.org/10.3390/rs10010082

Chicago/Turabian Style

Li, Yuanhao, Andrea Monti Guarnieri, Cheng Hu, and Fabio Rocca. 2018. "Performance and Requirements of GEO SAR Systems in the Presence of Radio Frequency Interferences" Remote Sensing 10, no. 1: 82. https://doi.org/10.3390/rs10010082

APA Style

Li, Y., Monti Guarnieri, A., Hu, C., & Rocca, F. (2018). Performance and Requirements of GEO SAR Systems in the Presence of Radio Frequency Interferences. Remote Sensing, 10(1), 82. https://doi.org/10.3390/rs10010082

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop