On the Impact and Mitigation of Signal Crosstalk in Ground-Based and Low Altitude Airborne GNSS-R

Global Navigation Satellite System Reflectometry (GNSS-R) is a powerful way to retrieve information from a reflecting surface by exploiting GNSS as signals of opportunity. In dual antenna conventional GNSS-R architectures, the reflected signal is correlated with a clean replica to obtain the specular reflection point delay and Doppler estimates, which are further processed to obtain the GNSS-R product of interest. An important problem that may appear for low elevation satellites is signal crosstalk, that is the direct line-of-sight signal leaks into the antenna dedicated to the reflected signal. Such crosstalk may degrade the overall system performance if both signals are very close in time, similar to multipath in standard GNSS receivers, the reason why mitigation strategies must be accounted for. In this article: (i) we first provide a geometrical analysis to justify that the estimation performance is only affected for low height receivers; (ii) then, we analyze the impact of crosstalk if not taken into account, by comparing the single source conditional maximum likelihood estimator (CMLE) performance in a dual source context with the corresponding Cramér–Rao bound (CRB); (iii) we discuss dual source estimators as a possible mitigation strategy; and (iv) we investigate the performance of the so-called variance estimator, which is designed to eliminate the coherent signal part, compared to both the CRB and non-coherent dual source estimators. Simulation results are provided for representative GNSS signals to support the discussion. From this analysis, it is found that: (i) for low enough reflected-to-direct signal amplitude ratios (RDR), the crosstalk has no impact on standard single source CMLEs; (ii) for high enough signal-to-noise ratios (SNR), the dual source estimators are efficient irrespective of the RDR, then being a promising solution for any reflected signal scenario; (iii) non-coherent dual source estimators are also efficient at high SNR; and (iv) the variance estimator is efficient as long as the non-coherent part of the signal is dominant.


Introduction
The use of Global Navigation Satellite System (GNSS) signals spans over a plethora of applications, from its original navigation purposes [1], to precise time synchronization, remote sensing of the ionosphere, radio occultation for atmospheric sounding, or reflectometry (GNSS-R) for Earth observation [2,3], to name a few. The concept of GNSS-R of interest in this contribution is roughly 30 years old [4]. Instead of considering reflected signals (i.e., multipath) as a nuisance or interference that the receiver must eliminate to obtain the best positioning solution (which is the case in standard GNSS), reflected signals are used as an additional source of information (i.e., signals of opportunity) for a different purpose. As stated in [2], by correctly exploiting the scattered signal, the receiver can be used as an altimeter or a scatterometer, able to estimate height, roughness, or other properties of the reflecting surface, as well as the corresponding surface wind, soil moisture, tenna cross-polarization, and receiver geometry. Because the GNSS-R product is obtained from the resulting DDM, crosstalk may significantly degrade the system performance as the reflected signal may be distorted (i.e., similar to multipath in GNSS receivers), and therefore, it is fundamental to have a theoretical analysis of such performance loss and solutions to counteract its impact. This is a problem that was already pointed out in [40], where the authors proposed an estimation technique, the so-called variance estimator (VE), which aims to cancel the coherent part of the signal. This approach may have some limitations; for instance, it may not be suitable for reflected signals with a strong coherent part such as specular reflections over ice. In addition, to the best of the authors' knowledge, such a VE has not been properly characterized in the literature from an estimation point of view. The correct analysis is obtained by comparing the estimator performance with the corresponding Cramér-Rao bound (CRB).
Considering the crosstalk problem in cGNSS-R architectures, this article provides the following contributions: • From an estimation point of view, we assess the impact that a possible crosstalk may have on the reflected signal time-delay estimation performance, if not correctly accounted for. This is done by comparing the single source conditional maximum likelihood estimator (CMLE) performance in a dual source context with the corresponding single source CRB. • By resorting to dual source estimators, we show the optimal time-delay estimation performance in a dual source crosstalk context, as a function of path separation and the reflected-to-direct signal amplitude ratio (RDR), compared to the corresponding dual source CRBs. • As a complementary analysis of practical importance, we assess the robustness of such dual source estimators under a misspecified number of sources, i.e., when two sources are estimated, but crosstalk is not present in the signal. • The performance of the VE, proposed in [40] for crosstalk mitigation, is compared to the CRB and non-coherent dual source estimators, which to the best of the authors' knowledge is an important missing estimation performance analysis in the literature.

Preliminaries: A Geometrical Analysis
First, we provide a geometrical analysis of the realistic time difference between reflected and direct paths: ∆τ = |τ 1 − τ 0 | (where τ 0 is the delay for the reflected path and τ 1 the delay for the crosstalk LOS signal) for different receiver heights and satellite local elevations (h and e in the simplified geometry presented in Figure 1). We compute such that ∆τ using the geometrical model in [4], and the corresponding results are depicted in Figure 2 (i.e., the right-hand side plot is a zoom in of the left-hand side plot, to clearly see the low receiver height cases). Notice that these plots are given w.r.t. the local elevation at the specular point. We also show in Figure 2 the values equivalent to ∆τ = 1, 1/2, 1/4 and 1/8 L1 C/A chips, ∆τ = 1 and 1/2 L5-I chips (i.e., the same for the Galileo E5a or E5b), and the ∆τ corresponding to the first zero-crossing of the E5 ACF, i.e., refer to the table in Figure 3 (left) for the corresponding zero-crossing values. For completeness, we show the ACF for representative GNSS signals in Figure 3 (right).
Notice that the widest ACF function is the one obtained for the legacy GPS L1 C/A signal, potentially leading to the worst case (i.e., the largest range of scenarios where a crosstalk may induce an estimation performance degradation). Therefore, we must take into account that (i) only values below ∆τ = 1 L1 C/A chip are of interest (i.e., otherwise, there is almost no intersection of the reflected and direct ACFs), and (ii) only satellite local elevations that are relatively low may induce a crosstalk (i.e., for high satellite local elevations, this is very unlikely if using a standard antenna geometry). From the results in Figure 2, we can conclude that the estimation performance analysis under crosstalk conditions makes sense only for low height receivers (h < 1 km). In addition, if large bandwidth signals such as GPS L5-I or Galileo E5a are used, then this restricts the impact to ground-based receivers (h < 100 m). In the case of exploiting a Galileo Alternative BOC (AltBOC) E5 signal [41,42], the first zero of the ACF is at ≈0.017 L1 C/A chips, and therefore, the geometrical analysis suggests that a possible parameter estimation performance degradation induced by crosstalk is very unlikely, which is a first interesting result for crosstalk mitigation.

Single and Dual Source GNSS Signal Model
We consider that a transmitter T and a receiver R (see Figure 1) have uniform linear motions such that the positions can be described as p T (t) = p T + v T t and p R (t) = p R + v R t, where p and v are the position and velocity vectors, respectively. Under such conditions, the distance between T and R at instant t can be approximated by a first order distancevelocity model: where d is the T-to-R absolute distance when t = 0, v is the T-to-R radial velocity, τ is the time-delay due to the propagation path, (1 − b) is the dilatation induced by the Doppler effect, and c is the speed of light in a vacuum. A band-limited signal s(t), with bandwidth B, is transmitted by T over a carrier frequency f c (λ c = c/ f c ). Notice that this signal model encompasses any GNSS signal. Using (1), the single source complex analytic signal at the output of the receiver's antenna is: with i = 0 for the reflected signal, of interest in this study, and i = 1 for the crosstalk signal, which is the LOS signal in this study. w R (t) is a zero-mean white complex circular Gaussian noise, i are the amplitude and phase of the complex coefficients induced by the propagation characteristics, polarization mismatches, antenna gains, and the reflecting surface for the reflected signal. Under the narrowband signal hypothesis, i.e., the time-frequency product B · T c (with T c the coherent integration time) is smaller than the inverse Doppler term c/v, the Doppler effect on the band-limited baseband signal s(t) is usually neglected, s((1 − b)(t − τ)) ≈ s(t − τ). Therefore, the baseband output of the receiver's Hilbert filter with bandwidth F s containing one signal (i.e., direct or reflected) is: where Considering the acquisition of N = N 2 − N 1 + 1 samples at a sampling rate T s = 1/F s , the discrete single source conditional signal model (CSM) is, Similarly, the dual source complex baseband output of the receiver's Hilbert filter with bandwidth F s containing a single specular reflection and a direct signal (e.g., when the antenna pointing downwards also receives the LOS signal) is: and the corresponding discrete dual source CSM is, with A(η 0 , η 1 ) = [a(η 0 ), a(η 1 )] and α T = ρ 0 e jφ 0 , ρ 1 e jφ 1 .

Single and Dual Source CRBs
Considering the optimal single source scenario or the dual source crosstalk case, the user must estimate the parameters of the signal of interest θ i . Therefore, it is fundamental to know the achievable estimation performance limits by computing the corresponding single and dual source CRBs, as well as to assess the estimators' performance w.r.t. such CRBs.
The CRB is given by the inverse of the Fisher information matrix (FIM). If we define , was recently derived in [43] (i.e., and given in [44] for the general wideband signal case). The inversion of F θ i | i ( i ) proposed in [43] leads to a closed-form CRB expression for the delay, Doppler, phase, and amplitude estimation, where a noteworthy feature is that the CRBs are expressed w.r.t. the baseband signal samples, therefore being valid and easy to compute for any GNSS signal. Such CRBs were studied for different GNSS signals in [41,42].
In the dual source case, the parameters to be estimated are T = [σ 2 n , θ T 0 , θ T 1 ]. The dual source FIM is: where: (i) F θ 0 | ( ) and F θ 1 | ( ) correspond to the FIMs of the signals when they are totally decoupled, derived in [43]; and (ii) F θ 0 ,θ 1 | ( ) = F θ 1 ,θ 0 | ( ) T characterizes the interference between both signals. A compact closed-form expression of F | ( ), which depends only on the baseband signal samples, in the vein of [43], was very recently derived in [45].

Single and Dual Source CMLEs
In the single source case, the CMLEs of delay, Doppler, and phase (1S-CMLE) for each individual satellite are expressed as: Notice that the phase CMLE is given by the argument of the cross-ambiguity function evaluated at the delay and Doppler CMLEs. Then, if the delay Doppler CMLE reaches its asymptotic performance, then so does the phase estimate.
The dual source CMLE (2S-CMLE) that maximizes the likelihood is such that ( η 0 , η 1 ) maximizes the projection of the received signal onto the subspace defined by the mixing matrix A [46]. As a consequence, the estimated parameters ( ρ 0 , ρ 1 , φ 0 , φ 1 ) depend on the delay/Doppler estimates, where the projectors are defined as to be computationally complex, it is asymptotically efficient for the dual source Gaussian CSM (i.e., in the large sample and/or high SNR regimes [47,48]).

CLEAN-RELAX Estimator
In order to deal with the 2S-CMLE heavy computational load, alternative estimators can be used. A possible solution is the RELAX version of the so-called CLEAN estimator algorithm (CRE) [49]. This estimator assumes that the two signals are decoupled. The method first estimates the parameters ( τ 0 , b 0 , ρ 0 , φ 0 ) of the strongest signal, and it subtracts this estimated signal from the measured input, then estimates the secondary signal parameters ( τ 1 , b 1 , ρ 1 , φ 1 ). These two estimation steps are iterated until the cost function, which is the resulting likelihood, is stabilized. In GNSS, the same algorithm is called the multipath estimating delay lock loop (MEDLL) and is typically used to mitigate multipath [50]. The CRE can be seen as a sub-optimal estimator because it works only under certain assumptions. Indeed, considering that the two signals can be decoupled stands only when ∆τ is large enough compared to the sharpness of the ACF. This estimator is not expected to perform well when the two signals are very close in time. On the other hand, when the signals are clearly apart, the CRE asymptotically behaves as the 2S-CMLE.

Non-Coherent and Variance Estimators
It is common practice in GNSS-R to perform the non-coherent integration of individual coherent 1S-CMLEs up to some tens of seconds [40]. Indeed, the non-coherent approach is useful to deal with long integration times and when the signal itself is non-coherent, which is typically the case when dealing with complex reflecting surfaces. Considering a total signal duration T = N nc T c with T c the individual CMLE coherent observation time in (10), we can express both estimators as follows (superscripts (·) C and (·) NC for coherent and non-coherent). First, let us define y k (η , the output of the coherent matched filter at k th time slot for a coherent observation time T c , then: The same applies for the dual source CMLEs, which can be formulated both in a coherent and non-coherent manner, as well as their approximation using a CRE approach. Regarding the non-coherent VE derived in [40], it is formulated as: where the coherent part of the signal is subtracted from the non-coherently averaged one.

Results and Discussions on Coherent Estimation
In this section, the crosstalk impact on coherent estimation is investigated. To do this, a prior study on the 1S-CMLE is used to estimate the reflected path in the presence of crosstalk; this aims at asserting the need for a mitigation technique based on either 2S-CMLE or CRE. Such estimators' performances are then presented in a properly specified scenario (presence of crosstalk) and then in a misspecified one (no crosstalk).

Crosstalk Impact on GNSS-R Time-Delay Estimation
Exploiting the single and dual source CRBs introduced in Section 2.3, we can properly characterize the impact that two sources (crosstalk) may have on the reflected or direct signal parameter estimation w.r.t. the single source case. We consider in this section a specular reflection and a receiver height that lead to crosstalk interference (i.e., this is the case for low height receivers, of interest in this article). The two parameters investigated in this study are the reflected vs. direct signal delay difference, ∆τ, and the ratio between signal amplitudes |Γ| where Γ is a complex number such that Γ = ρ 1 e jφ 1 /ρ 0 e jφ 0 . The first parameter, ∆τ, is determined by the system geometry, that is receiver height and satellite elevation (refer to Section 2.1). The second parameter, |Γ|, depends on the antenna radiation pattern, the polarization mismatch between the leaking signal and the antenna, and the incidence angle (again via the system geometry). Given the fact that we consider a system with two antennas, we assume that the modulus |Γ| is smaller than one, that is the reflected signal has more energy than the crosstalk signal from the LOS. The other parameters in the signal definition (7) are set to a fixed value for the rest of this study: the Doppler frequency of the reflected signal F d,0 = f c · b 0 = 30Hz, the Doppler frequency of the LOS signal F d,1 = f c · b 1 = 50Hz, the phase of the reflected complex amplitude φ 0 = π/3, and the phase of the LOS complex amplitude φ 1 = π/4.

Analysis Setup
As representative case studies, we consider |Γ| = {0.01, 0.1, 0.3, 0.5}, ∆τ = {1/8, 1/4, 1/2, 1} L1 C/A chips for a GPS L1 C/A signal, and ∆τ = {1/2, 1} L5-I chips for a GPS L5-I signal, covering a wide range of possible values. For each pair {∆τ, |Γ|}, we obtain the estimation performance for the GPS L1 C/A signal scenario, as well as for a GPS L5-I signal (equivalent to Galileo E5a or E5b). Notice that the L1 C/A chip is 10 times larger than the E5a/L5-I chip (i.e., BPSK(1) vs. BPSK (10)). For illustrative purposes, we show the squared modulus of the cross-correlation function (CCF) of a BPSK signal for different ∆τ and |Γ| = 0.5 in Figure 4.  In this section, two scenarios are of interest: • Case (1): suboptimal single source estimation, that is two sources are present, but the corresponding estimator considers only one source. • Case (2): optimal dual source estimation, that is it is known that two sources are present, and the corresponding estimator is matched to this.
In the sequel, we show the root mean squared error (RMSE) results for the time-delay estimation of the stronger signal, being for instance the main parameter of interest for GNSS-R altimetry, with respect to the signal-to-noise ratio (SNR) at the output of the matched filter.

Suboptimal Single Source Estimation
This analysis provides insights into the performance loss when the possible crosstalk is not taken into account. We compare the RMSE results obtained with the 1S-CMLE introduced Section 3.1 with the single source CRB given in Section 2.3. Notice that ∆τ = 1/2 L5-I chips = 0.05 L1 C/A chips and ∆τ = 1 L5-I chips = 0.1 L1 C/A chips. The results are shown in Figure 5.
In the sequel, in all the figures presented, the CMLE/CRE behavior can be divided into three distinct regions: (i) in the a priori region, at a low SNR, the estimate behaves as a random variable uniformly distributed on the search area and leads to a very large RMSE; in this case, the estimate is useless for the user; (ii) the threshold region (i.e., where the estimator RMSE drops significantly) is useful to know the minimum SNR at the output of the matched filter to be in optimal operation conditions; and (iii) the asymptotic region is where the estimator behaves correctly, given that the estimator is asymptotically efficient (i.e., convergence to the CRB).  From the GPS L1 C/A results, we can draw three main conclusions: (i) In a dual source context, the single source estimators only converge to the bound for very low values of |Γ| (i.e., 0.01 and 0.1 in the ∆τ = 1, 1/2, 1/4 L1 C/A chips results or only 0.01 in the case of ∆τ = 1/8 L1 C/A chips);. In this case, there is no impact of a possible crosstalk. (ii) For higher values of |Γ|, the 1S-CMLE never converges to the 1S-CRB, and therefore, the crosstalk induces a performance degradation. This significantly impacts the receiver performance for |Γ| = 0.5 and for ∆τ = 1/8 L1 C/A chips. It is worth pointing out that in such conditions, the estimator RMSE reaches a constant floor, which is typical of a biased estimator (easy to spot in Figure 5d-f). This is due to the fact that for a significant value of |Γ|, the crosstalk distorts the CCF and shifts its maximum. Consequently, even if, for larger |Γ|, the threshold seems to be earlier, the estimated value turns out to be biased. (iii) These results justify the need for dual source estimators or other robust solutions if a possible strong crosstalk is expected. For the GPS L5-I case, we can see the same behavior for ∆τ = 1 L5-I chips and a performance saturation for ∆τ = 1/2 L5-I chips (even for |Γ| = 0.1). This is because the achievable RMSE with an L5-I signal is far below the corresponding one with the L1 C/A signal (i.e., the L5-I CCF is sharper than the L1 C/A CCF); therefore, small deviations on the CCF have a larger impact on the final performance.

Optimal Dual Source Estimation
This analysis provides insights into the optimal achievable performance with respect to the single source case. Recall that with respect to the previous case, we know that two sources are present, and the parameters associated with both sources are estimated, either with a CRE or a 2S-CMLE when the CRE does not converge because both sources are too close in time. We compare these dual source estimators with the single and dual source CRBs. Results are shown in Figure 6.
From the results shown in Figure 6, we can draw different conclusions: (i) First, we notice that at a sufficiently high SNR, the estimators converge to the CRB, irrespective of the values of ∆τ and |Γ|, which is a known result for CMLE [48]. This implies that using a dual source estimator and under certain conditions, we can mitigate the impact of any possible crosstalk. (ii) It is also worth pointing out the difference between the 1S-CRB and 2S-CRB, which is limited for most of the range of ∆τ considered in this study. This implies that once the estimator converges to the bound, the crosstalk impact on the estimation performance is limited as well (with respect to the optimal single source case). (iii) The third point is the ability of the CRE to converge to the bound for values as low as ∆τ = 1/4 chips (F s = 8 MHz) for the GPS L1 C/A or 0.1 L1 C/A chips (F s = 20 MHz) for the GPS L5-I. For the GPS L1 C/A case with ∆τ = 1/8 chips or the GPS L5-I case with ∆τ = 0.05 L1 C/A chips, the CRE did not converge, the reason why the results were obtained with the 2S-CMLE. Notice that this threshold can be increased with a better resolution, that is a higher signal bandwidth. (iv) Notice also the strong impact that |Γ| has on the convergence threshold. While for low values such as |Γ| = 0.01 or 0.1, the convergence is very fast, for larger values, this may imply several dB of difference on the optimal operation point. (v) Lastly, notice that the CMLE performance cannot be below the CRB, as is the case in the plots. This is because of a reduced parameter space search in the implementation, to avoid a huge computational burden. In any case, we can conclude that (i) for low values of |Γ| and using a dual source estimator, crosstalk has almost no impact, and (ii) for large values of |Γ|, we need a higher SNR in order to mitigate the crosstalk.

On the Dual Source Estimators' Robustness: Misspecified Number of Sources
In the previous Sections 4.1.2 and 4.2, we provided a time-delay estimation performance analysis in a dual source context. First, we considered a single source estimator, in order to assess the impact of a possible crosstalk (i.e., not taken into account) w.r.t. the optimal single source case, and then, we considered dual source estimators to assess the impact w.r.t. the single source case in an optimal setup. A fundamental question of practical interest is: what is the performance of dual source estimators when only one source is present? In other words, what is the robustness of dual source estimators when the number of sources is misspecified? In practice, this is the case of interest, because the crosstalk is not always present.
We assess the performance obtained with a CRE considering two sources, but a single source is present in the signal. For the simulation, we used the same parameters of Doppler frequency and phase as before for the reflected signal, and we simply set the modulus |Γ| to zero. Results for both GPS L1 C/A and GPS L5-I signals are shown in Figure 7, where it is clear that the CRE is robust in the case of a misspecified number of sources. Indeed, it properly estimates the reflected signal parameters and sees the rest as noise. This implies that the CRE is a promising option to mitigate the possible crosstalk.

Results and Discussion on Non-Coherent Estimation
Finally, we want to link the variance method to the previous CRBs and maximum likelihood estimation, that is to the corresponding estimation performance, which, to the best of our knowledge, has not been done in the literature. In addition, we also compare the results with a non-coherent CRE (NC-CRE), which must be used in practice to perform long integration times.

Analysis Setup
Recall that the VE given in (17) is designed to eliminate the contribution of the coherent part of the signal. Therefore, to cover a set of representative case studies, we consider four different scenarios, where for all of them, there are two GPS L1 C/A signals, F s = 8 MHz, ∆τ = 1/4 L1 C/A chips, |Γ| = 0.5, and a total signal duration of 20 ms (i.e., T = 20 ms, T c = 1 ms, and N nc = 20; one PRN is 1 ms long). The direct LOS signal is considered to be coherent during the observation time, that is it has the same phase during T. For the reflected signal whose delay we want to estimate in this study and taking into account that we use 20 GPS L1 PRN codes, we consider the following: Notice that these four cases cover from the fully non-coherent reflected signal in Scenario #1 to the almost coherent reflected signal in Scenario #4. In these simulations, we arbitrarily set the complex amplitude phase of the reflected path equal to the LOS signal's one (one could set it to another fixed value without changing the general outputs of this study).

On the Estimation Performance of the Variance Estimator and Non-Coherent CRE
The RMSE results obtained for the four scenarios with both the VE and NC-CRE are shown in Figure 8. First, notice that the VE performance for a fully non-coherent scenario (#1) reaches the corresponding dual source CRB (i.e., which is almost equivalent to the single source CRB in these scenarios; refer to the middle-left plot of Figure 6) at around 17 dB; then, in this case the estimator is asymptotically efficient and behaves similar to the 1S-CMLE without crosstalk. In contrast, when the coherent part of the reflected signal increases, its performance drastically degrades, and the VE is not a valid approach anymore (i.e., Scenarios #3 and #4), as expected by construction. Therefore, the VE is only useful if considering reflected signals with a dominant non-coherent component such as ocean surfaces, but it is not for specular reflections over ice or for certain land surfaces with a significant coherent signal part. Regarding the NC-CRE, it is asymptotically efficient irrespective of the signal coherency, therefore being a very powerful seamless alternative to the VE. Notice that the NC-CRE for the estimation of the reflected delay does not fully converge to the corresponding 2S-CRB until the convergence of the estimation of the secondary delay (in this case, the LOS delay), at 30 dB, but the performance degradation between the first threshold and 30 dB is minor. Then, both methods have their pros and cons, the VE being easier to implement and the NC-CRE being more robust to the type of surface.

Conclusions
In this contribution, we analyzed, from an estimation point of view, the possible crosstalk in ground-based and low altitude airborne GNSS-R. Both the crosstalk impact on standard architectures and possible mitigation strategies were discussed. First, it was found that the single source CMLE performance is not affected by crosstalk if the reflected-to-direct signal amplitude ratio is low enough; otherwise, there is a performance degradation that must be accounted for. Second, it was shown that both CRE and dual source CMLE are asymptotically optimal irrespective of the amplitude ratio, then being a promising crosstalk mitigation strategy. Notice that using the CRE is limited by the time separation among both sources. However, such analysis was performed in a coherent signal context, and in practice, the reflected signal is non-coherently averaged over long integration times. That is the reason why a complementary analysis was provided to assess the performance of an NC-CRE, which was shown to be also asymptotically efficient irrespective of the signal coherency. In contrast, the VE used in the literature was shown to be a valid approach only for non-coherent signals, for which this estimator is asymptotically efficient. Otherwise, for instance in close-to-specular reflections, the VE suffers a significant performance degradation. Then, we can conclude that both the NC-CRE and VE are good ways to mitigate the crosstalk impact, with their pros and cons. Notice that the performance analysis provided in this study investigated the impact of the LOS signal on the reflected one, but the results are also valid to analyze the impact of the reflected signal on the LOS one. The only difference is that the LOS signal is stronger than the reflected signal, or vice versa.
As a side note, the robustness of the CRE for a misspecified number of sources was also discussed. It is also worth pointing out that to avoid the crosstalk, a possible alternative is to exploit large bandwidth Galileo E5 signals, for which a possible parameter estimation performance degradation induced by crosstalk is very unlikely. Finally, another possible way to avoid such crosstalk, not discussed in this contribution, may be to consider antenna arrays with adaptive beamforming capabilities, where a beam can be steered towards the reflection point and a null placed on the direct signal direction of arrival.

Conflicts of Interest:
The authors declare no conflict of interest.