GNSS Spoofing Detection Based on Signal Power Measurements: Statistical Analysis

A threat to GNSS receivers is posed by a spoofing transmitter that emulates authentic signals but with randomized code phase and Doppler values over a small range. Such spoofing signals can result in large navigational solution errors that are passed onto the unsuspecting user with potentially dire consequences. An effective spoofing detection technique is developed in this paper, based on signal power measurements and that can be readily applied to present consumer grade GNSS receivers with minimal firmware changes. An extensive statistical analysis is carried out based on formulating a multihypothesis detection problem. Expressions are developed to devise a set of thresholds required for signal detection and identification. The detection processing methods developed are further manipulated to exploit incidental antenna motion arising from user interaction with a GNSS handheld receiver to further enhance the detection performance of the proposed algorithm. The statistical analysis supports the effectiveness of the proposed spoofing detection technique under various multipath conditions.


Introduction
The received GNSS signal power at the output of a 3 dB gain hemispherical linearly polarized antenna at ground level is approximately −130 dBm [1].This makes GNSS receivers susceptible to nearby noise jammers and standoff spoofers (SS) that can easily transmit power levels well above −130 dBm.A high processing gain based on a long integration time is often the only option available to overcome a noise jammer.Nevertheless, if the GNSS receiver undergoes random motion, then the channel decorrelates quickly such that attaining such large processing gains to overcome jamming is neither feasible nor desirable from an operational perspective.Also a jammer is relatively easy to locate with radio direction finding and to potentially disable as its spectrum is significantly larger than the ambient noise [2,3].In addition, the noise jammer is at least detectable as the spectral power in the affected GNSS receiver band will be abnormally high.Hence the jammer can deny service but the user is aware of being jammed, limiting the damage potential of the jammer.A more insidious threat is the standoff spoofer that broadcasts a set of replicas of the authentic satellite vehicle (SV) signals visible to the mobile GNSS receiver [2].Disruption of GNSS services is achieved by randomly modulating the code phase over a small region of the overall Code Delay Space (CDS) that is commensurate with a target area.The spoofing attack is assumed to happen during the acquisition stage.Therefore, it is not possible to identify the SS signal based on the code phase as corresponding to an outlier navigation solution.The SS is assumed to remain synchronized with currently visible GNSS signals and then transmit a set of signals that would correspond to the typical GNSS signals observable by a receiver in the target area.Note that an effective SS does not necessarily synthesize a specific counterfeit location for a specific GNSS receiver but rather aims to disrupt GNSS services over a general target area by matching the Doppler offset of the replicated SV signals and adjusting the code phase such that it is commensurate with the intended target region.Hence the GNSS receiver cannot easily detect the contribution of these counterfeit signals as obvious outliers.An unaware receiver computes the navigation solution based on the SS generated counterfeit signals which are passed on to the user as being reliable with potentially damaging consequences.
GNSS receivers tethered to a wireless data service provider will typically provide the user with an aided-GNSS (AGNSS) service, significantly reducing the CDS corresponding to a physical area of several square kilometres [4].Hence there is a diminishing gain for the spoofer attempting to affect a larger target area than this.Hence the counterfeit SS navigation solutions will be construed as plausible.As such, receiver-autonomous integrity monitoring (RAIM) and fault detection and exclusion (FDE) are ineffective in discriminating signals sourced from the envisioned SS [5].
The typical handheld consumer GNSS receiver coherently integrates the signal for about 10 to 20 ms resulting in a correlation peak in the CDS that has a spread in Doppler of about 100 Hz, which is commensurate with the Doppler spread of typical urban traffic (<50 km/hr) [6].Even if the GNSS receiver is equipped with other ancillary sensors such that the receiver velocity vector is independently known, this cannot be used to discriminate the SS signal as multipath Doppler spreading is approximately equivalent for both the SS and the authentic SV signals.
Note that the receiver processing gain used for suppressing a jamming signal is not effective in the case of the spoofer signal.Consequently, the spoofer transmit power can be orders of magnitude less than that of the noise jammer, which makes the spoofer source much more difficult to locate and disable through radio direction finding and beam forming.
The objective of this paper is to address a computationally efficient processing technique that can be added to relatively unsophisticated consumer grade GNSS receivers to discriminate the spoofer signals transmitted by an SS.The proposed processing is based on estimating and comparing the receiver signal power with a set of thresholds to verify the authenticity of the signal.The detection problem is formulated based on a Rayleigh fading multipath scenario.Nevertheless, it is shown that although suboptimal, the deduced expressions can be utilized for spoofing detection in a generalized Rician multipath channel with minimal performance degradation.
The proposed technique is further extended to include incidental motions of the handheld receiver, instigated through the user interaction with the handset device, in the form of spatial translation and polarization rotation.User interaction with the handheld creates variability in the antenna response, which can be transformed into a diversity gain that adds to the general processing gain of the receiver [7][8][9][10].This processing gain enhances the estimation of the received signal power of the correlation peaks, that, is necessary information in spoofer discrimination.A case study based on GPS L1 C/A signals is developed to demonstrate the effectiveness of the proposed technique.Nevertheless, this technique can be directly extended for other GNSS signal formats such as GPS L2 C/A and GLONASS.
The rest of the paper is organized as follows.In Section 2 the system definition and the assumptions are given.Section 3 formulates a multihypothesis detection problem and focuses on the statistical evaluation of the proposed technique, with the conclusions provided in Section 4.

System Definition
This paper considers the analysis of individual GNSS satellite signals, while realizing that simultaneous processing of the available GNSS signals provide extra diversity that can be used to further improve the performance of the proposed spoofer detection technique.The received complex GNSS baseband signal is denoted here by where the signal component of g(t) is represented by s(t) = A(t)s o (t), where "t" is time, A(t) is the channel response to the incident signal at the antenna, and s o (t) is the complex baseband component of the satellite signal, which can be written as where d(t) is the navigation data modulation, c(t) is the Pseudo Random Noise (PRN) code, τ is the code phase, Δ f represents carrier frequency offset (due to the Doppler of the GNSS signal as well as any frequency offset of the receivers local oscillator), and ψ is the initial phase offset.s o (t) is known to the receiver except for the navigation data, the code phase, the carrier frequency offset, and the initial phase offset ψ.The received signal, g(t), is corrupted by additive noise (WGN) which has an equivalent complex baseband representation denoted by w(t).It is assumed that w(t) is a complex normal random process, independent of the signal, and has a Power Spectral Density (PSD) that is constant within the bandwidth of the received signal.
The GNSS receiver integrates a temporal snapshot of g(t) over the interval of t ∈ [0, T I ], where T I is typically smaller than the duration of one navigation data bit (20 ms).The signal snapshot of g(t) is collected by the receiver and then despread by a locally generated copy of s o (t) during the initial acquisition.The initial acquisition is typical of a multihypothesis detection in which the receiver searches the Code Doppler Space (CDS) for the frequency offset Δ f and the code delay, τ [11,12].Note that the initial phase offset ψ is not known to the receiver during the initial acquisition and as such the output of the despreading matched filter is a random complex variable.
The despread baseband signal samples at a correlator output are represented by where " * " is a complex conjugate, the subscript "n" denotes the nth signal sampling interval which extends over t ∈ [(n−1)T I , nT I ], and s n , w n are the postintegration signal and the WGN components, respectively.In addition, τ and Δ f represent the estimated code phase and Doppler based on the initial acquisition which consists of a maximum likelihood search over the CDS of a signal sample, x n;τ,Δ f , such that where { τ, Δ f } are the maximum likelihood (ML) estimates of the true code phase and Doppler frequency, respectively.The estimated code phase and Doppler are then passed on to the tracking loops to facilitate further receiver processing.Consequently, N signal samples, namely, can be collected and used for spoofer detection.

Theoretical Analysis of Spoofer Detection
A hypothetical scenario is considered based on an SS transmitting spoofing signals in an urban environment as shown in Figure 1.The authentic signal and the spoofer signal are affected by multipath fading and therefore, the received signal power is random in space and polarization.In other words, multipath fading results in signal power fluctuation when the receiver is spatially translated or undergoes polarization changes due to rotation.Unlike the authentic signal power which is insensitive to signal power variations arising from pathloss in the target area (this is due to the fact that the satellite-receiver separation is approximately unchanged over a period of several minutes), the spoofer signal power varies with variation in the spooferreceiver separation.An empirical model of order n can be utilized to model the spoofer signal power variation due to pathloss as where R1 is a reference range, d is the spoofer-receiver range, n is the pathloss exponent, ρ is the average spoofer SNR at d, and ρ (sp) R1 is the average received spoofer SNR at d = R1 in dBs.Note that, for the spoofer to be effective, the average received spoofer signal power needs to be higher than that of the authentic signal in the target area.Therefore, the received signal power from a standoff spoofer varies significantly with range due to pathloss, meaning that the spoofer signal power is abnormally higher than that of the authentic signal when the receiver is in the proximity of the standoff spoofer.This characteristic of the spoofer signal can be exploited to limit the effectiveness of the SS in its target area based on comparing the measured signal power against a preset threshold.
As stated earlier, a receiver records N signal samples with each of these n = 1, . . ., N signal samples belonging to one of the three hypotheses, namely, the noise hypothesis H0, the authentic signal hypothesis H1, and the spoofer signal hypothesis H2 as H0 : where the normalization t+T t |s o (t)| 2 dt = 1 is assumed and τ, Δ f are suppressed for notational convenience.A (a) (t) and A (sp) (t) represent channel gains associated with the authentic and the spoofer signals, respectively.Consequently, a detection variable, r = h(x), can be formulated to decide between the three hypotheses of (6) based on comparing "r" with a set of thresholds, ρ 1 , ρ 2 , as shown in Figure 2. Note that h(x) is a function that maps the measured signal samples x to a single variable, r, which is a sufficient statistic with respect to H0, H1, and H2.As will be shown in Section 3.1, r can be found from the probability density functions (PDF) of x [13] or alternatively from the PDFs of "r" which are denoted here by H0 : f r|H0 (r) H1 : f r|H1 (r) H2 : f r|H2 (r), (7) where f (•) denotes a PDF.

International Journal of Navigation and Observation
One optimization criteria for determining the thresholds (ρ 1 , ρ 2 ) is based on minimizing the probability of error, namely, where P(Hi) for i = 0, 1, 2 are the probabilities of H0, H1, and H2 states and P(Hi|Hi) denotes the conditional probability that indicates that of deciding Hi if Hi is correct.Consequently, where F r|Hi (•) denotes a cumulative distribution function (CDF) of the random variable "r" under H i .As can be seen from ( 9), P e is a function of the authentic signal, the spoofer, and the noise statistics.Therefore, any optimization based on minimizing the probability of error hinges on knowing the spoofer signal statistics, which is not available to an unsuspecting receiver given the capricious nature of a spoofer.
Alternatively, a second optimization can be made based on maximizing the probability of detection for a given probability of false alarm.Assuming ρ 2 > ρ 1 , the threshold ρ 1 can be determined based on selecting a probability of false alarm P FA1 as Therefore, As is evident from (11), ρ 1 depends on P FA1 and on the noise statistic, which is approximately known to the receiver.As stated earlier, the average spoofer SNR is not known and varies with varying spoofer-receiver separation due to pathloss, spoofer transmit power variations, and so forth.However, the average authentic line of sight (LOS) SNR is approximately known, given that the average LOS CNR of GNSS signals at the ground level is typically within [40-50] (dB-Hz), which maps into a postprocessing SNR of approximately [10][11][12][13][14][15][16][17][18][19][20] dB based on 1 ms of coherent integration.This a priori information can be used to determine a second threshold, ρ 2 , based on selecting a probability of false alarm associated with H2 as Given that the satellite geometry is not known to an acquiring receiver, it is reasonable to assume that SVs are approximately uniformly distributed in the sky.Consequently, the PDF of the average post processing SNR of the authentic GNSS signals, ρ (a) , can be approximated as where U(ρ L , ρ H ) denotes a uniform PDF and ρ L ≈ 10, ρ H ≈ 20 dB denotes the lower and the upper bounds of the uniform distribution.Consequently, (a) .( 14) ρ 2 can be numerically computed by inserting ( 14) into (12).Finally, the probability of detection associated with H2 can be computed as 3.1.Spoofer Detection Based on a Moving Antenna.As stated earlier, the typical usage mode of a handheld receiver includes incidental motion in the form of spatial translation, polarization rotation, and blocking of the receiver antenna.It is known that any temporal variation in the antenna response results in a temporal signal decorrelation in a multipath environment such that extra diversity branches can be made available for receiver processing [7][8][9][10].
To exploit the extra processing gain arising from antenna motion, the statistical properties of x need to be considered.Distribution of scatterers in many multipath environments such as indoors or urban areas approximately resembles a uniform sphere of scatterers [9,14].The correlation coefficient between signal samples, s = [s 1 , . . ., s n ], collected through spatially translating an antenna over an arbitrary trajectory in a Rayleigh fading environment that resembles a sphere of scatterers can be shown to be [7] [C s ] mn = ηsinc k 0 p mn , (16) where k 0 = 2π/λ is the propagation constant, p mn = |p m − p n | is the spatial separation between the antenna positions at which signal samples x m and x n are collected (see Figure 3), and η is the variance of s.Consequently, x m are statistically uncorrelated if the spatial separation between the antenna positions at which the samples are measured are greater than half a carrier wavelength (in GPS L1 frequency this maps into a spatial separation of 10 cm), resulting in the approximation of C s ≈ ηI N , where I N is an N × N identity matrix.
Rotation is another form of user interaction with a handheld receiver that results in variation in antenna's polarization.Variation in antenna's polarization is known to result in signal decorrelation.It can be shown that the covariance of signal samples measured through polarization rotation of a handheld antenna follows from [15] as where ψ mn is the angular separation of the polarization vectors at which signal samples x m and x n are collected (see Figure 3).Note that only three degrees of freedom are realizable based on a polarization rotation of a linearly polarized antenna [16].Therefore, N ≤ 3 uncorrelated signal samples are realizable based on a polarization rotation.A combination of polarization rotation and spatial translation can be utilized to further increase the number of diversity branches [9].The cross-covariance arising from a combined spatial-polarization translation of a GNSS handheld antenna can be shown to be [9] [C s ] mn = ηsinc k 0 p mn cos ψ mn . ( As can be seen from ( 18), the receiver motion in the form of a combined translation in space and rotation of polarization decorrelates the received signal and therefore can be utilized to synthesize several diversity branches useful for receiver processing.

Uncorrelated Rayleigh Fading Channels.
Assume that N uncorrelated signal samples are obtained based on a combined spatial and polarization translation of a GNSS handheld receiver in an uncorrelated Rayleigh fading channel such that C s = ηI N .Consequently, x = {x 1 , . . .x N } are jointly CN zero-mean RVs with x = s + w ∼ CN(0, C x ).C x = C s +C w denotes a covariance matrix of x with C w as the noise covariance matrix.To simplify the expressions to follow and without any loss of generality, the noise covariance is normalized such that C w = I N .Therefore, the SNR can be written as Consequently, the signal samples collected by a moving antenna in an uncorrelated Rayleigh fading channel are distributed according to x ∼ CN(0, (η + 1)I).It can be shown that is a sufficient statistics with respect to the hypotheses H0, H1, and H2 and as such is the detection variable [13].The thresholds (ρ 1 , ρ 2 ) can be found by determining the PDF of r and substituting in ( 10)-( 15) for any given P FA1 , P FA2 .Note that r is a measure of the received signal power.Therefore, the detection problem is based on comparing the received signal power, r, with a set of thresholds, (ρ 1 , ρ 2 ), to determine the authenticity of the received signal.For the spoofer to be effective, the spoofer signal power must be higher than that of the authentic signal in the target area such that the ML search in the CDS results in selecting the spoofer signal which has the largest correlation peaks.Therefore, r, which is a measure of the received signal power, can be utilized to discriminate the spoofer from the authentic signals.

Generalized Rician Channels.
In a generalized Rician channel, the channel gain, A(t), is a random variable distributed according to CN(μ, η/2) where μ = |μ| √ 2 exp( jα(t)) is the complex mean with α(t) denoting the phase of the complex mean and η/2 is the variance of the in-phase and the quadrature-phase Gaussian components of the channel gain.Consequently, x are jointly CN RVs and are distributed according to T is an N × 1 vector with α i denoting the phase, C x = C s + C w is a covariance matrix of x, and C w = I N is the normalized noise covariance.In a Rician channel, the average SNR, ρ, can be defined as ρ = 2|μ| 2 + η, and the magnitude of the mean, |μ|, and the variance, η, are related through the Rician K-factor, κ, such that κ = |μ| 2 /η.Since the angle of arrival (AoA) of the dominant signal component is not known to the receiver, μ cannot be estimated and therefore m and subsequently κ are unknown which makes it impossible to formulate a sufficient statistics based on a likelihood ratio test [13].Nevertheless, as will be shown here, the performance of the spoofer detection is approximately insensitive to the variation in the K-factor, κ, and to the cross-correlation of signal samples s as long as the cross-correlation remains moderately low, for example, <0.7.This is reasonable since diversity gain arising from combining equal-power diversity branches remains mostly unchanged for branch cross-correlations <0.7.Therefore, the suboptimal detection variable of (20) can be applied for spoofing detection in a generalized Rician channel with small performance degradation.Figure 4 shows ρ 2 being computed from (14) for various authentic and spoofer channel Kfactors (κ (a) and κ (sp) ) and based on ρ (a) = 15 (dB) and ρ (sp) = 20 (dB), two typical P FA2 = 0.01, 0.1, and for N = 1,3,5.As can be seen from Figure 4, smaller K-factors result in larger ρ 2 values.This is due to the increased uncertainty in the received signal power as the K-factor decreases.Nevertheless, the variation in ρ 2 is limited to a few dB and as such the K-factor does not play a major role in the optimization problem and may be ignored in the expense of slightly lower performance.Therefore, (20) can be applied to a generalized Rician channel as a suboptimal detector.In addition, as can be seen from this figure, a larger N results in a smaller ρ 2 for the same performance requirement of P FA2 .This is due to the diversity gain made available through the extra diversity branches for N > 1.
Figures 5 and 6 show the receiver operating characteristics (ROC) based on the detection variable of (20) and for ρ (a) = 15 dB, various N and ρ (sp) , and based on κ (a) = κ (sp) = 1 and κ (a) = 10, κ (sp) = 1.As can be seen in these figures, the detection performance improves with increasing the number of diversity branches, N. Also, larger ρ (sp) result in a better detection due to the further separation between the PDFs of the authentic and spoofing signals.When a stronger LOS signal component is present (κ (a) = 10 in Figure 6), a better detection performance is realized due to the reduced uncertainty in the authentic signal power.Note that setting P FA2 = 0 in (12) results in ρ 2 = ∞ and therefore P D2 = 0.This corresponds to a receiver not equipped with any spoofer detection.
To provide an alternative measure of performance improvement, the probability of error P e of (9) can be used.Figure 7 shows P e for various N, κ (a) , and ρ (a) = 15 and based on κ (sp) = 1 and ρ (sp) = 25 (dB).As can be seen from this figure, P e is approximately independent of the exact value of K-factor, which emphasizes the previous observations of Figure 4 where it was shown that the threshold ρ 2 is not very sensitive to the variations of the K-factor.P e decreases rapidly with increasing the number of diversity branches.The latter demonstrates the performance enhancement arising from the extra diversity branches made available through utilising a moving antenna. Figure 8 shows P e for ρ (a) = 15 (dB) and various ρ (sp) , κ (a) = κ (sp) = 1, and for N = 1, 3, 5, 10, 15, 20.Similarly, larger N and larger ρ (sp) , which provide a better separation between the authentic and the spoofing signal PDFs, result in a smaller P e .This is further demonstrated in Figure 9 where P e is plotted for N = 1-20, κ (a) = κ (sp) = 1, and ρ (sp) = 10-25 (dB).Note that, as the PDFs of the authentic and the spoofer signals become more alike, for example, ρ (a) = ρ (sp) , P e becomes larger.
As stated earlier, the spoofer signal power is affected by pathloss quantifiable by (5).As a result, the received spoofer SNR varies with the proximity to the spoofer transmitter.To provide an average measure of performance enhancement arising from utilizing the proposed technique, the average probability of error, P e , can be defined as P e = ρmax ρmin P e ρ (sp) dρ (sp) .(21) Figure 10 shows P e for various N based on κ (a) = κ (sp) = 1 and ρ (a) = 15 (dB), and for ρ min = 15, ρ max = 25 (dB).The effect of diversity is further emphasized in this figure where the average probability of error decreases by increasing N such that P e 0.19 for N = 20, implying that the proposed technique is very effective in reducing the spoofer effectiveness in the target area.

Conclusions
A multi-hypothesis detection problem was formulated based on a likelihood ratio test applicable to GNSS spoofing International Journal of Navigation and Observation detection.A straight forward spoofing detection technique based on signal power measurements was proposed and was shown to be effective for verifying the authenticity of the received GNSS signals in urban multipath environments, meaning that the spoofer signal power is abnormally higher than that of the authentic signal when the receiver is in the proximity of the standoff spoofer.
The proposed processing was further extended to exploit extra diversity branches made available based on a moving handheld receiver and was shown to further improve the spoofer detection performance.Unlike the previously proposed antispoofing techniques, the proposed technique does not require any hardware modification and can be readily applied to any handheld GNSS receiver with minimal firmware changes.It was shown that the proposed technique is largely insensitive to uncertainties in the statistical properties of the multipath channel as long as the collected signal samples are not strongly correlated.A suboptimal detector was proposed and effectively applied to a generalized Rician channel in which the channel parameters are not available to the receiver.An extensive statistical analysis was performed to assess the performance of the proposed technique.It was shown that the average probability of error can be reduced to less than 20% in a typical urban environment.

Figure 1 :
Figure 1: Hypothetical stand-off spoofer scenario in an urban canyon.The contours represent the random average spoofer signal power.The average authentic signal power is approximately constant over the entire area given that the receiver-satellite range is approximately unchanged and hence the pathloss.

Figure 2 :
Figure 2: A diagram of the PDFs of the detection variable "r" under H0, H1, and H2 hypotheses.

Figure 3 :
Figure 3: Spatial translation and polarization rotation of a GNSS handheld antenna.