Next Article in Journal
Semiconductor Quantum Dots in Chemical Sensors and Biosensors
Next Article in Special Issue
A Real-Time Smart Sensor for High-Resolution Frequency Estimation in Power Systems
Previous Article in Journal
Detection and Classification of Human Body Odor Using an Electronic Nose
Previous Article in Special Issue
A MEMS-Based Flow Rate and Flow Direction Sensing Platform with Integrated Temperature Compensation Scheme
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bistatic Radar Configuration for Soil Moisture Retrieval: Analysis of the Spatial Coverage

by
Nazzareno Pierdicca
1,
Ludovico De Titta
1,
Luca Pulvirenti
1,* and
Giuliano Della Pietra
2
1
Department of Electronic Engineering, Sapienza University of Rome, via Eudossiana 18, 00184 Rome, Italy
2
Space Engineering S.p.A., 91, Via dei Berio, Rome 00155, Italy
*
Author to whom correspondence should be addressed.
Sensors 2009, 9(9), 7250-7265; https://doi.org/10.3390/s90907250
Submission received: 21 August 2009 / Revised: 8 September 2009 / Accepted: 9 September 2009 / Published: 10 September 2009
(This article belongs to the Special Issue Sensor Configuration and Smart Sensors)

Abstract

:
Some outcomes of a feasibility analysis of a spaceborne bistatic radar mission for soil moisture retrieval are presented in this paper. The study starts from the orbital design of the configuration suitable for soil moisture estimation identified in a previous study. This configuration is refined according to the results of an analysis of the spatial resolution. The paper focuses on the assessment of the spatial coverage i.e., on the verification that an adequate overlap between the footprints of the antennas is ensured and on the duty cycle, that is the fraction of orbital period during which the bistatic data are acquired. A non-cooperating system is considered, in which the transmitter is the C-band Advanced Synthetic Aperture Radar aboard Envisat. The best performances in terms of duty cycle are achieved if the transmitter operates in Wide Swath Mode. The higher resolution Image Swath Modes that comply with the selected configuration have a duty cycle that is never less than 12% and can exceed 21%. When Envisat operates in Wide Swath Mode, the bistatic system covers a wide latitude range across the equator, while in some of the Image Swath Modes, the bistatic measurements, collected from the same orbit, cover mid-latitude areas. In the latter case, it might be possible to achieve full coverage in an Envisat orbit repeat cycle, while, for a very large latitude range such as that covered in Wide Swath Mode, bistatic acquisitions could be obtained over about 65% of the area.

1. Introduction

A spaceborne bistatic radar system is defined when antennas for reception and transmission are physically separated and located aboard two spacecraft. The system could be either cooperating (transmitter and receiver designed for the specific bistatic application), or non-cooperating (receiver designed independently of the transmitter) [1]. Non-cooperating systems can use already-existing satellite instruments, such as Synthetic Aperture Radars (SARs), as transmitters (or illuminators) [2]. Bistatic applications were investigated in the fields of target detection [3,4] and planetology [5]. Concerning Earth Observation, the specular reflection from navigation satellites was exploited for scatterometric applications over the oceans [6]. The possibility to envisage bistatic configurations for SAR interferometry was studied in [7,8]. Recently, experiments involving the X-band have been carried out (e.g., [911]) to assess the additional information contained in the bistatic reflectivity of targets.
As for land applications, Ceraldi et al. [12] demonstrated the potential to detect the signal scattered in the specular direction for retrieving bare soil parameters. However, they outlined that measuring the coherent component of the scattered signal by means of bistatic radars is problematic, because the spatial resolution becomes very poor around the specular direction, and ambiguities between surface points before and after the specular point occur [13]. Zavorotny and Voronovich stated, in a brief report [14], that considering a GPS signal impinging on a rough land surface, soil moisture is related to the ratio between horizontally and vertically polarized scattered waves.
In a previous work [2], we reported on a theoretical investigation aiming at identifying the best bistatic measurement configuration, in terms of incidence angle, observation direction, polarization, and frequency band, for soil moisture content (SMC) retrieval. In addition, we evaluated the improvement of the estimation accuracy with respect to the conventional backscattering measurements. While in [2] we did not tackle the problem of verifying the identified radar configurations from a technical point of view, this work deals with spatial resolution, spatial coverage and duty cycle (i.e., the fraction of the orbit in which suitable bistatic data collection is ensured). It aims at evaluating whether an adequate overlap between the footprints of the transmitting and receiving antennas is ensured for the bistatic observation geometry that guarantees both an improvement of the quality of the SMC retrieval and a quite good spatial resolution.
We point out that investigations on the design of a spaceborne mission implementing a configuration of bistatic radars devoted to specific environmental applications are generally lacking in the literature. Moreover, the literature generally deals with fixed bistatic configurations (i.e., the static design), without considering the orbital dynamics, while the feasibility of maintaining the selected configuration of sensors along the orbit of the spacecraft on which the receiver is installed, ensuring good performances in terms of operational parameters such as spatial coverage, is a crucial aspect.
In the present study, the nominal bistatic mission is considered, so that problems such as the orbit and attitude controls are not tackled, thus implicitly assuming that the passive system is periodically controlled in order to maintain the designed bistatic formation.
In Section 2, the findings regarding the bistatic configurations most suitable for SMC retrieval are summarized and the selection of one configuration, based on the evaluation of the spatial resolution, is described. In Section 3, the orbit design approach is depicted, while, in Section 4, the results of the analysis of the spatial coverage are presented and discussed. Section 5 draws the main conclusions.

2. Selection of a Bistatic Configuration

2.1. Configurations Suitable for Soil Moisture Retrieval

Bistatic configurations are defined in terms of frequency, polarization, and transmitter-target-receiver relative geometry. This geometry is shown in the left panel of Figure 1, in which the incident plane is assumed as a reference coordinate plane, and the configuration is identified through the incidence angle (θi) with respect to the vertical, the zenith (θs) and azimuth (φs) scattering angles, and their acceptable ranges to be compatible with the considered application.
In [2], we adopted a well-established electromagnetic surface scattering model, such as the Advanced Integral Equation Model (AIEM) [15], to numerically simulate the bistatic measurements, for the purpose of identifying the configurations suitable for soil moisture estimation. We underline that the AIEM model has a large range of validity, so that the study described in [2] was carried out under different conditions. L-, C- and X-bands were considered, both Gaussian and Exponential autocorrelation functions were adopted, and the parameters characterizing the surface roughness, i.e., the standard deviation of heights s and the correlation length l, varied within fairly standard intervals, encompassing most of the typical agricultural fields in bare soil conditions. The results we obtained can be therefore considered fairly general.
The findings in [2] indicate that, especially at C-band, the quality of soil moisture retrievals improves by complementing bistatic measurements with monostatic ones provided by already-existing SARs. This could be expected because, for bistatic systems, the detected field is determined by a different behavior of the scattering with respect to the monostatic case, and in the specular configuration also by a coherent scattering mechanism as opposed to a diffuse one, thus providing independent information on the target parameters. The need to take advantage of various sources of information for tackling the complex and usually ill-posed problem of SMC retrieval was indeed pointed out in numerous literature studies (e.g., [16], where polarimetric data were used).
The range of valuable scattering directions, for a C-band radar, singled out in [2] is sketched in the right panel of Figure 1. Observations around nadir, at θs below 4°–10°, form a cone (i.e., φs can assume any value) whose vertical axis corresponds to the z-axis (in a Cartesian reference system), while the vertex is in the origin. Furthermore, for φs less than 40°–50°, the zenith scattering angle can indefinitely be increased by including the specular direction and some grazing observations. The illuminator should observe the Earth with a quite small incidence angle, approximately between 15° and 35°. In [2] it was demonstrated that, at C-band, the standard deviation of the SMC retrieval error can be reduced up to a factor 3 with respect to that achievable with monostatic observations, by integrating backscattering and bistatic measurements.

2.2. Spatial Resolution Analysis

The ground range and azimuth resolutions could be very poor in some bistatic configurations, thus implying a bad radar image quality. For instance, the ground range resolution is critical in specular configuration [1,12]. Basing on these considerations, we have performed an evaluation of the bistatic spatial resolutions in order to put an additional constraint to the sensor configurations selected in [2]. It is worth noting that no adequate model has been found in the literature for such an exercise, because only few cases addressed the spaceborne receiver configuration. In [1], this matter was analyzed, but only for a two-dimensional (2-D) configuration restricted to the bistatic plane. General considerations and formulas valid for the 3-D case were provided in [13].
Starting from the relationships found in [13], we have derived the formulas for evaluating the ground range (ρgr) and azimuth (ρa) bistatic resolutions, by doing the following assumptions: i) flat and fixed Earth; ii) satellites orbiting at the same height (and therefore flying with the same velocity); iii) equal coherent integration times. Then, we have normalized the bistatic resolutions to the monostatic ones, obtaining equations depending on the relative geometry of the satellites carrying the sensors and independent from orbital and antenna parameters. In this section, we do not report all the mathematical manipulations, for the sake of conciseness, but only the result. More details can be found in the Appendix. The final equations are the following:
ρ gr ρ gr back = 2   sin   θ i sin 2   θ i + sin 2   θ s 2   sin   θ i   sin   θ s cos   φ s
ρ a ρ a back = 2   cos   θ i F ( θ i , θ s , φ s , δ T , δ R )
where:
F ( θ i , θ s , φ s , δ T , δ R ) = cos 2   θ i ( cos 4   θ i   cos 2   δ T + sin 2   δ T ) + cos 2   θ s [ 1 sin 2   θ s ( 1 + cos 2   θ s ) cos 2 ( φ s δ R ) ] + + 2   cos   θ i   cos   θ s [ cos 2   θ i   cos   δ T   cos   δ R + sin   δ T   sin   δ R sin 2   θ s ( cos 2   θ i   cos   δ T   cos   φ s + sin   δ T   sin   φ s )   cos   ( φ s δ R ) ]
In Equations (1) and (2a), the back superscript indicates the spatial resolutions computed for the conventional monostatic case (i.e., the backscattering measurement) and δT and δR are the angles complementary to those identified by the unit vector normal to the incidence plane and the velocity vectors (see Appendix).
Starting from the above formulas, and supposing that the satellites’ velocity vectors are directed normally to the incidence plane (i.e., δT = δR =90°), maps of ρgr/ρgrback and ρa/ρaback have been generated as function of θs (in the range [0°–60°]) and φs (in the range [0°–180°]), for a fixed value of θi. Figure 2 shows the maps for the minimum and maximum values of θi recommended for estimating soil moisture (i.e., 15° and 35°).
By observing Figure 2, it can be noted that the ground range resolution is more critical than the azimuth one. The latter does not exceed 2ρaback, while ρgr can be several times larger than ρgrback, especially near the specular direction (as expected). As an additional constraint to the bistatic sensors configurations selected in [2], we have chosen ρgr < ρgrback, thus finding that φs must be limited in the interval [90°−270°], that is in the backward quadrant (see Figure 1 right panel and the Appendix). We have therefore considered the following ranges for the zenith and azimuth scattering angles: θs ∈ [0°−8°] and φs ∈ [90°−270°].

3. Orbit Design

According to the chosen frequency (i.e., C-band), we have made reference to the ASAR instrument onboard Envisat. The real operative Image Swath Modes (ISMs) and Wide Swath Mode (WSM) have been taken into account for the analysis. ISMs guarantee a fairly high resolution, but have a smaller near and far range capability that implies worse monostatic (and thus bistatic) spatial coverage. WSM presents opposite characteristics. Table 1 reports some parameters of the ASAR illuminator in Image and Wide Swath modes. Note that, among the seven Envisat/ASAR Image Swath Modes, only four were selected, because the other ones observe the Earth with incidence angles larger than 35°, which are not compatible with our requirements for SMC retrieval.
We have firstly designed the orbit of the receiver by choosing its observation geometry at the equator based on the requirements for SMC estimation (i.e., the static design) and then the orbits of both platforms have been propagated to evaluate the antenna footprint superimposition, i.e., the spatial coverage of the bistatic acquisitions. Note that the time when the passive satellite is over the equator has been assumed as the initial one [17].
We have made a number of assumptions. Firstly, that the bistatic system is non-cooperative, so that orbit, attitude (including yaw-steering) and antenna pointing of the active system are not subordinated to the bistatic acquisition requirements. This is of course a worst case assumption, which aims to use existing radar as source of opportunity. Secondly, the transmitting antenna is right-looking. Moreover, the satellite carrying the receiver (hereafter denoted also as the passive satellite) flies in formation with that carrying the illuminator (the active one), on a parallel orbit, thus establishing a “parallel orbit pendulum” configuration, in which the platforms move along orbits with the same inclination, but different ascending nodes. In principle, also the “leader-follower” configuration, in which the spacecraft fly on the same orbit, but with different crossing times of ascending node, would be possible. However, it does not allow any of the bistatic observation geometries previously selected if a non-cooperative bistatic system and a side-looking illuminator are assumed.
Another hypothesis we have made for the static design of the orbits is that the receiver performs bistatic acquisitions only along the ascending pass, since the crossing of the orbits near the poles makes the transmitter to illuminate out of the footprint of the receiver. Note that an eventual left/right looking capability of the receiver would dramatically change the bistatic angles with respect to the required ones, whereas only a cooperative transmitter with left/right looking capability would enable valuable data acquisition in the descending pass.
The orbit of a spacecraft is described by means of the well-known six Keplerian parameters: semi-major-axis, eccentricity, inclination, right ascension of the ascending node (Ω), perigee argument, and mean anomaly (M). We have considered Envisat as the active satellite, so that its orbital parameters are fixed. The passive satellite semi-major-axis, eccentricity, and inclination have been chosen coincident with those of the active one, in order to reduce orbit maintenance operations. In addition, the same perigee argument has been considered, to minimize the instantaneous satellite velocity differences and in turn, the along-track relative displacements [17].
The remaining design parameters (Ω and M) define the difference between the orbits of the two satellites. The differences between the ascending node right ascensions (ΔΩ) and the mean anomalies (ΔM), i.e., the relative position of the satellites (shown in Figure 3), have been computed by means of the formulas provided in [17]. ΔM has been chosen in order to enable the passive satellite to be located in the incidence plane of the illuminator at the initial time, thus fixing the observation geometry when the receiver is over the equator. Such a condition corresponds to φs = 180°, in this case.
Figure 3 shows the position of the active and passive satellites when the latter is over the equator and the line joining the two positions (i.e., the baseline). This direction has to be comprised within the azimuth aperture of both the transmitting and receiving antennas. To calculate ΔΩ and ΔM, we have assumed θi equal to the maximum among those suggested in [2] (i.e., 35°), in order to increase the range of latitudes for bistatic acquisitions [17]. As for θs at the initial epoch, we have firstly considered that the receiver should be in the backward quadrant because of the constraint we have imposed for spatial resolution. Then we have made the hypothesis that, for a soil moisture application, the minimum bistatic swath should be 10 km.
By choosing θs = 0° at the initial epoch, the time interval during which θs is in a useful range (i.e., [0°–8°]) when the orbits are propagated is maximized (note that at the equator the baseline is at its maximum, as will be discussed in Section 4). However, looking at Figure 4, that shows the observation geometry at the initial time, it can be noted that such a choice would imply a very small bistatic swath at the equator (theoretically, only one point), because for θs = 0° the receiver is just above the boundary of the area illuminated by Envisat, with the maximum allowable zenith angle (i.e., 35°), and cannot accomplish a forward observation. For the same reason, the ISMs whose far range looking angle is less than 35° (see Table 1) do not perform bistatic acquisitions over the equator, as will be shown in Section 4. The receiving near range point (NRRX) has to be located within the area illuminated by the transmitting antenna (see Figure 4) and its distance to the transmitter far range point (FRTX) cannot be less than 10 km. Hence, the best choice for the initial θs turned out to be 1°, which is (approximately) the minimum zenith scattering angle for which the overlapping between the footprints of the antennas occurs with the selected minimum width of the bistatic swath. Table 2 reports the selected configuration and the orbital parameters at the initial epoch, that is the static design.
The Satellite Tool Kit (available at: http://www.agi.com) has been successively used to propagate the orbits, taking into account only the orbital perturbations due to the J2 geo-potential harmonic, to maintain a heliosynchronous orbit and assuming negligible the orbital decay (i.e., the effect of the differential aerodynamic drag), as done in [18].

4. Discussion of the Outcomes

To produce coverage data, an appropriate software tool, which accounts for spacecraft propagated orbits, yaw-steering maneuvers of both spacecraft, sensors pointing geometry, Earth rotation and estimates the targeted area, has been developed. At each time point of the simulation, whose overall period is 35 days, i.e., the orbit repeat cycle of Envisat, the software evaluates the area of the Earth surface observed simultaneously by the two antennas on the basis of their aperture angles and of the satellites’ positions. The portion of this area which is actually observed according to the observation geometry selected in section 2, i.e., θI < 35°, θs ∈ [0°−8°] and φs ∈ [90°−270°], is the area target. To perform the simulation, we have considered the transmitter able to illuminate according to the real access capabilities of the reference mission (i.e., Envisat, see Table 1). As for the receiving antenna, we have assumed a main lobe azimuth aperture of 6° (see [17]) which is capable to point toward the area illuminated by the transmitter by means of attitude maneuvers or electronic steering. The latter does not ensure the zero Doppler condition of bistatic measurements, but its purpose is the superposition of the antenna footprints.
Table 3 reports the results. It can be noted that WSM presents the best results in terms of duty cycle (∼40%) thanks to the greater access capability. Higher resolution modes achieve duty cycles not less than 12% (actually 11.9%) and that can exceed 21%, which is a still acceptable performance, considering that we are in the hypothesis of a fully non-cooperative system, i.e., the maximum duty cycle is 50% for a fixed looking side of the transmitting antenna. It is worth noting that we have decided to discard the swaths smaller than 10 km in the computation of the fraction of orbital period during which the bistatic data are acquired and this choice implies an overall decrease of the resulting duty cycle.
ISM1 and WSM reach the highest latitudes (almost 72° in both the geographic hemispheres). In addition, the far range capability of the considered transmitter operating mode limits the latitude range of bistatic acquisitions. As stated when discussing Figure 4, since the passive satellite orbit is defined according to the maximum incidence angle identified for the application (i.e., 35°) during static design, the operational modes ISM1, ISM2 and ISM3 do not perform bistatic acquisitions at low latitudes but they focus on mid-latitude belts.
Figure 5 shows the Earth maps with the acquired targeted areas (red dots) for WSM (upper panel), ISM4 (central panel) and ISM2 (lower panel). The ground tracks of the active and passive satellites are plotted in green and blue, respectively. Although our simulation considers the entire Envisat orbit repeat cycle (35 days), Figure 5 regards a time-window of one week, for the sake of clarity. The different latitude bands covered in WSM, ISM4 and ISM2 are clearly visible.
To quantify the spatial coverage of the bistatic measurements, we have evaluated the ratio between the sum of the areas of bistatic targets and the area of the Earth surface within the covered latitude range. Such a computation sums up all the areas imaged during a complete Envisat cycle of 35 days, including also zones already observed in previous orbits. This result is therefore an estimate of the fractional coverage that does not ensure that every point within a specific latitude range is actually covered by a bistatic acquisition, but that can indicate whether full coverage could be potentially obtained in such a range. The ratio is minimum for ISM4 and is in the order of 65% for WSM, while for ISM1 and ISM2 our evaluation indicates that full coverage might be achieved within 35 days.
Finally, Figure 6 shows the trends of the baseline absolute value |B| (distance between the satellites) and of its components (Bx,By,Bz) throughout one orbit, that we obtained as a consequence of our design. The baseline is evaluated on the orbital reference frame centered in the centre of mass of the active satellite, with x-axis towards its instantaneous velocity vector, y-axis normal to the orbital plane and z-axis in direction of the Earth center. The distance between the satellites achieves its maximum value at the equator and the minimum one near the poles. By represents the main component of B and changes its sign at the orbit interceptions near the poles, when the satellites change their relative orientation, while it achieves the maximum value at the equator. Also Bz goes to zero at the orbits interceptions and the minimum baseline at the orbit crossing, i.e., the safety distance between the satellites, has turned out to be equal to 35.5 km, a result similar to that found in [18] in which Envisat was used as the active satellite. Observing Figure 6, it can be noted that the maximum value of the baseline (at the equator), is slightly larger than 500 km.

5. Conclusions

Some outcomes of a coverage analysis of a spaceborne bistatic mission for SMC estimation have been presented. The study has started from the identification of the bistatic configurations suitable for SMC retrieval, in terms of frequency and ideal transmitter-target-receiver relative geometry accomplished in a previous study. An evaluation of the spatial resolution of the bistatic system has been firstly carried out. It has led us to restrict our analysis to observations around nadir and in the backward quadrant. Then, the study has assessed the feasibility of the nominal mission in terms of spatial coverage and duty cycle, assuming a non-cooperative system, with Envisat/ASAR as illuminator.
We have shown that the Wide Swath Mode has the best duty cycle (∼40%). Higher resolution image modes yield acceptable results in terms of this parameter (never less than 12%), considering that, in our hypothesis of fully non-cooperative system, the maximum achievable value is 50%. It has been found, that, for high resolution modes, the minimum and maximum latitudes do not necessarily identify one latitude belt, so that their setting may allow focusing on a particular geographic region, such as mid-latitude areas. In a selected latitude band, it might be possible to achieve full coverage within the Envisat orbit repeat cycle (35 days), while for a larger latitude range, almost covering the entire planet, it is unfeasible.
The improvement of soil moisture retrieval that we have demonstrated, in a previous work, to be achievable through a bistatic mission, may be actually obtained by designing a simple and cheap (even in terms of power requirements) receiver that flies in formation with a standard C-band SAR. Combining the active monostatic with the bistatic measurements can strengthen the contribution of single frequency radar observations in many disciplines needing an evaluation of soil moisture, such as hydrology, agriculture and climate change monitoring.
The methodology based on a study of the sensitivity to a geophysical target parameter of bistatic radars, followed by a system performance analysis can be useful for other applications (e.g., vegetation biomass retrieval) and for future missions. For instance, this approach can provide support to a possible employment of a passive receiver flying in formation with the forthcomong European Space Agency (ESA) Sentinel satellites. In particular, the wide swath capability of the Sentinel-1 C-band SAR is expected to improve the coverage performances of the bistatic mission with respect to that computed for Envisat.

6. Appendix Spatial Resolution

In this section, Equations (1) and (2) are derived starting from the relationships found in [13]. The hypotheses of flat Earth and same satellite heights hold true. We make reference to the geometry shown in Figure 7.
Figure 7. Geometry for evaluating the bistatic spatial resolution. T, R and P represent the positions of transmitter, receiver and target, respectively.
Figure 7. Geometry for evaluating the bistatic spatial resolution. T, R and P represent the positions of transmitter, receiver and target, respectively.
Sensors 09 07250f7
By assuming the incidence plane as a reference, let us consider a three-dimensional Cartesian coordinates system in which the (x,z)-plane is coincident with the incidence one and the (x,y)-plane is considered to be the (flat) Earth surface. Let us fix: i) the transmitter position T (xT, yT, H) in such a reference frame; ii) the incidence angle θi; iii) the elevation θs and azimuth φs scattering angles. Then, we can determine, as function of these parameters, the following ones: iv) the ground target position P (xT + θi, 0, 0); v) the receiver position R (xT + Htanθi + Htanθscosφs, yT + Htgθssinφs, H).
Denoting by x0, y0 and z0 the versors of the considered Cartesian system, the vectors directed from the transmitting (uT) and receiving (uR) antennas to the target are (see Figure 7):
u T = H   tan   θ i x 0 H z 0
u R = H   tan   θ s   cos   φ s x 0 H   tan   θ s   sin   φ s y 0 H z 0
The transmitter-to-target and receiver-to-target ranges are the modules of uT and uR, respectively:
| u T | = H / cos   θ i
| u R | = H / cos   θ s
Now we can consider the relationship for the ground range resolution proposed in [13]:
ρ gr = c 2 B W | Γ xy   u eff |
where c denotes the velocity of light, Bw is the pulse bandwidth and Γxy is the projector onto the horizontal plane. 2ueff is defined as [13]:
2 u eff = ( u T | u T | + u R | u R | )
Accounting for Equations (3)(6) it turns out:
2 u eff = ( sin   θ s   cos   φ s sin   θ i ) x 0 + sin   θ s   sin   φ s y 0 + ( cos   θ i + cos   θ s ) z 0
The application of the Γxy operator leads to the following relationship:
2 Γ xy   u eff = ( sin   θ s   cos   φ s sin   θ i ) x 0 + sin   θ s   sin   φ s y 0
By substituting the norm of 2Γxy ueff in (7), we can derive the final expression of ρgr:
ρ gr = c / B W sin 2   θ i + sin 2   θ s 2   sin   θ i   sin   θ s   cos   φ s
Considering that the ground range resolution of a conventional monostatic radar is c/(2Bwsinθi), the ρgr/ρgrback ratio is given by (1).
It is worth noting that, if θs is equal to 0°, the ρgr/ρgrback ratio does not depend on φs and is equal to 2 for every θi. Moreover, if we consider φs equal to 0° (or around it), when θs increases starting from 0°, the ρgr/ρgrback ratio grows too much and is always larger than 2. On the other hand, if we assume φs equal to 180° (or around it), when θs increases starting from 0°, the ρgr/ρgrback ratio decreases and is always smaller than 2. These considerations can be useful to choose the sensors observation geometry in order to achieve a satisfactory ground range resolution.
To derive the azimuth resolution, we have to consider the transmitting (vT) and receiving (vR) satellite velocity vectors:
v T = v T   ( cos   δ T   x 0 + sin   δ T   y 0 )
v R = v R   ( cos   δ R   x 0 + sin   δ R   y 0 )
where vT and vR are the amplitudes of vT and vR, respectively, and δT and δR indicate the angles between these vectors and the incidence plane (see Figure 7).
The expression for the azimuth resolution proposed in [13] is:
ρ a = λ 2 ξ int | Γ xy   ω eff |
where ξint is the time interval during which the target is observed, λ is the electromagnetic wavelength and the vector ωeff is given by [13]:
2 ω eff = v Tn | u T | + v Rn | u R |
In (15), vTn represents the component of vT normal to uT and vRn denotes the component of vR normal to uR. vTn can be expressed as:
v Tn = v T ( v T u T | u T | ) u T | u T | = v T v T   sin 2   θ i   cos   δ T   x 0 + v T   cos   θ i   sin   θ i   cos   δ T z 0
By replacing (12) in (16) it turns out:
v Tn = v T   cos 2   θ i   cos   δ T x 0 + v T   sin   δ T y 0 + v T   cos   θ i   sin   θ i   cos   δ T z 0
The projection of vTn on the (x,y)-plane is:
Γ xy   v Tn = v T   cos 2   θ i   cos   δ T x 0 + v T   sin   δ T y 0
The calculation of vRn is similar to that accomplished to derive vTn:
v Rn = v R ( v R u R | u R | ) u R | u R | = v R v R   sin 2   θ s   cos   φ s   cos   ( φ s δ R ) x 0 v R   sin 2   θ s   sin   φ s   cos   ( φ s δ R ) y 0 v R   sin   θ s   cos   θ s   cos   ( φ s δ R ) z 0
By replacing (13) in (19) and by applying the operator that project onto the horizontal plane, it turns out:
Γ xy   v Rn = v R [ cos   δ R sin 2   θ s   cos   φ s   cos   ( φ s δ R ) ] x 0 + v R [ sin   δ R sin 2   θ s   sin   φ s   cos   ( φ s δ R ) ] y 0
By taking into account the relationships (5), (6), (15), (18) and (20), we can determine the projection of the vector 2ωeff, on the horizontal plane:
2 Γ xy   ω eff = 1 H { v T   cos 3   θ i   cos   δ T + v R   cos   θ s [ cos   δ R sin 2   θ s   cos   φ s   cos   ( φ s δ R ) ] } x 0 + + 1 H { [ v T   cos   θ i   sin   δ T + v R   cos   θ s [ sin   δ R sin 2   θ s   sin   φ s   cos   ( φ s δ R ) ] ] } y 0
Having assumed that the two spacecraft fly at the same height, we can make the simplifying hypothesis that the velocity vectors have the same amplitude (i.e., vT = vR), thus obtaining:
2 | Γ xy   ω eff | = v T F ( θ i , θ s , φ s , δ T , δ R )
where:
F ( θ i , θ s , φ s , δ T , δ R ) = cos 2   θ i ( cos 4   θ i   cos 2   δ T + sin 2   δ T ) + cos 2   θ s [ 1 sin 2   θ s ( 1 + cos 2   θ s )   cos 2   ( φ s δ R ) ] + + 2   cos   θ i   cos   θ s [ cos 2   θ i   cos   δ T   cos   δ R + sin   δ T   sin   δ R sin 2   θ s ( cos 2   θ i   cos   δ T   cos   φ s + sin   δ T   sin   φ s ) cos ( φ s δ R ) ]
The final expression for the azimuth resolution of a bistatic system is derived by replacing (22) in (14):
ρ a = H   λ / ξ int v T F ( θ i , θ s , φ s , δ T , δ R )
The azimuth resolution of a conventional monostatic SAR is given by:
ρ a back = D 2 H   λ / ξ int 2 v T   cos   θ i
where D is the dimension of the antenna. Therefore, the ρa/ρaback ratio turns out to be given by (2a).

Acknowledgments

This work was realized under ESA/ESTEC contract 19173/05/NL/GLC, “Use of Bistatic Microwave Measurements for Earth Observation”. We thank Dr. N. Floury from ESA for his useful advice.

References and Notes

  1. Willis, N.J. Bistatic Radar; Artech House: New York, NY, USA, 1991. [Google Scholar]
  2. Pierdicca, N.; Pulvirenti, L.; Ticconi, F.; Brogioni, M. Radar bistatic configurations for soil moisture retrieval: a simulation study. IEEE Trans. Geosci. Rem. Sens 2008, 46, 3252–3264. [Google Scholar]
  3. Burkholder, R.J.; Gupta, J.; Johnson, J.T. Comparison of monostatic and bistatic radar images. IEEE Antenn. Propag. Mag 2003, 45, 41–50. [Google Scholar]
  4. Soumekh, M. Moving target detection in foliage using along track monopulse synthetic aperture radar imaging. IEEE Trans. Image Proc 1997, 6, 1148–1163. [Google Scholar]
  5. Simpson, R.A. Spacecraft studies of planetary surfaces using bistatic radar. IEEE Trans. Geosci. Rem. Sens 1993, 31, 465–482. [Google Scholar]
  6. Garrison, J.L.; Komjathy, A.; Zavorotny, V.U.; Katzberg, S.J. Wind speed measurement using forward scattered GPS signals. IEEE Trans. Geosci. Rem. Sens 2002, 40, 50–65. [Google Scholar]
  7. Moccia, A.; Fasano, G. Analysis of spaceborne tandem configurations for complementing COSMO with SAR interferometry. EURASIP J. Appl. Signal Process 2005, 20, 3304–3315. [Google Scholar]
  8. Sanz-Marcos, J.; Lopez-Dekker, P.; Mallorqui, J.J.; Aguasca, A.; Prats, P. SABRINA: a SAR bistatic receiver for interferometric applications. IEEE Geosci. Rem. Sens. Lett 2007, 4, 307–311. [Google Scholar]
  9. Walterscheid, I.; Ender, J.H.G.; Brenner, A.R.; Loffeld, O. Bistatic SAR processing and experiments. IEEE Trans. Geosci. Rem. Sens 2006, 44, 2710–2717. [Google Scholar]
  10. Klare, J.; Walterscheid, I.; Brenner, A.R.; Ender, J.H.G. Evaluation and optimisation of configurations of a hybrid bistatic SAR experiment between terraSAR-X and PAMIR. Proc. IEEE/IGARSS 2006, 1208–1211. [Google Scholar]
  11. Rodriguez-Cassola, M.; Baumgartner, S.V.; Krieger, G.; Nottensteiner, A.; Horn, R.; Steinbrecher, U.; Metzig, R.; Limbach, M.; Prats, P.; Fischer, J.; Schwerdt, M.; Moreira, A. Bistatic spaceborne-airborne experiment TerraSAR-X/F-SAR: data processing and results. Proc. IEEE/IGARSS 2008, 3, III-451–III-454. [Google Scholar]
  12. Ceraldi, E.; Franceschetti, G.; Iodice, A.; Riccio, D. Estimating the soil dielectric constant via scattering measurements along the specular direction. IEEE Trans. Geosci. Rem. Sens 2005, 43, 295–305. [Google Scholar]
  13. Walterscheid, I.; Brenner, A.R.; Ender, J.H.G. Results on bistatic synthetic aperture radar. Electron. Lett 2004, 40, 1–2. [Google Scholar]
  14. Zavorotny, V.U.; Voronovich, A.G. Bistatic GPS signal reflections at various polarizations from rough land surface with moisture content. Proc. IEEE/IGARSS 2000, 7, 2852–2854. [Google Scholar]
  15. Chen, K.S.; Wu, T.D.; Tsang, L.; Li, Q.; Shi, J.; Fung, A.K. Emission of rough surfaces calculated by the integral equation method with comparison to three-dimensional moment method simulations. IEEE Trans. Geosci. Rem. Sens 2003, 41, 90–101. [Google Scholar]
  16. Pierdicca, N.; Castracane, P.; Pulvirenti, L. Inversion of electromagnetic models for bare soil parameter estimation from multifrequency polarimetric SAR data. Sensors 2008, 8, 8181–8200. [Google Scholar]
  17. D’Errico, M.; Moccia, A. The BISSAT mission: a bistatic SAR operating in formation with COSMO/SkyMed X-band radar. Proc. IEEE Aerospace Conference, Big Sky Montana, USA, March 9–16, 2002; 2, pp. 809–818.
  18. Moccia, A.; Vetrella, S.; Bertoni, R. Mission analysis and design of a bistatic synthetic aperture radar on board a small satellite. Acta Astronaut 2000, 47, 819–829. [Google Scholar]
Figure 1. Left panel: geometric elements that identify the transmitter-target-receiver (Tx-TG-Rx) bistatic configuration. Right panel: sketch of observing configurations suitable for SMC retrieval.
Figure 1. Left panel: geometric elements that identify the transmitter-target-receiver (Tx-TG-Rx) bistatic configuration. Right panel: sketch of observing configurations suitable for SMC retrieval.
Sensors 09 07250f1
Figure 2. Maps of ρgr/ρgrback (left panels) and ρa/ρaback (right panels), for θi = 15° (upper panels) and θi = 35° (lower panels) in the (θs,φs)-plane. Ground range resolution has been upper limited to 4 times the backscattering value for the sake of figure clarity.
Figure 2. Maps of ρgr/ρgrback (left panels) and ρa/ρaback (right panels), for θi = 15° (upper panels) and θi = 35° (lower panels) in the (θs,φs)-plane. Ground range resolution has been upper limited to 4 times the backscattering value for the sake of figure clarity.
Sensors 09 07250f2
Figure 3. Relative position of the satellites at the initial epoch (γ is the yaw-steering angle of the active satellite).
Figure 3. Relative position of the satellites at the initial epoch (γ is the yaw-steering angle of the active satellite).
Sensors 09 07250f3
Figure 4. Bistatic observation geometry at the initial time in terms of zenith incidence (θi) and scattering (θs) angles, receiver (Rx) and transmitter (Tx) near range (NR) and far range (FR) points, and heights (H) of the satellites.
Figure 4. Bistatic observation geometry at the initial time in terms of zenith incidence (θi) and scattering (θs) angles, receiver (Rx) and transmitter (Tx) near range (NR) and far range (FR) points, and heights (H) of the satellites.
Sensors 09 07250f4
Figure 5. Data acquisition (red dots) for a one week scenario considering a passive system acquiring ASAR signal in WSM (upper panel), ISM4 (central panel) and ISM2 (lower panel). Green/blue lines are the Envisat/passive satellite ground tracks.
Figure 5. Data acquisition (red dots) for a one week scenario considering a passive system acquiring ASAR signal in WSM (upper panel), ISM4 (central panel) and ISM2 (lower panel). Green/blue lines are the Envisat/passive satellite ground tracks.
Sensors 09 07250f5
Figure 6. Baseline length (absolute value) and components in the active satellite orbital reference frame versus the orbital period percentage. One orbit is considered for the sake of figure clarity.
Figure 6. Baseline length (absolute value) and components in the active satellite orbital reference frame versus the orbital period percentage. One orbit is considered for the sake of figure clarity.
Sensors 09 07250f6
Table 1. Some parameters of the ASAR illuminator in Image Swath and Wide Swath Modes. Note that both the swath and the ground range resolution refer to the backscattering acquisition; regarding the latter, its value is related to a nominal incidence angle (e.g., 23° for ISM2).
Table 1. Some parameters of the ASAR illuminator in Image Swath and Wide Swath Modes. Note that both the swath and the ground range resolution refer to the backscattering acquisition; regarding the latter, its value is related to a nominal incidence angle (e.g., 23° for ISM2).
ModeSwath [km]Near Range [deg]Far Range [deg]Range Resolution [m]

ISM110515.022.930
ISM210519.226.730
ISM38226.031.430
ISM48831.036.330
WSM4051537100
Table 2. Selected bistatic configuration (in terms of zenith incidence angle, and zenith and azimuth scattering angles), and orbital parameters at the initial epoch (static design).
Table 2. Selected bistatic configuration (in terms of zenith incidence angle, and zenith and azimuth scattering angles), and orbital parameters at the initial epoch (static design).
zenith incidence angle θi [deg]35
zenith scattering angle θs [deg]1
azimuth scattering angle φs [deg]180
Semi-major-axis [km]7,159.48
Eccentricity0.00115
Inclination [deg]98.5
Perigee argument [deg]90
Ascending nodes difference (ΔΩ) [deg]4.20
Mean anomaly difference (ΔM) [deg]0.91
Table 3. Coverage statistics and duty cycle, assuming Envisat in different Swath Modes as illuminator. Minimum and maximum acquisition latitudes, and mean width of the bistatic swath (minimum is 10 km) are reported. Note that for ISM1, ISM2 and ISM3, two latitude belts are identified.
Table 3. Coverage statistics and duty cycle, assuming Envisat in different Swath Modes as illuminator. Minimum and maximum acquisition latitudes, and mean width of the bistatic swath (minimum is 10 km) are reported. Note that for ISM1, ISM2 and ISM3, two latitude belts are identified.
ModeISM1ISM2ISM3ISM4WSM

Duty Cycle [%]13.012.911.921.540.1
Lat. Min. 1[deg]−71.7−66.2−53.2−43.5−71.7
Lat. Max. 1[deg]−48.5−39.6−24.242.971.6
Lat. Min. 2[deg]48.940.525.6
Lat. Max. 2[deg]71.766.251.3
Swath [km]39.937.229.224.139.1

Share and Cite

MDPI and ACS Style

Pierdicca, N.; De Titta, L.; Pulvirenti, L.; Della Pietra, G. Bistatic Radar Configuration for Soil Moisture Retrieval: Analysis of the Spatial Coverage. Sensors 2009, 9, 7250-7265. https://doi.org/10.3390/s90907250

AMA Style

Pierdicca N, De Titta L, Pulvirenti L, Della Pietra G. Bistatic Radar Configuration for Soil Moisture Retrieval: Analysis of the Spatial Coverage. Sensors. 2009; 9(9):7250-7265. https://doi.org/10.3390/s90907250

Chicago/Turabian Style

Pierdicca, Nazzareno, Ludovico De Titta, Luca Pulvirenti, and Giuliano Della Pietra. 2009. "Bistatic Radar Configuration for Soil Moisture Retrieval: Analysis of the Spatial Coverage" Sensors 9, no. 9: 7250-7265. https://doi.org/10.3390/s90907250

APA Style

Pierdicca, N., De Titta, L., Pulvirenti, L., & Della Pietra, G. (2009). Bistatic Radar Configuration for Soil Moisture Retrieval: Analysis of the Spatial Coverage. Sensors, 9(9), 7250-7265. https://doi.org/10.3390/s90907250

Article Metrics

Back to TopTop