A feasibility study of 5G positioning with current cellular network deployment

This research examines the feasibility of using synchronization signals broadcasted by currently deployed fifth generation (5G) cellular networks to determine the position of a static receiver. The main focus lies on the analysis of synchronization among the base stations of a real 5G network in Milan, Italy, as this has a major impact on the accuracy of localization based on time of arrival measurements. Understanding such properties, indeed, is fundamental to characterize the clock drifts and implement compensation strategies as well as to identify the direct communication beam. The paper shows how the clock errors, i.e., inaccurate synchronization, among 5G base stations exhibit a significant bias, which is detrimental for precise cellular positioning. By compensating the synchronization errors of devices’ clocks, we demonstrate that it is in principle possible to localize a static user with an accuracy of approximately 8–10 m in non-obstructed visibility conditions, for urban and rural scenarios, using the deployed 5G network operating at 3.68 GHz and relying on broadcast signals as defined by 5G Release 15 standard. This work has been funded by the European Space Agency (ESA) Navigation Innovation and Support Program (NAVISP) Element 2 pillar which aims at improving the competitiveness of the industry of the participating States in the global Positioning, Navigation and Timing (PNT) market.

to hardware limitation.Specifically, cellular positioning capabilities are strongly dependent on precise network synchronization.If the local clocks of cellular base stations are not accurately synchronized to a common network timing reference, a bias error will be present in measuring time-related parameters, thus impairing precise positioning in Time of Arrival (ToA) and Time Difference of Arrival (TDoA) methods.It is worth noting that this issue does not affect Round Trip Time (RTT) methods, in which range errors are mainly due to the propagation delays and clock drifts accumulated during the RTT.
The synchronization problem has been widely investigated in the literature [24][25][26][27] , but it still remains one of the biggest sources of error for precise positioning in cellular networks.For 5G communications, uplink and downlink transmissions between base stations and the UEs are scheduled in temporal slots, thus synchronization among the clocks is required to have a univocal reference of timing.Given the unfeasibility of providing atomic clocks at the UE side synchronized with Navigation Satellite System (GNSS), a Precision Time Protocol version 2 (PTPv2) 28 was created to meet the synchronization requirements.The PTPv2 protocol behaves according to a primary-and-secondary clock paradigm.The primary reference clock is deployed in the backhaul of the cellular network and provides timestamps to the secondary clock which, on the contrary, is placed at the TRPs (i.e., at the cellular base stations).Then, the secondary clock adjusts its local clock to keep precise time alignment with the primary clock.Although the PTPv2 protocol attracted significant interest for its capabilities of jamming resistance and indoor deployment, actual implementations could only achieve a clock synchronization that is accurate up to ±1.5 µs, as recommended by the International Telecommunication Union (ITU) 29 .Converting such time error into the distance, an error of up to ±450 m is introduced, confirming that the cellular network has never been conceived for accurate positioning purposes so far.Therefore, the clock synchronization issue emerges as a key limiting factor for the roll-out of precise 5G LCSs.A definite solution has yet to come, but a possible option to mitigate the problem could be including more precise oscillators such as the Temperature-Controlled Crystal Oscillators (TCXOs), with accuracy of approximately ±100 ns, or the Oven-Controlled Crystal Oscillators (OCXOs), with accuracy of ±50 ns 30 .Moreover, additional solutions based on the already available 5G Integrated Access and Backhaul (IAB) 31 could calculate the synchronization errors across TRPs, store them in a database and share them with the UE as assistance information.
Since 3GPP Release 16 positioning features are foreseen to be operational not earlier than late 2023, and PRSs are not available in the current release, we here assess the positioning capability of current 5G networks using the SSBs 32 .SSBs represent a special category of 5G signals that a TRP periodically broadcasts over spatial communication beams, enabling a UE to establish a cellular connection.The details of the initial access phase are here omitted; the reader can refer to 33 .
Recommended periodicity of SSBs is of 20 ms, although 3GPP standard foresees also other periodicities, i.e., 5 ms, 10 ms, 40 ms, 80 ms and 160 ms 34 .Different time-domain patterns for SSB transmission exist, differing according to the signal frequency f (distinguishing among frequencies lower than 3 GHz, between 3 GHz and 6 GHz, and higher than 6 GHz) and Sub-Carrier Spacing (SCS).The patterns are referred to as Case A, B, C, D, and E 35 , as detailed in Table 1.Each pattern is also characterized by the number of consecutive SSBs ( N SSB ) that are sent in a burst.Intuitively, a pattern specifies how many spatial directional beams a TRP uses to enable its discovery by the UEs. Figure 1 reports all the possible SSB patterns highlighting that a higher number of beams is used for high frequencies (e.g., millimeter waves).Higher values of N SSB indicate a TRP able to densely scan the spatial domain through more directional beams.
In this paper, we measure the ToA of SSB in real urban and rural scenarios with the final aim of assessing the impact of the synchronization error on 5G cellular positioning with the current 5G network deployment.This is done by knowing the periodicity of the SSB pattern and the exact position of the TRPs.The main scientific contribution relies upon a thorough study of the synchronization error (or, equivalently, range bias) and its impact on the positioning accuracy, analyzing urban and rural scenarios with TRPs in visibility and non-visibility conditions.We analyze the decorrelation period and the time drifts of the internal clocks at the base stations.Finally, we investigate the positioning capabilities of the real 5G system assessing both raw and synchronizationcompensated ranging measurements.
The paper is organized as follows: Section Results reports the main outcome of the measurement campaigns, i.e., the range bias distributions, the clock decorrelation period and the positioning performances.We further examine the reported analysis in Section Discussion, where we derive conclusions and suggest future research directions.Lastly, in Section Methods, we provide information about the data gathering and data filtering processes, together with the main adopted performance metrics and positioning algorithms.

Frequency
Frequency Frequency The analysis of the clock error is reported in Fig. 4 in terms of the histogram of the range bias, i.e., the synchronization error converted in meters, for each of the eight SSB beams of TRPs in visibility.The measured distribution reports a twofold behavior: for PCIs 252 and 411 the range bias distribution is mono-modal (note the presence of one clear peak), while for PCI 383 it presents a wider dispersion without a clearly defined shape.Considering that the visibility condition was well-respected, the mono-modal distribution for PCIs 252 and 411 is somehow understandable as the LOS distribution of the time of flight error is little affected by multipath components.However, we remark that the range error has a mean value (i.e., bias) of 73 m and 79 m for PCIs 411 and 252, respectively.As discussed later in Section Positioning performance, compensation of such bias is a key operation to unlock precise positioning services in 5G networks.Moving to the analysis of the TRP with PCI 383, we notice a completely different distribution.Given that the scanner simultaneously measured the 5G signals received from all TRPs, we attribute such an anomalous behavior to the TRP with PCI 383 itself.In this particular case, TRP with PCIs 383 adopted the PTPv2 standard which provides a high-quality frequency  www.nature.com/scientificreports/synchronization irrespective of the network load through physical layer (Ethernet interfaces).This protocol can provide a timing accuracy of at least ±100 ns (i.e., ±30 m), which can be improved with more precise local oscillators.This aspect is further investigated in Section Clock offset coherence time.
The analysis for non-visibility conditions is reported in Fig. 5. Compared to the LOS case in Fig. 4, the major difference is in the shape of the distribution, which manifests different peaks mostly likely generated by the multipath phenomenon.Due to the NLOS, in fact, the observed range error is generated not only by the clock error, but also by the multipath contributions.Moreover, wider supports, large mean errors (the overall mean for non-visibility TRPs is of 150 m), and absence of received signal is also experienced.Due to the higher emitted power of TRPs, i.e., around 10 − 20 W, the LOS contributions are likely to be present, but hard to isolate.A joint tracking over time of the range biases that accounts for the different visibility conditions is advisable 36 .
Clock offset coherence time.The second analysis focuses on the temporal variability of the synchronization error, with the objective of quantifying the coherence time defined as the time over which the clock error auto-correlation coefficient drops below a predefined threshold.This is extremely useful to understand how often information on clock correction should be provided to the UEs (in case of UE-based localization).The normalized autocorrelation function 37 has been computed and compared to the threshold obtained using two standard errors of the sample autocorrelation as the confidence bounds.We refer to the Section Methods for more details.With this purpose, we consider the PCI 252 (which is in visibility) and plot in Fig. 6 the (normalized) autocorrelation function of the range error for each SSB beam.From the plots in the figure, we first notice that the range error of SSB 3 is more correlated over time (i.e., more stable over time), compared to the other SSBs.This behavior was expected as the beam for SSB 3 is directed towards the location of the scanner, i.e., along the geometrical LOS.On the contrary, the range error extracted from other beams have shorter coherence time (e.g., the autocorrelation is even impulsive for SSB 7), demonstrating the uncorrelation property of a range bias which is frequently subject to multipath phenomenon.Multipath main effects are random spikes in the range bias, which increase the signal decorrelation.It follows that the analysis on the autocorrelation properties can be used to detect the LOS beam.This is verified by the analysis on the received power in Section Methods where we clearly identify SSB 3 as the beam along the direct path.As concerns the coherence time, by setting a threshold of 0.2 to the autocorrelation function, we get a coherence time of 30 − 35 seconds.This value is indicative to suggest an update of clock synchronization with a period not larger than 30 − 35 seconds.Note that after 75 seconds, even if in visibility, the range bias becomes completely uncorrelated since the autocorrelation drops below the confidence bounds ( ±2σ ), which are computed considering the range bias as a Gaussian white noise process with a standard deviation σ of about 1/ √ N , with N the number of range bias samples.To complete the analysis, in Table 2, we show the coherence time, indicated with T COH , and the associated mean power for each PCI (i.e., 252, 411 and 383) and SSB.Note that PCI 411 holds a similar behavior with respect to PCI 252 and that the SSB with higher coherence time corresponds to the higher received mean power, verifying the detection of LOS beams through the autocorrelation property.
On the contrary, while the TRPs with PCI 252 and 411 hold a much narrow range bias distribution, the TRP with PCI 383 has a much higher drift in time.To highlight this aspect, in Fig. 7 we report the range biases, and associated received power, over the measuring time.To identify the LOS beam, we also report the corresponding power per SSB together with the associated range bias values in the radial axes.Considering the LOS beam of PCI 383, i.e., SSB 4, we measure a rate of change in range bias over time, i.e., clock drift, of 0.3 m/s which causes the wide support in Fig. 4.This is compatible with the OCXOs which are typically adopted in cellular TRPs 30 and hold a drift between 0.1 and 10 part per billion (ppb), equivalent to 0.03 m/s and 3 m/s.In order to discuss the idea that the drift was caused by the near site distance, we point out that the drift was observed in all the three PCIs belonging to the site in discussion, i.e., PCI 382, 383 and 402, and in all five reception points.Given the severe observed offset, compensating for these types of errors becomes of paramount importance.Moreover, while stable range biases, as the ones belonging to PCI 252 and 411, can be compensated with an a-priori clock-offset estimation, the drift has to be constantly tracked and communicated as assistance data to the UEs.
Positioning performance.This last assessment has the objective of verifying the impact of synchronization errors on the positioning performance of current 5G network deployments.In particular, we perform two experiments to verify the effect of estimating and compensating the range bias in UE localization.
We consider two different positioning scenarios, an urban area in Milan (Fig. 2) in the Leonardo Campus of the Politecnico di Milano and a rural area in the suburbs of Milan (see Fig. 8).In the former, the multipath is expected to be relevant and play a crucial role, often creating multi-modal distributions of range errors.In the latter, on the other hand, LOS condition is mostly available, but the larger distance between the TRPs and UE might cause a reduction of received power.For the evaluation of 5G positioning in the urban scenario, we considered three TRPs in LOS, i.e., PCIs 383, 411 and 252, while for the rural scenario we selected four TRPs in LOS, i.e., PCIs 733, 427, 613 and 750.Regarding the synchronization technology, PCI 733 adopts the Precision Time (PTPv2) protocol, while PCIs 427, 613 and 750 use SFP via GNSS.The selection of LOS TRPs is based on both the analysis of clock offset distributions, as described in Section Clock offset statistics, and from geographical positions with respect to the UE.
Three different localization methods are compared: a snapshot (i.e., without tracking over time) Non-linear Least Squares (NLS) algorithm without range bias correction (blue scatter points); an NLS with range bias correction (green scatter points); a tracking over time by an Extended Kalman Filter (EKF) with range bias correction (purple scatter points).For each method, the 95% confidence ellipse (95% of estimates fall inside the ellipse) is reported.The CramérRao Bound (CRB) ellipse (in yellow), considering the standard deviation on the range bias and achievable accuracy with the given bandwidth, is shown as reference.Lastly, the blue icon indicates the UE ground truth position; while the red icons identify the TRPs.The mean of the position estimate and the error bias are denoted with a black star marker and solid black lines, respectively.For the assessment, we considered the following performance metrics: Mean Absolute Error (MAE), Distance Root Mean Square (DRMS), and error ellipse with 95% confidence.Moreover, for each method, we plot the Cumulative Density Function (CDF) of UE positioning error, while a complete list of results, i.e., MAE, Circular Error Probable (CEP) 95, bias norm and DRMS, for both urban and rural scenarios, is reported in Table 3.
The results of static UE positioning for the urban scenario are presented in Fig. 9. Comparing the NLS approaches, we remark the beneficial effect of range bias correction, which reduces the MAE of UE localization from 111 m to 10 m.Analyzing the variance of UE location estimates (removing the bias) we can use the CRB as benchmark.Specifically, the NLS with range bias correction has a CEP 95 of 16.4 m which is higher than the CRB of 10.3 m.Adopting a tracking filter for UE positioning and by taking into account the correlation of the UE Moving to the rural scenario, we have an increased number of TRPs in visibility but with a much higher distance compared to the urban area (see the scale in the bottom-right corner of Figs. 2 and 8).The results of the analyses are reported in Fig. 10, where the same methods adopted in the urban area are used.The UE localization with NLS without range bias correction achieves an MAE of 141 m, i.e., 30 m bigger than the urban case.The range correction and use of EKF lead to increased positioning performance, with MAEs of 24.6 m and 10.7 m, respectively.
Discussion.According to the analyses and results reported in Section Results, we can conclude that synchronization constitutes a major issue for precise cellular positioning.Considering the current network deployments characterized by 3GPP Release 15 standard, relative-narrow bandwidth of SSBs (7.2 MHz) and high clock drifts between TRPs, a positioning accuracy of 8.2 m and 10.7 m has been achieved in urban and rural environments, respectively, upon compensating for synchronization errors.We remark that this is done without any alteration of the network and gathering only passive signals.To achieve such performance, a tracking filter together with almost perfect correction of synchronization errors has been used.
While the deployment of future 5G releases should improve cellular coverage and use larger signal bandwidth (e.g., by the PRSs in 3GPP Release 16), network synchronization will still remain a great challenge.In contexts where the relative time error between nearby base stations is the most important factor, novel synchronization algorithms are being defined.Promising solutions are the over-the-air synchronization (OAS) 38 or even the optical fiber synchronization 39 .As explained in Section Methods, we adopted a single 5G scanner for the experiments, gathering separately, i.e., at different time, the range bias statistics and the measurements adopted for positioning.This led to an unavoidable positioning error quantified with the standard deviation of the range-biases distributions ( 13 − 15 m in our case).In order to obtain cm-level accuracy, as required for challenging applications (e.g., cellular positioning for automated driving 40,41 ), range bias should be measured and shared with the UEs.
Future works could extend our study in this direction, exploiting not only higher signal bandwidths, but also precise beam management and joint use of multi-constellation GNSS and 5G observations.Moreover, advanced synchronization techniques to track in real-time the range corrections are also of interest, as well as the use of filtering techniques to jointly estimate the UE position and clock offsets of TRPs.

Methods
To study the effect of synchronization in the currently deployed 5G networks, we conducted a number of experiments with a R&S TSMA6 scanner 42 , which is capable of simultaneously measuring signals between 350 MHz and 6 GHz.The TSMA6 scanner makes use of an internal pulse per second (PPS) signal generator to determine the time of arrival of a 5G signal.The scanner requires a configuration procedure to select the frequency bands, load modules (e.g., 5G scanner and GNSS), and acquire GNSS signal.When it is fully configured, it is able to simultaneously measure multiple 5G signals (from different TRPs), outputting a measurement every 20 ms (the sampling frequency of the instrument is 30.72 MHz).Among them, the ones of interest for our study are the www.nature.com/scientificreports/SSBs.Specifically, the scanner outputs the ToA of each detected SSB: this information is available under the technical name of "ToA PPS" 42 since it uses the GNSS PPS for time reference.The "ToA PPS" refers to the time of arrival of the PSS within the SSB.We shall notice that a post-processing is performed by the scanner in order to retrieve the final Time of Arrival (ToA) value.The scanner's internal processing involves the computation of the cross-correlation of the received signals, which is used to determine the top-N multipath components (i.e., the N components with the highest cross-correlation value).The decision value for this ranking is set to the average Signal-to-Interference-plus-Noise Ratio (SINR), which provides a balance between signal strength and quality.Among these top-N multipath components, the one with the lowest "ToA PPS" is selected for ToA estimation.This approach ensures robustness in mitigating the ambiguities due to early arrivals.It is important to note that the R&S scanner does not permit full access to physical level parameters, as it automatically performs TOA estimation processing.Therefore, the analyses we carried out are based on the processed output provided by the scanner.The datasheet of the R&S scanner specifies a granularity of 5 ns with good GNSS visibility.This level of precision allows for accurate measurements of the ToA of 5G signals, which is crucial for our study.For measurement uncertainties, the scanner reports the chosen average SINR, which can be used as quality metrics of the measurements, and also the Synchronization Signal (SS)-Received Signal Received Power (RSRP) 7 , which we referred to as received power.
In Fig. 11, we show a concise representation of beams, i.e., SSBs, sent by the TRPs.Since the SSBs are transmitted over different spatial directions, the scanner acquires signals coming from both LOS and NLOS beams.Knowing the 5G frame structure, the employed numerology, and SSB transmission pattern and measuring the ToA of each SSB (if any), we estimate the distance between the transmitting TRP and receiving UE (i.e., the scanner).The measured delay of the received signal is composed of the Time of Flight (ToF) of the propagation, multipath effects, the additional contribution of desynchronization, and other sources of unknown errors inside the hardware.
We define with t (Start) tx the time of transmission of the beginning of the 5G frame at the TRP, with t tx,i and t the instants of transmission of the SSB k at TRP i and its reception at the scanner, respectively.The clock offsets of the receiver (i.e., the scanner) and of TRP i are indicated with t UE and t i , respectively.Defining with p TRP i and x the position of the TRP i and of the scanner, respectively, our goal is to estimate the true ToF τ from TRP i and SSB k to the scanner.The estimate is defined as: where t mp and t n are the range errors due to multipath and noise at the radio receiver, respectively.Given that our experiments were conducted in controlled (to the best of our possibilities) LOS conditions, where multipath effects are minimized and LOS is mainly guaranteed (at least for one beam), we assume the ToA measurements gathered by the R&S equipment refer to the direct path.As concerns the receiver noise, we observe that the standard deviation of the ranging error due to this source of error is nearly an order of magnitude smaller than the standard deviation of synchronization error.For this reason, we focus on the range biases due to synchronization error only, neglecting the characterization of error contributions t mp and t n , leading to the approximation in (1).The ToA measured by the receiving scanner for TRP i and SSB k is:  for SSB 0 with numerology 1, f = 3.68 GHz and SSB Case C (see Fig. 1) is about 2 OFDM symbol intervals, i.e., 71.875 µ s.We can then transform (1) into: This is possible if both the R&S and the TRPs have the common reference GNSS time and, therefore, they are able to nominally compute the same start of the frame t (Start) tx (except from the local clock errors).After having measured the ToF with (3), by knowing the exact distance between the transmitting TRPs and the receiving UE it is possible to compute the synchronization error e  43 , by placing the RTK receiver right below the 5G antenna of the R&S scanner.Since RTK technology guarantees an accuracy of few centimeters 44 , considering the final error on 5G positioning in the order of several meters, we can safely use the RTK output as the exact UE coordinate.On the other hand, the coordinates of TRPs were measured by merging RTK surveys of the site boundaries (we remark that accessing the 5G sites and reaching the heights of the antennas is extremely challenging) with architectural drawings of the sites, resulting into highly-reliable information in the order of decimeters and still much lower than the 5G positioning error.We highlight that the clock difference t i − t UE cannot be separated into the two offset components as they are related to perfectly synchronized time-axes which cannot be measured.Therefore, we can just measure the difference in offset between two asynchronous clocks.
The ranging error will only consist of the offset between the clock of the scanner and the clock of the TRPs, under the assumption of LOS propagation.To guarantee such assumption, we need to understand which beam is in LOS, i.e., which beam is pointing towards the scanner.To this aim, knowing the azimuth and elevation of the panel array of the TRP, in Fig. 12 we report the associated power per SSB, i.e., per beam, together with the associated range bias values in the radial axes.By knowing the technical implementation features of the Vodafone 5G TRPs, we are able to state that the represented TRP in LOS (PCI 252) covers a spatial sector of 120 degrees in azimuth with 6 high-beams of roughly 20 deg beamwidth each and 2 low-beams of roughly 60 deg beamwidth each.The high-beams are referred to SSB numbers 0, 1, 2, 3, 4 and 5, while the low-beams with SSB numbers 6 and 7.For this specific TRP with PCI 252 in visibility condition, the beam with higher received power is the SSB number 3, with an average received power of about −90 dBm.On the contrary, the received power for the other beams is as low as −120 dBm.A peculiar fact is that, analyzing the spreads of the range bias samples, they appear to hold the same support, i.e., about 60 m as also shown in Fig. 4.This is due to the fact that the scanner is able to receive even very attenuated signals and that the low-power received signals are obtained from a direct path.Indeed, the radiation pattern of the antenna panel, even if highly directional, does not retain a null gain in the undesired directions of propagation.This aspect can be non-optimal for communication aspects (where   high received power enables higher throughput), but it is still useful for positioning (where the key information is about the shortest path).For the study in Section Clock offset coherence time, we computed the normalized autocorrelation function of the clock error time series y t as 37 : where r k is the autocorrelation at lag k, ȳ is the mean of the time series and T is the total number of observations.The coherence time is defined as the lag over which the autocorrelation drops below a given threshold.The threshold is obtained using two standard errors of the sample autocorrelation in the confidence bounds, where the formula for the standard error using Bartlett's approximation is: where q is the lag beyond which the autocorrelation function is 0.
Given the availability of a single 5G scanner, we carried out the range bias estimation and measurement collection for UE positioning on two separate datasets.In the urban scenario, we considered four reception points for range error computation at the beginning of the experiment and one reception point at the end for testing.On the contrary, in the rural area, we did the opposite procedure.Both the testing and measurements phases last about 250 s, while the intervals between each measurement point in both scenarios were approximately of 10 minutes.Given the unfeasibility of tracking in real-time the range biases, we used the median values of each range bias distribution as a-priori compensation, which were then averaged across the measurement points.
We should also note that, for real applications, the computed coherence time should serve as a reference time interval for estimating the clock offsets.In other words, range biases should be estimated and compensated at least once within the coherence time, or even more such as by continuously tracking.
UE positioning is performed every T seconds by gathering all ToF measurements from all the measured SSBs within the temporal window.Denoting with τ i = { τ i (ℓ)} L i ℓ=1 the set of ToF measurements from TRP measurements from TRP i at position p TRP i in the interval T, we compute the median value over the measured ToF for each TRP such that the aggregated vector ρ of all ToF measurements is created as: where med(•) is the median operator and N is the number of detected TRPs.Then, the UE bidi- mensional position estimate x , is computed according to the NLS algorithm with Jacobian matrix [H] i−th row = [H(x)] i−th row = ∂h i (x) ∂x , where h i (x) = �p TRP i − x� 2 is the ToF measurement model.The NLS algorithm is implemented with the iterative search of the estimate according to the Gauss-Newton methodology 45 or EKF 46 with random walk motion model 47 .
The 95% confidence ellipse are obtained from the error covariance matrix C defined as: . The 95% confidence level is computed considering an horizontal accuracy of 2σ H = 2 √ Tr (C) , where Tr (•) denotes the trace of the argument matrix.Finally, the MAE can be extracted as In order to have a lower bound on the positioning error, we computed the CRB by deriving the Fisher information matrix of the range-based localization problem 48 , which accounts for both the TOF measurement error and the geometrical condition of TRPs and UE.For CRB computation, we assumed the TOF measurements as Gaussian, unbiased and uncorrelated, and each link as single-path LOS.It follows that the standard deviation of the ranging measurement of TRP i is: where β is the effective bandwidth (here assumed to be equal to 7.2 MHz), and SNR i E s,i /N 0 is the signal-to- noise ratio which is directly computed by the R&S scanner.E s,i is the average received energy from TRP i, while N 0 /2 represents the spectral density of the additive white Gaussian noise.
The final CRB covariance matrix C CRB of the position estimate is obtained as: where G = (H T H) −1 is the Geometric Dilution of Precision (GDoP) 49 of the positioning problem and σ τ = E i σ 2 τ ,i .

Figure 2 .
Figure 2. Representation of the urban scenario selected for the experiments and composed of three 5G sites (red squares), each with 3 TRPs (red icons) with indicated PCI identifiers, and 5 UE measurement points (blue icons).The power map is computed with MATLAB ray-tracing software at ground level.

Figure 3 .
Figure 3. Three sites of the urban scenario in Fig.2, with related TRPs.For each TRP, the color of PCI number indicates the visibility condition: green color is for LOS, red color is for NLOS.

Figure 4 .Figure 5 .
Figure 4. Distribution of the synchronization error in meters, i.e., range bias, for the TRPs in visibility (identified by their PCIs) for each of the 8 SSB beams.

Figure 6 .
Figure 6.Normalized autocorrelation function of the range biases of TRP with PCI 252 for each SSBs.

Figure 7 .
Figure 7. Evolution of range bias (in m) over time for PCI 252, 411 and 383 with corresponding 2D representation of received power per beam.Colors indicate the value of measured received power in dBm.

Figure 8 .
Figure 8. Simulated representation of the rural scenario composed of 12 TRPs (red icons) identified by the indicated PCI number and a static measurement point (blue icon).The power map is computed with MATLAB ray-tracing software at ground level.

Figure 9 .Figure 10 .
Figure 9. Urban scenario: (a) scatterplot of the position estimates in case of range correction (green squares), without range correction (blue squares) and EKF.The ellipses are obtained with 95% of confidence level.Black squares and solid black lines represent the mean position estimates and the error bias, respectively, while the CRB is highlighted with orange line.(b) CDF of the UE positioning error.
thus it is possible to compute the ToF by subtracting the time interval between the start of the 5G frame t (Start) tx and the expected transmission of the SSB t (k) tx,i to the measured arrival time ToA (k) i .As an example, the interval (1)

Figure 11 .
Figure 11.SSB beam sweeping along different directions.The receiver, i.e., R&S TSMA6 scanner, gathers both direct and reflected signals, experimenting a multipath delay.t i and t UE are the clock offset of the TRP i and UE, respectively.On the contrary, the distance between the TRP i and UE are indicated with d i−UE .
between TRP i and UE as: where d i−UE = �p TRP i − x� 2 is the distance between TRP i and the receiving scanner, while c is the propaga- tion speed.The positions of the measurement points have been measured with Real-Time Kinematic (RTK) technology, and positioning services provided by SPIN3 GNSS

7 (
a) 3D representation of received power per beam/SSB.
2D representation of received power per beam/SSB.

Figure 12 .
Figure 12.The radial axes represents the range bias in meters while the color indicates the received power for each measure.(a) and (b) show the beams in 3D and 2D, respectively.

Table 2 .
Clock offset coherence time T COH (Table2a) and mean received power (Table2b) per PCI for each SSB.The SSBs with higher coherence time are highlighted in bold.