Tag: Per Enge

  • Per Enge memorial webcast this Saturday

    Per Enge memorial webcast this Saturday

    On Saturday, Nov. 10, Stanford colleagues of Professor Per Enge will host a celebration of his life. A live webcast of the event will be available here at  at 1 p.m. Pacific Standard Time.

    The video will be available afterwards on this site.

    GPS World extends this invitation to join Per’s family and friends in celebrating the wonderful life that he led and the extraordinary impact he had on those around him. Guests at the live event will be invited to share their favorite memories.

    The following statement was recently read into the minutes of the Stanford Faculty Senate:


    Per K. Enge, the Vance D. and Arlene C. Coffman Professor in the School of Engineering and one of the world’s foremost experts in GPS technologies, passed away on April 22, 2018, at his home in Mountain View, California. He was 64.

    Per Enge, Professor and Director, Stanford university Center for Position Navigation and Time
    Per Enge, Professor and Director, Stanford University Center for Position Navigation and Time. (Photo: Stanford University)

    Per was born Oct. 29, 1953, in Bergen, Norway. He immigrated at the age of 2 to the United States with his father and mother.

    He earned his B.S. in electrical engineering at the University of Massachusetts at Amherst in 1975 and his MS and PhD at the University of Illinois Urbana-Champaign in 1979 and 1983, respectively. He met his wife of 38 years, Elaine, while at UMass. His son, Nick, a Stanford graduate and now a lecturer at the University of Texas at Austin, said that despite his academic upbringing, his dad was a middling student until Elaine introduced him to the library at UMass, where she was most often found.

    While pursuing his advanced degrees, Per worked in industry, where he first gained experience using radio signals for terrestrial navigation. He then took a position as assistant professor at Worcester Polytechnic Institute, where he designed a radio-navigation positioning system that today has more than 1.5 million marine and land users. He also started to work on augmenting GPS so that it could be safely used for aeronautical navigation. Due to this effort he was recruited by Stanford University for a one-year visiting professorship in 1993 that was ultimately extended to full professorship.

    While at Stanford, Per became one of the FAA’s most trusted advisors on the provision of aircraft guidance. He oversaw the development of two different systems that today allow airplanes to safely determine their positions within meters regardless of the weather conditions.

    Per was a member of the National Academy of Engineers, a Fellow of the Institute of Navigation, and a Fellow of the Institute of Electrical and Electronics Engineers.

    Per is particularly remembered as a teacher and mentor. He designed a freshman course in electric cars and aircraft and helped launch a popular massive open online course (MOOC) for the GPS community outside Stanford. He leaves behind a strong legacy of more than 40 Ph.D. students, co-workers and colleagues who have been inspired by his genuine joy in being able to work in such an exciting field.

    Per is survived by his wife, Elaine, of Mountain View, and a son, Nick, of Austin. In lieu of flowers, donations in Enge’s memory can be made to a new graduate student scholarship fund under the Stanford department of Aeronautics and Astronautics. Donations can be made at https://gps.stanford.edu/resources/giving

     

     

  • In memoriam Per Enge

    With great sadness we must report that Per Enge passed away on April 22, at home and surrounded by family. Per was a genial friend and colleague to many, and a pillar of the PNT community. He is greatly missed by all.

    At the culmination of his long, fruitful career he served as the Vance and Arlene Coffman Professor of Engineering at Stanford University, where he also directed the Stanford Center for Position Navigation and Time.

    For many years he conducted research funded by the Federal Aviation Administration, directed at safe and secure air navigation and leading to development of the Wide Area Augentation System (WAAS) and Local Area Augmentation Systems (LAAS). WAAS became fully operational for aviation in the United States in 2003 and is currently carried by more than 110,000 aircraft; similar systems have been deployed in Europe, Japan and India.

    Per Enge at National Cheng Kung University (courtesy Shau Shiu Jan).

    He received the Kepler Award from the Institute of Navigation in 2000 and was inducted into the GPS Hall of Fame by the U.S. Air Force in 2012. He served as a member of the Space-Based Position Navigation and Time Federal Advisory Committee since 2007. In 2013 he received the GNSS Leadership Award for Signals from this magazine, for signal design including national differential GPS, satellite-based augmentation systems, and alternative positioning, navigation and timing sources. He co-wrote Global Positioning System: Signals, Measurements, and Performance.

    Always an educator, Per served as instructor, mentor and gentle encourager of many, many Ph.D. and other graduate-level students at Stanford who have gone on to distinguished careers of their own. In a lifetime marked by great achievements, this is perhaps his greatest and ultimately will be the most far-reaching.

    Born in Norway and brought to the U.S. at age 2, he received a B.S.E.E. from the University of Massachusetts and M.S.E.E. and Ph.D. degrees from the University of Illinois.

    Speaking at GPS World dinner, accepting Signals Leadership 2013 award. (Photo: GPS World file)

    In remarks on accepting the GNSS Leadership Award for Signals, Per cited Faflick’s theorem, “that you will never ever work on any projects that are both interesting and important.” After calling out both GNSS and WAAS as exceptions to the theorem, he identified a third outlier: spoofing.

    “Today’s e-security is based on three security factors: what we know (passwords), what we carry (key fob), and what we are (fingerprints, iris scan). And it is not enough. To meet this challenge, we need to rejuvenate the original security factor: location. In the past, transactions were secured by our presence. In the world of e-commerce, this factor has disappeared, and we must use GNSS to approximate this ancient and effective security factor.

    “All of this will require the best effort of this precious community of ours.”

    Further biographical details are available in an article published by the Stanford News. Among the tributes included there is this one by Brad Parkinson, who recruited Enge to Stanford in the early 1990s. “Anyone who works in GPS is aware of Per and his influence. He was just an intellectually talented person who could understand many scientific nuances and integrate them in ways others could not.”

    Teaching the massive online open course.

    The article also reminds us that he co-originated and co-taught, with Frank van Diggelen, a massive open online course to share GPS knowledge with a worldwide audience, far beyond Stanford’s walls. Titled “GPS: An Introduction to Satellite Navigation, with an interactive Worldwide Laboratory using Smartphones,” it enrolled 31,000 people from 192 countries.  It is available here.

    Per’s Stanford colleagues Sherman Lo, Todd Walter and Sam Pullen assisted GPS World with this article and provided these photos from personal archives. The Stanford group is working on setting up a scholarship in Per’s name.  More information on it and how to support it will be on the SCPNT website once it becomes available.

    Frank van Diggelen has sent further photos, below.

    At the Stanford GPS Lab with colleagues from Stanford and DLR (German aerospace agency).

     

     

     

     

     

     

     

    Dinner discussions with US-EU bilateral group.
    With Alan Chen, Sherman Lo and an early spectrum image of GIOVE-A (or Galileo).
    Visiting Neuschwanstein Castle in Bavaria after 2005 European Navigation Conference.
    At the Stanford Center Position, Navigation and Time, which he co-founded in 2005.
    Team China Consumers at GPS World dinner 2010. The winning team in the Grand Game of GNSS.
    Fierce “opponents” (examiners) for Ph.D. defense of Ignacio Fernández Hernández of EC/Galileo. Aalborg University, Denmark.
    Prepared to come aboard in Kobenhavn.
    The co-authors of Global Positioning System: Signals, Measurements, and Performance (with Pratap Misra).

     

    A road warrior for GNSS.

     

     

    Faculty of the GNSS Summer School at Svalbard, Norway (Arctic Ocean, 78.7° N).

     

     

     

     

  • Per Enge appointed to Satelles board of directors

    Per Enge appointed to Satelles board of directors

    Per Enge, Professor and Director, Stanford university Center for Position Navigation and Time

    Satelles, a secure time and location solutions company, has appointed Per Enge to its board of directors. Satelles provides a time and location solutions delivered over the Iridium constellation of 66 low-earth-orbiting satellites.

    Enge is the Vance and Arlene Coffman Professor of Aeronautics and Astronautics for Stanford University, where he is also the director of the Stanford Center for Position Navigation and Time.

    “I am eager to join the Satelles Board of Directors and look forward to supporting the management team,” Enge said. “I am encouraged by the progress Satelles has made and continue to have confidence in the leadership team and future growth of the business.”

    Enge’s laboratory has worked with the U.S. Coast Guard to design a medium frequency radio system to broadcast differential GPS corrections to maritime users, and this system has been implemented as a worldwide standard.

    His laboratory also worked with the U.S. Federal Aviation Administration to develop WAAS, the Wide-Area Augmentation System that provides GPS integrity data to airborne users. Today, WAAS is carried by more than 100,000 aircraft, and similar systems have been implemented in Europe, India and Japan.

    Enge also serves on the board of directors of Amida Technologies, and he serves as a technical advisor to Polaris Wireless.

    He has received the Kepler, Thurlow and Burka Awards from the Institute of Navigation for his work. He is a Fellow of the Institute of Electrical and Electronics Engineers. He is a member of the National Academy of Engineering and a fellow of the Institute of Navigation.

    Enge received his Ph.D. in electrical engineering from the University of Illinois in 1983. In 2012, the U.S. Air Force inducted Enge into the GPS Hall of Fame.

    “It is with great pleasure that we welcome Per to Satelles Board of Directors,” said Michael O’Connor, Satelles CEO. “Per has distinguished himself as a technology innovator and brings to our board of directors deep expertise in global navigation satellite systems. His wealth of experience and expertise in GPS and other technologies adds new depth to our board as we continue to deliver Satellite Time and Location  to users around the world. We look forward to working with Per on our mission is to deliver trusted time and location solutions that augment and enhance existing solutions — including GPS.”

  • Expert Opinions: The effect of LEO constellations on GNSS services

    Expert Opinions: The effect of LEO constellations on GNSS services

    Q: What is the potential for low-Earth orbit constellations to augment services provided by the four medium-Earth orbit GNSS?

    Doug Taggart, President, Overlook Systems Technologies, Inc.
    Doug Taggart, President, Overlook Systems Technologies, Inc.

    A: With more than one hundred GNSS satellites broadcasting on three or more frequencies, our international constellation of medium-Earth orbiting (MEO) satellites will provide a combination of path diversity and frequency diversity. However, satellites in low-Earth orbit (LEO) should be added to our MEO mélange to provide orbital diversity and thus cyber safety. The LEO satellites would have 20 dB less path loss and compel jammers and spoofers to become conspicuous. Even with only one LEO in view, we would be able to use the LEO signal as a hot clock to improve the robustness of GNSS signal acquisition by our users. For timing applications, a solitary LEO satellite would enable time transfer to fixed locations worldwide.


    Per Enge, Professor and Director, Stanford university Center for Position Navigation and Time
    Per Enge, Professor and Director, Stanford university Center for Position Navigation and Time

    A: While it is prudent to take advantage of multiple PNT sources, the devil is in the details. Are users seeking more availability, accuracy, integrity and/or resilience to fill gaps? What is the complexity and cost for integration in user equipment, the reliability compared to other augmentations, the applications to be supported, vulnerability to interference, and so on? Additionally, all things from space may not be the best solution when all user needs and vulnerabilities are factored in.

  • 2013 Leadership Awards: The Honorees Speak

    2013 Leadership Awards: The Honorees Speak

    The GPS World 2013 Leadership Awards.
    The GPS World 2013 Leadership Awards.

    This article reproduces the acceptance speeches given by the winners of GPS World’s 2013 Leadership Awards, at the Leadership Dinner in Nashville in September.

    The Leadership Dinner was sponsored by Lockheed Martin and Exelis.

    Nominations for the 2013 Leadership Awards came from the Editorial Advisory Board of GPS World, from the 2012 Award winners, and from a select handful of industry executives. A similar group of GNSS community members, roughly double in size, voted for the final Award winners presented here.


    Success with a New System, QZSS

    Japanese Government Approves a Four-Satellite Constellation

    Headshot: Satoshi KogureIt is my great honor to receive this GPS World Leadership Award 2013 on behalf of my team. I’d like to express my gratitude to my family, my colleagues, and especially my team members. We jointly developed the first satellite of the Quasi-Zenith Satellite System in Japan.

    The project was established after long discussion between the Japanese government, industries, and user communities, as well as other GNSS providers. QZSS is now providing GPS availability and capability enhancement to Asia and the Oceania region as a regional augmentation system. We are very proud of a great achievement on the first satellite technical validation. Our first satellite can provide the same signal-in-space user range error performance as the latest GPS Block IIR-M and IIF space vehicles.

    As a result, considering our achievement, the Japanese government has decided to proceed with the QZSS program to the second phase. So we will have three additional satellites by 2017 and we’ll start the augmentation service surrounding Japan, East Asia, and the Pacific region from 2018, using the four-satellite constellation.

    The Japanese government is conducting the system design for the operational QZSS, and the Japanese Aerospace Exploration Agency is continually supporting the government’s activity. We expect that QZSS will bring great social benefits to life in the Asia-Pacific region in a wide variety of fields.

    Again, on behalf of all the people who have been involved in the QZSS project, I give my sincere appreciation for this reward. I would also like to acknowledge everybody in this room,  because we have learned so many things from this wonderful society. So it is my great pleasure to share my gratitude with you. Thank you very much.

    The One Constant: New Challenges

    Jamming, Spoofing, Squeezing, Security, and More

    Headshot: Per EngeFirst and foremost, thank you for this wonderful award. It caused me to reflect on all of the wonderful people that I have worked with over the years. This morning, I tried to list these people for NDGPS, WAAS, and APNT; and the lists were too long to repeat in five minutes. For WAAS alone, the list includes people from the Federal Aviation Administration, Stanford University, AJ Systems, RTCA, Raytheon, Zeta, Inc., and the Mitre Corporation. This beautiful award reflects on the efforts of everyone on these long lists.

    These lists also remind me of the remarkable community that is GNSS today — all of us. Apparently, we still have our work cut out for us. After all, I understand from this very conference that:

    • GNSS is blocked;
    • GNSS is jammed;
    • GNSS spectrum will be squeezed;
    • GNSS is spoofed;
    • GNSS is scintillated;
    • GNSS is too expensive;
    • GNSS is Swiss cheese (thank you, Nunzio).

    Honestly, I was not really aware that GNSS is all that fragile. It seems to work fine to me. But let us take the point: these threats are worth our worry, and mitigation will require the full strength of our far-flung community.
    Spoofing is worth some special mention. It is certainly interesting, and the theory will be gorgeous. However, I remind you of Faflick’s theorem:

    “In your professional life, you will work on many interesting projects. With luck, you will work on some important projects; those will be important to your company, your nation to the global community. However, (Faflick asserts) you will never ever work on any projects that are both interesting and important.”

    Well, it seems that GNSS and WAAS are exceptions. But two such small exceptions should not cause us to turn our back on a theorem as powerful or persuasive as the one provided by Faflick. If we apply the theorem, then spoofing cannot be important, because it is interesting. However, I think spoofing is also an exception to Faflick’s theorem.
    After all, today’s e-security is based on three security factors: what we know (passwords), what we carry (key fob), and what we are (fingerprints, iris scan). And it is not enough. For example, our health records will soon be online, and the damage caused by losing control of these health records would be great. To meet this challenge, we need to rejuvenate the original security factor: location. In the past, transactions were secured by our presence. In the world of e-commerce, this factor has disappeared, and we must use GNSS to approximate this ancient and effective security factor.

    All of this will require the best effort of this precious community of ours.

    Science Powers New Applications

    Radio Occultation Techniques May Warn of Natural Disasters

    Photo: Attila KomjathyI am deeply honored to receive GPS World’s GNSS Leadership Award this year. Many thanks to the magazine and the international group of GNSS experts for recognizing my contributions to the GNSS community with this prestigious distinction.

    I have been privileged to receive support and guidance from a great number of colleagues and friends over the years — too many to list, but I thank them all. If I am a good researcher, it is because I had excellent mentors. Chief among them is Professor Richard Langley of the University of New Brunswick, sitting in the audience, who believed in me and taught me to leave no stones unturned in my research and to continue questioning myself along the way.

    Professors George Born and Penina Axelrad also had immeasurable influence on me at the University of Colorado in becoming a GNSS researcher. I am greatly indebted to Tony Mannucci and Brian Wilson at NASA’s Jet Propulsion Laboratory for creating a stimulating, conducive, and rewarding research environment in my group at JPL. I would like to thank my wife Katie and sons David and Adam who have been most patient and who have supported me all these years.

    At the early stages of my career, working with GPS data at the University of New Brunswick, we routinely had access to about 40–50 GPS receivers worldwide, and this was after a laborious process of hunting down individual RINEX files. Now, we process data from about 1,200 GPS and GLONASS stations daily in an automated fashion including daily, hourly, and streaming GPS sites. Data centers including ours at JPL generate real-time products to retrieve three-dimensional electron densities using global assimilative ionospheric modeling and other techniques. As computational power increases, I envision that we will see new scientific discoveries evolving with real-life applications and tangible benefits to society. These emerging applications already include the detection of small perturbations in the ionosphere using GNSS real-time measurements that may be indicative of natural hazards generating acoustic and gravity waves in the atmosphere.

    Why is this important to us, other than being of scientific interest?

    The advances in GNSS ionospheric measurements may soon be capable of augmenting tsunami early-warning systems, for instance. Other potential applications include detecting volcanic eruptions in remote corners of our planet using GNSS data. Less innocent yet equally important may be detecting nuclear tests using GNSS ionospheric measurements.

    Another scientific application in the works is the creation of three-dimensional coupled ocean-atmosphere-ionosphere modeling capabilities. These models are at the initial phases of development, running on supercomputers including one at JPL. In a few years, we expect to see that coupled physics-based modeling will help detect and confirm natural hazards and artificial explosions on the Earth’s surface.

    In a not too distant future, I anticipate that we will see natural hazard-generated signatures detected using low-Earth orbiting satellites. Initial results suggest that some of the largest events including large earthquakes, tsunamis, and meteor events may be detected in radio occultation measurements. I expect to see a fleet of low-Earth orbiters monitoring Earth’s ionosphere in real time, looking for signatures of tsunamis, earthquakes, volcanic eruptions, meteor events, large forest fires, or nuclear tests.

    It is also tempting to use this technology and associated algorithms to claim ionospheric precursors of earthquakes. While I am not convinced about the technical feasibility of ionospheric earthquake precursors, it is a fascinating research area, which has generated heated debates in the community over the years.

    With the abundance of different constellations including GPS, GLONASS, BeiDou, QZSS, and others using ground and space-based observations, we are clearly entering a new era when mid-ocean 10-centimeter-level tsunami waves may be detected in the ionosphere at 400 kilometers altitude in real time! Could you have imaged even only 10 years ago that by looking at real-time ionospheric measurements we could infer 10-centimeter-level tsunami waves well before they reach a shoreline? This might help decision makers to save human lives and potentially billions in material damages.

    In closing, let me repeat that I am deeply honored to have received the GPS World Leadership Award. I look forward to continue pushing the boundaries of scientific and technical discoveries using GNSS technologies.

    Thank you very much.

    Getting to the Plus in GNSS

    Multi-Constellation, Environments, Multi-Sensors

    Headshot: Peter GrognardDear friends, dear colleagues, it is with pride that I receive this award tonight, in this very town where I attended ION GPS for the first time in 1998, and where I met some of you for the very first time.

    Septentrio didn’t formally exist then, but while at the Interuniversity Microelectronics Center IMEC in Belgium, and just returning from a four-year assignment in the United States, where I witnessed the beginning of civil GPS, I was preparing the creation of the company that I am so proud of, and that is first and foremost a very international family of unique, talented people that I enjoy working with every day!

    Hence, my first words of gratitude go to my world-class team that for 15 years now has been at the forefront of advanced GNSS receiver development.

    I also take the opportunity to thank some of our long-standing customers and partners, who all those years have given us their business and confidence. In particular, let me thank here tonight the European Space Agency, with whom we have worked hand in hand for more than 10 years. Thank you for your confidence and business — it is a tremendous honor to work with you, and we have enjoyed the truly historic milestones on the Galileo program as much as you have.

    As the motivation for tonight’s award, a number of important new developments were mentioned, as achieved by Septentrio throughout an already very busy year:

    • the successful position/velocity/time (PVT) calculation using BeiDou satellites in a three-constellation solution on a commercial PolaRx4 platform just days following the official release of BeiDou’s Interface Control Document;
    • the first autonomous real-time Galileo PVT calculation. The standalone position was calculated from in-orbit navigation messages using a standard PolaRx4 GNSS receiver equipped with commercially released firmware;
    • the first four-constellation PVT performed by a standard commercial receiver. The 3D-position fix happened shortly after the Galileo IOV satellites began transmitting, for the first time ever, a fully usable navigation message as part of an ESA experiment;
    • and last, but certainly not least, the first Galileo PRS-only positioning for European Secured Navigation using only the encrypted Galileo Public Regulated Service (PRS) signals broadcast by the four Galileo In-Orbit Validation satellites.

    I’ll share with you some future directions for us, outlining three major themes running across our developments.
    When we first met here 15 years ago, the focus of our industry was all about GPS. Today, our favorite yearly rendezvous has been rightly renamed ION GNSS PLUS. GPS is not alone anymore. Not alone, because who today would give up the great signal availability that multi-constellation reception brings us? Not alone, because who would claim that GPS can survive in its own bubble of safety without considering the environmental conditions it lives in?

    Finally, not alone because modern positioning and navigation require multi-sensor fusion to perform at their best.

    Multi-Constellation Reception. Septentrio has, from its very beginning, focused on developing multi-frequency, multi-constellation receivers, at a time when this was far from common. Our first proprietary ASIC was a highly integrated GPS-GLONASS-SBAS chip, which we used as the core of the first generation of Septentrio multi-system/multi-frequency receivers.

    “All the signals in the sky” was, is, and will always remain one of our key leitmotifs. We follow with great interest the new developments around the world, for instance for future-generation EGNOS and Galileo systems. As a practical example, I’d like to mention the work that we conduct with European partners on EGNOSv3.
    Also, it’s my pleasure to announce here tonight the development of our fourth-generation correlator ASIC, GReCo4: all–in-view, all-signals, all-Septentrio!

    The consistent focus on multi-system, multi-frequency receivers was obviously driven by the desire to make very powerful, robust GNSS receivers that would provide high-precision PVTs under the widest range of circumstances.

    Our unwavering support for Galileo from the early days was also based on the fact that we could, at the turn of the century, easily receive 15–20 satellites: GPS+GLONASS+EGNOS.

    Environments. Obviously, no matter how many independent systems are combined, GNSS remains vulnerable and sensitive to all kinds of interfering sources, human-made and made by Mother Nature. That’s why in recent years, we have been investing a lot of effort in understanding and mitigating the environmental disturbance that affects GNSS receivers.

    While we must be humble for, and cannot (completely) tame the forces of Mother Nature, Septentrio has spent a lot of effort on ionospherical receivers. These developments have been conducted by a Septentrio-led consortium of European and Brazilian partners, with as a tangible outcome the Septentrio PolaRxS receiver, aimed at measuring the ionosphere. This unique state-of-the-art receiver is being deployed worldwide by our customers and partners, which helps us improve our understanding of ionopheric behaviour and its impact on GNSS receivers. Our goal: making the most robust professional receivers in the world.

    With respect to interference, both intentional and unintentional, Septentrio will continue to improve and expand its solutions. Septentrio’s state-of-the-art Advanced Interference Mitigation (AIM) is AIMed at providing our customers the most robust reception under adverse and hostile conditions.

    Multi-Sensors. “All of the signals in the sky” only works if your receiver sees the sky. If there are no satellite signals — for instance, because of blockage of the signals in a container terminal — one obviously needs other sensors to take over and guarantee the continued availability of a PVT-solution.

    For several years, we have been offering combined GNSS-INS solutions for a variety of applications, such as management of container terminals. One of our very first customers was the Port of Singapore Authority.
    Septentrio will continue to work along these three exciting axes. We look forward to communicating new trends and developments to you wonderful, global GNSS+ community!

     

  • Innovation: Getting Control

    Innovation: Getting Control

    Off-the-Shelf Antennas for Controlled-Reception-Pattern Antenna Arrays

    By Yu-Hsuan Chen, Sherman Lo, Dennis M. Akos, David S. De Lorenzo, and Per Enge

    GPS World photo
    INNOVATION INSIGHTS by Richard Langley

    THE ANTENNA IS A CRITICAL COMPONENT OF ANY GNSS RECEIVING EQUIPMENT. It must be carefully designed for the frequencies and structures of the signals to be acquired and tracked. Important antenna properties include polarization, frequency coverage, phase-center stability, multipath suppression, the antenna’s impact on receiver sensitivity, reception or gain pattern, and interference handling. While all of these affect an antenna’s performance, let’s just look at the last two here.

    The gain pattern of an antenna is the spatial variation of the gain, or ratio of the power delivered by the antenna for a signal arriving from a particular direction compared to that delivered by a hypothetical isotropic reference antenna. Typically, for GNSS antennas, the reference antenna is also circularly polarized and the gain is then expressed in dBic units.

    An antenna may have a gain pattern with a narrow central lobe or beam if it is used for communications between two fixed locations or if the antenna can be physically steered to point in the direction of a particular transmitter. GNSS signals, however, arrive from many directions simultaneously, and so most GNSS receiving antennas tend to be omni-directional in azimuth with a gain roll-off from the antenna boresight to the horizon.

    While such an antenna is satisfactory for many applications, it is susceptible to accidental or deliberate interference from signals arriving from directions other than those of GNSS signals. Interference effects could be minimized if the gain pattern could be adjusted to null-out the interfering signals or to peak the gain in the directions of all legitimate signals. Such a controlled-reception-pattern antenna (CRPA) can be constructed using an array of antenna elements, each one being a patch antenna, say, with the signals from the elements combined before feeding them to the receiver. The gain pattern of the array can then be manipulated by electronically adjusting the phase relationship between the elements before the signals are combined. However, an alternative approach is to feed the signals from each element to separate banks of tracking channels in the receiver and form a beam-steering vector based on the double-difference carrier-phase measurements from pairs of elements that is subsequently used to weight the signals from the elements before they are processed to obtain a position solution. In this month’s column, we learn how commercial off-the-shelf antennas and a software-defined receiver can be used to design and test such CRPA arrays.


    “Innovation” features discussions about advances in GPS technology, its applications, and the fundamentals of GPS positioning. The column is coordinated by Richard Langley, Department of Geodesy and Geomatics Engineering, University of New Brunswick. To contact him with topic ideas, email him at lang @ unb.ca.


    Signals from global navigation satellite systems are relatively weak and thus vulnerable to deliberate or unintentional interference. An electronically steered antenna array system provides an effective approach to mitigate interference by controlling the reception pattern and steering the system’s beams or nulls. As a result, so-called controlled-reception-pattern-antenna (CRPA) arrays have been deployed by organizations such as the U.S. Department of Defense, which seeks high levels of interference rejection.

    Our efforts have focused on developing a commercially viable CRPA system using commercial off-the-shelf (COTS) components to support the needs of Federal Aviation Administration (FAA) alternative position navigation and timing (APNT) efforts. In 2010, we implemented a seven-element, two-bit-resolution, single-beam and real-time CRPA software receiver. In 2011, the receiver was upgraded to support all-in-view, 16-bit-resolution with four elements.

    Even though we can implement these CRPA software receivers in real time, the performance of anti-interference is highly dependent on the antenna array layout and characteristics of the antenna elements. Our beamforming approach allows us to use several COTS antennas as an array rather than a custom-designed and fully calibrated antenna. The use of COTS antennas is important, as the goal of our effort is to develop a CRPA for commercial endeavors — specifically for robust timing for the national airspace. Hence, it is important to study the geometry layout of the individual antennas of the array to assess the layouts and to determine how antenna performance affects the array’s use.

    In our work, we have developed a procedure for calculating the electrical layouts of an antenna array by differential carrier-phase positioning. When compared to the physical layout, the results of electrical layouts can be used to determine the mutual coupling effect of each combination. Using the electrical layout, the resultant gain patterns can be calculated and used to see the beamwidth and the side-lobe issue. This is important as these factors have significant effects on anti-interference performance. This study focuses on understanding the performance effects of geometry and developing a method for describing the best geometry.

    We adopted three models of COTS antenna and two possible layouts for a four-element array. Then, signal collection hardware consisting of four Universal Software Radio Peripheral (USRP) software-defined radios and one host personal computer was assembled to collect array data sets for each layout/antenna combination. Our developed CRPA software receiver was used to process all data sets and output carrier-phase measurements.

    In this article, we will present the pattern analysis for the two selected layouts and describe how we collected the experimental data. We’ll then show the results of calculating the electrical spacing for the layouts are compare them to the physical layouts. Lastly, we’ll show the resulting patterns, discuss the antenna mutual coupling effects, and give our conclusions.

    Antenna Array Pattern Analysis

    Pattern is defined as the directional strength of a radio-frequency signal viewed from the antenna. The pattern of an antenna array is the product of the isotropic array factor and the isolated element pattern. We assume that the pattern of each element is identical and only consider the isotropic array factor. FIGURE 1 shows the coordination of an antenna array. The first element is set as a reference position. The x-axis is the east direction, the y-axis is the north direction, and the z-axis is the up direction. The baseline vector of the ith antenna is given by I-pi and I-r is the unit vector to the satellite.

    I-Fig1
    Figure 1. Antenna array geometry and direction of satellite. Array elements are identified as E#1, E#2, E#3, and E#4.

    The isotropic array factor is given by

    I-Eq1   (1)

    where λ is wavelength, and Ai is a complex constant. Currently, we only implement a four-element-array CRPA software receiver in real time. Hence, we analyze two kinds of layout of half-wavelength four-element arrays: a symmetrical Y array and a square array. Each antenna is separated from its nearest neighbor by a half wavelength. FIGURE 2 shows photos of the two layouts. FIGURE 3 shows the physical layouts.

    I-Fig2
    Figure 2. Photos of antenna arrays (left: Y array; right: square array).
    I-Fig3top
    Figure 3A. Physical layout of antenna arrays (Y array).
    I-Fig3bottom
    Figure 3B. Physical layout of antenna arrays (square array).

    The antenna patterns towards an elevation angle of 90 degrees, computed using equation 1 and the design layouts, are shown in FIGURE 4. One of the key characteristics of a pattern is the beamwidth, which is defined as the angle with 3-dB loss. FIGURE 5 shows the patterns in elevation angle where the beamwidth of the Y layout is 74 degrees and 86 degrees for the square layout. A narrow beamwidth will benefit anti-interference performance particularly if the interference is close to the direction of a target satellite.

    I-Fig4
    Figure 4. Patterns of antenna arrays (left: Y array; right: square array).
    FIGURE 5 Pattern beamwidths of Y and square arrays (3 dB beamwidth shown).
    Figure 5. Pattern beamwidths of Y and square arrays (3 dB beamwidth shown).
    Specifications of COTS Antennas

    Typically, the COTS antenna selection is determined by high gain and great out-of-band rejection. TABLE 1 shows the specifications of the three antenna models used in this article. These antennas are all patch antennas. The antennas are equipped with surface-acoustic-wave filters for rejecting out-of-band signals. A three-stage low noise amplifier with over 30 dB gain is also embedded in each antenna.

    I-T1
    Table 1. Specifications of COTS antennas used.
    Signal Collection Hardware and Experimental Setup

    The hardware used to collect the antenna array datasets is shown in FIGURE 6 with block-diagram representation in FIGURE 7. The hardware includes a four-element antenna array, four USRP2 software radio systems and one host computer. The signal received from the COTS antenna passes to a USRP2 board equipped with a 800–2300 MHz DBSRX2 programmable mixing and down-conversion daughterboard. The individual USRP2 boards are synchronized by a 10-MHz external common clock generator and a pulse-per-second (PPS) signal. The USRP2s are controlled by the host computer running the Ubuntu distribution of Linux. The open-source GNU Radio software-defined radio block is used to configure USRP2s and collect datasets. All USRP2s are configured to collect the L1 (1575.42 MHz) signal. The signals are converted to near zero intermediate frequency (IF) and digitized to 14-bit complex outputs (I and Q).

    I-Fig7
    Figure 6. Photo of the signal collection hardware.
    I-Fig6
    Figure 7. Block diagram of the signal collection hardware.

    The sampling rate is set as 4 MHz. The host computer uses two solid state drives for storing data sets. For our study, a 64-megabytes per second data transfer rate is needed. The fast solid state drives are especially useful when using high bandwidth signals such as L5, which will require an even higher data streaming rate (80 megabytes per second per channel).

    To compare the physical and electrical layouts of the antenna arrays, we set up the signal collection hardware to record six data sets for the two layouts and the three antenna models as shown in TABLE 2. All of the data sets were five minutes long to obtain enough carrier-phase measurements for positioning.

    I-T2
    Table 2. Experimental setups.
    Logging Carrier-Phase Measurements

    To calculate the precise spacing between the antenna elements, hundreds of seconds of carrier-phase measurements from each element are needed. The collected data sets were provided by our in-house-developed CRPA software receiver. The receiver was developed using Visual Studio under Windows. Most of source code is programmed using C++. Assembly language is used to program the functions with high computational complexity such as correlation operations. The software architecture of the receiver is depicted in FIGURE 8. This architecture exploits four sets of 12 tracking channels in parallel to process each IF signal from an antenna element. Each channel is dedicated to tracking the signal of a single satellite. The tracking channels output carrier-phase measurements to build the steering vectors for each satellite. The Minimum Variance Distortionless Response (MVDR) algorithm was adopted for adaptively calculating the weights for beamforming. Here, there are 12 weight sets, one for each satellite in a tracking channel, for the desired directions of satellites.

    Figure 8. Block diagram of the software architecture.
    Figure 8. Block diagram of the software architecture.

    Using the pre-correlation beamforming approach, the weights are multiplied with IF data and summed over all elements to form 12 composite signals. These signals are then processed by composite tracking channels. Finally, positioning is performed if pseudoranges and navigation messages are obtained from these channels. FIGURE 9 is the graphical user interface (GUI) of the CRPA software receiver. It consists of the channel status of all channels, carrier-phase differences, positioning results, an east-north (EN) plot, a sky plot, a carrier-to-noise-density (C/N0) plot and the gain patterns of the array for each tracked satellite. In the figure, the CRPA software receiver is tracking 10 satellites and its positioning history is shown in the EN plot. The beamforming channels have about 6 dB more gain in C/N0 than the channels of a single element. In each pattern, the direction with highest gain corresponds to the direction of the satellite. While the CRPA software receiver is running, the carrier-phase measurements of all elements and the azimuth and elevation angle of the satellites are logged every 100 milliseconds. Each data set in Table 2 was processed by the software receiver to log the data.

    FIGURE 9 Screenshot of the controlled-reception-pattern-antenna software-receiver graphical user interface.
    Figure 9. Screenshot of the controlled-reception-pattern-antenna software-receiver graphical user interface.
    Electrical Layout of Antenna Array – Procedure

    The procedure of calculating the electrical layout of an antenna array is depicted in FIGURE 10. The single-difference integrated carrier phase (ICP) between the signals of an element, i, and a reference element, j, is represented as:

    I-Eq2   (2)

    where rkij is differential range toward the kth satellite between the ith and jth antenna elements (a function of the baseline vector between the ith and jth elements), δLij is the cable-length difference between the ith and jth antenna elements, Nkij is the integer associated with Φkij , εkij and  is the phase error. The double-difference ICP between the kth satellite and reference satellite l is represented as:
    I-Eq3   (3)

    The cable-length difference term is subtracted in the double difference. Since the distances between the antenna elements are close to one wavelength, equation (3) can be written as:
    I-Eq4   (4)

    where i-rk is the unit vector to satellite k, pij is the baseline vector between the ith and jth elements. By combining all the double-difference measurements of the ijth pair of elements, the observations equation can be represented as:
    I-Eq5      (5)

    From the positioning results of composite channels, the azimuth and elevation angle of satellites are used to manipulate matrix G. To solve equation (5), the LAMBDA method was adopted to give the integer vector N. Then, pij  is solved by substituting N into equation (5). Finally, the cable-length differences are obtained by substituting the solutions of N and pij into equation (2).

    This approach averages the array pattern across all satellite measurements observed during the calibration period.

    FIGURE 10 Procedure for calculating antenna-array electrical spacing.
    Figure 10. Procedure for calculating antenna-array electrical spacing.
    Electrical Layout of Antenna Array – Results

    Using the procedure in the previous section, all electrical layouts of the antenna array were calculated and are shown in FIGURES 11 and 12. We aligned the vectors from element #1 to element #2 for all layouts. TABLE 3 lists the total differences between the physical and electrical layouts. For the same model of antenna, the Y layout has less difference than the square layout. And, in terms of antenna model, antenna #1 has the least difference for both Y and square layouts. We could conclude that the mutual coupling effect of the Y layout is less than that of the square layout, and that antenna #1 has the smallest mutual coupling effect among all three models of antenna for these particular elements and observations utilized.

    FIGURE 11 Results of electrical layout using three models of antenna compared to the physical layout for the Y array.
    Figure 11. Results of electrical layout using three models of antenna compared to the physical layout for the Y array.
    I-Fig12
    Figure 12. Results of electrical layout using three models of antenna compared to physical layout for the square array.
    Table 3. Total differences between physical and electrical layouts.
    Table 3. Total differences between physical and electrical layouts.

    To compare the patterns of all calculated electrical layouts, we selected two specific directions: an elevation angle of 90 degrees and a target satellite, WAAS GEO PRN138, which was available for all data sets. The results are shown in FIGURES 13 and 14, respectively. From Figure 13, the beamwidth of the Y layout is narrower than that of the square layout for all antenna models. When compared to Figure 5, this result confirms the validity of our analysis approach. But, in Figure 14, a strong sidelobe appears at azimuth -60º in the pattern of Y layout for antenna #2. If there is some interference located in this direction, the anti-interference performance of the array will be limited. This is due to a high mutual coupling effect of antenna #2 and only can be seen after calculating the electrical layout.

    I-Fig13
    Figure 13. Patterns of three models of antenna and two layouts toward an elevation angle of 90 degrees.
    I-Fig14
    Figure 14. Patterns of three models of antenna and two layouts toward the WAAS GEO satellite PRN138.
    Conclusions

    The results of our electrical layout experiment show that the Y layout has a smaller difference with respect to the physical layout than the square layout. That implies that the elements of the Y layout have less mutual coupling. For the antenna selection, arrays based on antenna model #1 showed the least difference between electrical and physical layout. And its pattern does not have a high grating lobe in a direction other than to the target satellite.

    The hardware and methods used in this article can serve as a testing tool for any antenna array. Specifically, our methodology, which can be used to collect data, compare physical and electrical layouts, and assess resultant antenna gain patterns, allows us to compare the performances of different options and select the best antenna and layout combination. Results can be used to model mutual coupling and the overall effect of layout and antenna type on array gain pattern and overall CRPA capabilities. This procedure is especially important when using COTS antennas to assemble an antenna array and as we increase the number of antenna elements and the geometry possibilities of the array.

    Acknowledgments

    The authors gratefully acknowledge the work of Dr. Jiwon Seo in building the signal collection hardware. The authors also gratefully acknowledge the Federal Aviation Administration Cooperative Research and Development Agreement 08-G-007 for supporting this research. This article is based on the paper “A Study of Geometry and Commercial Off-The-Shelf (COTS) Antennas for Controlled Reception Pattern Antenna (CRPA) Arrays” presented at ION GNSS 2012, the 25th International Technical Meeting of the Satellite Division of The Institute of Navigation, held in Nashville, Tennessee, September 17–21, 2012.

    Manufacturers

    The antennas used to construct the arrays are Wi-Sys Communications Inc., now PCTEL, Inc. models WS3978 and WS3997 and PCTEL, Inc. model 3978D-HR. The equipment used to collect data sets includes Ettus Research LLC model USRP2 software-defined radios and associated DBSRX2 daughterboards.


    Yu-Hsuan Chen is a postdoctoral scholar in the GNSS Research Laboratory at Stanford University, Stanford, California.

    Sherman Lo is a senior research engineer at the Stanford GNSS Research Laboratory.

    Dennis M. Akos is an associate professor with the Aerospace Engineering Science Department in the University of Colorado at Boulder with visiting appointments at Luleå Technical University, Sweden, and Stanford University.

    David S. De Lorenzo is a principal research engineer at Polaris Wireless, Mountain View, California, and a consulting research associate to the Stanford GNSS Research Laboratory.

    Per Enge is a professor of aeronautics and astronautics at Stanford University, where he is the Kleiner-Perkins Professor in the School of Engineering. He directs the GNSS Research Laboratory.

    FURTHER READING

    • Authors’ Publications

    “A Study of Geometry and Commercial Off-The-Shelf (COTS) Antennas for Controlled Reception Pattern Antenna (CRPA) Arrays” by Y.-H. Chen in Proceedings of ION GNSS 2012, the 25th International Technical Meeting of The Institute of Navigation, Nashville, Tennessee, September 17–21, 2012, pp. 907–914 (ION Student Paper Award winner).

    “A Real-Time Capable Software-Defined Receiver Using GPU for Adaptive Anti-Jam GPS Sensors” by J. Seo, Y.-H. Chen, D.S. De Lorenzo, S. Lo, P. Enge, D. Akos, and J. Lee in Sensors, Vol. 11, No. 9, 2011, pp. 8966–8991, doi: 10.3390/s110908966.

    “Real-Time Software Receiver for GPS Controlled Reception Pattern Array Processing” by Y.-H. Chen, D.S. De Lorenzo, J. Seo, S. Lo, J.-C. Juang, P. Enge, and D.M. Akos in Proceedings of ION GNSS 2010, the 23rd International Technical Meeting of The Institute of Navigation, Portland, Oregon, September 21–24, 2010, pp. 1932–1941.

    “A GNSS Software Receiver Approach for the Processing of Intermittent Data” by Y.-H. Chen and J.-C. Juang in Proceedings of ION GNSS 2007, the 20th International Technical Meeting of The Institute of Navigation, Fort Worth, Texas, September 25–28, 2007, pp. 2772–2777.

    • Controlled-Reception-Pattern Antenna Arrays

    “Anti-Jam Protection by Antenna: Conception, Realization, Evaluation of a Seven-Element GNSS CRPA” by F. Leveau, S. Boucher, E. Goron, and H. Lattard in GPS World, Vol. 24, No. 2, February 2013, pp. 30–33.

    “Development of Robust Safety-of-Life Navigation Receivers” by M.V.T. Heckler, M. Cuntz, A. Konovaltsev, L.A. Greda, A. Dreher, and M. Meurer in IEEE Transactions on Microwave Theory and Techniques, Vol. 59, No. 4, April 2011, pp. 998–1005, doi: 10.1109/TMTT.2010.2103090.

    Phased Array Antennas, 2nd Edition, by R. C. Hansen, published by John Wiley & Sons, Inc., Hoboken, New Jersey, 2009.

    • Antenna Principles

    “Selecting the Right GNSS Antenna” by G. Ryley in GPS World, Vol. 24, No. 2, February 2013, pp. 40–41 (in PDF of 2013 Antenna Survey.)

    GNSS Antennas: An Introduction to Bandwidth, Gain Pattern, Polarization, and All That” by G.J.K. Moernaut and D. Orban in GPS World, Vol. 20, No. 2, February 2009, pp. 42–48.

    A Primer on GPS Antennas” by R.B. Langley in GPS World, Vol. 9, No. 7, July 1998, pp. 50-54.

    • Software-Defined Radios for GNSS

    “A USRP2-based Reconfigurable Multi-constellation Multi-frequency GNSS Software Receiver Front End” by S. Peng and Y. Morton in GPS Solutions, Vol. 17, No. 1, January 2013, pp. 89-102.

    Software GNSS Receiver: An Answer for Precise Positioning Research” by T. Pany, N. Falk, B. Riedl, T. Hartmann, G. Stangl, and C. Stöber in GPS World, Vol. 23, No. 9, September 2012, pp. 60–66.

    Simulating GPS Signals: It Doesn’t Have to Be Expensive” by A. Brown, J. Redd, and M.-A. Hutton in GPS World, Vol. 23, No. 5, May 2012, pp. 44–50.

    Digital Satellite Navigation and Geophysics: A Practical Guide with GNSS Signal Simulator and Receiver Laboratory by I.G. Petrovski and T. Tsujii with foreword by R.B. Langley, published by Cambridge University Press, Cambridge, U.K., 2012.

    “A Real-Time Software Receiver for the GPS and Galileo L1 Signals” by B.M. Ledvina, M.L. Psiaki, T.E. Humphreys, S.P. Powell, and P.M. Kintner, Jr. in Proceedings of ION GNSS 2006, the 19th International Technical Meeting of The Institute of Navigation, Fort Worth, Texas, September 26–29, 2006, pp. 2321–2333.

  • Expert Advice: Location by Database

    Expert Advice: Location by Database

    Tarun Bhattacharrya, Hassan El-Sallabi, Jian Zhu, Jeff Wu, and Per Enge.

    Radio-Frequency Pattern Matching

    By Tarun Bhattacharrya, Hassan El-Sallabi, Jian Zhu, Jeff Wu, and Per Enge

    Radio-frequency pattern matching (RFPM) is the engine that enables the use of mobile-phone signals to locate wireless devices in any environment, including dense downtown areas and indoors. This exciting technology leverages the power of the database to improve location accuracy to within 50 meters in even the toughest signal environments. Significant advances in RFPM technology have been made over the last 10 years. The system described here is deployed in more than 24 wireless networks to provide the location of E-911 callers and help save lives. For simplicity, we focus on the RFPM using signal strengths even though the technology also works with arrival times, signal-to-noise ratios, differential signal strengths and any signal parameter that varies in a predictable fashion over the coverage area.

    Like GPS, RFPM is based on correlation. However, it does not correlate a received spread-spectrum code with a replica code stored in the receiver. Rather, it correlates the signal strength of cell-phone signals measured by the roving phone to a database that contains a map of those signal strengths for the covered area. Consider Figure 1. It shows this key correlation operation. As shown, the database contains a k-vector for each location within the covered area, where the k elements give the estimated strength for the k mobile phone signals that can be received at the given grid point. These k-vectors are typically stored over a 10- or 30-meter grid. This grid of predicted signal strengths is built in advance and is updated only when the topography of the wireless network changes. Thankfully, base stations do not generally move!

    Figure 1. Radio-frequency pattern matching of n-vector from mobile user to k-vectors within database.

    The mobile phone provides the network measurement report (NMR) in real time. This report does not require any network hardware or on-phone software beyond that required by the 2G, 3G and LTE standards for all mobile phones. Thus, the Polaris Wireless solution is capable of locating any mobile phone over any air interface. The NMR is also shown in Figure 1. It contains an n-vector of received signal strengths, where k ≥ n. A multiplicity of n-vectors are backhauled to the server that contains the database. They are correlated with the k-vectors, and the estimated location of the mobile phone is the location associated with the maximum correlation.

    For Example, San Francisco

    Figures 2, 3, and 4 explode the RFPM database for the financial district of San Francisco. Figure 2 is the top view, and the Bay Bridge is shown heading northwest across the Bay. The numbered black dots are some of the base stations in action for this area. Figure 3 digs down one level. It shows the individual k-vectors contained within the database. As shown, this database is based on a 30-meter grid. Figure 4 is a super-zoom that explodes the individual k-vectors. As shown, each of these vectors contains an element for each base station that can be received at the given location. In Figure 4, each element is color coded to correspond to the strength for the signal from the given base station.

    Figure 2. Coverage area of an RFPM database within San Francisco.
    Figure 3. Zoomed view of San Francisco database showing a multiplicity of k-tuples.
    Figure 4. Radio-frequency pattern matching of n-vector from mobile user to k-vectors within database.

    Building the Database

    RFPM accuracy depends strongly on the quality of the database, which needs to be built with great care. In fact, signal propagation depends on the network topology including:
    ◾    antenna location, heights, patterns, effective radiated power, tilt, and azimuth
    ◾    cell type, such as micro-cell, macro-cell, indoor or distributed antenna systems.

    Signal propagation also depends on information available     from geographical information systems such as:
    ◾    tree canopy
    ◾    height of buildings and terrain
    ◾    topography (water, open area, suburban, urban)
    ◾    roads.

    With this data, the signal strength radiating from a base station can be estimated. This is not a simple business. For example, the calculation must identify the points where terrain or buildings interrupts the ray from the transmitter to the receiver. It must also identify the points where these obstacles break the Fresnel zone that surrounds the ray.

    Finally, these open-loop predictions are tuned based on a sparse set of measurements. Once tuned, the database is time invariant or nearly so. If minor changes are made to the network topography, the open loop predictions alone are sufficient to accommodate the changes. If network changes are significant, such as the building of many new base stations, then the open-loop predictions must be updated, and a new set of measurements used to tune the predictions.

    Figure 5 shows a typical map of signal strengths surrounding one mobile phone in a completely open area. Absent terrain and buildings, the signal strengths vary rather smoothly. Figure 6 is for one of the transmitters in the San Francisco financial district, which is a much more complicated urban environment due to the dense concentration of high-rise buildings and uneven terrain. In this case, the signal-strength signature has a gratifying abundance of detail. This detail enables RFPM to work very well in the complicated signal environments that we find in downtown areas and also indoors. In short, RFPM benefits from the buildings and terrain that hinder satellite measurements.

    Figure 5. Predicted signal strength for a transmitter surrounded by open ground.
    Figure 6. Predicted signal strength for one transmitter in the San Francisco financial district.

    Performance and Summary

    RFPM works well. It provides high accuracy in a in a wide variety of environments. Polaris Wireless routinely tests the accuracy of its solution in urban settings. Table 1 shows the results of such evaluations, based on measurement sets that are not used to tune the database.

    Table 1. Evaluations based on routine accuracy tests of RFPM in urban settings.

    These days, robust navigation for downtown and indoors is based on an expanding suite of location technologies. These include: assisted GPS, new satellite constellations (Galileo, GLONASS, Compass, and so on), inertial measurements, Wi-Fi ranging, and signals from low-Earth orbit. RFPM, and its unique reliance on database-derived location, should remain an important part of this mix.


    Tarun Bhattacharrya is vice president of research at Polaris Wireless. He earned his Ph.D. in electrical engineering from the Indian Institute of Science.

    Hassan El-Sallabi received his D.Sc. in electrical and communications engineering from Helsinki University of Technology, Finland. At Polaris he works on RF propagation modeling.

    Jian (JET) Zhu received his Ph.D. in electrical engineering from Georgia Institute of Technology; he is a research engineer at Polaris.

    Jeff Wu focuses on algorithm development for propagation modeling at Polaris, and is a Ph.D. candidate in electrical engineering at Stanford.

    Per Enge is the Kleiner Perkins professor of engineering at Stanford University, where he directs the Stanford Center for Position, Navigation, and Time. He is also a technical advisor to Polaris Wireless.

  • Enge Pilots Polaris Wireless/GPS Integration

    Polaris Wireless, a provider of software-based wireless location solutions, announced that Per K. Enge of Stanford University has joined its executive team as the chief technical advisor. In an interview with GPS World, Enge describes how combined wireless location signatures and GPS show an exciting way forward for location in dense urban environments, where the wireless solution actually improves as GPS degrades.

    Polaris Wireless, a provider of software-based wireless location solutions, based in Mountain View, California, announced that Per K. Enge has joined its executive team as the chief technical advisor to the company’s CEO, Manlio Allegra.

    Enge is a professor of aeronautics and astronautics at Stanford University, where he directs the GPS Research Laboratory. He is also is co-author of the textbook Global Positioning System: Signals, Measurements, and Performance and has received the Kepler, Thurlow, and Burka Awards from the Institute of Navigation (ION) for his work. He received his Ph.D. from the University of Illinois in 1983, where he designed a direct-sequence multiple-access communication system that provided an orthogonal signal set to each user.

    Polaris Wireless’s core product is a Wireless Location Signature (WLS) for wireless network operators to identify the location of a wireless device to within 50 meters, at low cost, even in urban and indoor environments.  Polaris WLS is designed to scale as an operator’s network grows, providing location capabilities in 2G (GSM/CDMA), 3G (UMTS/CDMA 2000), and emerging 4G (LTE) air interfaces, as well as indoor technologies such as W-iFi, DAS, and Femtocells.  Founded in 2003, Polaris also counts among its customers law enforcement/government agencies and location-based application companies.

    In an interview with GPS World, Enge described how the Polaris location solution is actually enhanced by buildings, in direct opposition to the way GPS operates, and how he sees integration of the two technologies as leading the way forward for location-based services.

    “A cell phone general receives signals from 3 to 10 base stations, and is required to report the signal to noise ratio of those, even beyond base station used.  We backhaul that information into the Polaris server, where we do multilateration based on signal strength measurements. Imagery shows how RF signal stregnth varies in the city. Polaris uses the base station signals (similar to how Skyhook using Wifi access points), but the cell infrastructure use is better because it is always there and used by every phone, while Wifi receiver is not on all phones yet, just smartphones.

    “Marriage to GPS is very cool, because the Polaris technology loves buildings, they introduce the variation, the grain in the signal-strength map. It’s jagged, it’s got nooks and crannies, which are well-predicted by our propagation models. So those buildings enhance our technology.

    “We’ve got it pretty widely deployed, on 24 E911 networks, 24 wireless carriers that we have commercial relationships with in the United States, from quite big to more boutique providers.  Some are smaller, regional, other larger, like Verizon though not the total network.

    Michael Doherty, Polaris Wireless product marketing and communications director, added that “We also work outside the United States for location surveillance for governments who are also looking to deploy for commercial services.  The solution  is software-based, in the control plane of the handset, it is already within the handset, using network measurement reports (NMRs). Emergency caller location for cell phones, E911, is just the first of applications to come, it’s where we cut our teeth.”

    Polaris Wireless has been a research partner and sponsor at the Stanford GPS Research Laboratory that Enge directs. “Now I’ve stepped up my relationship with them, to start talking about getting the next level of accuracy.  They joined about 4 or 5 years ago, and funded the work of David De Lorenzo , a Ph.D. student who did first big set of measurements that we did on indoor navigation, on GPS penetration into buildings. On campus.  We’re aviation people, so this was a new area, good for us, we expanded into that area and got a lot smarter about how to navigate indoors.”

    De Lorenzo’s 2009 Institute of Navigation International Technical Meeting paper, “Design and Performance of a Minimum-Variance Hybrid Location Algorithm Utilizing GPS and Cellular Received Signal Strength for Positioning in Dense Urban Environments,” can be accessed here.

    “In New York, three-quarters of cell calls get 50-meter accuracy,” stated Doherty. “99 percent get 150-meter. Polaris solutions will improve that by a factor of 2. That’s a forward-looking  statement.”

    “We are entering a significant growth period for the company. From 2007–08 we saw our business start to take off, from core E911, take off in demand from overseas governments, for anti-crime, anti-terrorism. From 2007 to 2010, Polaris Wireless had an average 66 per increase per year in bookings, in contracts signed. That’s a healthy indicator of future revenues

    “Year over year 2010 to 2011, we saw a quadrupling of bookings, and  a revenue increase of 30 percent, as we started to see those deployments taking place, and revenues coming in as a result. we expect that trajectory to continue through 2012 and 2013. Again, forward-looking statements.”

    “During the deep recession, when U.S. governments were cutting back, we went overseas and were able to ride that wave.”

    As to particulars on directions that Polaris research and development might take under his advisement, Enge volunteered “One of the really big initiatives is the inclusion of map-matching, knowing that we’re probably on a street, especially if our velocity is appreciable, going beyond a snapshot solution fix, taking advantage of knowing where streets our.  We have cached that as a particle-filter problem, a Kalman filter associated with a set of hypotheses. When you have a sequence of snapshots, the particle filter is a very powerful way of putting that together.

    “How do we know if we’re in a building? By looking at the NMRs. Now, how high are we in that builidng, so that’s aided by the database.

    “The big core enterprise of Polaris, stripping away PNT, the strength is the underlying [cellular network] database.  Polaris is where the database meets navigation.  To make E911, law enformcement, navigation, all better.

    “One application I’m particularly keen on is for campuses.The risk for muggings is greatest when walking to a car in a distant parking lot. We’re thinking about how we can help that student when assaulted. You can’t pull out your phone and call E911 nor run to a blue tower [emergency phone].  But you might have your keys in your hand, so if we make it so you can just push a button on it, the NMR goes to an emergency center on campus, then we can be more helpful in helping that student.”

    Doherty inserted, “The beauty of our technology is that no one has to do anything. Our technology is capable of locating any device connected to the cellular network, as long as the device is powered on. It’s a software deployment, nothing to go out and deploy equipment cell tower by cell tower.  Nor can bad guys disable anything in the phone to void this service. Polaris’s solution resides in the network, not in the handset.  The confluence of signals in the network drives the solution.  There is no chip in the handset that can be disabled or turned off.  Whether GPS is enabled or not.”

    In a final comment on concerns for the future, Enge said, “Personal privacy devices, PPDs, the GPS jammers. I’m pretty worried about that. Unlike the Polaris technology, GPS technology can be slapped on someone’s car without an organization making sure the rule of law is preserved.  We’re hoping to make the PW technology complementary to GPS, so that when GPS is jammed, the PW solution will still provide a location.

    Doherty added, “It is complementary today. Cellular networks covers a large area. Take for example, downtown San Francisco: The wireless signature solution is absolutely the best to locate there. In wide-open country, GPS has an advantage. That complementarity is happening today.

    “Having Per onboard, we’re going to make that complementarity work better and better.”

  • Innovation: Digging into GPS Integrity

    Innovation: Digging into GPS Integrity

    Charting the Evolution of Signal-in-Space Performance by Data Mining 400,000,000 Navigation Messages

    By Liang Heng, Grace Xingxin Gao, Todd Walter, and Per Enge

    There are four important requirements of any navigation system: accuracy, availability, continuity, and integrity. In this month’s column we take a look at one particular aspect of GPS integrity: that of the signal in space and find out how trustworthy is the satellite ephemeris and clock information in the broadcast navigation message.

    GPS World photo
    INNOVATION INSIGHTS by Richard Langley

    BUT THE GREATEST OF THESE IS INTEGRITY. There are four important requirements of any navigation system: accuracy, availability, continuity, and integrity.

    Perhaps the most obvious navigation system requirement, accuracy describes how well a measured value agrees with a reference value, typically the true value. In the case of GPS, we might talk about the accuracy of a range measurement. A receiver actually measures a pseudorange — a biased and noisy measure of the geometric range between the receiver and the satellite. After correcting for satellite ephemeris and satellite clock errors (the primary so-called signal-in-space errors), receiver clock errors, and atmospheric effects, we can get an estimate of the geometric range. How well we account for these errors or biases, will determine the accuracy of the corrected pseudorange measurement and ultimately, the accuracy of a derived position.

    A navigation system’s availability refers to its ability to provide the required function and performance within the specified coverage area at the start of an intended operation. In many cases, system availability implies signal availability, which is expressed as the percentage of time that the system’s transmitted signals are accessible for use. In addition to transmitter capability, environmental factors such as signal attenuation or blockage or the presence of interfering signals might affect availability.

    Ideally, any navigation system should be continuously available to users. But, because of scheduled maintenance or unpredictable outages, a particular system may be unavailable at a certain time. Continuity, accordingly, is the ability of a navigation system to function without interruption during an intended period of operation. More specifically, it indicates the probability that the system will maintain its specified performance level for the duration of an operation, presuming system availability at the beginning of that process.

    The integrity of a navigation system refers to its trustworthiness. A system might be available at the start of an operation, and we might predict its continuity at an advertised accuracy during the operation.

    But what if something unexpectedly goes wrong? If some system anomaly results in unacceptable navigation accuracy, the system should detect this and warn the user. Integrity characterizes a navigation system’s ability to provide this timely warning when it fails to meet its stated accuracy. If it does not, we have an integrity failure and the possibility of conveying hazardously misleading information. GPS has built into it various checks and balances to ensure a fairly high level of integrity. However, GPS integrity failures have occasionally occurred.

    In this month’s column we take a look at one particular aspect of GPS integrity: that of the signal in space and find out how trustworthy is the satellite ephemeris and clock information in the broadcast navigation message.


    The Navstar Global Positioning System is so far the most widely used space-based positioning, navigation, and timing system. GPS works on the principle of trilateration, in which the measured distances from a user receiver to at least four GPS satellites in view, as well as the position and clock data for these satellites, are the prerequisites for the user receiver to fix its exact position. For most GPS Standard Positioning Service (SPS) users, real-time satellite positions and clocks are derived from ephemeris parameters and clock correction terms in navigation messages broadcast by GPS satellites. The GPS Control Segment routinely generates navigation message data on the basis of a prediction model and the measurements at more than a dozen monitor stations. The differences between the broadcast ephemerides/clocks and the truth account for signal-in-space (SIS) errors. SIS errors are usually undetectable and uncorrectable for stand-alone SPS users, and hence directly affect the positioning accuracy and integrity. Nominally, SPS users can assume that each broadcast navigation message is reliable and the user range error (URE) derived from a healthy SIS is at the meter level or even sub-meter level. In practice, unfortunately, SIS anomalies have happened occasionally and UREs of tens of meters or even more have been observed, which can result in an SPS receiver outputting a hazardously misleading position solution. Receiver autonomous integrity monitoring (RAIM) or advanced RAIM is a promising tool to protect stand-alone users from such hazards; however, most RAIM algorithms assume at most one satellite fault at a time. Knowledge about the SIS anomalies in history is very important not only for assessing the GPS SIS integrity performance but also for validating the fundamental assumption of RAIM.

    A typical method for calculating SIS UREs is to compare the broadcast ephemerides/clocks with the precise, post-processed ones. Although this method is very effective in assessing the GPS SIS accuracy performance, few attempts have been made to use it to assess the GPS SIS integrity performance because broadcast ephemeris/clock data obtained from a global tracking network sometimes contain errors caused by receivers or data conversion processes and these errors usually result in false SIS anomalies. In this article, we introduce a systematic methodology to cope with this problem and screen out all the potential SIS anomalies in the past decade from when Selective Availability (SA) was turned off.

    GPS SIS Integrity

    The integrity of a navigation system refers — just as it does to a person — to its honesty, veracity, and trustworthiness. In the case of GPS, this includes the integrity of the ephemeris and clock data in the broadcast navigation messages. We refer to this as signal-in-space integrity.

    GPS SIS URE. As indicated by the name, GPS SIS URE is the pseudorange modeling inaccuracy due to operations of the GPS ground control and the space vehicles. Specifically, SIS URE includes satellite ephemeris and clock errors, satellite antenna performance variations, and signal imperfections, but not ionospheric or tropospheric delay, multipath, or any errors due to user receivers. SIS URE is dominated by ephemeris and clock errors because antenna variations and signal imperfections are at a level of millimeters or centimeters.

    In broadcast navigation messages, there is a parameter called user range accuracy (URA) that is intended to be a conservative representation of the standard deviation (1-sigma) of the URE at the worst-case location on the Earth. For example, a URA index value of 0 means that the 1-sigma URE is expected to be less than 2.4 meters, and a URA index value of 1 means that the 1-sigma URE is expected to be greater than 2.4 meters but less than 3.4 meters, and so on. In the past several years, most GPS satellites have a URA index value of 0. A nominal URA value, in meters, can be computed as X = 2(1+N/2), where N is the index value, for index values of 6 or less. For 6 < N < 15, X = 2(N-2).

    GPS SPS SIS Integrity. In the SPS Performance Standard (PS), as well as the latest version of the Interface Specification (IS-GPS-200E), the GPS SPS SIS URE integrity standard assures that for any healthy SIS, there is an up-to-10−5 probability over any hour of the URE exceeding the not-to-exceed (NTE) tolerance without a timely alert during normal operation. The NTE tolerance is currently defined to be 4.42 times the upper bound (UB) on the URA value broadcast by the satellite. Before September 2008, the NTE tolerance was defined differently, as the maximum of 30 meters and 4.42 times URA UB. The reason for the “magic” number 4.42 here is the Gaussian assumption of the URE, although this assumption may be questionable. (4.42 sigma corresponds to a probability level of 99.999 percent (1 – 10–5)).

    In this article, a GPS SPS SIS anomaly is defined as a threat of an SIS integrity failure; that is, a condition during which an SPS SIS marked healthy results in a URE exceeding the NTE tolerance. Because the definition of the NTE tolerance is different before and after September 2008, we consider both of the two NTE tolerances for the sake of completeness and consistency.

    Methodology

    The SIS anomalies are screened out by comparing broadcast ephemerides/clocks with precise ones. As shown in Figure 1, the whole process consists of three steps: data collecting, data cleansing, and anomaly screening.

    Inn-Fig1 Source: Richard Langley
    Figure 1. Framework of the whole process. XYZB values refer to the coordinates of satellite position and satellite clock bias.

    In the first step, the navigation message data files are downloaded from the International GNSS Service (IGS). In addition, two different kinds of precise ephemeris/clock data are downloaded from IGS and the National Geospatial-Intelligence Agency (NGA), respectively. The details about these data sources will be discussed in the next section.

    Since each GPS satellite can be observed by many IGS stations at any instant, each navigation message is recorded redundantly. In the second step, a data-cleansing algorithm exploits the redundancy to remove the errors caused on the ground. This step distinguishes our work from that of most other researchers because the false anomalies due to corrupted data can be mostly precluded.

    The last step is computing worst-case SIS UREs as well as determining potential SIS anomalies. The validated navigation messages prepared in the second step are used to propagate broadcast orbits/clocks at 15-minute intervals that coincide with the precise ones. A potential SIS anomaly is claimed when the navigation message is healthy and in its fit interval with the worst-case SIS URE exceeding the SIS URE NTE tolerance.

    Data Sources

    We obtained broadcast navigation message data and precise ephemeris and clock data from publicly available sources.

    Broadcast Navigation Message Data. Broadcast GPS navigation message data files are available at IGS Internet sites. All the data are archived in Receiver Independent Exchange (RINEX) navigation file format, which includes not only the ephemeris/clock parameters broadcast by the satellites but also some information produced by the ground receivers, such as the pseudorandom noise (PRN) signal number and the transmission time of message (TTOM).

    The IGS tracking network is made up of more than 300 volunteer stations all over the world (a map is shown in Table 1) ensuring seamless, redundant data logging. Since broadcast navigation messages are usually updated every two hours, no single station can record all navigation messages. For the ease of users, two IGS archive sites, the Crustal Dynamics Data Information System (CDDIS) and the Scripps Orbit and Permanent Array Center (SOPAC), provide two kinds of ready-to-use daily global combined broadcast navigation message data files, brdcddd0.yyn and autoddd0.yyn, respectively, where ddd is the day of year yy. Unfortunately, these files sometimes contain errors that can cause false anomalies.

     Table 1. Comparison of IGS and NGA precise ephemeris/clock data. Source: Richard Langley
    Table 1. Comparison of IGS and NGA precise ephemeris/clock data.

    Therefore, we devised and implemented a data-cleansing algorithm to generate the daily global combined navigation messages, which are as close as possible to the navigation messages that the satellites actually broadcast, from all available navigation message data files of all IGS stations. The data-cleansing algorithm is based on majority vote, and hence all values in our data are cross validated. Accordingly, we name our daily global combined navigation messages “validated navigation messages,” as shown in Figure 1.

    Precise Ephemeris and Clock Data. Precise GPS ephemerides/clocks are generated by some organizations such as IGS and NGA that routinely post-process observation data. Precise ephemerides/clocks are regarded as “truth” because of their centimeter-level accuracy.

    Table 1 shows a side-by-side comparison between IGS and NGA precise ephemeris/clock data, in which the green- and red-colored text implies pros and cons, respectively. For NGA data, the only con is that the data have been publicly available only since January 4, 2004. As a result, for the broadcast ephemerides/clocks before this date, IGS precise ephemerides/clocks are the only references. Nevertheless, care must be taken when using IGS precise ephemerides/clocks due to the following three issues.

    The first issue with the IGS precise ephemerides/clocks is the relatively high rate of bad/absent data, as shown in the third row of Table 1. For a GPS constellation of 27 healthy satellites, 1.5 percent bad/absent data means no precise ephemerides or clocks for approximately 10 satellite-hours per day. This issue can result in undetected anomalies (false negatives).

    The second issue is that, as shown in the fourth row of Table 1, IGS switched to IGS Time for its precise ephemeris/clock data on 22 February, 2004. The IGS clock is not synchronized to GPS Time, and the differences between the two time references may be as large as 3 meters. Fortunately, the time offsets can be extracted from the IGS clock data files. Moreover, a similar problem is that IGS precise ephemerides use a frame aligned to the International Terrestrial Reference Frame (ITRF) whereas broadcast GPS ephemerides are based on the World Geodetic System 1984 (WGS 84). The differences between ITRF and the versions of WGS 84 used since 1994 are on the order of a few centimeters, and hence a transformation is not considered necessary for the purpose of our work.

    The last, but not the least important, issue with the IGS precise ephemerides is that the data are provided only for the center of mass (CoM) of the satellite. Since the broadcast ephemerides are based on the satellite antenna phase center (APC), the CoM data must be converted to the APC before being used. Both IGS and NGA provide antenna corrections for every GPS satellite. Although the IGS and the NGA CoM data highly agree with each other, the IGS satellite antenna corrections are quite different from the NGA’s, and the differences in z-offsets can be as large as 1.6 meters for some GPS satellites. The reason for these differences is mainly due to the different methods in producing the antenna corrections: the IGS antenna corrections are based on the statistics from more than 10 years of IGS data, whereas the NGA’s are probably from the calibration measurements on the ground. In order to know whose satellite antenna corrections are better, the broadcast orbits for all GPS satellites in 2009 were computed and compared with three different precise ephemerides: IGS CoM + IGS antenna corrections, IGS CoM + NGA antenna corrections, and NGA APC. Generally, the radial ephemeris error is expected to have a zero mean. However, the combination “IGS CoM + IGS antenna corrections” results in radial ephemeris errors with a non-zero mean for more than half of the GPS satellites. Therefore, the NGA antenna corrections were selected to convert the IGS CoM data to the APC.

    Data Cleansing

    Figure 2 shows a scenario of data cleansing. Owing to accidental bad receiver data and various hardware/software bugs, a small proportion of the navigation data files from the IGS stations have defects such as losses, duplications, inconsistencies, discrepancies, and errors. Therefore, more than just removing duplications, the generation of validated navigation messages is actually composed of two complicated steps.

     Figure 2. A scenario of data cleansing: In the figure, the GPS satellite PRN32 started to transmit a new navigation message at 14:00. Receiver 1 had not observed the satellite until 14:36, and hence the TTOM in its record was 14:36. Additionally, Receiver 1 made a one-bit error in ∆n (4.22267589140 × 10-9 11823 × 2−43 π). Receiver 2 perhaps had some problems in its software: the IODC was unreported and both the toc and ∆n were written weirdly. Receiver n used an incorrect ranging code, PRN01, to despread and decode the signal of PRN32; fortunately, all the parameters except TTOM were perfectly recorded. Moreover, the three receivers interpreted URA (SV accuracy) differently. A computer equipped with our data cleansing algorithms is used to process all the data from the receivers. The receiver-caused errors are removed and the original navigation message is recovered. Source: Richard Langley
    Figure 2. A scenario of data cleansing: In the figure, the GPS satellite PRN32 started to transmit a new navigation message at 14:00. Receiver 1 had not observed the satellite until 14:36, and hence the TTOM in its record was 14:36. Additionally, Receiver 1 made a one-bit error in ∆n (4.22267589140 × 10-9 11823 × 2−43 π). Receiver 2 perhaps had some problems in its software: the IODC was unreported and both the toc and ∆n were written weirdly. Receiver n used an incorrect ranging code, PRN01, to despread and decode the signal of PRN32; fortunately, all the parameters except TTOM were perfectly recorded. Moreover, the three receivers interpreted URA (SV accuracy) differently. A computer equipped with our data cleansing algorithms is used to process all the data from the receivers. The receiver-caused errors are removed and the original navigation message is recovered.

    First step. Suppose that we want to generate the validated navigation messages for day n. In the first step, we apply the following operations sequentially to each RINEX navigation data file from day n − 1 to day n + 1:

    1) Parse the RINEX navigation file;

    2) Recover least significant bit (LSB);

    3) Classify URA values;

    4) Remove the navigation messages not on day n;

    5) Remove duplications;

    6) Add all remaining navigation messages into the set O.

    The reason why the data files from day n − 1 to day n + 1 are considered is that a few navigation messages around 00:00 can be included in some data files on day n − 1, and a few navigation messages around 23:59 can be included in some data files on day n + 1. The LSB recovery is used here to cope with the discrepant representations of floating-point numbers in RINEX navigation files. The URA classifier is employed to recognize and unify various representations of URA in the files. The duplication removal is applied because some stations write the same navigation messages repeatedly in one data file, which is unfavorable to the vote in the second step.

    Second Step. At the end of the first step, we have a set O that includes all the navigation messages on day n. The set O still has duplications because a broadcast navigation message can be reported by many IGS stations. However, as shown in Figure 2, duplications of a broadcast navigation message may come with different errors and are not necessarily identical. Several other examples of such problems can be found in our journal paper listed in Further Reading. Fortunately, most orbital and clock parameters are seldom reported incorrectly, and even when errors happen, few stations agree on the same incorrect value. In our work, these parameters are referred to as robust parameters. On the contrary, some parameters, such as TTOM, PRN, URA and issue of data clock (IODC), are more likely to be erroneous and when errors happen, several stations may make the same mistake. These parameters are referred to as fragile parameters. The cause of the fragility is either the physical nature (for example, TTOM, PRN) or the carelessness in hardware/software implementations (for example, URA, IODC).

    Majority vote is applied to all fragile parameters (except TTOM, which is determined by another algorithm described in our journal paper) under the principle that the majority is usually correct. Meanwhile, the robust parameters are utilized to identify the equivalence of two navigation messages — two navigation messages are deemed identical if and only if they agree on all the robust parameters, although their fragile parameters could be different. Therefore, the goal of duplication removal and majority vote is a set P, in which any navigation message must have at least one robust parameter different from any other and has all fragile parameters confirmed by the largest number of stations that report this navigation message.

    After the operations above, we have a set P in which there are no duplicated navigation messages in terms of robust parameters and all fragile parameters are as correct as possible. A few navigation messages in P still have errors in their robust parameters. These unwanted navigation messages feature a small number of reporting stations. Finally, the navigation messages confirmed by only a few stations being discarded and the survivors are the validated broadcast navigation messages, stored in files sugldddm.yyn. For further details of our algorithms, see our journal paper.

    Anomaly Screening

    The validated broadcast navigation messages prepared using the algorithm described in the previous section were employed to propagate broadcast satellite orbits and clocks. For each 15-miniute epoch, t, that coincides with precise ephemerides/clocks, the latest transmitted broadcast ephemeris/clock is chosen to calculate the worst-case SIS URE – the maximum SIS URE that a user on Earth can experience.

    Finally, a potential GPS SIS anomaly is claimed when all of the following conditions are fulfilled.

    • The worst-case SIS URE exceeds the NTE tolerance;
    • The broadcast navigation message is healthy; that is,
      • The RINEX field SV health is 0, and
      • The URA UB ≤ 48 meters;
    • The broadcast navigation message is in its fit interval; that is, ∆t = t − TTOM ≤ 4 hours;
    • The precise ephemeris/clock is available and healthy.

    Results

    A total of 397,044,414 GPS navigation messages collected by an average of 410 IGS stations from June 1, 2000 (one month after turning off SA), to August 31, 2010, have been screened. The NGA APC precise ephemerides/clocks and the IGS CoM precise ephemerides/clocks with the NGA antenna corrections were employed as the truth references. Both old and new NTE tolerances were used for determining anomalies.

    Before interpreting the results, it should be noted that there are some limitations due to the data sources and the anomaly-determination criteria. First, false anomalies may be claimed because there may be some errors in the precise ephemerides/clocks or the validated navigation messages. Second, some short-lived anomalies may not show up if they happen to fall into the 15-minute gaps of the precise ephemerides/clocks. Third, some true anomalies may not be detected if the precise ephemerides/clocks are temporarily missing. The third limitation is especially significant for the results before January 3, 2004, because only the IGS precise ephemerides/clocks are available, which feature a high rate of bad/absent data. (For example, the clock anomaly of Space Vehicle Number (SVN) 23/PRN23 that occurred on January 1, 2004 is missed by our process because the IGS precise clocks for PRN23 on that day were absent.) Last but not least, users might not experience some anomalies because a satellite was not trackable at that time, or the users were notified via a Notice Advisory to Navstar Users (NANU). (A satellite may indicate that it is unhealthy through the use of non-standard code or data. The authors’ future work will include using observation data to verify the potential anomalies found in the results presented here.) Therefore, all the SIS anomalies claimed in this article are considered to be potential and under further investigation.

    Potential SIS Anomalies. A total of 1,256 potential SIS anomalies were screened out under SPS PS 2008 (or 374 potential SIS anomalies under SPS PS 2001). Figure 3 shows all these anomalies in a Year-SVN plot. It can be seen that during the first year after SA was turned off, SIS anomalies occurred frequently for the whole constellation.

     Figure 3. Potential SIS anomalies from June 1, 2000, to August 31, 2010. The horizontal lines depict the periods when the satellites were active (not necessarily healthy). The color of the lines indicates the satellites' block type, as explained by the top left legend. Source: Richard Langley
    Figure 3. Potential SIS anomalies from June 1, 2000, to August 31, 2010. The horizontal lines depict the periods when the satellites were active (not necessarily healthy). The color of the lines indicates the satellites’ block type, as explained by the top left legend.

    Moreover, 2004 is apparently a watershed: before 2004, anomalies occurred for all GPS satellites (except two satellites launched in 2003, SVN45/PRN21 and SVN56/PRN16) whereas after 2004, anomalies occurred much less frequently and more than 10 satellites have never been anomalous. Figure 4 further confirms the improving GPS SIS integrity performance in the past decade, no matter which SPS PS is considered.

     Figure 4. Number of potential SIS anomalies per year. The SIS performance was improved during the past decade. There were 0 anomalies in 2009 according to SPS PS 2001 and this number is represented by 0.1 in the figure. Source: Richard Langley
    Figure 4. Number of potential SIS anomalies per year. The SIS performance was improved during the past decade. There were 0 anomalies in 2009 according to SPS PS 2001 and this number is represented by 0.1 in the figure.

    Therefore, it is possible to list all potential SIS anomalies from January 4, 2004, to August 31, 2010, in a compact table: Table 2. Most anomalies in the table have been confirmed by NANUs and other literature. The table reveals an important and exciting piece of information: never have two or more SIS anomalies occurred simultaneously since 2004. Accordingly, in the sense of historical GPS SIS integrity performance, it is valid for RAIM to assume at most one satellite fault at a time.

     Table 2. List of potential anomalies from January 4, 2004, to August 31, 2010. Source: Richard Langley
    Table 2. List of potential anomalies from January 4, 2004, to August 31, 2010.

    Validated Navigation Messages. For the purpose of comparison and verification, the IGS daily global combined broadcast navigation message data files brdcddd0.yyn and autoddd0.yyn were used to propagate broadcast satellite orbits and clocks as well. The NGA APC precise ephemerides/clocks were employed for the truth references. The SPS PS 2008 NTE tolerance was used for determining anomalies. The other criteria for anomaly screening that are the same as in the previous section were still applied.

    All the potential SIS anomalies for 2006–2009 were found based on the three kinds of daily combined broadcast navigation messages. Table 3 shows a comparison of the total hours of the anomalies per year. It can be seen that brdcddd0.yyn and autoddd0.yyn result in approximately 11 times more false anomalies than true ones. Moreover, all potential anomalies derived from sugldddm.yyn are confirmed by brdcddd0.yyn and autoddd0.yyn, which indicates that our sugldddm.yyn does not introduce any more false anomalies than brdcddd0.yyn and autoddd0.yyn.

     Table 3. Total hours of anomalies per year computed from three different kinds of daily global combined broadcast navigation messages. Source: Richard Langley
    Table 3. Total hours of anomalies per year computed from three different kinds of daily global combined broadcast navigation messages.

    Conclusion

    In this article, the GPS SIS integrity performance in the past decade was assessed by comparing the broadcast ephemerides/clocks with the precise ones. Thirty potential anomalies were found. The fundamental assumption of RAIM is valid based on a review of the GPS SIS integrity performance in the past seven years.

    Acknowledgments

    The authors gratefully acknowledge the support of the Federal Aviation Administration. This article contains the personal comments and beliefs of the authors, and does not necessarily represent the opinion of any other person or organization.

    The authors would like to thank Mr. Tom McHugh, William J. Hughes FAA Technical Center, for his valuable input to the data-cleansing algorithm. This article is based on the paper “GPS Signal-in-Space Integrity Performance Evolution in the Last Decade: Data Mining 400,000,000 Navigation Messages from a Global Network of 400 Receivers” to appear in the Institute of Electrical and Electronics Engineers (IEEE) Transactions on Aerospace and Electronic Systems..


    Liang Heng is a Ph.D. candidate under the guidance of Professor Per Enge in the Department of Electrical Engineering at Stanford University.

    Grace Xingxin Gao is a research associate in the GPS Research Laboratory of Stanford University.

    Todd Walter is a senior research engineer in the Department of Aeronautics and Astronautics at Stanford University.

    Per Enge is a professor of Aeronautics and Astronautics at Stanford University, where he is the Kleiner-Perkins, Mayfield, Sequoia Capital Professor in the School of Engineering. He directs the GPS Research Laboratory, which develops satellite navigation systems based on GPS.


    FURTHER READING

    • Authors’ Research Papers

    “GPS Signal-in-Space Integrity Performance Evolution in the Last Decade: Data Mining 400,000,000 Navigation Messages from a Global Network of 400 Receivers” by L. Heng, G.X. Gao, T. Walter, and P. Enge in Transactions on Aerospace and Electronic Systems, the Institute of Electrical and Electronics Engineers, accepted for publication.

    “GPS Signal-in-Space Anomalies in the Last Decade: Data Mining of 400,000,000 GPS Navigation messages” by L. Heng, G.X. Gao, T. Walter, and P. Enge in Proceedings of ION GNSS 2010, the 23rd International Technical Meeting of The Satellite Division of the Institute of Navigation, Portland, Oregon, September 21–24, 2010, pp. 3115–3122.

    “GPS Ephemeris Error Screening and Results for 2006–2009” by L. Heng, G.X. Gao, T. Walter, and P. Enge in Proceedings of ION ITM 2010, the 2010 International Technical Meeting of the Institute of Navigation, San Diego, California, January 24–26, 2010, pp. 1014–1022.

    • Earlier Work on Assessing GPS Broadcast Ephemerides and Clocks

    “GPS Orbit and Clock Error Distributions” by C. Cohenour and F. van Graas in Navigation, Vol. 58, No. 1, Spring 2011, pp. 17–28.

    “Statistical Characterization of GPS Signal-in-Space Errors” by L. Heng, G.X. Gao, T. Walter, and P. Enge in Proceedings of ION ITM 2011, the 2011 International Technical Meeting of the Institute of Navigation, San Diego, California, January 24–26, 2011, pp. 312–319.

    “Broadcast vs. Precise GPS Ephemerides: A Historical Perspective” by D.L.M. Warren and J.F. Raquet in GPS Solutions, Vol. 7, No. 3, 2003, pp. 151–156, doi: 10.1007/s10291-003-0065-3.

    “Accuracy and Consistency of Broadcast GPS Ephemeris Data” by D.C. Jefferson and Y.E. Bar-Sever in Proceedings of ION GPS-2000, the 13th International Technical Meeting of the Satellite Division of The Institute of Navigation, Salt Lake City, Utah, September 19–22, 2000, pp. 391–395.

    “The GPS Broadcast Orbits: An Accuracy Analysis” by R.B. Langley, H. Jannasch, B. Peeters, and S. Bisnath, presented in Session B2.1-PSD1, New Trends in Space Geodesy at the 33rd COSPAR Scientific Assembly, Warsaw, July 16–23, 2000.

    • Signal-in-Space Anomalies

    “GNSS: The Present Imperfect” by D. Last in Inside GNSS, Vol. 5, No. 3, May 2010, pp. 60–64.

    “Investigation of Upload Anomalies Affecting IIR Satellites in October 2007” by K. Kovach, J. Berg, and V. Lin in Proceedings of ION GNSS 2008, the 21st International Technical Meeting of the Satellite Division of The Institute of Navigation, Savannah, Georgia, September 16–19, 2008, pp. 1679–1687.

    Global Positioning System (GPS) Standard Positioning Service (SPS) Performance Analysis Report No. 58, July 31, 2007, Reporting Period: 1 April – 30 June 2007.

    Discrepancy Report, DR No. 55, “GPS Satellite PRN18 Anomaly Affecting SPS Performance” by N. Vary, FAA William J. Hughes Technical Center, Pomona, New Jersey, April 11, 2007.

    “GPS Receiver Responses to Satellite Anomalies” by J.W. Lavrakas and D. Knezha in Proceedings of the 1999 National Technical Meeting of The Institute of Navigation, San Diego, California, January 25–27, 1999, pp. 621–626.

    • GPS Integrity and Receiver Autonomous Integrity Monitoring

    “Prototyping Advanced RAIM for Vertical Guidance” by J. Blanch, M.J. Choi, T. Walter, P. Enge, and K. Suzuki in Proceedings of ION GNSS 2010, the 23rd International Technical Meeting of the Satellite Division of The Institute of Navigation, Portland, Oregon, September 21–24, 2010, pp. 285–291.

    “The Integrity of GPS” by R.B. Langley in GPS World, Vol. 10, No. 3, March 1999, pp. 60–63.

    • International GNSS Service Ephemerides and Clocks

    “On the Precision and Accuracy of IGS Orbits” by J. Griffiths and J.R. Ray in Journal of Geodesy, Vol. 83, 2009, pp. 277–287, doi: 10.1007/s00190-008-0237-6.

    “The International GNSS Service: Any Questions?” by A.W. Moore in GPS World, Vol. 18, No. 1, January 2007, pp. 58–64.

    International GNSS Service Central Bureau website.

    • National Geospatial-Intelligence Agency Ephemerides and Clocks

    “NGA’s Role in GPS” by B. Wiley, D. Craig, D. Manning, J. Novak, R. Taylor, and L. Weingarth in Proceedings of ION GNSS 2006, the 19th International Technical Meeting of the Satellite Division of The Institute of Navigation, Fort Worth, Texas, September 26–29, 2006, pp. 2111–2119.

    National Geospatial-Intelligence Agency, Geoint Sciences Office, Global Positioning System Division website.

    • Antenna Phase Center Corrections

    “Generation of a Consistent Absolute Phase-center Correction Model for GPS Receiver and Satellite Antennas” by R. Schmid, P. Steigenberger, G. Gendt, M. Ge, and M. Rothacher in Journal of Geodesy, Vol. 81, No. 12, 2007, pp. 781–798, doi: 10.1007/s00190-007-0148-y.

    “The Block IIA Satellite: Calibrating Antenna Phase Centers” by G.L. Mader and F.M. Czopek in GPS World, Vol. 13, No. 5, May 2002, pp. 40–46.

    • GPS Interface and Performance Specifications

    Navstar GPS Space Segment / Navigation User Interfaces, Interface Specification, IS-GPS-200 Revision E, prepared by Science Applications International Corporation, El Segundo, California, for Global Positioning System Wing, June 2010.

    Global Positioning System Standard Positioning Service Performance Standard, 4th edition, by the U.S. Department of Defense, Washington, D.C., September 2008.

    Global Positioning System Standard Positioning Service Performance Standard, 3rd edition, by the U.S. Department of Defense, Washington, D.C., October 2001.

  • Integrity for Non-Aviation Users: Moving Away from Specific Risk

    Non-aviation users of satellite- and ground-based augmentation systems do not require the conservative level of integrity built into these systems for aviation users. Removing it can produce substantial benefits in terms of smaller error bounds and improved availability.

    By Sam Pullen, Todd Walter, and Per Enge

    Both space-based and ground-based augmentation systems (SBAS and GBAS, respectively) are designed to enhance standalone GNSS navigation to meet the requirements of civil aviation. SBAS and GBAS corrections and integrity information are also available to the non-aviation user population, such as automobiles, buses, and trains on land as well as ships near shore. This much larger user base can benefit as much from the integrity components of SBAS and GBAS as from the increased accuracy obtained from applying SBAS and GBAS pseudorange corrections. However, there are significant differences between the aviation interpretation of navigation integrity and the interpretation that would be natural to most users.

    SBAS and GBAS provide integrity in a multi-step procedure that is laid out in the RTCA Minimum Operational Performance Standards (MOPS) for the FAA versions of both systems: DO-229D for the Wide Area Augmentation System (WAAS) and DO-253C for the Local Area Augmentation System (LAAS). These systems indicate which ranging measurements should be excluded as unsafe to use and provide bounding error standard deviations, or sigmas, for the remaining usable measurements. Each aircraft uses this information to compute vertical and horizontal protection levels that define position-domain error bounds at desired probabilities. This process is straightforward, logical, and is not limited to aviation users. However, the requirements and assumptions underlying it make it very conservative.

    SBAS and GBAS are designed to meet integrity requirements defined in terms of what is known as specific risk. Briefly, this means that all safety requirements must be met for the worst combination of knowable or potentially foreseeable circumstances under which an operation may be conducted. Some variable factors important to safety, such as the user’s satellite geometry, are known by definition. Others, such as receiver thermal noise, are random and unpredictable. But several factors that are critical to GNSS performance, such as multipath and ionospheric errors, are neither completely random nor deterministic. Specific risk typically treats all error sources that are not completely random in a worst-case manner. SBAS and GBAS are designed to mitigate specific risk to support civil aviation, and the resulting conservatism makes SBAS and GBAS less attractive to non-aviation users who expect tighter protection levels relative to nominal system accuracy.

    Fortunately, non-aviation users need not apply all MOPS procedures required of aviation users if their own safety requirements differ. Most users define integrity in average or ensemble terms, meaning that everything not known in practice is treated as random and is probabilistically mixed (or convolved) together. The protection levels valid for these users would be much lower than for aviation users, even though the stated bounding probability is the same. This contrast is illustrated in Figure 1, which shows example bounds on 2-D vertical errors at a probability of 0.95 (the 95th percentile, or 95 percent) for accuracy and a probability of 1–10-7 for integrity. The term VPE stands for vertical position error, while VPL stands for vertical protection level. Analogous terms (HPE and HPL) and a similar picture exist in two dimensions for horizontal errors.

    Only one 95 percent error bound is shown in Figure 1 because this probability can be observed, estimated, and modeled with theory and reasonable amounts of data (hundreds or thousands of independent samples). This is not at all the case at the very small probability of 10-7 that applies to aviation precision approach: it is roughly equivalent to one event in 47.5 years per 150-second precision-approach interval. Both theory and data fall far short of being able to predict such rare-event errors. Extrapolating from available data to 1–10-7 using Gaussian distributions is perilous because the Gaussian distribution almost never applies at such small probabilities. Mixed-Gaussian models, other so-called fat-tailed distributions, and inflation of Gaussian parameters help address this, but the uncertainty regarding the true error distribution results in significantly different error bounds depending on the assumptions that are made. The same is true regarding the effects of faults and anomalies that are more probable than 10-7 but are still rare and poorly understood.

    In the end, different means of assessing these uncertainties and various degrees of user risk aversion result in different 1–10-7 protection levels, as shown in Figure 1. It is this difference that we wish to quantify and exploit in this article.

    Average versus Specific Risk

    The concept of average or ensemble risk is intuitive to those with a background in probability and is one of the key principles of probabilistic risk assessment (PRA). Thus, it helps to examine it first.

    Average risk is the probability of unsafe conditions based upon the convolved (averaged) estimated probabilities of all unknown events. More specifically, probability distributions are derived (based on the best available knowledge) for all unknown parameters relevant to user safety, and these are combined (by probabilistic convolution) to create an overall distribution that represents safety risk as a function of the known parameters. This straightforward, natural interpretation of probability and uncertainty has a major advantage in that it cleanly separates the probabilistic calculation of safety risk from users’ aversion to risk. By keeping risk probability and risk aversion (or severity) separate, a final risk consequence measure can be derived that supports apples-to-apples comparisons of alternatives. One useful result of this is known as the value of information (VOI). By comparing the risk outcomes of two scenarios in which the latter case has additional information (for example, from an additional sensor or integrity monitor), the risk-reduction benefit of the added information can be traded off against the cost and complexity that it introduces to the system. Similar comparisons can be made for any definition of risk, but the definition and use of VOI in an average-risk framework makes the most sense in both theory and practice.

    Turning to specific risk, no single definition exists within the aviation safety community, to our knowledge. This is partially because of the uniqueness and complexity of the concept and partially because multiple inconsistent interpretations appear to exist. Therefore, we provide our own definition: Specific risk is the probability of unsafe conditions subject to the assumption that all credible unknown events that could be known occur with a probability of one (on a risk-by-risk basis).

    To understand how specific risk differs from average risk, it helps to start with a fault-tree representation of risk in which loss of integrity (LOI) can result from any of the nodes of the tree. Figure 2 shows a simplified example of a fault tree for CAT I GBAS. It shows the allocation of the CAT I total integrity risk requirement of 2 × 10-7 per approach to the various possible causes of integrity loss. In specific-risk analysis, each type of failure shown in the tree, if deemed to be a credible failure (meaning, in practice, that its assumed prior probability is larger than compared to its allocation in the fault tree), is assessed that the failure is guaranteed to occur in a worst-case fashion. This means that the variables that describe this particular failure scenario take the values that maximize the hazard to users. In an average-risk analysis, these variables would take many values according to their own probability distributions, and these distributions would be convolved together to provide an overall representation of risk under that scenario. Instead, one scenario drives the specific risk assessment for a particular user class, and it is the worst one possible from that user’s standpoint. (Another user class would be evaluated under a different set of parameters corresponding to the separate worst case for that user.) The improbability of the worst-case combination of parameters is not considered as long as the probability of the failure scenario as a whole is deemed high enough to be of concern.

    Figure 2. Fault tree for CAT I GBAS integrity. Source: Sam Pullen, Todd Walter, and Per Enge
    Figure 2. Fault tree for CAT I GBAS integrity.

    Since GNSS augmentation systems contain multiple levels of health monitoring, the worst-case scenario is usually the one that maximizes the probability of an undetected hazardous error for a particular user class. Hazardous error is typically defined as any error that exceeds a pre-defined safety zone known as an alert limit (AL) or any error that exceeds the computed protection level (PL), which allows integrity to be defined separately from the intended application. Both definitions are conservative in that all errors exceeding AL or PL are treated as equally hazardous. In other words, an error just above AL is treated as just as dangerous as an error of 10 × AL. They are also misleading when used in specific-risk analyses because the resulting worst-case conditions are those that give errors just above AL or PL, as these are the generally hardest for monitoring algorithms to detect.

    The use of specific risk in aviation is an evolution of deterministic guidelines for tolerable risk that date back to an earlier era when flying was more dangerous. It remains dominant in aviation safety assessment because it is partly responsible for the development of safer and more reliable air transportation. However, it has several important weaknesses compared to average risk. The first is that the degree of risk aversion preferred for aviation is buried within the hazard probabilities generated by specific risk — it cannot be separated out. This means that specific-risk results do not translate well to other classes of users, as very few users would happen to have the same risk preferences that have evolved within aviation over several decades. In addition, specific risk makes a distinction between unknown events that could be known and those that are both rare and completely unknowable. A very risk-averse value of information is much different than the risk-neutral one built into PRA, as it severely penalizes systems that do not include all potentially-informative sensors. Since each sensor added to a system provides less benefit than the last, almost all cost-effective systems choose to include less than the maximum possible number of sensors.

    The conservatism implicit in specific-risk assessment severely penalizes users. Although PRA would show that the combination of factors (shown in an example induced by extreme ionospheric spatial decorrelation) needed to produce a 40-meter error in a CAT I GBAS system is exceedingly improbable (almost certainly below 10-10 per approach), specific risk forces a significant part of the GBAS risk-mitigation effort to be targeted at this scenario. In this case, since monitoring is not guaranteed to detect the anomaly in time, the only recourse is geometry screening, a cumbersome technique in which the ground system continually evaluates the worst-case error and, if it exceeds a 28-meter tolerable limit at the CAT I decision height, determines which broadcast parameters to inflate such that all satellite geometries causing worst-case errors exceeding 28 meters are made unavailable (because the inflated VPL is larger than the 10-meter CAT I VAL). The result of this procedure is much lower user availability than would be achieved without inflation. SBAS pays a similar penalty, as we will see later. The broadcast grid ionospheric vertical error values that bound worst-case ionospheric errors (and thus the resulting protection levels) are much higher than they would be if the unusual combination of factors needed to create the worst-case error scenario were not the dominant concern.

    To the extent that loss of availability represents a safety issue at the airspace level, the worst-case focus that results from specific risk is not optimal even from a safety standpoint. But this is not the only concern. Specific risk requires a great deal of development and testing to identify and mitigate a handful of very peculiar, non-representative conditions. When schedule and resources are limited, other potential threats that are easier to foresee but seem extremely improbable are often neglected. One example is the treatment of multiple hardware failures. If individual failures are assumed to be statistically independent, the probability of multiple simultaneous failures is very small. However, while statistical independence is a common assumption in probability classes because it makes calculations easier, it rarely applies in the real world. Because satellites and ground receivers are similar, if not identical, the presence of a failure in one unit may suggest a common cause or at least a common vulnerability, meaning that the probability of additional failures is much higher than independence would suggest. Thus, assuming independence by default could lead to neglecting entire categories of risk that are more threatening than the worst-case events that dominate specific risk.

    Maximum WAAS Errors, Protection

    To investigate the conservatism built into SBAS and GBAS specific risk assessment, maximum WAAS horizontal and vertical position errors over time (as measured by the Performance Analysis Network (PAN) maintained by the William J. Hughes FAA Technical Center) have been examined and compared to the protection levels computed when the maximum errors occurred. This study begins with PAN Report #8 (covering January to March 2004 — shortly after WAAS commissioning in mid-2003) and extends through PAN Report #34 (covering July to September 2010). Each PAN report covers three months of observed WAAS performance.

    Figure 3 shows the 38 WAAS reference stations (WRSs) used by the PAN to collect position error and protection level information (some of these stations were not active in 2004 and thus were not used in earlier PAN reports). While measurements from these stations are used to generate WAAS corrections and error bounds, they are also used by the PAN as static pseudo-users that compute WAAS-corrected positions and protection levels according to the aircraft user algorithms specified in the WAAS MOPS. The resulting positions are compared to the known, pre-surveyed positions of each station to derive estimates of vertical and horizontal position errors (VPE and HPE) once per second.

    Figure3 Source: Sam Pullen, Todd Walter, and Per Enge
    Figure 3. WAAS PAN reference station network.

    Figure 3 groups these stations into three sets of stations based on their presumed quality of WAAS coverage. These sets are unofficial and were created for the purposes of this study. The seven stations in the inner set are expected to have good WAAS coverage at all times because they are surrounded by other stations. The 13 stations in the outer set are expected to only have acceptable coverage because s
    ome of them are at the edges of CONUS. The remote stations provide coverage to the inner and outer regions as well as the best possible coverage of their own regions. Because the remote stations extend beyond the primary coverage region of WAAS in CONUS, errors at these stations are not considered here.

    Figure 4 is a 2-D plot of position error versus protection level in the vertical axis (that is, VPE versus VPL) for all epochs and stations during the three months covered by the recent WAAS PAN Report #34 (July 1–September 30, 2010). These results are typical of the entire period since WAAS commissioning in 2003, particularly the last several years. The vertical lines on the plot indicate the 95th-percentile, 99th percentile, and maximum VPEs in this period (1.2, 1.8, and 7 meters, respectively). The maximum VPE occurred at Barrow, AK, which is one of the most remote stations in the WAAS network (see Figure 3). In comparison, the lowest VPLs (intended to be 1–10-7 bounds on VPE) are in the range of 10–15 meters, and values as high as 40 meters are not uncommon. The most demanding approach operation that WAAS supports, LPV-200, allows approaches to a 200-foot minimum decision height and requires that VPL be below a vertical alert limit (VAL) of 35 meters. HPL must also be below a horizontal alert limit (HAL) of 45 meters. When this is not the case, the approach operation is not available; thus these higher VPLs extract a significant cost.

    Figure 4. WAAS vertical protection level versus vertical position error (June–September 2010). Source: Sam Pullen, Todd Walter, and Per Enge
    Figure 4. WAAS vertical protection level versus vertical position error (June–September 2010).

    Figure 5 and Figure 6 (for vertical and horizontal errors, respectively) span the entire period of WAAS PAN Reports used in this study. VPL represents the VPL at the station and time of the maximum VPE; it is not the largest VPL recorded at a particular station. The horizontal errors shown in Figure 6 are defined analogously. Note that the station that observes the largest horizontal error in a given PAN report may differ from the one that observes the largest vertical error.

    Figures 5 and 6 demonstrate that, while both 95 percent and maximum errors are quite low and are within the expected range of each other, the protection levels associated with the maximum errors greatly exceed them. This pattern is clearer in Figure 5 for vertical errors because maximum VPL tends to be more consistent across PAN reports, but it is true for horizontal errors as well.

    Figure 5. WAAS vertical errors and protection levels from 2004–2010. Source: Sam Pullen, Todd Walter, and Per Enge
    Figure 5. WAAS vertical errors and protection levels from 2004–2010.
    Figure 6. WAAS horizontal errors and protection levels from 2004–2010. Source: Sam Pullen, Todd Walter, and Per Enge
    Figure 6. WAAS horizontal errors and protection levels from 2004–2010.

    Figures 7 and 8 clarify this relationship by plotting the ratio of VPL to VPE and HPL to HPE for the station and time of the maximum error. The mean of this ratio is very high and is about the same in both cases: 5.38 for vertical and 5.21 for horizontal. Figure 7 shows a steady upward trend in the ratio that is mostly due to WRS improvements that resulted in maximum VPE being reduced over time. This trend is clearly visible in Figure 5 and appears to exceed the weaker trend of lowering VPL due to WAAS algorithm enhancements. The same trend is visible in the horizontal Figures 6 and 8 but is weaker due to the greater variability of maximum HPL over time.

    To evaluate the significance of the large PL-to-max-PE ratios in the WAAS PAN database, we need to approximate the number of independent samples from which the maximum errors were derived. As noted before, WAAS protection levels represent error bounds at the 1–10-7 probability level based on specific risk. With one measurement being collected at each operational station every second, a total of about 4.25 billion samples were collected in the PAN reports from January 2004 to September 2010. Note that measurements from remote stations are included in this count, but they are also represented in the conclusions because their PL-to-max-PE ratios are very similar to the ones shown in Figures 7 and 8. Translating this number into the number of statistically independent samples depends on the interval between independent measurements. Because both nominal and rare-event errors affect this interval, it is hard to estimate. Our best guess is a range between roughly 30 and 150 seconds, suggesting that the PAN database contains between 2.8 × 107 and 1.4 × 108 independent samples. Both of these numbers suggest that WAAS protection levels are very conservative from the perspective of average risk.

    Figure 7. Ratio of VPL to VPE from 2004–2010. Source: Sam Pullen, Todd Walter, and Per Enge
    Figure 7. Ratio of VPL to VPE from 2004–2010.
    Figure 8. Ratio of HPL to HPE from 2004–2010. Source: Sam Pullen, Todd Walter, and Per Enge
    Figure 8. Ratio of HPL to HPE from 2004–2010.

    Adjusting for Average-Risk Users

    Using the above results, a preliminary estimate of the reduced WAAS protection levels that would apply to average-risk users can be made. Figure 9 shows a comparison between the actual 95 percent WAAS VPL and HPL and the adjusted VPL and HPL potentially achievable with WAAS (for the same 1–10-7 bounding probability) for average-risk users. The actual WAAS VPLs are taken from the more recent WAAS PAN Reports starting from #24 (covering January to March 2008) as the period from 2008 to 2010 includes most of the WAAS algorithm improvements introduced since commissioning in 2003. The actual 95 percent VPLs and HPLs represent the largest reported 95th-percentile values among the stations within CONUS for each quarterly period. The lower adjusted VPLs and HPLs are derived by dividing each VPL by a factor of 4.0 and each HPL by a factor of 2.5. These two reduction factors are derived from Figures 7 and 8, respectively, as conservative estimates of the ratio between protection levels and maximum position errors. Note that the factor of 2.5 for horizontal errors does not include the 12-meter error in Cleveland from PAN Report #13, as this is thought to be spurious (that is, not representative of actual WAAS behavior).

    Figure 9. Projected WAAS protection level reductions for average-risk users. Source: Sam Pullen, Todd Walter, and Per Enge
    Figure 9. Projected WAAS protection level reductions for average-risk users.

    While projections based on these reduction factors are imprecise, they demonstrate the much lower error bounds that non-aviation users with an average-risk safety perspective could achieve. Most non-aviation users operate on land or sea and will be primarily concerned with horizontal error bounds. Figure 9 suggests that the typical 95th percentile WAAS HPLs of 15–20 meters (for the worst location in CONUS) can be reduced to 6–8 meters and still provide a confident 1–10-7 error bound.

    It is important to emphasize that these preliminary projections for average-risk users are just that. In order to formally establish new integrity requirements and protection levels for existing systems, the hazardously misleading information (HMI) analyses previously done for these systems need to be redone using the principles of PRA and average risk. While the original development of the WAAS and LAAS HMI analyses was lengthy and resource-intensive, almost all of the detailed work is already complete. As long as the original analyses are available, it is a much smaller task to take these results and create PRAs out of them by extracting the original specific-risk assumptions and applying average-risk principles instead.

    LAAS Users. Since the first GBAS ground station design (the Honeywell SLS-4000 LAAS Ground Facility) was certified for CAT I use in 2009 and has not yet been approved for operations at a specific airport, much less data is available to do a preliminary analysis for GBAS similar to the one done for WAAS above. However, the degree of sigma inflation in the parameters broadcast by CAT I LAAS is approximately known, meaning that it can be more-precisely removed from the current LAAS protection levels to estimate what they would be for average-risk users.

    Figure 10 shows the degree of inflation applied to the broadcast σvertical_iono_gradient (or σvig) parameter in order to protect against the worst-case ionospheric anomaly described previously. This result is for the SPS-standard 24-satellite constellation over a 24-hour period at the LAAS installation at Newark Airport, New Jersey (the method used by the Honeywell SLS-4000 is somewhat different). While not all epochs require inflation, a majority cause the nominal σvig value to be increased by a factor of 2 or more, which significantly decreases CAT I availability and currently makes it impossible to take advantage of the Differentially Corrected Positioning Service (DCPS) for non-CAT-I operations.

    Figure 10. Typical σvig inflation factors for CAT I LAAS. Source: Sam Pullen, Todd Walter, and Per Enge
    Figure 10. Typical σvig inflation factors for CAT I LAAS.

    Because of the extreme rarity of the worst-case event that dictates this inflation, it would likely not be needed for average-risk users. Figure 11 shows how much the σvig inflation in Figure 10 increases the LAAS VPL at Newark for the standard 24-satellite constellation. The VPL reduction from removing the inflation is not as dramatic as the potential reductions shown for WAAS in Figure 9, but they are significant relative to the 10-meter VAL for LAAS CAT I approaches. Furthermore, the pre-inflated nominal value of σvig for LAAS is 6.4 millimeters/kilometer, which is much higher than the actual one-sigma nominal gradient value of 1–2 mm/km because, under specific risk, the very worst nominal data must be bounded (also, worst-case tropospheric gradients must also be bounded by σvig). Other broadcast parameters that affect VPL, such as σpr_gnd and the ephemeris P-value that bounds worst-case ephemeris failures, would also be reduced significantly by switching to average risk. Overall, it is likely that LAAS protection levels based on average risk would be reduced from the current specific-risk PLs by about the same range of factors (2–5) observed from WAAS data.

    Figure 11. Impact of σvig inflation on LAAS VPL. Source: Sam Pullen, Todd Walter, and Per Enge
    Figure 11. Impact of σvig inflation on LAAS VPL.

    User Performance Improvements

    This discussion assumes that most non-aviation users who are not encumbered by the history of aviation standards development will prefer to quantify risk using PRA and the average-risk approach. As noted earlier, average risk better matches most users’ intuitive understanding of uncertainty and has the enormous advantage of separating risk quantification from risk aversion. Regardless of how risk-averse or conservative a given operator is, his or her model of risk aversion can be applied most efficiently to a risk-neutral calculation of risk that fairly represents all aspects of uncertainty. Inserting risk aversion into the calculation of risk, as done in the specific-risk approach, is both inefficient and non-optimal from a safety perspective because extensive focus on a few extreme worst-case events drives attention away from other events.

    The HPL reductions for average-risk users illustrated here would be significant for many classes of ground and marine transportation users. They would allow operations with tighter physical safety margins to be supported. Users who gain no particular benefit from tighter protection levels would still obtain much higher availability of integrity, as a 25-meter HPL could be supported by much poorer satellite geometries than would otherwise be the case. In other words, users that can tolerate 25-meter horizontal error bounds would be able to operate safely a much higher percentage of the time, because the degree of GNSS constellation deterioration needed to exceed this limit would occur much less often. These benefits do not only apply at the 1–10-7 probability level, as they would scale to the higher probabilities (1–10-4 to 1–10-6) that many non-aviation applications would be most concerned with.

    While very few non-aviation users of GNSS today have real-time safety requirements similar to those of civil aviation, the number of such users will likely increase as the coverage of augmented GNSS (and the availability of integrity from standalone receiver-autonomous integrity monitoring, or RAIM) expands. The evolution of standalone civil GPS usage provides a precedent: as basic GPS accuracy improved from tens of meters to several meters, and the cost of user equipment dropped, more and more uses were discovered. A similar, although smaller-scale, trend is likely to occur as the advantages of augmented GNSS become more available and better understood. The primary beneficiaries are likely to be intelligent road-transport systems, train services, and marine transportation in restricted waters.

    One application where tight real-time integrity bounds would be useful is in harbor and marina entry and exit; see Figure 12, taken from a Google map of a marina in San Diego, California. Based on the earlier analysis, two typical 1−10-7 horizontal protection levels are shown: 18 meters using the unchanged WAAS MOPS approach, and 7 meters based upon modifying the broadcast bounding parameters to represent average risk (these HPLs are bounds on error in either direction, positive or negative; thus the 2-D error bounding circle has a diameter of twice the HPL).

    Figure 12. Example of reduced protection levels for harbor/marina access. Source: Sam Pullen, Todd Walter, and Per Enge
    Figure 12. Example of reduced protection levels for harbor/marina access.

    When the resulting error bounds are compared, the relative advantage of the smaller bound for this application is immediately apparent. In general, when HPL is significant compared to potential obstacles, its significance varies with the square of HPL rather than HPL itself, as the area being protected matters more than either linear direction. In this example, the ratio of HPLs being compared is 18/7, or 2.57, but the ratio of HPL-squared is much larger: 182/72 = 6.61.

    When real-time integrity is not needed, augmented GNSS provides an easy means to guarantee or certify vehicle locations after the fact with great precision and reliability, without the need for post-processing. Vehicle and cargo tracking based on standalone GPS is common today, a certification of the correctness of the tracking data to probabilities suitable for legal or commercial guarantees is lacking. For this, error bounds at 1–10-4 to 1– 10-6 probabilities are likely sufficient, and would allow HPLs of below 5 meters from WAAS and below 3 meters from LAAS. In some scenarios, the difference between a 5-meter and a 15-meter guarantee would be minor, but in others, it could make a substantial difference.

    As noted earlier, even for uses where the required HPL (as represented by the safe error limit, or HAL, for a particular application) is satisfied by the existing WAAS and LAAS protection levels, the use of modified average-risk protection levels increases the availability of integrity, which is most often expressed as the probability or percentage of time (over all satellite geometries and othe
    r variable system states) that the integrity requirement is met throughout an operation (in simple terms, that HPL ≤ HAL). For user locations within good WAAS or LAAS coverage, the most variable element over time is satellite geometry. Decreasing HPL by a factor of 2.5 or more substantially increases the margin between HPL and HAL and makes it far less likely that the satellite geometry will degrade to the point where HPL exceeds HAL. For example, if the unmodified WAAS HPL equals HAL at an (un-weighted) HDOP of about 1.5, the resulting satellite availability (an upper bound on overall availability) for the SPS-standard 24-satellite GPS constellation would be roughly 98.5 percent. This means that the satellites in view (in this case, all satellites above 5 degrees elevation at a location in CONUS) would provide HDOP ≤ 1.5 about 98.5 percent of the time. However, the modified average-risk HPL (using the factor-of-2.5 reduction) would roughly translate into a limiting HDOP of about 3.75. This allows the required integrity bound to be satisfied by much poorer GPS geometries and gives a satellite availability of greater than 99.9 percent. Thus, when integrity is needed, this much greater availability of integrity is a major advantage.

    Summary

    SBAS and GBAS broadcasts are freely available to all GNSS users, most of whom will have different definitions of acceptable risk. These users are not optimally served at present and may hesitate to take advantage of SBAS and GBAS as a result.

    Using years of collected data for the FAA WAAS system and analysis of the inflation factors built into the CAT I version of the FAA LAAS system, it appears that average-risk users of WAAS and LAAS would be adequately supported by protection levels that are 2 to 5 times lower than those currently derived by aviation users. The fact that two different approaches used to examine WAAS and LAAS suggest similar levels of over-conservatism lends credence to these estimates. While further validation by full-scale probabilistic risk assessments is necessary, we conclude that non-aviation users willing to accept average risk would obtain much better performance and availability from simple modifications to the existing SBAS and GBAS protection level calculations specified for aviation users.

    Acknowledgments

    We thank the FAA Satellite Navigation Program Office for its support of our research on WAAS and LAAS. However, the opinions expressed here are solely our own. We thank Jim Kelly and Tim Murphy for their explanations of the evolution of today’s SBAS and GBAS integrity requirements. We also thank the FAA Technical Center for its efforts in collecting and publishing WAAS error data over the last decade using its Performance Analysis Network (PAN).


    Sam Pullen is a senior research engineer at Stanford University, where he is the director of the Local Area Augmentation System (LAAS) research effort. He has supported the FAA and others in developing GNSS system concepts, requirements, integrity algorithms, and performance models since obtaining his Ph.D. from Stanford in Aeronautics and Astronautics.

    Todd Walter is a senior research engineer in the Department of Aeronautics and Astronautics at Stanford University. He received his Ph.D. from Stanford and is currently working on the Wide Area Augmentation System (WAAS), defining future architectures to provide aircraft guidance, and on assuring integrity on GPS III.

    Per Enge is a professor of aeronautics and astronautics at Stanford University, where he is the Kleiner Perkins, Mayfield, Sequoia Capital Professor in the School of Engineering. He directs the GPS Research Laboratory and received his Ph.D. from the University of Illinois.