Next Article in Journal
Deep Fully Convolutional Networks for Cadastral Boundary Detection from UAV Images
Next Article in Special Issue
Impact of Signal Quantization on the Performance of RFI Mitigation Algorithms
Previous Article in Journal
Spatiotemporal Analysis of MODIS NDVI in the Semi-Arid Region of Kurdistan (Iran)
Previous Article in Special Issue
Interference Mitigation for Synthetic Aperture Radar Based on Deep Residual Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

LEO to GEO-SAR Interferences: Modelling and Performance Evaluation

by
Antonio Leanza
1,
Marco Manzoni
1,
Andrea Monti-Guarnieri
1,* and
Marco di Clemente
2
1
Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milan, Italy
2
Agenzia Spaziale Italiana, Via del Politecnico, 00133 Roma, Italy
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(14), 1720; https://doi.org/10.3390/rs11141720
Submission received: 6 June 2019 / Revised: 12 July 2019 / Accepted: 17 July 2019 / Published: 20 July 2019
(This article belongs to the Special Issue Radio Frequency Interference (RFI) in Microwave Remote Sensing)

Abstract

:
This paper proposes a statistical model to evaluate the impact of the signal backscattered by low Earth orbiting (LEO) synthetic aperture radar (SAR) and received by GEO-stationary orbiting SAR. The model properly accounts for the bistatic backscatter, the number of LEO-SAR satellites and their duty cycles. The presence of many sun-synchronous, dawn-dusk satellites creates a 24 h periodic pattern in interference that should be considered in the acquisition plan of future geostationary SAR. The model, implemented by a numerical simulator, allows also the prediction of performance in future scenarios of many LEO-SAR. Examples and evaluations are made here for X band.

Graphical Abstract

1. Introduction

Geosynchronous synthetic aperture radar (SAR) are acquiring momentum due to their capability to perform imaging over wide areas with a revisit otherwise impossible with low Earth orbiting (LEO) or medium Earth orbiting (MEO) SAR [1,2,3,4]. In particular, the concept of a geostationary orbiting SAR, that should fulfill ITU requirements allowing a slight, non-zero eccentricity orbit to provide a synthetic aperture, is quite attractive for their capability of near-continuous imaging. Images are obtained by integrating over long intervals, up to hours, to compensate for the huge spread losses, while enabling the use of current technologies in terms of power and antenna reflectors [4,5,6].
Nearly half of the Earth’s surface is seen from Geostationary (GEO) orbit under a small angular aperture of ~19°: this is the key to fast, super-continental access, but at the same time, it exposes these systems to many potential radio frequency interferences (RFI) coming from anywhere. In [7], the impact of disturbances coming from on-ground, or due to the signal back-scattered from LEO-SAR, has been analyzed under very general assumptions. It has been shown that C and X band systems are more promising than L band due to the lower population of RFI and higher antenna directivity.
The balance between ground and LEO RFI is different in the two frequencies: the C-band systems are more subject to on-ground RFI, since the band is open to ground services like weather radar [8,9], whereby X band is primarily allocated to ESS-Active (Exploration Satellite Service), that should, in principle, exclude on-ground RFI. However, the number of present and future X band SARs, like the forthcoming constellations of light SARs: ICEYE, Capella Space or UmbraLab, is so large to motivate the in-depth analysis performed in the paper, and the need for a performance model capable of adapting to present and future scenarios. While the numeric evaluations are here performed in X band, the methodology developed can be extended to any frequency.
The paper is organized as follows: Section 2 provides the model for estimating the mean power scattered by the terrain, illuminated by the LEO-SAR, and captured by the GEO-SAR. Section 3 provides a statistical characterization of that backscatter, made by analyzing the current LEO-SAR scenario. In Section 4 a performance evaluation is made as function of the hours of the day and the integration interval, i.e., accounting for the impact of many dawn-dusk satellites. Section 5 discusses possible mitigation approaches. The last two sections are about overall discussions and conclusions.

2. LEO-to-GEO RFI Evaluation

The geostationary radar geometry is sketched in Figure 1. The GEO-SAR, orbiting at about 42,000 km from Earth center, receives the contribution of its own backscatter, superposed to the backscatter originated by other N LEO satellites, L1…LN. As the GEO-SAR range is about 43 times that of LEO-SAR, its power density on the ground is expected to be ~432 = 1800 times lower than that of a LEO-SAR, by assuming the same transmitted power and antenna area. Therefore, it is unlikely that a GEO-SAR would disturb a LEO-SAR, but the opposite is a challenge for the GEO-SAR.
The evaluation of the signal to interference ratio (SIR) can carried out by comparing the backscattered echo power, as from Radar Equation [1,2,3,4,5,6,7], to the total RFI power. We assume that the GEO and LEO-SARs have uniform illumination and that transmit for all time the mean power, the product of peak power and the duty cycle. Notice that while the GEO-SAR is nearly stationary, the LEO-SAR moves by km/s, generating a fast time-varying disturbance.
The mean signal power received from GEO-SAR can be approximately expressed by assuming that all its mean power, P T g , is distributed to the ground beam surface (the wide green spot in Figure 1), then a fraction of this power, proportional to the mean monostatic backscatter coefficient σ G 0 , is scattered back to the sensor, and partly gathered by the GEO-SAR antenna [4,5,6]:
P G = P T G · σ G 0 · 1 4 π R G 2 · λ 2 4 π G G 0 · η T
where R G is the GEO-SAR range, η T the total losses, λ is the wavelength and G G 0 is the mean antenna directivity at the beam center (see Figure 2).
The mean RFI power received by the GEO-SAR and originated by the i-th LEO-SAR can be computed similarly to (1), but accounting for the fast time-variation of the LEO-SAR along its orbit:
P i ( t ) = v i ( t ) · r i · d i · P T i · σ 0 ( θ i L ( t ) , φ i L ( t ) , θ i G ( t ) , φ i G ( t ) ) 4 π R i G 2 ( t ) · λ 2 4 π G G ( θ i G ( t ) , φ i G ( t ) ) · η T
The terms have the following meanings:
  • vi(t) is a binary (1/0), “visibility” index, set to one if the LEO-SAR is in the part of Earth accessible by GEO-SAR, that is, the field of view assumed in Figure 2;
  • ri is the ratio between GEO-SAR and LEO-SAR bandwidth;
  • di is the LEO-SAR duty cycle along orbit;
  • P T i is the total LEO-SAR transmitted power;
  • σ 0 ( θ i L , φ i L , θ i G ( t ) , φ i G ( t ) ) is the bistatic scatter coefficient from the i-th LEO-SAR direction, defined by incidence and azimuth angles θ i L , φ i L to the GEO-SAR one, defined by θ i G ( t ) , φ i G ( t ) ;
  • R i G is the range from LEO-SAR scene center to GEO-SAR;
  • G G ( θ i G ( t ) , φ i G ( t ) ) is the GEO-SAR antenna directivity to the LEO-SAR beam center, see Figure 2.
Note that we do assume the same losses in both GEO, in (1) and LEO, in (2), that is equivalent to ignore them as for the signal-to-interference ratio (SIR).

Bistatic Scattering

The LEO-to-GEO bistatic scattering, σ 0 in (2) is to be evaluated for the proper geometry, where the illumination from LEO-SAR, most of them being polar orbiting, is nearly orthogonal to the backscatter in the direction of the geostationary SAR, for which range is almost south-north. The change in backscatter plane, shown in Figure 3. on the left is a major difference respective to [7].
The actual σ 0 has been calculated by means of the Gilbert–Johnson model [10] for category ‘soil’, and assuming X band, see also [11], and mean incidence θ i L = 30 ° . Examples of evaluation are shown in Figure 4 for different polarizations, scatter directions (back and forward, see Figure 3. on the right), and two LEO orbits: a 95° inclined, nearly polar one, and a 50° inclined nearly equatorial orbit.
We remark that the assumption of the ‘soil’ backscatter model is a consequence of the GEO-SAR longitude, here assumed located in central Europe. If the GEO-SAR were placed in mid ocean, evaluation of ocean clutter will be needed, which is strongly dependent on the wind-speed. The long term average wind speed in the central and south Atlantic is in the order of 5–10 m/s, for which X-band backscatter ranges from −15 to −10 dB, depending on the wind angle and illumination [12]. Therefore, the evaluation so far performed on the land, in Figure 4, provides conservative results.

3. Daily Averaged RFI Estimation

The present scenario of LEO-SARs missions suggests to evaluate the mean level of RFI, by averaging all the items in (2):
P ¯ ( t ) = v ¯ ( t ) · r ¯ · d ¯ · P ¯ T · σ 0 ´ ( t ) 4 π R G 2 ¯ ( t ) · λ 2 4 π G G ¯ ( t ) ,
but retaining the dependence with time, that is the hour of the day, following the daily periodicity of both GEO and sun synchronous (SS) LEO orbits.
In order to perform the time average, we have explored 30 X-band LEO-SAR presently operating, listed in Table 1. Parameters have been taken from the EOPortal directory and WMO Oscar web servers, while information not available were fitted with averages, particularly for the military satellites. For all satellites, we assumed incidence angle of 30°, right looking, and transmissions in both H and V polarization.
In order to evaluate those terms in (3) that depend on the LEO-SAR position, i.e., the backscatter, the range to GEO-SAR and the GEO-SAR antenna gain, we averaged over 10 days by taking orbits from NORAD Space-Objects Tracking Service [13]. We have limited all the analysis to a square box of 140°, in both longitudes and latitudes, centered on the GEO longitude, outside which the bistatic backscatter is negligible.

3.1. Mean Number of Satellites

The histogram of visibility with time, and the consequent mean number of satellites, plotted in Figure 6, confirms the daily periodic trend that is at the base of our choice to evaluate (3), with an average number that is roughly 1/3 of the satellites.

3.2. LEO Sensor Daily RFI Power

The daily distribution of RFI power is the information complementary to the mean number of satellites, used in (3) To express the mean power. In Table 1 we have classified three types of orbits: Sun synchronous (SS), dawn dusk, DD, and all the others. The DD is a special case of the SS orbit, whereby the LEO-SAR rides the terminator between day and night. In such conditions, solar panels can always see the Sun, without being shadowed by the Earth. This simplifies the Radar payload thermal design, and minimizes the power needed from batteries. Therefore, DD is the orbit chosen by many LEO-SAR, particularly those using high powers and high orbit duty cycle, that are likely to be those most disturbing to GEO-SAR (Figure 5).
In order to appreciate the contributions of the three different orbits, we have represented in Figure 6 the daily RFI for the DD, the SS, and the generic non-SS missions. Noticeably, all DD LEO-SARs interferences are concentrated in two timeslots during the 24 h, of which hours depends on the GEO-SAR longitudes. In addition, for SS satellites the interferences are concentrated on precise time slots, but different for each satellite. Finally, the non-SS missions, that have near-equatorial orbit, contribute with an approximately uniform disturbance, yet of minor extent.
In the RFI plotted in Figure 6 left and right, a sharp peak was removed by nulling about 10 s of acquisitions, where LEO-SAR illumination intersected the main lobe of GEO-SAR antenna beam. Such an interval is a minimal loss of data compared to the long GEO-SAR integration time, typically hours.

3.3. Total RFI Power

The resulting hourly scenario of total RFI power from the 30 LEO satellites listed in Table 1, as seen by a GEO-SAR placed at 10° E longitude, is plotted in Figure 7. This result can be interpreted as the product between the mean number of sensors, as shown in Figure 5, and the weighted average of RFI like those in Figure 6. Quite noticeable is the time-variation of the disturbances, with a dynamic of 18 dB, that confirm the usefulness of our daily estimate.

4. Performance Evaluation and Future Trend

The interference of LEO-SAR would depend on the GEO-SAR parameters in terms of antenna, power, and scene center, but also on the hour of the day in which GEO-SAR acquires, and the image time. Here we evaluate performance by assuming the GEO-SAR parameters as from Table 2.
The noise equivalent sigma zero (NESZ) can be computed by equating the GEO and LEO power (1), (2) and integrating in the image time, Δt:
σ Z 0 = 1 Δ t t 0 t 0 + Δ t v ¯ ( t ) · r ¯ · d ¯ · P ¯ T P T G · σ 0 ¯ ( t ) σ G 0 · R G 2 R G 2 ¯ ( t ) · G ¯ G ( t ) G G 0 d t
and adding the thermal noise contribution [2,4,5,6].
The NESZ is reported in Figure 8 on the left, as function of the synthetic aperture and the time of the day. We notice the two red bands where NESZ achieves very bad values, up to 5 dB, preventing imaging for up to three hours in a day. However, for the rest of the time, acceptable values in the range from –14 to –22 dB are achieved. Furthermore, the expected variation of NESZ due to variation in the number of satellites has been plotted in Figure 8 On the right. This allows a prediction of GEO-SAR performance in the future, under the assumption that the time-averaged values in (3) will not change too much.

5. Mitigation

The mitigation strategies for RFI can be based either on coordination, or diversity in space, frequency, Doppler, code or time domains, or a mixture of all these [8].
The space domain mitigation is suggested after Figure 4: one should choose the GEO-SAR longitude so that the regions on the Earth, in yellow, where GEO-SAR is mostly sensitive to LEO-SAR backscatter are not frequently imaged by LEO-SARs, which may be those in open ocean.
The frequency domain mitigation is already accounted for by the term ri in (3) A fine-tuning needs knowledge of LEO-SAR frequency usage, that could be done by on-ground measures.
As for Doppler domain, it is easy to show that, for distributed targets, no separation is possible since LEO-SAR backscatter observed by GEO-SAR is white, then sample-to-sample independent. Let us evaluate the LEO-to-GEO bistatic Doppler of target P:
f L E O G E O ( t ) = 1 λ R b ( P ) t = 1 λ t ( | S L E O ( t ) P | + | S G E O ( t ) P | ) 1 λ t ( | S L E O ( t ) P | ) = f L E O ( t ) 2
where the bistatic range, Rb, has been computed by assuming that the GEO-SAR position, SGEO, does not change appreciably with time when compared with the LEO-SAR one, SLEO(t). Then, the total Doppler bandwidth is half of the LEO-SAR monostatic:
B L E O G E O = B L E O 2 = v L L E O
This results in B L E O G E O = 700 7000 Hz for a LEO-SAR with antenna length in the range of 1–10 m, whereas the pulse repetition frequency (PRF) of a typical GEO-SAR is at most 350 Hz for 400 km ground coverage [6]. Then, GEO-SAR always under-samples the bistatic backscatter from LEO-SAR, that would result in purely white process, making it indistinguishable from thermal noise.
While any separation in Doppler or even code domain would fail for distributed targets, the case of backscatter from isolated point targets is rather different. It is quite reasonable to assume that the pulse transmitted by LEO-SAR would be totally defocused when demodulated with a maybe different central frequency phase (also considering the different clocks) and correlated with a different pulse. Therefore, we have not to worry too much about the peak RFI, but only its mean value.
The time domain mitigation is the one suggested by our approach, and noticeable in Figure 8. The figure shows that there are two intervals over 24 h when NESZ is very bad. However, the GEO-SAR orbit is such that there are also two intervals when no imaging can be performed. Therefore, the best solution would be to match those intervals. To better understand how it can be done, let us evaluate the resolution achievable by a sub-aperture ΔTs, at different hours of the day. The GEO orbit in the Earth fixed frame is approximately an ellipse in the equatorial plane that can be expressed in polar format, as in [5]:
r G ( t ) = R g
ψ G ( t ) = ψ 0 + 2 e · s i n ( Ω E ( t t p ) )
where e is the eccentricity, ΩE = 2π/86,400 the Earth angular speed, tp the time the satellite crosses the perigee and ψ0 the GEO-SAR longitude. The achievable synthetic aperture in an interval ΔTs can be approximated, for short intervals, as, where the instantaneous velocity is:
v = t R g · ψ ( t ) = 2 e · R g · Ω E · s i n ( Ω E ( t t p ) )
The GEO-SAR resolution is then computed for each time of the day, t, as:
ρ a z = λ 2 L s = λ 2 v Δ T s = λ 4 e · R g · Ω E · s i n ( Ω E ( t t p ) ) · Δ T s
and plotted in Figure 9 or different times since the perigee pass, tp, and different sub-apertures. We assumed an eccentricity e < 0.002, the maximum allowed by ITU for geostationary satellites [14], then Figure 9 plots the best resolution available in X band. Notice that the resolution diverges for t = t p + N · 12 h (N integer), which corresponds to the stationary points in the orbit where v = 0 in (8).
The natural mitigation of RFI would then be to synchronize the perigee time, when no imaging is possible, see Figure 9, with the passage of DD LEO-SAR below the GEO-SAR, as suggested by Figure 8 and in the graphical abstract of this paper. This would improve the NESZ from 5 to less than 15 dB.
As a final remark, we notice that as bistatic backscatter in the end behaves like white noise, an additional mitigation, on the top of those so far discussed, is achieved just by increasing the integration time. This would be equally effective for any noise source. However, is does not appear by looking at NESZ profiles in Figure 8 (left). In fact, for a GEO-SAR the noise improvement provided by longer integration is perfectly balanced by the reduced signal that results from the better resolution [6]. Eventually one could trade the resolution to increase the number of looks, and then both the interferometric and radiometric quality.

6. Discussions

The 43× longer distance from Earth of GEO-SAR in respect to a LEO impacts in two destructive ways: (1) the power density on the ground is 422 = 1600 times lower at same antenna size, due to the bigger swath, and (2) the GEO-SAR antenna has poor directivity since the whole Earth is seen under a very small angle, just 19°. This means that RFI from LEO-SAR may be a potential killer for the mission. In fact, we see that in the present condition the RFI that is mostly driving the NESZ, and not the thermal noise.
The benefit and innovation of this paper is to underline two elements: (1) a proper backscatter model, that accounts for the almost orthogonal planes between illumination (from LEO-SAR) and backscatter (to GEO-SAR), and (2) the periodic illumination provided by dawn-dusk SAR, currently the majority. Both results in non-stationarity of the disturbances, respectively in space and time, and can be exploited to mitigate the impact of interferences.
It is suggested that the orbit longitude and time to perigee should be critically designed to minimize the RFI impact.
This analysis has been conducted with information about the present LEO-SAR scenario taken from public available sources and by filling many incomplete data. The real scenario could be quite different. However, this paper is intended as a framework where more precise information, maybe by on-ground measures, could be provided to get better predictions.
The theory exposed here is useful when most of the disturbances come from DD satellites, that are time synchronized and polar orbiting, and therefore causing scattering on the orthogonal plane. Conversely, it is totally useless for near equatorial-orbiting satellites. In the present scenario, there are seven non SS LEO-SAR (RISAT-2, TecSAR, Topaz1), over a total of 30, moreover, their mean backscatter was 5 dB lower that the other SS satellites. This seems to support the present framework.

7. Conclusions

The paper proposed a statistical model to evaluate the impact of the signal scattered from LEO-SAR illumination in the direction of a geostationary SAR. The analysis is carried out on an hourly basis, taking into consideration both the daily cycle of both GEO-SAR and sun synchronous satellites. Performance has been estimated by assuming an X band system, and the present scenario of 30 LEO-SARs. Statistics were evaluated by a computer simulation that averaged the interference from many orbits. The analyses show a dynamic range of 17 dB in noise equivalent sigma zero, depending on the hour of the day and the image time. This result, mostly due to the prevalence of dawn-dusk satellites, suggests an efficient mitigation that is implemented by synchronizing the GEO-SAR orbit with the DD passage, whereby other mitigation strategies have been analyzed.
The model is suited to be extended to future scenarios with different numbers of satellites, thanks to the statistics computed based on the present scenario.

Author Contributions

A.M.G. for conceptualization and writing, A.L. for investigation, implementation of the model into the software code and performance analysis, M.M. for support to formal development and revision, M.d.C. for project administration.

Funding

The paper is the result of an activity performed in cooperation between Agenzia Spaziale Italiana, ASI and Politecnico di Milano, in the frame of GEOSAR project (CUP F43C1700010005).

Acknowledgments

The authors would like to thank Davide Giudici and Daniele Mapelli at Aresys for providing the support for the evaluation X-band bistatic backscatter.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. 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]
  2. Long, T.; Hu, C.; Ding, Z.; Dong, X.; Tian, W.; Zeng, T. Geosynchronous SAR: System and Signal Processing; Springer: Berlin/Heidelberg, Germany, 2018. [Google Scholar]
  3. Madsen, S.N.; Chen, C.; Edelstein, W. Radar options for global earthquake monitoring. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Toronto, ON, Canada, 24–28 June 2002; Volume 3, pp. 1483–1485. [Google Scholar]
  4. 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] [Green Version]
  5. Prati, C.; Rocca, F.; Giancola, D.; Guarnieri, A.M. Passive Geosynchronous SAR System Reusing Backscattered Digital Audio Broadcasting Signals. IEEE Trans. Geosci. Remote Sens. Piscataway 1998, 36, 1973–1976. [Google Scholar] [CrossRef]
  6. 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]
  7. Li, Y.; Monti Guarnieri, A.; Hu, C.; Rocca, F. Performance and Requirements of GEO-SAR Systems in the Presence of Radio Frequency Interferences. J. Remote Sens. 2018, 10, 82. [Google Scholar] [CrossRef]
  8. National Academy of Sciences. A Strategy for Active Remote Sensing Amid Increased Demand for Radio Spectrum; National Academies Press: Washington, DC, USA, 2015. [Google Scholar]
  9. 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]
  10. Gilbert, M.S.; Johnson, J.T. A study of the higher-order small-slope approximation for scattering from a Gaussian rough surface. Waves Random Media 2003, 13, 137–149. [Google Scholar] [CrossRef]
  11. Comite, D.; Pierdicca, N.; Ticconi, F.; Guerriero, L. On the Modeling of the Bistatic Coherent Scattering from a Rough Surface. In Proceedings of the 2018 IEEE International Symposium on Antennas and Propagation & USNC/URSI National Radio Science Meeting, Boston, MA, USA, 8–13 July 2018; pp. 2473–2474. [Google Scholar]
  12. Romeiser, R.; Breit, H.; Eineder, M.; Runge, H.; Flament, P.; De Jong, K.; Vogelzang, J. Current measurements by SAR along-track interferometry from a Space Shuttle. IEEE Trans. Geosci. Remote Sens. 2005, 43, 2315–2324. [Google Scholar] [CrossRef]
  13. Kelso, T.S. CelesTack. Available online: www.celestrack.com (accessed on 31 January 2019).
  14. ITU. Radio Regulations Articles. 2016. Available online: http://search.itu.int/history/HistoryDigitalCollectionDocLibrary/1.43.48.en.101.pdf (accessed on 30 September 2018).
Figure 1. Geometry of a geostationary synthetic aperture radar (SAR), “Geo”, affected by interferences due to the signal backscattered from a set of low Earth orbiting (LEO) SARs, L1 and Li.
Figure 1. Geometry of a geostationary synthetic aperture radar (SAR), “Geo”, affected by interferences due to the signal backscattered from a set of low Earth orbiting (LEO) SARs, L1 and Li.
Remotesensing 11 01720 g001
Figure 2. Directivity (in dB) for a 5 m diameter, X band antenna, paced on a GEO-SAR satellite and pointing at Central Italy.
Figure 2. Directivity (in dB) for a 5 m diameter, X band antenna, paced on a GEO-SAR satellite and pointing at Central Italy.
Remotesensing 11 01720 g002
Figure 3. Models for the bistatic backscatter from the i-th LEO-SAR to the GEO-SAR. (a) Geometry, as seen from the target position: Notice the change of backscatter plane between LEO and GEO directions. (b) Cases of forward and back-scattering. S i L and S i G represent the incoming signal from LEO-SAR and outgoing to the GEO-SAR.
Figure 3. Models for the bistatic backscatter from the i-th LEO-SAR to the GEO-SAR. (a) Geometry, as seen from the target position: Notice the change of backscatter plane between LEO and GEO directions. (b) Cases of forward and back-scattering. S i L and S i G represent the incoming signal from LEO-SAR and outgoing to the GEO-SAR.
Remotesensing 11 01720 g003
Figure 4. Bistatic backscatter evaluated as function of the geographic position of the LEO-SAR and observed by a GEO-SAR located at 10° E longitude (red cross) according to [10,11], for X-band, ‘soil’. Upper four panels: Polar orbiting LEO-SAR. Lower panels: nearly equatorial orbiting LEO-SAR.
Figure 4. Bistatic backscatter evaluated as function of the geographic position of the LEO-SAR and observed by a GEO-SAR located at 10° E longitude (red cross) according to [10,11], for X-band, ‘soil’. Upper four panels: Polar orbiting LEO-SAR. Lower panels: nearly equatorial orbiting LEO-SAR.
Remotesensing 11 01720 g004
Figure 5. Estimation of the probability density of LEO-SAR visibility in GEO-SAR: Earth view (left) and mean number of satellites, (right), as function of time, derived from analysis of Table 1.
Figure 5. Estimation of the probability density of LEO-SAR visibility in GEO-SAR: Earth view (left) and mean number of satellites, (right), as function of time, derived from analysis of Table 1.
Remotesensing 11 01720 g005
Figure 6. Estimation of the radio frequency interferences (RFI) power, as a function of the hour of the day, for three different cases: a SS/DD satellite with near polar orbit (left), two satellites with low inclined orbit (center), and a set of five satellites with SS orbits in three different planes (right).
Figure 6. Estimation of the radio frequency interferences (RFI) power, as a function of the hour of the day, for three different cases: a SS/DD satellite with near polar orbit (left), two satellites with low inclined orbit (center), and a set of five satellites with SS orbits in three different planes (right).
Remotesensing 11 01720 g006
Figure 7. Hourly estimation of the mean LEO-to-GEO RFI power, due to the 30 LEO-SAR in Table 1.
Figure 7. Hourly estimation of the mean LEO-to-GEO RFI power, due to the 30 LEO-SAR in Table 1.
Remotesensing 11 01720 g007
Figure 8. (Left) Mean noise equivalent sigma zero (NESZ), computed by assuming as noise the bistatic backscatter from the present scenario of 30 LEO-SARs in Table 1, as a function of the hour of the day and the integration time. (Right) variation of the NESZ predicted as function of the number of X-band LEO-SAR.
Figure 8. (Left) Mean noise equivalent sigma zero (NESZ), computed by assuming as noise the bistatic backscatter from the present scenario of 30 LEO-SARs in Table 1, as a function of the hour of the day and the integration time. (Right) variation of the NESZ predicted as function of the number of X-band LEO-SAR.
Remotesensing 11 01720 g008
Figure 9. Azimuth resolution achievable by an X band GEO-SAR with the maximum allowed eccentricity of 0.002 [14], as function of the hour of the day since the perigee pass (t = 0), computed for different subapertures. The shaded areas mark times when bad performance is achieved: It is recommended that they match the periodical peaks in Figure 7.
Figure 9. Azimuth resolution achievable by an X band GEO-SAR with the maximum allowed eccentricity of 0.002 [14], as function of the hour of the day since the perigee pass (t = 0), computed for different subapertures. The shaded areas mark times when bad performance is achieved: It is recommended that they match the periodical peaks in Figure 7.
Remotesensing 11 01720 g009
Table 1. Parameters of the 30 LEO-SAR satellites assumed in the study.
Table 1. Parameters of the 30 LEO-SAR satellites assumed in the study.
System# SatsAverage Power [W]Bandwidth [MHz]Orbit 1Orbit Duty CycleOrbit Inclination 2
Cosmo-SkyMED484090SS/DD15%97.9°
IGS-Radar2 + 1 + 2700 22002SS10% 297.5° + 97.3° + 97.4°
Kompsat-51400100SS/DD15%97.6°
Meteora-M N21500100SS15%98.8°
Risat-21400100-15%41°
SAR-Lupe5250150SS10% 298.2°
TSX/TDX2400150SS/DD15%97.44°
TecSAR1750200-15%41°
Topaz4 + 1700 22002-10% 2123° + 106°
SeoSAR-PAZ1400100SS/DD10% 297.44°
Yaogan4700 22002SS15%97.8° + 97.5° + 97.3° + 97.9°
1 SS = sun synchronous, DD = dawn dusk. 2 True values for many military SARs are unknown, we assumed those of SAR-Lupe.
Table 2. Parameters for the X-band GEO-SAR.
Table 2. Parameters for the X-band GEO-SAR.
SymbolDefinitionValue
λwavelength3 cm
AG0antenna equiv. area 20 m2
RGrange37,000 km
BG bandwidth7.5 MHz
θGincidence angle50°
PGmean tx power400 W
beam center (lat, long)42°N, 10°E
ψ0mean longitude10°E

Share and Cite

MDPI and ACS Style

Leanza, A.; Manzoni, M.; Monti-Guarnieri, A.; di Clemente, M. LEO to GEO-SAR Interferences: Modelling and Performance Evaluation. Remote Sens. 2019, 11, 1720. https://doi.org/10.3390/rs11141720

AMA Style

Leanza A, Manzoni M, Monti-Guarnieri A, di Clemente M. LEO to GEO-SAR Interferences: Modelling and Performance Evaluation. Remote Sensing. 2019; 11(14):1720. https://doi.org/10.3390/rs11141720

Chicago/Turabian Style

Leanza, Antonio, Marco Manzoni, Andrea Monti-Guarnieri, and Marco di Clemente. 2019. "LEO to GEO-SAR Interferences: Modelling and Performance Evaluation" Remote Sensing 11, no. 14: 1720. https://doi.org/10.3390/rs11141720

APA Style

Leanza, A., Manzoni, M., Monti-Guarnieri, A., & di Clemente, M. (2019). LEO to GEO-SAR Interferences: Modelling and Performance Evaluation. Remote Sensing, 11(14), 1720. https://doi.org/10.3390/rs11141720

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