LEO-PNT Performance Metrics: An Extensive Comparison Between Different Constellations Kaan Γ‡elikbilek1,* , Elena Simona Lohan1 1 Tampere University, Korkeakoulunkatu 7, 33720, Tampere, Finland Abstract Recent years have shown that Low-Earth Orbit (LEO) satellites are becoming the leading idea for the future of the space industry, gathering heavy investment from technology giants, as seen from several providers of mega-constellations, such as SpaceX (Starlink), Eutelsat (OneWeb), Iridium, or Amazon (Kuiper). LEO satellites are suitable not only for communication purposes, but they also hold a strong potential for Positioning, Navigation and Timing (PNT) applications, as their proximity to Earth results in fast satellite movements in orbit as well as in high received signal strengths, which may translate into better PNT signals compared to already available alternatives. In this work, we show the viability of LEO-PNT constellations by providing a comprehensive performance comparison based on coverage, Dilution of Precision (DOP) and Carrier-to-Noise Ratio (𝐢/𝑁0 ) metrics between several existing and upcoming constellations, as well as two theoretical LEO-PNT constellations, by considering them as dedicated LEO-PNT systems. Our results show that, among the existing and upcoming LEO constellations, Starlink, OneWeb, Xona and Centispace show great promise for future PNT solutions, and that alternative designs that are on par, or perhaps even better, are still possible. Keywords LEO Constellations, GNSS Positioning, LEO-PNT, Constellation Comparison, Link Budget Simulation 1. Introduction 1.1. Motivation for LEO-PNT Systems Technological advances in satellite technologies have opened up the possibility of cheaper and bulk satellite production, i.e., in the form of CubeSats [1, 2, 3, 4], which makes the concept of Low-Earth Orbit (LEO)-based applications a highly promising and exciting new area of business and research. From academia to industry, many players are interested in designing, launching, or innovating the existing and upcoming LEO satellites, with promises to transform our everyday habits. Possible applications of LEO systems are to offer broadband connectivity globally, to enable leading scientific research in fields within the space industry, to offer Earth-sensing solutions, and, more recently, to complement the existing Global Navigation Satellite Systems (GNSS) for robust and seamless navigation. There are already thousands of satellites launched in LEO orbits placed between 200 and 2000 km above the Earth [5, 6] and these satellites support various communication [7, 8, 9, 10], remote sensing [11] and Earth sensing applications [12, 13]. It is expected that the number of LEO satellites will continue to increase at a fast pace in the near future due to the fast time-to-market and lower launching costs compared to Medium Earth Orbit (MEO) and Geosynchronous Equatorial Orbit (GEO) satellite launches. The suitability of LEO satellites as Positioning, Navigation and Timing (PNT) solutions however, has not yet been addressed extensively in the current literature and it has started to gain momentum only in the past few years. Currently, global positioning solutions rely heavily on GNSS, which are susceptible to malicious interferences. Therefore, there is a timely need of complementary, or even alternative, global, and preferably low-cost, positioning methods. The Low-Earth Orbit-based Positioning, Navigation, and Timing (LEO-PNT) concept is the perfect candidate due to LEO satellite’s proximity to Earth, which translates to stronger signals compared to MEO and GEO signals at similar carrier frequencies, their WIPHAL 2024: Work-in-Progress in Hardware and Software for Location Computation, June 25-27, 2024, Antwerp, Belgium * Corresponding author. $ kaan.celikbilek@tuni.fi (K. Γ‡elikbilek); elena-simona.lohan@tuni.fi (E. S. Lohan)  0000-0001-5170-8656 (K. Γ‡elikbilek); 0000-0003-1718-6924 (E. S. Lohan) Β© 2024 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). CEUR ceur-ws.org Workshop ISSN 1613-0073 Proceedings capability for global coverage, and their potential for flexible, scenario-specific designs, e.g. through optimization methods [14]. These new LEO signals could be exploited for PNT in the inevitable event that GNSS signals become either unavailable (e.g., in deep urban canyons, under dense foliage, during long periods of interferences) or untrustworthy (e.g., under malicious spoofing attacks). 1.2. Current LEO Landscape The current LEO landscape includes new constellations such as Centispace [15, 16] and Xona [17], as well as older constellations such as Orbcomm [18], Iridium [19], Globalstar [20], Starlink [8], or Kuiper [10]. These constellations have potential for PNT services, either as signals of opportunity (e.g., Iridium, Starlink) or as dedicated PNT systems (e.g., Centispace, Xona). Although many aspects of these constellations are not public knowledge, either due to legal concerns or due to design uncertainties, research is being conducted for the different parts of LEO satellite systems in a speedy manner. For LEO-PNT applications, research work have mostly been focused on integrated LEO and GNSS solutions [21, 22], as well as on meta-signals and opportunistic signal frameworks [23, 24] and on alternative positioning based on Doppler integration [25, 26]. In [27], the authors combine these three focuses and argue that the unknown nature of the LEO satellite signals –due to private operators tendency to not share technical information– make the opportunistic approach a necessity and present their signal model and estimation procedure for multiple scenarios involving multiple-constellation PNT using OneWeb, Iridium NEXT, Starlink, and Orbcomm constellations. Their results show the feasibility of LEO and GNSS integration, as well as Doppler-based positioning implementations together with pseudorange based methods. In addition to PNT solutions, research work on LEO satellites expand to other fields as well, such as generic constellation designs for multi-purpose applications [14, 28], LEO-based applications for autonomous vehicles [29, 30], and LEO network design [31, 32]. The variety of research in possible applications of LEO satellites further shows the potential benefits of LEO-PNT solutions in the relatively near future. 1.3. Paper Goal and Contributions Motivated by the fact that very few performance comparisons among LEO-PNT constellations have been published so far, we present a simulation-based extensive comparison of several LEO-PNT performance metrics for eight LEO satellite constellations (one as the GNSS benchmark, three selected among existing mega-constellations, two selected among on-going LEO-PNT designs and two experimental ones, based on prior work by the Authors). The comparisons are made with a MATLAB-based [33] in-house developed constellation simulator (see section 2.3), under several indoor/outdoor scenarios. Our main contributions are: β€’ Providing an extensive performance comparison, based on coverage, Dilution of Precision (DOP), and Carrier-to-Noise Ratio (𝐢/𝑁0 ) metrics between nine constellations, based on three types of models: models relying on existing mega-constellations (Kuiper, OneWeb, Starlink), models relying on existing smaller-sized LEO-PNT constellations (Centispace and Xona) and Authors’ derived models, based on direct parameter optimization [14, 4]. The European GNSS constellation Galileo is included as a benchmark. β€’ Discussing the meaning of the obtained results under indoor and outdoor scenarios and empha- sizing the open challenges in designing a LEO-PNT system. β€’ Providing system recommendations on constellation design aspects of possible future LEO-PNT constellations. Table 1 Parameters of Included Satellite Constellations (as of April 2024) Constellation # of Satellites Altitude [km] Inclination [deg] EIRP Bandwidth Carrier Frequency * Kuiper[10] 7774 590-650 33-80 76 dBm 400 MHz 18 GHz * Starlink Gen-1[7] 4408 540-570 53-97.6 69.1 dBm 250 MHz 12 GHz * Starlink Gen-2[8] 29988 340-615 33-148 69.1 dBm 250 MHz 12 GHz * OneWeb[9] 7808 1200 40-88 65 dBm 250 MHz 12 GHz *CentiSpace[34] 190 975-1100 55-88 ** 65 dBm ** 4.092 MHz ** 1561.098 MHz * Xona[35] 492 ** 1080 53-90 ** 65 dBm ** 4.092 MHz ** 1561.098 MHz Experimental 1 [14] 280 825 72 69.1 dBm 10 MHz 12 GHz Experimental 2 [4] 420 600 76 69.1 dBm 10 MHz 12 GHz Galileo[36] 27 23222 56 59 dBm 24.552 MHz 1575.42 MHz * : These constellations have multiple-shells, and the indicated altitude and inclination ranges mean that individual shells have particular values for that parameter which falls between the given ranges. ** : These parameters are based on our assumptions as exact values could not be found in public resources. 2. Methodology and Target Performance Metrics 2.1. Relevant Constellations As mentioned in subsection 1.3, our study uses the Galileo constellation as the GNSS benchmark, and includes 7 known LEO constellations: 3 mega constellations, 2 smaller-scale constellations, and 2 experimental ones. Even if possibly outdated, the constellation parameters (i.e., orbital altitude, number of satellites, number of orbital planes, phasing angle between orbital planes, inclination angle of the orbital plane, constellation topology, Right Ascension of the Ascending Node (RAAN) and eccentricity) for the existing constellations can be obtained from public sources, with some exceptions. Most of the relevant parameters for Xona were published in their patent application [35], however a few parameters, i.e. the altitude and satellite’s exact operating Effective Isotropic Radiated Power (EIRP), are not given exactly. In a similar manner, to the best of the Authors’ knowledge, the exact channel parameters for Centispace are also not available publicly. Therefore, some assumptions were made as seen in Table 1 in order to fill the gaps. Since Centispace has been designed as an augmentation system for Beidou, we assume that the missing parameters are similar to Beidou’s B1 band. As for Xona, the orbital plane altitude is taken as the mentioned upper limit from [35], and the channel parameters are again assumed to be similar to Beidou’s for comparison purposes with Centispace’s constellation design. In addition to the mentioned known constellations, we provide 2 experimental single-shell LEO-PNT constellation designs, that have been obtained in our earlier studies [4, 14] that we name; i) "Experimental 1", and ii) "Experimental 2". The important parameters for all the relevant constellations that we selected for our comparisons are seen in Table 1. The Starlink constellation includes two generations: Gen-1 refers to the satellites launched according to the initial constellation design from 2018 and Gen-2 refers to the satellites launched during 2020-2022 with an alternate design. 2.2. LEO-PNT Metrics In order to evaluate the performance of LEO-PNT constellations we selected metrics that are related to the geometry between the users and the satellites, as well as to the reception quality of the received signal. The metrics that we selected for our comparisons are: coverage, Geometric Dilution of Precision (GDOP), Position (3D) Dilution of Precision (PDOP) and 𝐢/𝑁0 . The coverage reflects the signal-reception percentage; for most PNT solutions, a good reception requires having at least 4 satellites in view, i.e., 4-fold coverage. However, not every constellation we use in our comparison is PNT focused, and may not be optimized for 4-fold coverage. Thus, we compute both 4-fold and 1-fold coverage: the coverage here is computed as the percentage of the number of users that have at least 4 and 1 satellites in view, respectively. The GDOP and PDOP are DOP metrics that reflect the geometry of the user and the constellation [37]; the PDOP reflects the 3D position accuracy and the GDOP reflects the joint 3D position-and-timing accuracy. Similarly, the 𝐢/𝑁0 reflects the quality of the received signal, and it is calculated via eq. (1), where 𝐡𝑀 is the bandwidth of the channel, and Signal-to-Noise Ratio (SNR) is calculated for each satellite-user pair. 𝐢/𝑁0 π‘‘π΅βˆ’π»π‘§ = 𝑆𝑁 𝑅 + 10π‘™π‘œπ‘”10 (𝐡𝑀 ) (1) A high PNT performance with respect to these metrics is reached for high 4-fold coverage, low values for DOP metrics, and high values for 𝐢/𝑁0 . 2.3. Simulation Environment In our previous works [4, 14], a detailed LEO constellation simulator has been developed in MATLAB for LEO-PNT performance analysis, combining MATLAB libraries with the external QuaDRiGa channel library [38], used for the link-budget modelling. Our simulator mimics a satellite constellation from a set of inputs: the constellation parameters, the start time and the duration of the simulation, and the user information (i.e., position and velocity vectors at each time instant and number of users placed according to a uniform distribution on Earth). In addition, a secondary set of input parameters are provided for QuaDRiGa models, which includes: satellite EIRP, receiver sensitivity, atmospheric attenuation effects, and scenario information. We performed simulations for 2 different receiver sensitivity values; i) -125 dBm, representing the low-sensitivity case, and ii) -185 dBm, representing the high-sensitivity case, under 6 different scenarios provided by the QuaDRiGa library that correspond to different Line of Sight (LOS) and Non-Line of Sight (NLOS) conditions: β€’ "Indoor 1 (I-1)"; Indoor, Rural and NLOS β€’ "Indoor 2 (I-2)"; Indoor, Urban and NLOS β€’ "Outdoor 1 (O-1)"; Outdoor, Urban and NLOS β€’ "Outdoor 2 (O-2)"; Outdoor, Urban and LOS β€’ "Outdoor 3 (O-3)"; Outdoor, Rural and LOS β€’ "Outdoor 4 (O-4)"; Outdoor, Rural and NLOS The simulations consider users uniformly spread on Earth and stationary, and rural/urban choice changes the number of channel clusters and paths, as well as parameters related to large-scale fading decorrelation distances and inter-parameter correlations within QuaDRiGa. Indoor scenarios are considered with 50 meter penetration. Constellations are initialized according to their own parameters as in Table 1. Each simulation has a duration of 1 hour with 1 minute samples, meaning 600 Monte-Carlo runs per satellite in the constellation. Scenarios assume summer conditions for the atmospheric models. 3 attenuation models are taken into consideration in the link budget; i) atmospheric absorption, ii) rain, and iii) fog, all calculated via MATLAB’s internal functions. We assume light rain with 2.5 mm/h rate and a cloud liquid water density of 0.5 g/m3. Temperature (𝑇 ), dry air pressure (𝑝) and water-vapor density (𝜌) are modeled from the satellite altitude (β„Ž), as given in equations (2), (3) and (4) respectively [39]. 𝑇 (β„Ž) = 286.8374 βˆ’ 4.7805β„Ž βˆ’ 0.1402β„Ž2 [𝐾] (2) 𝑝(β„Ž) = 1008.0278 βˆ’ 113.2494β„Ž + 3.9408β„Ž2 [β„Žπ‘ƒ π‘Ž] (3) 2 𝜌(β„Ž) = 8.988 exp(βˆ’03614β„Ž βˆ’ 0.005402β„Ž βˆ’ 0.001955β„Ž ) 3 [𝑔/π‘š3] (4) 3. Comparative Results and Discussion Tables 2 and 3 show the results of the simulations for high and low sensitivity cases respectively, providing LEO-PNT metrics averaged over the Monte-Carlo runs. The immediate thing to notice between the two tables is the coverage difference; as expected, the coverage is significantly higher in high-sensitivity than in low-sensitivity case for all constellations. Further examining the coverage values, we see that: i In Table 3, the cases with 0% 4-fold coverage show the GDOP, PDOP, 𝐢/𝑁0 and the received power values as N/A, meaning not applicable, as DOP metrics have no physical meaning outside of instances where PNT solutions exist. ii In instances where the coverage does not change between Tables 2 and 3, the average 𝐢/𝑁0 and received power values are the same (i.e., Starlink O-2), but for instances that change (i.e., Starlink O-1), the better coverage instances include more conditions with weak signals, thus both the average 𝐢/𝑁0 and received power values decrease. iii It can be seen from Table 3 that Starlink, Xona and Centispace are the constellations that provide resilience to indoor and low sensitivity conditions, and are still able to provide acceptable DOP metrics given the four-fold coverage is achieved. In the less challenging scenarios (O-2, O-3 and O-4), every LEO constellation is able to achieve acceptable metrics in low sensitivity conditions, with Starlink, OneWeb and Xona being able to achieve full coverage. In comparison, Galileo is able to either operate fully, as seen in O-2 and O-3 scenarios, or not able to operate at all due to the low sensitivity of the receivers. iv Table 2 shows a different picture, that even in high-sensitivity conditions, it can be challenging to achieve a four-fold coverage. Centispace and the experimental designs show a similar picture; they are able to achieve acceptable LEO-PNT metrics with a relatively low number of LEO satellites, yet they struggle to achieve a complete four-fold coverage; while Xona’s design seems to be better performing but still cannot guarantee complete four-fold coverage in all scenarios. Among the mega-constellations, Starlink and OneWeb achieve the best performance and are very similar, with Kuiper falling behind; providing a better LEO-PNT performance compared to the smaller LEO constellations, but failing to achieve a complete four-fold coverage in any scenario. Comparing the rest of Table 2 to the Galileo entry, which serves as the GNSS benchmark, it is clear that similar, if not better, PNT performance can be achieved with the considered LEO constellation designs in comparison. Galileo values serve as a benchmark and show what metrics are obtained with the current GNSS systems in the considered scenarios, as well as what potential a LEO-PNT constellation has. Values in Table 2 in particular, as the sensitivity allows for very high coverage for almost all designs, show what the main appeal of a dedicated LEO-PNT constellation is over the GNSS constellations; a significant improvement in 𝐢/𝑁0 for the same frequency band, and therefore, the possibility of shifting operations to higher frequency bands, which are less crowded. While the best performing LEO constellation in this study is indeed Starlink, we would like to emphasize that this does not strictly mean that a mega-constellation is necessary to achieve good LEO-PNT services. In fact, the four-fold coverage difference seen between Kuiper, Xona and Centispace shows that simply increasing the number of satellites within a constellation does not directly translate into improved PNT performance. The smaller scale LEO-PNT constellations and the author team’s experimental designs provide a more balanced approach between cost/complexity of the constellation and its PNT performance. Their 𝐢/𝑁0 is a direct improvement compared to Galileo, and the DOP metrics, while slightly worse, are still within the same performance range of < 10. We would like to further demonstrate two aspects via Fig. 1. The first aspect is that, Similar trends are seen for both I-1 and O-2 scenarios. Indoor and outdoor conditions affect the constellations in the same way by lowering their 𝐢/𝑁0 values drastically; which can impact the coverage depending on receiver sensitivity levels. Fig. 1a and 1b show the LEO-PNT metrics of the considered constellations with respect to each other for the high-sensitivity case, for two selected scenarios; I-1 and O-2. Looking at Fig. 1a and 1b, Starlink and Oneweb yields drastically better 𝐢/𝑁0 and GDOP than Galileo while equally maintaining the complete coverage, and Kuiper showing slightly worse 𝐢/𝑁0 and GDOP compared to the other mega-constellations, yet fails to keep up with the complete coverage provided by Galileo. Aside from the mega-constellations, Xona, Centispace and the experimental designs also provide 𝐢/𝑁0 improvements over Galileo. Xona yields a slight improvement with respect to the GDOP, while Centispace and the experimental designs have slightly worse GDOP compared to Galileo, but still remain within good GDOP value ranges. We can also see that the average number of visible satellites are higher compared to Galileo for the mega-constellations, Xona, and the experimental designs, which (a) I-1 for -185 dBm (c) GDOP in I-1 for -185 dBm (b) O-2 for -185 dBm (d) 𝐢/𝑁0 in I-1 for -185 dBm Figure 1: Left: LEO-PNT metrics for all constellations; Right: GDOP and 𝐢/𝑁0 versus carrier frequency in I-1 imply better stability in cases of satellite failures. The second aspect is the carrier frequency effect on LEO-PNT. Fig. 1c and 1d presents GDOP and 𝐢/𝑁0 respectively for a subset of four constellations; Galileo, Oneweb, Centispace, and Experimental 1; representing the GNSS benchmark, mega-constellations, small-scale LEO constellations, and potential LEO-PNT constellations respectively. Scenario I-1 is shown as an example, with the observation that similar trends have been noticed for O-2 scenario as well. The GDOP is not influenced much by the carrier frequency, with the exception of Galileo, whose signals fail to penetrate the indoor scenario with enough strength to be received by the receiver after around 20 GHz. On the other hand, 𝐢/𝑁0 drops in a consistent manner as for all constellations as the carrier frequency increases. Together, this shows us that going above X-band or Ku-band for the carrier frequency would risk operational failure for GNSS systems, while LEO-PNT systems would be able to support up to Ka-band carrier frequencies. We note that LEO constellations operating in K-band can reach to the same 𝐢/𝑁0 values that the GNSS constellations have in L-band, and that if the same carrier frequency is used, the improvement is significant in favor of the LEO constellations versus Galileo, and GNSS by extension. 4. Conclusion This article has shown an extensive comparison between seven LEO satellite constellations and Galileo constellation as the GNSS benchmark (provided in Table 1) in terms of LEO-PNT metrics. We selected the coverage, the GDOP, and the 𝐢/𝑁0 as relevant LEO-PNT performance metrics to be in used comparisons. Using a MATLAB simulation created in-house for analysis, Tables 2 and 3 detail the impact of different scenarios and receiver sensitivities on these performance metrics. Indoor/outdoor, rural/urban, and LOS/NLOS scenarios are considered in the link budget, and compared for low (i.e., -125 dBm) versus high (i.e., -185 dBm) receiver sensitivity cases. We show that, among the available and upcoming LEO constellations, Starlink, OneWeb, and Xona are the most promising for future LEO-PNT applications. We also show that experimental designs that can reach similar GDOP and coverage performance as Galileo (and GNSS by extension) are possible, which provide a direct improvement on 𝐢/𝑁0 . This fact also hints at the possibility that alternative LEO-PNT constellation designs can be further optimized for even better LEO-PNT solutions. Acknowledgments This work was supported by the Jane and Aatos Erko Foundation and by Teknologiateollisuus 100-year Foundation, under the project INCUBATE. The Authors also thank Prof. B. Eissfeller, from University of the Bundeswehr Munich for his constructive feedback on our LEO-PNT-related research in the team, and Dr. R. Morales-Ferre, for developing parts of the custom simulator used in this work. References [1] A. Zeedan, T. Khattab, Cubesat communication subsystems: A review of on-board transceiver architectures, protocols, and performance, IEEE Access 11 (2023) 88161–88183. doi:10.1109/ ACCESS.2023.3304419. [2] M. Centenaro, C. E. Costa, F. Granelli, C. Sacchi, L. Vangelista, A survey on technologies, standards and open challenges in satellite iot, IEEE Comm. Surveys & Tutorials 23 (2021) 1693–1720. doi:10.1109/COMST.2021.3078433. [3] I. F. Akyildiz, A. Kak, S. Nie, 6G and beyond: The future of wireless communications systems, IEEE Access 8 (2020). doi:10.1109/ACCESS.2020.3010896. [4] R. Morales Ferre, J. Praks, G. Seco-Granados, E. S. Lohan, A feasibility study for signal-in-space design for LEO-PNT solutions with miniaturized satellites, IEEE Journal on Miniaturization for Air and Space Systems 3 (2022) 171–183. doi:10.1109/JMASS.2022.3206023. [5] U. of Concerned Scientists, UCS satellite database, 2023. URL: https://www.ucsusa.org/resources/ satellite-database. [6] J. C. McDowell, The low earth orbit satellite population and impacts of the SpaceX Starlink constellation, The Astrophysical Journal Letters 892 (2020) L36. doi:10.3847/2041-8213/ ab8016. [7] Starlink, SpaceX non-geostationary satellite system - technical information to supplement schedule, 2020. URL: https://fcc.report/IBFS/SAT-MOD-20200417-00037/2274316.pdf. [8] Starlink, Amendment to pending application for the SpaceX gen2 NGSO satellite system, 2021. URL: https://fcc.report/IBFS/SAT-AMD-20210818-00105/12943361.pdf. [9] Y. Henri, The OneWeb Satellite System, Handbook of Small Satellites: Technology, Design, Manu- facture, Applications, Economics and Regulation, Springer, Cham, 2020, pp. 1–10. doi:10.1007/ 978-3-030-20707-6_67-1. [10] K. S. LLC, Application of Kuiper systems LLC for authority to launch and operate a non- geostationary satellite orbit system in ka-band frequencies: Technical appendix, 2021. URL: https://apps.fcc.gov/els/GetAtt.html?id=285359&x=. [11] FCC, Exhibit 1 - FCC form 442 / HawkEye 360 Pathfinder Cluster / item 7 and 8: Purpose of experiment and duration, 2020. URL: https://apps.fcc.gov/els/GetAtt.html?id=186545&x=. [12] ESA, ICEYE ESA archive, about ICEYE, 2023. URL: https://earth.esa.int/eogateway/missions/iceye. [13] H. J. Krame, Blacksky constellation, 2016. URL: https://www.eoportal.org/satellite-missions/ blacksky-constellation. [14] K. Γ‡elikbilek, Z. Saleem, R. Morales Ferre, J. Praks, E. S. Lohan, Survey on optimization methods for LEO-satellite-based networks with applications in future autonomous transportation, Sensors 22 (2022). doi:10.3390/s22041421. [15] L. Yang, The Centispace-1: A LEO satellite-based augmentation system, in: 14th Meeting of the International Committee on Global Navigation Satellite Systems, International Committee on GNSS (ICG), 2019, pp. 342–352. Accessed on: June 25, 2024. [16] Z. Yang, H. Liu, C. Qian, B. Shu, L. Zhang, X. Xu, Y. Zhang, Y. Lou, Real-time estimation of low earth orbit (LEO) satellite clock based on ground tracking stations, Remote Sensing 12 (2020). doi:10.3390/rs12122050. [17] X. S. Systems, Xona Space PULSAR constellation description, 2022. URL: https://www.xonaspace. com/pulsar. [18] J. Harms, The Orbcomm experience, 2004. URL: https://artes.esa.int/sites/default/files/1_The_ Orbcomm_Experience.pdf. [19] C. Fossa, R. Raines, G. Gunsch, M. Temple, An overview of the IRIDIUM (R) low earth orbit (LEO) satellite system, in: Proceedings of the IEEE 1998 National Aerospace and Electronics Conference. NAECON 1998. Celebrating 50 Years (Cat. No.98CH36185), IEEE, IEEE, 1998, pp. 152–159. doi:10.1109/NAECON.1998.710110. [20] F. J. Dietrich, P. Metzen, P. Monte, The Globalstar cellular satellite system, IEEE Transactions on Antennas and Propagation 46 (1998) 935–942. doi:10.1109/8.686783. [21] T. Reid, T. Walter, P. Enge, D. Lawrence, H. Cobb, G. Gutt, M. O’Connor, D. Whelan, Navigation from low earth orbit–part 1: Concept, current capability, and future promise, in: J. Morton, F. van Diggelen, J. J. Spilker, B. Parkinson (Eds.), Position, Navigation, and Timing Technologies in the 21st Century, volume 2, John Wiley & Sons, Ltd, 2021, pp. 1359–1379. doi:10.1002/9781119458555. ch43a. [22] M. Hartnett, Performance Assessment of Navigation Using Carrier Doppler Measurements from Multiple LEO Constellations, Master’s thesis, Air Force Institute of Technology, Ohio, USA, 2022. [23] Z. Kassas, Navigation from low earth orbit–part 2: Models, implementation, and performance, in: J. Morton, F. van Diggelen, J. J. Spilker, B. Parkinson (Eds.), Position, Navigation, and Timing Technologies in the 21st Century, volume 2, John Wiley & Sons, Ltd, 2021, pp. 1381–1412. doi:10. 1002/9781119458555.ch43b. [24] A. Nardin, F. Dovis, J. A. Fraire, Empowering the tracking performance of LEO-based positioning by means of meta-signals, IEEE Journal of RF Identification 5 (2021) 244–253. doi:10.1109/JRFID. 2021.3077082. [25] M. Neinavaie, J. J. Khalife, Z. Kassas, Blind doppler tracking and beacon detection for opportunistic navigation with LEO satellite signals, 2021 IEEE Aerospace Conference (50100) (2021) 1–8. doi:10. 1109/AERO50100.2021.9438258. [26] C. Zhao, H. Qin, Z. Li, Doppler measurements from multiconstellations in opportunistic navigation, IEEE Transactions on Instrumentation and Measurement 71 (2022) 1–9. doi:10.1109/TIM.2022. 3147315. [27] Z. Kassas, S. Kozhaya, J. Saroufim, H. Kanj, S. Hayek, A look at the stars: Navigation with multi-constellation LEO satellite signals of opportunity, Inside GNSS (2023). [28] J. Saroufim, S. W. Hayek, Z. M. Kassas, Simultaneous LEO satellite tracking and differential LEO-aided IMU navigation, in: 2023 IEEE/ION Position, Location and Navigation Symposium, 2023, pp. 179–188. doi:10.1109/PLANS53410.2023.10140087. [29] Y. Li, H. Li, W. Liu, L. Liu, W. Zhao, Y. Chen, J. Wu, Q. Wu, J. Liu, Z. Lai, H. Qiu, A Networking Perspective on Starlink’s Self-Driving LEO Mega-Constellation, 17, Association for Computing Machinery, New York, NY, USA, 2023, pp. 1–16. doi:10.1145/3570361.3592519. [30] T. G. Reid, B. Chan, A. Goel, K. Gunning, B. Manning, J. Martin, A. Neish, A. Perkins, P. Tarantino, Satellite navigation for the age of autonomy, in: 2020 IEEE/ION Position, Location and Navigation Symposium (PLANS), IEEE, IEEE, 2020, pp. 342–352. [31] K. Zhang, L. Tang, C. Zhao, S. Zhong, H. Luo, Perfect hash-based routing lookup for LEO constellation backbone network, IEEE Transactions on Aerospace and Electronic Systems 59 (2023) 4844–4857. doi:10.1109/TAES.2023.3244897. [32] S. Ma, Y. C. Chou, H. Zhao, L. Chen, X. Ma, J. Liu, Network characteristics of LEO satellite constellations: A Starlink-based measurement from end users, in: IEEE INFOCOM 2023 - IEEE Conference on Computer Communications, 2023, pp. 1–10. doi:10.1109/INFOCOM53939.2023. 10228912. [33] T. M. Inc., MATLAB version: 9.13.0 (r2022b), 2022. URL: https://www.mathworks.com. [34] M. Xu, Centispace system experimental satellites PPP application test, in: Proceedings of the 2023 UN/Finland Workshop on GNSS, 2023, pp. 1–21. Accessed on: June 25, 2024. [35] X. S. S. Inc., Satellite for transmitting a navigation signal in a satellite constellation system, 2021. URL: https://patents.justia.com/patent/11668834. [36] E. S. A. (ESA), Galileo satellites, 2022. URL: https://www.esa.int/Applications/Navigation/Galileo/ Galileo_satellites. [37] E. S. Lohan, K. Borre, Accuracy limits in multi-GNSS, IEEE Transactions on Aerospace and Electronic Systems 52 (2016) 2477–2494. doi:10.1109/TAES.2016.150241. [38] S. Jaeckel, L. Raschkowski, K. BΓΆrner, L. Thiele, QuaDRiGa: A 3-D multi-cell channel model with time evolution for enabling virtual field trials, IEEE Transactions on Antennas and Propagation 62 (2014) 3242–3256. doi:10.1109/TAP.2014.2310220. [39] W. Wang, E. S. Lohan, J. Talvitie, MmWave solutions to cope with the data link challenges in aviation, Deliverable D4.3 within H2020-SESAR-2016-2, EMPHASIS, 2019. Table 2 LEO-PNT Performance Metrics for βˆ’185 dBm Receiver Sensitivity (high sensitivity) Coverage [%] Constellation Scenario GDOP [-] PDOP [-] 𝐢/𝑁0 [dB-Hz] Received Power [dB] Path-Loss [dB] 4-Fold 1-Fold I-1 100 100 3.88 3.34 8.67 βˆ’160.30 226.88 I-2 94.77 96.94 4.04 3.48 0.64 βˆ’168.34 236.43 O-1 100 100 3.86 3.31 24.80 βˆ’144.18 210.68 Galileo O-2 100 100 3.85 3.31 50.66 βˆ’118.32 184.82 O-3 100 100 3.85 3.31 51.07 βˆ’117.91 184.41 O-4 100 100 3.85 3.31 33.93 βˆ’135.04 201.54 I-1 99.23 99.25 1.18 1.05 14.39 βˆ’154.59 238.53 I-2 96.48 96.88 1.32 1.17 3.27 βˆ’165.70 252.01 O-1 99.27 99.29 1.18 1.05 35.27 βˆ’133.71 217.21 Kuiper O-2 99.27 99.29 1.17 1.04 64.05 βˆ’104.93 188.43 O-3 99.27 99.29 1.18 1.05 64.34 βˆ’104.63 188.13 O-4 99.27 99.29 1.18 1.04 48.37 βˆ’120.61 204.11 I-1 99.58 99.63 0.61 0.54 9.55 βˆ’159.43 233.00 I-2 93.15 93.92 0.63 0.56 0.45 βˆ’168.53 244.89 O-1 100 100 0.60 0.53 28.85 βˆ’140.12 212.62 OneWeb O-2 100 100 0.60 0.53 56.00 βˆ’112.97 185.47 O-3 100 100 0.60 0.53 56.32 βˆ’112.66 185.16 O-4 100 100 0.60 0.53 40.96 βˆ’128.01 200.51 I-1 100 100 0.89 0.80 19.16 βˆ’149.81 226.59 I-2 99.38 99.46 0.90 0.81 7.58 βˆ’161.40 239.21 O-1 100 100 0.86 0.77 39.40 βˆ’129.58 206.18 Starlink O-2 100 100 0.84 0.75 66.43 βˆ’102.55 179.15 O-3 100 100 0.83 0.73 66.73 βˆ’102.25 178.85 O-4 100 100 0.84 0.74 51.53 βˆ’117.45 194.05 I-1 96.77 100 5.94 5.31 40.70 βˆ’128.27 200.77 I-2 96.77 100 5.93 5.29 32.28 βˆ’136.70 209.20 O-1 96.77 100 5.93 5.31 57.57 βˆ’111.41 183.91 Centispace O-2 96.77 100 5.95 5.31 78.89 βˆ’90.09 162.59 O-3 96.77 100 5.94 5.31 79.24 βˆ’89.73 162.23 O-4 96.77 100 5.93 5.31 65.39 βˆ’103.59 176.09 I-1 100 100 3.03 2.71 40.49 βˆ’128.49 200.99 I-2 100 100 3.05 2.71 32.32 βˆ’136.65 209.16 O-1 100 100 3.03 2.71 57.31 βˆ’111.67 184.17 Xona O-2 100 100 3.05 2.73 78.53 βˆ’90.44 162.94 O-3 96.9 100 3.04 2.72 78.89 βˆ’90.09 162.59 O-4 100 100 3.04 2.72 65.11 βˆ’103.86 176.36 I-1 99.71 100 5.36 4.77 14.94 βˆ’154.03 230.73 I-2 93.10 95.48 5.27 4.67 3.78 βˆ’165.19 242.98 O-1 99.98 100 5.37 4.76 34.84 βˆ’134.13 210.73 Experimental 1 O-2 100 100 5.37 4.78 62.09 βˆ’106.89 183.49 O-3 100 100 5.37 4.78 62.40 βˆ’106.57 183.17 O-4 99.98 100 5.36 4.76 47.07 βˆ’121.90 198.50 I-1 99.19 100 5.98 5.36 18.34 βˆ’150.64 227.32 I-2 92.5 97.15 5.89 5.27 7.20 βˆ’161.77 240.10 O-1 99.71 100 5.91 5.32 37.06 βˆ’131.91 208.51 Experimental 2 O-2 99.75 100 5.92 5.33 63.95 βˆ’105.03 181.63 O-3 99.75 100 5.93 5.33 64.25 βˆ’104.72 181.32 O-4 99.75 100 5.94 5.32 49.19 βˆ’119.78 196.38 Table 3 LEO-PNT Performance Metrics for βˆ’125 dBm Receiver Sensitivity (low sensitivity) Coverage [%] Constellation Scenario GDOP [-] PDOP [-] 𝐢/𝑁0 [dB-Hz] Received Power [dB] Path-Loss [dB] 4-Fold 1-Fold I-1 0 0.38 N/A N/A N/A N/A 226.88 I-2 0 0.38 N/A N/A N/A N/A 236.43 O-1 0 0.38 N/A N/A N/A N/A 210.68 Galileo O-2 100 100 3.86 3.31 50.66 βˆ’118.32 184.82 O-3 100 100 3.85 3.30 51.07 βˆ’117.91 184.41 O-4 0 0.60 N/A N/A N/A N/A 201.54 I-1 0 3.21 N/A N/A N/A N/A 238.53 I-2 0 1.19 N/A N/A N/A N/A 252.01 O-1 0 5.02 N/A N/A N/A N/A 217.21 Kuiper O-2 99.27 99.29 1.18 1.05 64.05 βˆ’104.93 188.43 O-3 99.27 99.29 1.18 1.05 64.34 βˆ’104.63 188.13 O-4 99.27 99.29 1.28 1.14 48.38 βˆ’120.60 204.11 I-1 0 2.15 N/A N/A N/A N/A 233.00 I-2 0 1.67 N/A N/A N/A N/A 244.89 O-1 0 2.23 N/A N/A N/A N/A 212.62 OneWeb O-2 100 100 0.60 0.53 56.00 βˆ’112.97 185.47 O-3 100 100 0.60 0.53 56.32 βˆ’112.66 185.16 O-4 100 100 3.46 2.72 45.39 βˆ’123.58 200.51 I-1 5.81 9.96 4.17 3.34 46.80 βˆ’122.18 226.59 I-2 0 1.94 N/A N/A N/A N/A 239.11 O-1 80 81.52 3.47 2.79 45.75 βˆ’123.22 206.18 Starlink O-2 100 100 0.87 0.78 66.43 βˆ’102.55 179.15 O-3 100 100 0.86 0.77 66.73 βˆ’102.25 178.85 O-4 100 100 0.91 0.82 51.55 βˆ’117.43 194.05 I-1 33.04 70.35 7.67 6.83 51.08 βˆ’117.90 200.77 I-2 8.06 34.27 5.68 4.86 48.06 βˆ’120.92 209.20 O-1 96.52 100 5.97 5.28 57.59 βˆ’111.39 183.91 Centispace O-2 96.77 100 5.93 5.31 78.89 βˆ’90.09 162.59 O-3 96.77 100 5.93 5.31 79.24 βˆ’89.73 162.23 O-4 96.77 100 5.88 5.27 65.39 βˆ’103.59 176.09 I-1 55.65 72.92 5.91 5.12 49.99 βˆ’118.98 200.99 I-2 19 40.77 8.33 7.26 48.16 βˆ’120.81 209.16 O-1 100 100 3.14 2.78 57.34 βˆ’111.64 184.17 Xona O-2 100 100 3.03 2.71 78.53 βˆ’90.44 162.94 O-3 100 100 3.03 2.71 78.89 βˆ’90.09 162.59 O-4 100 100 3.05 2.72 65.11 βˆ’103.86 176.36 I-1 0 0 N/A N/A N/A N/A 230.73 I-2 0 0 N/A N/A N/A N/A 242.98 O-1 0 0 N/A N/A N/A N/A 210.73 Experimental 1 O-2 99.98 100 5.37 4.78 62.09 βˆ’106.89 183.49 O-3 99.98 100 5.37 4.78 63.12 βˆ’106.57 183.17 O-4 87.67 100 5.27 4.72 47.05 βˆ’121.92 198.50 I-1 0 2.56 N/A N/A N/A N/A 227.32 I-2 0 0.04 N/A N/A N/A N/A 241.10 O-1 0 8.23 N/A N/A N/A N/A 208.51 Experimental 2 O-2 99.67 100 5.94 5.35 63.95 βˆ’105.03 181.63 O-3 99.67 100 5.93 5.35 64.25 βˆ’104.72 181.32 O-4 91.52 100 6.08 5.40 49.19 βˆ’119.78 196.38