<?xml version="1.0" encoding="UTF-8"?>
<rss version="2.0"
	xmlns:content="http://purl.org/rss/1.0/modules/content/"
	xmlns:wfw="http://wellformedweb.org/CommentAPI/"
	xmlns:dc="http://purl.org/dc/elements/1.1/"
	xmlns:atom="http://www.w3.org/2005/Atom"
	xmlns:sy="http://purl.org/rss/1.0/modules/syndication/"
	xmlns:slash="http://purl.org/rss/1.0/modules/slash/"
	>

<channel>
	<title>GPS World &#187; GNSS</title>
	<atom:link href="http://www.gpsworld.com/category/gnss-system/feed/" rel="self" type="application/rss+xml" />
	<link>http://www.gpsworld.com</link>
	<description>The Business and Technology of Global Navigation and Positioning</description>
	<lastBuildDate>Mon, 10 Jun 2013 18:54:37 +0000</lastBuildDate>
	<language>en-US</language>
	<sy:updatePeriod>hourly</sy:updatePeriod>
	<sy:updateFrequency>1</sy:updateFrequency>
	<generator>http://wordpress.org/?v=3.5.1</generator>
		<item>
		<title>Air Force Video Explains GPS Role in Daily Life</title>
		<link>http://www.gpsworld.com/air-force-video-explains-gps-role-in-daily-life/?utm_source=rss&#038;utm_medium=rss&#038;utm_campaign=air-force-video-explains-gps-role-in-daily-life</link>
		<comments>http://www.gpsworld.com/air-force-video-explains-gps-role-in-daily-life/#comments</comments>
		<pubDate>Wed, 05 Jun 2013 20:15:52 +0000</pubDate>
		<dc:creator>GPS World staff</dc:creator>
				<category><![CDATA[Aviation & Space]]></category>
		<category><![CDATA[Defense News]]></category>
		<category><![CDATA[GNSS News]]></category>
		<category><![CDATA[Government News]]></category>
		<category><![CDATA[Latest News]]></category>
		<category><![CDATA[Air Force]]></category>
		<category><![CDATA[video]]></category>

		<guid isPermaLink="false">http://www.gpsworld.com/?p=21705</guid>
		<description><![CDATA[All of us in the GPS industry know someone who only thinks of GPS as a feature of their smartphone. You might direct them to a new YouTube video presented by the U.S. Air Force, which summarizes the worldwide role of GPS. It also touches on the GPS modernization program and new signals. The seven-minute [...]]]></description>
				<content:encoded><![CDATA[<p>All of us in the GPS industry know someone who only thinks of GPS as a feature of their smartphone. You might direct them to <a href="http://www.youtube.com/watch?v=chNQW22vVNI" target="_blank">a new YouTube video</a> presented by the U.S. Air Force, which summarizes the worldwide role of GPS. It also touches on the GPS modernization program and new signals.</p>
<p>The seven-minute video explains in simple terms how important GPS has become to everyday life — for aircraft and ship navigation, global financial transactions, precision agriculture, weather forecasting, disaster relief, and, of course, smartphones.</p>
<p><iframe src="http://www.youtube.com/embed/chNQW22vVNI" height="315" width="420" allowfullscreen="" frameborder="0"></iframe></p>
]]></content:encoded>
			<wfw:commentRss>http://www.gpsworld.com/air-force-video-explains-gps-role-in-daily-life/feed/</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Lockheed Martin Completes Functional Testing of GPS III Electronic Systems</title>
		<link>http://www.gpsworld.com/lockheed-martin-completes-functional-testing-of-gps-iii-electronic-systems/?utm_source=rss&#038;utm_medium=rss&#038;utm_campaign=lockheed-martin-completes-functional-testing-of-gps-iii-electronic-systems</link>
		<comments>http://www.gpsworld.com/lockheed-martin-completes-functional-testing-of-gps-iii-electronic-systems/#comments</comments>
		<pubDate>Wed, 05 Jun 2013 17:20:34 +0000</pubDate>
		<dc:creator>GPS World staff</dc:creator>
				<category><![CDATA[Aviation & Space]]></category>
		<category><![CDATA[Defense News]]></category>
		<category><![CDATA[GNSS News]]></category>
		<category><![CDATA[GPS Modernization]]></category>
		<category><![CDATA[Latest News]]></category>

		<guid isPermaLink="false">http://www.gpsworld.com/?p=21694</guid>
		<description><![CDATA[A Lockheed Martin-led industry team has completed successful functional integration tests of the spacecraft bus and network communications equipment on the first satellite of the next generation Global Positioning System, known as GPS III. The recent testing of GPS III space vehicle 1 (SV 1) bus — the portion of the space vehicle that carries [...]]]></description>
				<content:encoded><![CDATA[<p>A <a href="http://www.lockheedmartin.com/" target="_blank">Lockheed Martin</a>-led industry team has completed successful functional integration tests of the spacecraft bus and network communications equipment on the first satellite of the next generation Global Positioning System, known as GPS III.</p>
<p>The recent testing of GPS III space vehicle 1 (SV 1) bus — the portion of the space vehicle that carries mission payloads and hosts them in orbit — assured that all bus subsystems are functioning normally and ready for final integration with the satellite&#8217;s navigation payload. Systems tested included: guidance, navigation and control; command and data handling; on-board computer and flight software; environmental controls; and electrical power regulation. The SV 1 satellite&#8217;s network communication equipment subsystem that interfaces with the ground control segment and distributes data throughout the space vehicle also passed all tests as expected.</p>
<p>This milestone follows <a href="http://www.gpsworld.com/lockheed-martin-powers-on-first-gps-iii-satellite/" target="_blank">February&#8217;s successful initial power-on of SV 1</a>, which demonstrated the electrical-mechanical integration, validated the satellite&#8217;s interfaces, and led the way for functional and hardware-software integration testing.</p>
<p>&#8220;The successful completion of the SV 1 bus functional check out validates that the spacecraft is now ready to begin the next sequence of payload integration and environmental testing, prior to delivery,&#8221; explained Keoki Jackson, vice president of Lockheed Martin&#8217;s Navigation Systems mission area.</p>
<p>GPS III SV 1&#8242;s navigation payload, which is being produced by <a href="http://www.lockheedmartin.com/gps" target="_blank">ITT Exelis</a>, will be delivered to Lockheed Martin&#8217;s GPS Processing Facility (GPF) near Denver later in 2013. The hosted nuclear detection system payload has already been delivered and mechanically integrated. The satellite remains on schedule for flight-ready delivery to the U.S. Air Force in 2014.</p>
<p>GPS III is a critically important program for the Air Force, affordably replacing aging GPS satellites in orbit, while improving capability to meet the evolving demands of military, commercial and civilian users. GPS III satellites will deliver three times better accuracy and — to outpace growing global threats that could disrupt GPS service — up to eight times improved anti-jamming signal power for additional resiliency. The GPS III will also include enhancements adding to the spacecraft&#8217;s design life and a new civil signal designed to be interoperable with international global navigation satellite systems.</p>
<p>The U.S. Air Force has produced a video about the GPS satellite modernization program:<br />
<iframe src="http://www.youtube.com/embed/chNQW22vVNI" height="315" width="420" allowfullscreen="" frameborder="0"></iframe></p>
<p>Lockheed Martin is under contract for production of the first four GPS III satellites (SV 1-4), and <a href="http://www.gpsworld.com/air-force-awards-lockheed-martin-contracts-for-next-set-of-gps-iii-satellites/" target="_blank">has received advanced procurement funding</a> for long-lead components for the fifth, sixth, seventh and eighth satellites (SV 5-8).</p>
<p>The GPS III team is led by the <a href="http://www.losangeles.af.mil/library/factsheets/factsheet.asp?id=18830" target="_blank">Global Positioning Systems Directorate</a> at the U.S. Air Force Space and Missile Systems Center. Lockheed Martin is the GPS III prime contractor with teammates ITT Exelis, General Dynamics, Infinity Systems Engineering, Honeywell, ATK and other subcontractors. <a href="http://www.schriever.af.mil/library/factsheets/factsheet.asp?id=4045" target="_blank">Air Force Space Command&#8217;s 2nd Space Operations Squadron (2SOPS)</a>, based at Schriever Air Force Base, Colo., manages and operates the GPS constellation for both civil and military users.</p>
<p>Headquartered in Bethesda, Maryland, Lockheed Martin is a global security and aerospace company that employs about 118,000 people worldwide and is principally engaged in the research, design, development, manufacture, integration, and sustainment of advanced technology systems, products, and services. The corporation&#8217;s net sales for 2012 were $47.2 billion.</p>
]]></content:encoded>
			<wfw:commentRss>http://www.gpsworld.com/lockheed-martin-completes-functional-testing-of-gps-iii-electronic-systems/feed/</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Navtech Offers Condensed GNSS Signals and Systems Course</title>
		<link>http://www.gpsworld.com/navtech-offers-condensed-gnss-signals-and-systems-course/?utm_source=rss&#038;utm_medium=rss&#038;utm_campaign=navtech-offers-condensed-gnss-signals-and-systems-course</link>
		<comments>http://www.gpsworld.com/navtech-offers-condensed-gnss-signals-and-systems-course/#comments</comments>
		<pubDate>Tue, 04 Jun 2013 17:29:43 +0000</pubDate>
		<dc:creator>GPS World staff</dc:creator>
				<category><![CDATA[GNSS News]]></category>
		<category><![CDATA[course]]></category>
		<category><![CDATA[ION GNSS]]></category>
		<category><![CDATA[Navtech]]></category>

		<guid isPermaLink="false">http://www.gpsworld.com/?p=21590</guid>
		<description><![CDATA[Navtech is offering a four-day version of Course 551, &#8220;Using Advanced GPS/GNSS Signals and Systems,&#8221; customized for those attending the ION GNSS+ 2013 conference. This course will help attendees develop proficiency with advanced receiver processing of current, modernized, and new signals from GPS, GLONASS, Galileo, BeiDou, and QZSS. It teaches systems engineering skills, along with [...]]]></description>
				<content:encoded><![CDATA[<p>Navtech is offering a four-day version of Course 551, &#8220;Using Advanced GPS/GNSS Signals and Systems,&#8221; customized for those attending the ION GNSS+ 2013 conference.</p>
<p>This course will help attendees develop proficiency with advanced receiver processing of current, modernized, and new signals from GPS, GLONASS, Galileo, BeiDou, and QZSS. It teaches systems engineering skills, along with techniques for receiver processing and for assessing processing performance. Review problems, worked in class, help students understand and apply the key concepts.</p>
<p>Those who attend will become proficient with the essential aspects of using GPS and GNSS signals.</p>
<p><strong>Course days:</strong><br />
Friday, Saturday, September 13-14<br />
Monday, Tuesday, September 16-17</p>
<p><strong>Instructor:</strong> Dr. John Betz, MITRE</p>
<p>For more information, <a href="http://www.navtechgps.com/events/course_541_details/" target="_blank">visit the Navtech website</a>.</p>
]]></content:encoded>
			<wfw:commentRss>http://www.gpsworld.com/navtech-offers-condensed-gnss-signals-and-systems-course/feed/</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Embezzlement of GLONASS Funds Investigated</title>
		<link>http://www.gpsworld.com/embezzlement-glonass-funds-investgated/?utm_source=rss&#038;utm_medium=rss&#038;utm_campaign=embezzlement-glonass-funds-investgated</link>
		<comments>http://www.gpsworld.com/embezzlement-glonass-funds-investgated/#comments</comments>
		<pubDate>Mon, 03 Jun 2013 22:30:15 +0000</pubDate>
		<dc:creator>GPS World staff</dc:creator>
				<category><![CDATA[GLONASS]]></category>
		<category><![CDATA[GNSS News]]></category>
		<category><![CDATA[Latest News]]></category>

		<guid isPermaLink="false">http://www.gpsworld.com/?p=21556</guid>
		<description><![CDATA[The Russian Federal Security Service is investigation the embezzlement of billions of rubles from the construction of the GLONASS center in Korolyov, a town outside Moscow, Izvestia daily reports. According to information shared by the Russian Legal Information Agency, the Investigative Committee&#8217;s department for the Moscow Region has launched a preliminary probe into the case. [...]]]></description>
				<content:encoded><![CDATA[<p>The Russian Federal Security Service is investigation the embezzlement of billions of rubles from the construction of the GLONASS center in Korolyov, a town outside Moscow, Izvestia daily reports.</p>
<p>According to information shared by the <a href="http://rapsinews.com/news/20130530/267615850.html" target="_blank">Russian Legal Information Agency</a>, the Investigative Committee&#8217;s department for the Moscow Region has launched a preliminary probe into the case.</p>
<p>Construction of the GLONASS satellite navigation system control and support center began in June 2010 on the site used by TsNIImash, the head research company of Russia&#8217;s federal space agency. The center was supposed to hold equipment for collecting and processing the data supplied by the GLONASS global network.</p>
<p>The construction was financed by a federal program, with 1.050 billion ($33.22 million) allocated for the project. By the end of 2010, it came to light that construction costs had been overstated, Izvestia reports. An expert appraisal revealed that the contractor had rigged the costs. The government did not allocate additional funds, so construction was suspended in December 2011 when the Federal GLONASS Program for 2002-2011 ended. The construction of the building has never been completed.</p>
<p>In November 2012, the <a href="http://www.gpsworld.com/designer-of-glonass-navigation-system-fired-amid-embezzlement-scandal/" target="_blank">general designer of GLONASS, Yuri Urlichich, was dismissed</a> from his post as a result of the scandal.</p>
]]></content:encoded>
			<wfw:commentRss>http://www.gpsworld.com/embezzlement-glonass-funds-investgated/feed/</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Navigation Center for India&#8217;s SatNav System Inaugurated</title>
		<link>http://www.gpsworld.com/navigation-center-for-indias-sat-nav-system-inaugurated/?utm_source=rss&#038;utm_medium=rss&#038;utm_campaign=navigation-center-for-indias-sat-nav-system-inaugurated</link>
		<comments>http://www.gpsworld.com/navigation-center-for-indias-sat-nav-system-inaugurated/#comments</comments>
		<pubDate>Mon, 03 Jun 2013 22:11:35 +0000</pubDate>
		<dc:creator>GPS World staff</dc:creator>
				<category><![CDATA[Augmentation & Assistance]]></category>
		<category><![CDATA[GNSS News]]></category>
		<category><![CDATA[Latest News]]></category>
		<category><![CDATA[India]]></category>
		<category><![CDATA[INRSS]]></category>

		<guid isPermaLink="false">http://www.gpsworld.com/?p=21549</guid>
		<description><![CDATA[The Indian Space Research Organization (ISRO) Navigation Centre, an important element of the Indian Regional Navigation Satellite System (IRNSS), was inaugurated May 28. The INC has been established at the Indian Deep Space Network complex at Byalalu, about 40 kilometers from Bangalore, India. IRNSS, an independent navigation satellite system being developed by India, will have [...]]]></description>
				<content:encoded><![CDATA[<p><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/isroi.jpg"><img class="size-full wp-image-21552 alignright" alt="isroi" src="http://www.gpsworld.com/wp-content/uploads/2013/06/isroi.jpg" width="275" height="193" /></a>The <a href="http://www.isro.org/index.aspx" target="_blank">Indian Space Research Organization (ISRO) Navigation Centre</a>, an important element of the Indian Regional Navigation Satellite System (IRNSS), was inaugurated May 28. The INC has been established at the Indian Deep Space Network complex at Byalalu, about 40 kilometers from Bangalore, India.</p>
<p>IRNSS, an independent navigation satellite system being developed by India, will have a constellation of seven satellites that enables its users to determine their location and time accurately. These satellites will be positioned in geostationary and inclined geosynchronous orbits 36,000 kilometers above the Earth&#8217;s surface. IRNSS coverage will extend over India and the southeast Asia region. The satellites are equipped with high-precision atomic clocks and continuously transmit navigation signals to users.</p>
<p>As the focal point of many critical operations of IRNSS, the ISRO Navigation Centre (INC) is responsible for providing the time reference, generation of navigation messages, and monitoring and control of ground facilities including ranging stations of IRNSS. It hosts several key technical facilities for supporting various navigation functions.</p>
<p>Key to the navigation support is the time reference to which all ground systems and the satellite clocks are synchronized. This time reference is generated by the high-precision timing facility located at INC. This timing facility is equipped with high-stability, high-precision atomic clocks to provide stable and continuous time reference to the navigation system.</p>
<p>IRNSS will have a network of 21 ranging stations geographically distributed primarily across India. They provide data for the orbit determination of IRNSS satellites and monitoring of the navigation signals. The data from the ranging/monitoring stations is sent to the data processing facility at INC where it is processed to generate the navigation messages. The navigation messages are then transmitted from INC to IRNSS satellites through the spacecraft control facility at Hassan/Bhopal. The data processing and storage facilities at INC enable swift processing of data and support its systematic storage.</p>
<p>INC is connected to the ranging stations and to the satellite control facilities through two highly reliable dedicated communication networks consisting of satellite and terrestrial links. The hub for the satellite communication links is hosted at INC.</p>
<p>The INC was inaugurated by V. Narayanasamy, minister of state in the Indian prime minister&#8217;s office. Speaking on the occasion, Narayanasamy said he appreciated the commitment and dedication of Indian space scientists in realizing the objectives of the country&#8217;s space programme. The minister also gave away various awards instituted by Astronautical Society of India (ASI) and ISRO.</p>
]]></content:encoded>
			<wfw:commentRss>http://www.gpsworld.com/navigation-center-for-indias-sat-nav-system-inaugurated/feed/</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Innovation: GNSS Spoofing Detection: Correlating Carrier Phase with Rapid Antenna Motion</title>
		<link>http://www.gpsworld.com/innovation-gnss-spoofing-detection-correlating-carrier-phase-with-rapid-antenna-motion/?utm_source=rss&#038;utm_medium=rss&#038;utm_campaign=innovation-gnss-spoofing-detection-correlating-carrier-phase-with-rapid-antenna-motion</link>
		<comments>http://www.gpsworld.com/innovation-gnss-spoofing-detection-correlating-carrier-phase-with-rapid-antenna-motion/#comments</comments>
		<pubDate>Sat, 01 Jun 2013 17:40:47 +0000</pubDate>
		<dc:creator>Richard Langley</dc:creator>
				<category><![CDATA[GNSS]]></category>
		<category><![CDATA[Innovation]]></category>
		<category><![CDATA[Richard B. Langley]]></category>
		<category><![CDATA[meaconing]]></category>
		<category><![CDATA[spoofing]]></category>

		<guid isPermaLink="false">http://www.gpsworld.com/?p=21495</guid>
		<description><![CDATA[By Mark L. Psiaki with Steven P. Powell and Brady W. O’Hanlon IT’S A HOSTILE (ELECTRONIC) WORLD OUT THERE, PEOPLE. Our wired and radio-based communication systems are constantly under attack from evil doers. We are all familiar with computer viruses and worms hiding in malicious software or malware distributed over the Internet or by infected [...]]]></description>
				<content:encoded><![CDATA[<p><em>By Mark L. Psiaki with Steven P. Powell and Brady W. O’Hanlon</em></p>
<div id="attachment_730" class="wp-caption alignright" style="width: 129px"><a href="http://www.gpsworld.com/wp-content/uploads/2012/09/Langley-INTRO-T.jpg"><img class="size-full wp-image-730" alt="INNOVATION INSIGHTS with Richard Langley" src="http://www.gpsworld.com/wp-content/uploads/2012/09/Langley-INTRO-T.jpg" width="119" height="150" /></a><p class="wp-caption-text">INNOVATION INSIGHTS with Richard Langley</p></div>
<p><strong>IT’S A HOSTILE (ELECTRONIC) WORLD OUT THERE, PEOPLE.</strong> Our wired and radio-based communication systems are constantly under attack from evil doers. We are all familiar with computer viruses and worms hiding in malicious software or malware distributed over the Internet or by infected USB flash drives. Trojan horses are particularly insidious. These are programs concealing harmful code that can lead to many undesirable effects such as deleting a user’s files or installing additional harmful software. Such programs pass themselves off as benign, just like the “gift” the Greeks delivered to the Trojans as reported in Virgil’s Aeneid. This was a very early example of spoofing. Spoofing of Internet Protocol (IP) datagrams is particularly prevalent. They contain forged source IP addresses with the purpose of concealing the identity of the sender or impersonating another computing system.</p>
<p>To spoof someone or something is to deceive or hoax, passing off a deliberately fabricated falsehood made to masquerade as truth. The word “spoof” was introduced by the English stage comedian Arthur Roberts in the late 19th century. He invented a game of that name, which involved trickery and nonsense. Now, the most common use of the word is as a synonym for parody or satirize — rather benign actions. But it is the malicious use of spoofing that concerns users of electronic communications.</p>
<p>And it is not just wired communications that are susceptible to spoofing. Communications and other services using radio waves are, in principle, also spoofable. One of the first uses of radio-signal spoofing was in World War I when British naval shore stations sent transmissions using German ship call signs. In World War II, spoofing became an established military tactic and was extended to radar and navigation signals. For example, German bomber aircraft navigated using radio signals transmitted from ground stations in occupied Europe, which the British spoofed by transmitting similar signals on the same frequencies. They coined the term “meaconing” for the interception and rebroadcast of navigation signals (meacon = m(islead)+(b)eacon).</p>
<p>Fast forward to today. GPS and other GNSS are also susceptible to meaconing. From the outset, the GPS P code, intended for use by military and other so-called authorized users, was designed to be encrypted to prevent straightforward spoofing. The anti-spoofing is implemented using a secret “W” encryption code, resulting in the P(Y) code. The C/A code and the newer L2C and L5 codes do not have such protection; nor, for the most part, do the civil codes of other GNSS. But, it turns out, even the P(Y) code is not fully protected from sophisticated meaconing attacks.</p>
<p>So, is there anything that military or civil GNSS users can do, then, to guard against their receivers being spoofed by sophisticated false signals? In this month’s column, we take a look at a novel, yet relatively easily implemented technique that enables users to detect and sequester spoofed signals. It just might help make it a safer world for GNSS positioning, navigation, and timing.</p>
<hr />
<h6>“Innovation” is a regular feature that discusses advances in GPS technology andits applications as well as the fundamentals of GPS positioning. The column is coordinated by Richard Langley of the Department of Geodesy and Geomatics Engineering, University of New Brunswick. He welcomes comments and topic ideas. To contact him, see the “Contributing Editors” section on page 4.</h6>
<hr />
<p>The radionavigation community has known about the dangers of GNSS spoofing for a long time, as highlighted in the 2001 Volpe Report (see Further Reading). Traditional receiver autonomous integrity monitoring (RAIM) had been considered a good spoofing defense. It assumes a dumb spoofer whose false signal produces a random pseudorange and large navigation solution residuals. The large errors are easy to detect, and given enough authentic signals, the spoofed signal(s) can be identified and ignored.</p>
<p>That spoofing model became obsolete at The Institute of Navigation’s GNSS 2008 meeting. Dr. Todd Humphreys introduced a new receiver/spoofer that could simultaneously spoof all signals in a self-consistent way undetectable to standard RAIM techniques. Furthermore, it could use its GNSS reception capabilities and its known geometry relative to the victim to overlay the false signals initially on top of the true ones. Slowly it could capture the receiver tracking loops by raising the spoofer power to be slightly larger than that of the true signals, and then it could drag the victim receiver off to false, but believable, estimates of its position, time, or both.</p>
<p>Two of the authors of this article contributed to Humphreys’ initial developments. There was no intention to help bad actors deceive GNSS user equipment (UE). Rather, our goal was to field a formidable “Red Team” as part of a “Red Team/Blue Team” (foe/friend) strategy for developing advanced “Blue Team” spoofing defenses.</p>
<p>This seemed like a fun academic game until mid-December 2011, when news broke that the Iranians had captured a highly classified Central Intelligence Agency drone, a stealth Lockheed Martin RQ-170 Sentinel, purportedly by spoofing its GPS equipment. Given our work in spoofing and detection, this event caused quite a stir in our Cornell University research group, in Humphreys’ University of Texas at Austin group, and in other places. The editor of this column even got involved in our extensive e-mail correspondence. Two key questions were: Wouldn’t a classified spy drone be equipped with a Selective Availability Anti-Spoofing Module (SAASM) receiver and, therefore, not be spoofable? Isn’t it difficult to knit together a whole sequence of false GPS position fixes that will guide a drone to land in a wrong location? These issues, when coupled with apparent inconsistencies in the Iranians’ story and visible damage to the drone, led us to discount the spoofing claim.</p>
<p><b>Developing a New Spoofing Defense</b></p>
<p>My views about the Iranian claims changed abruptly in mid-April 2012. Todd Humphreys phoned me about an upcoming test of GPS jammers, slated for June 2012 at White Sands Missile Range (WSMR), New Mexico. The Department of Homeland Security (DHS) had already spent months arranging these tests, but Todd revealed something new in that call: He had convinced the DHS to include a spoofing test that would use his latest “Red Team” device. The goal would be to induce a small GPS-guided unmanned aerial vehicle (UAV), in this case a helicopter, to land when it was trying to hover. “Wow”, I thought. “This will be a mini-replication of what the Iranians claimed to have done to our spy drone, and I’m sure that Todd will pull it off. I want to be there and see it.” Cornell already had plans to attend to test jammer tracking and geolocation, but we would have to come a day early to see the spoofing “fun” — if we could get permission from U.S. Air Force 746th Test Squadron personnel at White Sands.</p>
<p>The implications of the UAV test bounced around in my head that evening and the next morning on my seven-mile bike commute to work. During that ride, I thought of a scenario in which the Iranians might have mounted a meaconing attack against a SAASM-equipped drone. That is, they might possibly have received and re-broadcast the wide-band P(Y) code in a clever way that could have nudged the drone off course and into a relatively soft landing on Iranian territory.</p>
<p>In almost the next moment, I conceived a defense against such an attack. It involves small antenna motions at a high frequency, the measurement of corresponding carrier-phase oscillations, and the evaluation of whether the motions and phase oscillations are more consistent with spoofed signals or true signals. This approach would yield a good defense for civilian and military receivers against both spoofing and meaconing attacks. The remainder of this article describes this defense and our efforts to develop and test it.</p>
<p>It is one thing to conceive an idea, maybe a good idea. It is quite another thing to bring it to fruition. This idea seemed good enough and important enough to “birth” the conception. The needed follow-up efforts included two parts, one theoretical and the other experimental.</p>
<p>The theoretical work involved the development of signal models, hypothesis tests, analyses, and software. It culminated in analysis and truth-model simulation results, which showed that the system could be very practical, using only centimeters of motion and a fraction of a second of data to reliably differentiate between spoofing attacks and normal GNSS operation.</p>
<p>Theories and analyses can contain fundamental errors, or overlooked real-world effects can swamp the main theoretical effect. Therefore, an experimental prototype was quickly conceived, developed, and tested. It consisted of a very simple antenna-motion system, an RF data-recording device, and after-the-fact signal processing. The signal processing used Matlab to perform the spoofing detection calculations after using a C-language software radio to perform standard GPS acquisition and tracking.</p>
<p>Tests of the non-spoofed case could be conducted anywhere outdoors. Our initial tests occurred on a Cornell rooftop in Ithaca, New York. Tests of the spoofed case are harder. One cannot transmit live spoofing signals except with special permission at special times and in special places, for example, at WSMR in the upcoming June tests. Fortunately, the important geometric properties of spoofed signals can be simulated by using GPS signal reception at an outdoor antenna and re-radiation in an anechoic chamber from a single antenna. Such a system was made available to us by the NASA facility at Wallops Island, Virginia, and our simulated spoofed-case testing occurred in late April of last year. All of our data were processed before mid-May, and they provided experimental confirmation of our system’s efficacy. The final results were available exactly three busy weeks after the initial conception.</p>
<p>Although we were convinced about our new system, we felt that the wider GNSS community would like to see successful tests against live-signal attacks by a real spoofer. Therefore, we wanted very much to bring our system to WSMR for the June 2012 spoofing attack on the drone. We could set up our system near the drone so that it would be subject to the same malicious signals, but without the need to mount our clumsy prototype on a compact UAV helicopter. We were concerned, however, about the possibility of revealing our technology before we had been able to apply for patent protection. After some hesitation and discussions with our licensing and technology experts, we decided to bring our system to the WSMR test, but with a physical cover to keep it secret. The cover consisted of a large cardboard box, large enough to accommodate the needed antenna motions. The WSMR data were successfully collected using this method. Post-processing of the data demonstrated very reliable differentiation between spoofed and non-spoofed cases under live-signal conditions, as will be described in subsequent sections of this article.</p>
<p><b>System Architecture and Prototype</b></p>
<p>The components and geometry of one possible version of this system are shown in FIGURE 1. The figure shows three of the GNSS satellites whose signals would be tracked in the non-spoofed case: satellites <i>j</i>-1, <i>j</i>, and <i>j</i>+1. It also shows the potential location of a spoofer that could send false versions of the signals from these same satellites. The spoofer has a single transmission antenna. Satellites <i>j</i>-1, <i>j</i>, and <i>j</i>+1 are visible to the receiver antenna, but the spoofer could “hijack” the receiver’s tracking loops for these signals so that only the false spoofed versions of these signals would be tracked by the receiver.</p>
<div id="attachment_21505" class="wp-caption alignnone" style="width: 586px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/In-Fig1.jpg"><img class=" wp-image-21505 " alt="Figure 1. Spoofing detection antenna articulation system geometry relative to base mount, GNSS satellites, and potential spoofer." src="http://www.gpsworld.com/wp-content/uploads/2013/06/In-Fig1.jpg" width="576" height="352" /></a><p class="wp-caption-text">Figure 1. Spoofing detection antenna articulation system geometry relative to base mount, GNSS satellites, and potential spoofer.</p></div>
<p>The receiver antenna mount enables its phase center to be moved with respect to the mounting base. In Figure 1, this motion system is depicted as an open kinematic chain consisting of three links with ball joints. This is just one example of how a system can be configured to allow antenna motion. Spoofing detection can work well with just one translational degree of freedom, such as a piston-like up-and-down motion that could be provided by a solenoid operating along the <i>z<sub>a</sub></i> articulation axis. It would be wise to cover the motion system with an optically opaque radome, if possible, to prevent a spoofer from defeating this system by sensing the high-frequency antenna motions and spoofing their effects on carrier phase.</p>
<p>Suppose that the antenna articulation time history in its local body-fixed (<i>x<sub>a</sub></i>, <i>y<sub>a</sub></i>, <i>z<sub>a</sub></i>) coordinate system is <b><i>b</i></b><sub>a</sub>(<i>t</i>). Then the received carrier phases are sensitive to the projections of this motion onto the line-of-sight (LOS) directions of the received signals. These projections are along  <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj1.jpg"><img class="alignnone  wp-image-21612" alt="Eq-rj1" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj1.jpg" width="22" height="16" /></a>, <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg"><img class="alignnone  wp-image-21611" alt="Eq-rj" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg" width="14" height="16" /></a>, and  <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-r-j+1.jpg"><img class="alignnone  wp-image-21610" alt="Eq-r-j+1" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-r-j+1.jpg" width="25" height="17" /></a> in the non-spoofed case, with <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg"><img alt="Eq-rj" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg" width="14" height="16" /></a>  being the known unit direction vector from the <i>j</i>th GNSS satellite to the nominal antenna location. In the spoofed case, the projections are all along <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rsp.jpg"><img class="alignnone  wp-image-21613" alt="Eq-rsp" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rsp.jpg" width="22" height="15" /></a>, regardless of which signal is being spoofed, with <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rsp.jpg"><img class="alignnone  wp-image-21613" alt="Eq-rsp" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rsp.jpg" width="22" height="15" /></a> being the unknown unit direction vector from the spoofer to the victim antenna. Thus, there will be differences between the carrier-phase responses of the different satellites in the non-spoofed case, but these differences will vanish in the spoofed case. This distinction lies at the heart of the new spoofing detection method. Given that a good GNSS receiver can easily distinguish quarter-cycle carrier-phase variations, it is expected that this system will be able to detect spoofing using antenna motions as small as 4.8 centimeters, that is, a quarter wavelength of the GPS L1 signal.</p>
<p>The UE receiver and spoofing detection block in Figure 1 consists of a standard GNSS receiver, a means of inputting the antenna motion sensor data, and additional signal processing downstream of the standard GNSS receiver operations. The latter algorithms use as inputs the beat carrier-phase measurements from a standard phase-locked loop (PLL).</p>
<p>It may be necessary to articulate the antenna at a frequency nearly equal to the bandwidth of the PLL (say, at 1 Hz or higher). In this case, special post-processing calculations might be required to reconstruct the high-frequency phase variations accurately before they can be used to detect spoofing. The needed post-processing uses the in-phase and quadrature accumulations of a phase discriminator to reconstruct the noisy phase differences between the true signal and the PLL numerically controlled oscillator (NCO) signal. These differences are added to the NCO phases to yield the full high-bandwidth variations.</p>
<p>We implemented the first prototype of this system with one-dimensional antenna motion by mounting its patch antenna on a cantilevered beam. It is shown in FIGURE 2. Motion is initiated by pulling on the string shown in the upper left-hand part of the figure. Release of the string gives rise to decaying sinusoidal oscillations that have a frequency of about 2 Hz.</p>
<div id="attachment_21506" class="wp-caption alignnone" style="width: 586px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/In-Fig2.jpg"><img class=" wp-image-21506 " alt="Figure 2. Antenna articulation system for first prototype spoofing detector tests: a cantilevered beam that allows single-degree-of-freedom antenna phase-center vibration along a horizontal axis." src="http://www.gpsworld.com/wp-content/uploads/2013/06/In-Fig2.jpg" width="576" height="325" /></a><p class="wp-caption-text">Figure 2. Antenna articulation system for first prototype spoofing detector tests: a cantilevered beam that allows single-degree-of-freedom antenna phase-center vibration along a horizontal axis.</p></div>
<p>The remainder of the prototype system consisted of a commercial-off-the-shelf RF data recording device, off-line software receiver code, and off-line spoofing detection software. The prototype system lacked an antenna motion sensor. We compensated for this omission by implementing additional signal-processing calculations. They included off-line parameter identification of the decaying sinusoidal motions coupled with estimation of the oscillations’ initial amplitude and phase for any given detection.</p>
<p>This spoofing detection system is not the first to propose the use of antenna motion to uncover spoofing, and it is related to techniques that rely on multiple antennas. The present system makes three new contributions to the art of spoofing detection: First, it clearly explains why the measured carrier phases from a rapidly oscillating antenna provide a good means to detect spoofing. Second, it develops a precise spoofing detection hypothesis test for a moving-antenna system. Third, it demonstrates successful spoofing detection against live-signal attacks by a “Humphreys-class” spoofer.</p>
<p><b>Signal Model Theory and Verification</b></p>
<p>The spoofing detection test relies on mathematical models of the response of beat carrier phase to antenna motion. Reasonable models for the non-spoofed and spoofed cases are, respectively:</p>
<p><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-1b1.jpg"><img class="alignnone  wp-image-21615" alt="Eq-1b" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-1b1.jpg" width="569" height="47" /></a>  (1a)</p>
<p><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-1a.jpg"><img alt="Eq-1a" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-1a.jpg" width="571" height="46" /></a>(1b)</p>
<p>where <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-0jk.jpg"><img class="alignnone  wp-image-21607" alt="Eq-0jk" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-0jk.jpg" width="18" height="25" /></a> is the received (negative) beat carrier phase of the authentic or spoofed satellite-<i>j</i> signal at the <i>k</i>th sample time <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tjmk.jpg"><img class="alignnone  wp-image-21614" alt="Eq-tjmk" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tjmk.jpg" width="29" height="20" /></a> . The three-by-three direction cosines matrix <b><i>A</i></b> is the transformation from the reference system, in which the direction vectors <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg"><img class="alignnone  wp-image-21611" alt="Eq-rj" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg" width="14" height="16" /></a>  and <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rsp.jpg"><img class="alignnone  wp-image-21613" alt="Eq-rsp" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rsp.jpg" width="22" height="15" /></a> are defined, to the local body-axis system, in which the antenna motion <b><i>b</i></b><i><sub>a</sub></i>(<i>t</i>) is defined. <em>λ</em> is the nominal carrier wavelength. The terms involving the unknown polynomial coefficients <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg"><img class="alignnone  wp-image-21619" alt="Eq-Bj0" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg" width="18" height="21" /></a>, <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg"><img class="alignnone  wp-image-21618" alt="Eq-Bj1" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg" width="19" height="22" /></a> , and <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg"><img class="alignnone  wp-image-21617" alt="Eq-Bj2" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg" width="22" height="22" /></a> model other low-frequency effects on carrier phase, including satellite motion, UE motion if its antenna articulation system is mounted on a vehicle, and receiver clock drift. The term <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-nj0k.jpg"><img class="alignnone  wp-image-21621" alt="Eq-nj0k" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-nj0k.jpg" width="23" height="26" /></a> is the receiver phase noise. It is assumed to be a zero-mean, Gaussian, white-noise process whose variance depends on the receiver carrier-to-noise-density ratio and the sample/accumulation frequency.</p>
<p>If the motion of the antenna is one-dimensional, then <b><i>b</i></b><i><sub>a</sub></i>(<i>t</i>) takes the form <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-ba1.jpg"><img class="alignnone  wp-image-21623" alt="Eq-ba1" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-ba1.jpg" width="92" height="22" /></a>, with <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-ba.jpg"><img class="alignnone  wp-image-21622" alt="Eq-ba" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-ba.jpg" width="16" height="22" /></a> being the articulation direction in body-axis coordinates and <i>r</i><i><sub>a</sub></i>(<i>t</i>) being a known scalar antenna deflection amplitude time history. If one defines the articulation direction in reference coordinates as <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-ra.jpg"><img class="alignnone  wp-image-21624" alt="Eq-ra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-ra.jpg" width="73" height="29" /></a> , then the carrier-phase models in Equations (1a) and (1b) become</p>
<p><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-2a.jpg"><img class="alignnone  wp-image-21625" alt="Eq-2a" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-2a.jpg" width="564" height="47" /></a>   (2a)</p>
<p><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-2b.jpg"><img class="alignnone  wp-image-21626" alt="Eq-2b" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-2b.jpg" width="569" height="44" /></a>  (2b)</p>
<p>There is one important feature of these models for purposes of spoofing detection. In the non-spoofed case, the term that models the effects of antenna motion varies between GPS satellites because the <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg"><img class="alignnone  wp-image-21611" alt="Eq-rj" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg" width="14" height="16" /></a> direction vector varies with <i>j</i>. The spoofed case lacks variation between the satellites because the one spoofer direction <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rsp.jpg"><img alt="Eq-rsp" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rsp.jpg" width="22" height="15" /></a> replaces <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg"><img alt="Eq-rj" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg" width="14" height="16" /></a> for all of the spoofed satellites. This becomes clear when one compares the first terms on the right-hand sides of Eqsuations (1a) and (1b) for the 3-D motion case and on the right-hand sides of Equations (2a) and (2b) for the 1-D case.</p>
<p>The carrier-phase time histories in FIGURES 3 and 4 illustrate this principle. These data were collected at WSMR using the prototype antenna motion system of Figure 2. The carrier-phase time histories have been detrended by estimating the <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg"><img alt="Eq-Bj0" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg" width="18" height="21" /></a>, <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg"><img alt="Eq-Bj1" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg" width="19" height="22" /></a> , and <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg"><img alt="Eq-Bj2" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg" width="22" height="22" /></a> coefficients in Equations (2a) and (2b) and subtracting off their effects prior to plotting. In Figure 3, all eight satellite signals exhibit similar decaying sinusoid time histories, but with differing amplitudes and some of them with sign changes. This is exactly what is predicted by the 1-D non-spoofed model in Equation (2a). All seven spoofed signals in Figure 4, however, exhibit identical decaying sinusoidal oscillations because the <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rsp-tra.jpg"><img alt="Eq-rsp-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rsp-tra.jpg" width="53" height="22" /></a> term in Equation (2b) is the same for all of them.</p>
<div id="attachment_21507" class="wp-caption alignnone" style="width: 610px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/In-Fig3.jpg"><img class=" wp-image-21507 " alt="Figure 3. Detrended carrier-phase data from multiple satellites for a typical non-spoofed case using a 1-D antenna articulation system." src="http://www.gpsworld.com/wp-content/uploads/2013/06/In-Fig3.jpg" width="600" height="386" /></a><p class="wp-caption-text">Figure 3. Detrended carrier-phase data from multiple satellites for a typical non-spoofed case using a 1-D antenna articulation system.</p></div>
<p>&nbsp;</p>
<div id="attachment_21508" class="wp-caption alignnone" style="width: 610px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/In-Fig4.jpg"><img class=" wp-image-21508 " alt="Figure 4. Multiple satellites’ detrended carrier-phase data for a typical spoofed case using a 1-D antenna articulation system." src="http://www.gpsworld.com/wp-content/uploads/2013/06/In-Fig4.jpg" width="600" height="396" /></a><p class="wp-caption-text">Figure 4. Multiple satellites’ detrended carrier-phase data for a typical spoofed case using a 1-D antenna articulation system.</p></div>
<p>As an aside, an interesting feature of Figure 3 is its evidence of the workings of the prototype system. The ramping phases of all the signals from <i>t</i> = 0.4 seconds to <i>t</i> = 1.4 seconds correspond to the initial pull on the string shown in Figure 2, and the steady portion from <i>t</i> = 1.4 seconds to <i>t</i> = 2.25 seconds represents a period when the string was held fixed prior to release.<b><br />
</b></p>
<p><b>Spoofing Detection Hypothesis Test</b></p>
<p>A hypothesis test can precisely answer the question of which model best fits the observed data: Does carrier-phase sameness describe the data, as in Figure 4? Then the receiver is being spoofed. Alternatively, is carrier-phase differentness more reasonable, as per Figure 3? Then the signals are trustworthy.</p>
<p>A hypothesis test can be developed for any batch of carrier-phase data that spans a sufficiently rich antenna motion profile <b><i>b</i></b><i><sub>a</sub></i>(<i>t</i>) or <em>ρ</em><i><sub>a</sub></i>(<i>t</i>). The profile must include high-frequency motions that cannot be modeled by the  <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg"><img alt="Eq-Bj0" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg" width="18" height="21" /></a>, <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg"><img alt="Eq-Bj1" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg" width="19" height="22" /></a> , and <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg"><img alt="Eq-Bj2" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg" width="22" height="22" /></a>quadratic polynomial terms in Equations (1a)-(2b); otherwise the detection test will lose all of its power. A motion profile equal to one complete period of a sine wave has the needed richness.</p>
<p>Suppose one starts with a data batch that is comprised of carrier-phase time histories for <i>L</i> different GNSS satellites: <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-0jk.jpg"><img class="alignnone  wp-image-21607" alt="Eq-0jk" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-0jk.jpg" width="18" height="25" /></a> for samples <i>k</i> = 1, &#8230;, <i>M</i><i><sub>j</sub></i> and for satellites <i>j</i> = 1,&#8230;, <i>L</i>. A standard hypothesis test develops two probability density functions for these data, one conditioned on the null hypothesis of no spoofing, <i>H</i><sub>0</sub>, and the other conditioned on the hypothesis of spoofing, <i>H</i><sub>1</sub>.  The Neyman-Pearson lemma (see Further Reading) proves that the optimal hypothesis test statistic equals the ratio of these two probability densities. Unfortunately, the required probability densities depend on additional unknown quantities. In the 1-D motion case, these unknowns include the <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg"><img alt="Eq-Bj0" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg" width="18" height="21" /></a>, <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg"><img alt="Eq-Bj1" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg" width="19" height="22" /></a> , and <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg"><img alt="Eq-Bj2" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg" width="22" height="22" /></a> coefficients, the dot product <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rsp-tra.jpg"><img class="alignnone  wp-image-21627" alt="Eq-rsp-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rsp-tra.jpg" width="53" height="22" /></a>, and the direction <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img class="alignnone  wp-image-21628" alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="15" height="21" /></a>  if one assumes that the UE attitude is unknown. A true Neyman-Pearson test would hypothesize <i>a priori</i> distributions for these unknown quantities and integrate their dependencies out of the two joint probability distributions. Our sub-optimum test optimally estimates relevant unknowns for each hypothesis based on the carrier-phase data, and it uses these estimates in the Neyman-Pearson probability density ratio. Although sub-optimal as a hypothesis test, this approach is usually effective, and it is easier to implement than the integration approach in the present case.</p>
<p>Consider the case of 1-D antenna articulation and unknown UE attitude. Maximum-likelihood calculations optimally estimate the nuisance parameters  <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg"><img alt="Eq-Bj0" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg" width="18" height="21" /></a>, <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg"><img alt="Eq-Bj1" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg" width="19" height="22" /></a> , and <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg"><img alt="Eq-Bj2" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg" width="22" height="22" /></a>  for <i>j</i> = 1, &#8230;, <i>L</i> for both hypotheses along with the unit vector <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="15" height="21" /></a> for the non-spoofed hypothesis, or the scalar dot product <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-nsix.jpg"><img class="alignnone  wp-image-21629" alt="Eq-nsix" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-nsix.jpg" width="79" height="23" /></a> for the spoofed hypothesis. The estimation calculations for each hypothesis minimize the negative natural logarithm of the corresponding conditional probability density. Because  <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg"><img alt="Eq-Bj0" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg" width="18" height="21" /></a>, <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg"><img alt="Eq-Bj1" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg" width="19" height="22" /></a> , and <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg"><img alt="Eq-Bj2" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg" width="22" height="22" /></a> enter the resulting cost functions quadratically, their optimized values can be computed as functions of the other unknowns, and they can be substituted back into the costs. This part of the calculation amounts to a batch high-pass filter of both the antenna motion and the carrier-phase response.</p>
<p>The remaining optimization problems take, under the non-spoofed hypothesis, the form:</p>
<p>find:      <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img class="alignnone  wp-image-21628" alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="17" height="24" /></a>    (3a)</p>
<p>to minimize:       <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jnonsp.jpg"><img class="alignnone  wp-image-21630" alt="Eq-Jnonsp" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jnonsp.jpg" width="272" height="54" /></a>  (3b)</p>
<p>subject to:             <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rasmall.jpg"><img class="alignnone  wp-image-21634" alt="Eq-rasmall" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rasmall.jpg" width="89" height="26" /></a>   (3c)</p>
<p>and, under the spoofed hypothesis, the form:</p>
<p>find:     <strong> <em>η</em></strong>    (4a)</p>
<p>to minimize:   <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jspn.jpg"><img class="alignnone  wp-image-21632" alt="Eq-Jspn" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jspn.jpg" width="199" height="53" /></a>      (4b)</p>
<p>subject to:     <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-111.jpg"><img class="alignnone  wp-image-21631" alt="Eq-111" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-111.jpg" width="83" height="20" /></a> .   (4c)</p>
<p>The coefficient <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj44.jpg"><img class="alignnone  wp-image-21637" alt="Eq-rj44" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj44.jpg" width="27" height="26" /></a> is a function of the deflections <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Pat.jpg"><img class="alignnone  wp-image-21638" alt="Eq-Pat" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Pat.jpg" width="74" height="27" /></a> for <i>k</i> = 1, &#8230;, <i>M</i><i><sub>j</sub></i>, and the non-homogenous term <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-zj4.jpg"><img class="alignnone  wp-image-21639" alt="Eq-zj4" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-zj4.jpg" width="20" height="27" /></a> is derived from the <i>j</i>th phase time history <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-0jk.jpg"><img class="alignnone  wp-image-21607" alt="Eq-0jk" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-0jk.jpg" width="18" height="25" /></a> for <i>k</i> = 1, &#8230;, <i>M</i><i><sub>j</sub></i>. These two quantities are calculated during the  <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg"><img alt="Eq-Bj0" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg" width="18" height="21" /></a>, <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg"><img alt="Eq-Bj1" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg" width="19" height="22" /></a>, <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg"><img alt="Eq-Bj2" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg" width="22" height="22" /></a> optimization. The constraint in Equation (3c) forces the estimate of the antenna articulation direction to be unit-normalized. The constraint in Eq. (4c) ensures that <em>η</em> is a physically reasonable dot product.</p>
<p>The optimization problems in Equations (3a)-(3c) and (4a)-(4c) can be solved in closed form using techniques from the literature on constrained optimization, linear algebra, and matrix factorization. The optimal estimates of <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img class="alignnone  wp-image-21628" alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="17" height="24" /></a> and <em>η</em> can be used to define a spoofing detection statistic that equals the natural logarithm of the Neyman-Pearson ratio:</p>
<p><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-y-small.jpg"><img class="alignnone  wp-image-21640" alt="Eq-y-small" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-y-small.jpg" width="196" height="29" /></a>(5)</p>
<p>It is readily apparent that <em>γ</em> constitutes a reasonable test statistic: If the signal is being spoofed so that carrier-phase sameness is the best model, then <em>η</em><i><sub>opt</sub></i> will produce a small value of  <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jsp-n.jpg"><img class="alignnone  wp-image-21642" alt="Eq-Jsp-n" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jsp-n.jpg" width="59" height="22" /></a>because the spoofed-case cost function in Equation (4b) is consistent with carrier-phase sameness. The value of <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jnonsp-r.jpg"><img class="alignnone  wp-image-21641" alt="Eq-Jnonsp-r" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jnonsp-r.jpg" width="76" height="21" /></a>, however, will not be small because the plurality of  <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg"><img class="alignnone  wp-image-21611" alt="Eq-rj" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg" width="14" height="16" /></a> directions in Equation (3b) precludes the possibility that any <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img class="alignnone  wp-image-21628" alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="15" height="21" /></a> estimate will yield a small non-spoofed cost. Therefore, <em>γ</em> will tend to be a large negative number in the event of spoofing because <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jnonsp-r.jpg"><img alt="Eq-Jnonsp-r" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jnonsp-r.jpg" width="76" height="21" /></a> &gt;&gt; <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jsp-n.jpg"><img alt="Eq-Jsp-n" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jsp-n.jpg" width="59" height="22" /></a> is likely. In the non-spoofed case, the opposite holds true: <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-ropt.jpg"><img class="alignnone  wp-image-21645" alt="Eq-ropt" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-ropt.jpg" width="26" height="20" /></a>  will yield a small value of <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jnonsp-r.jpg"><img alt="Eq-Jnonsp-r" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jnonsp-r.jpg" width="76" height="21" /></a>, but no estimate of <em>η</em> will yield a small <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-jspn2.jpg"><img class="alignnone  wp-image-21646" alt="Eq-jspn2" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-jspn2.jpg" width="46" height="24" /></a>, and <em>γ</em> will be a large positive number because  <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jnonsp-r.jpg"><img alt="Eq-Jnonsp-r" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jnonsp-r.jpg" width="76" height="21" /></a>&lt;&lt; <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jsp-n.jpg"><img alt="Eq-Jsp-n" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Jsp-n.jpg" width="59" height="22" /></a>.</p>
<p>Therefore, a sensible spoofing detection test employs a detection threshold <em>γ</em><i><sub>th</sub></i> somewhere in the neighborhood of zero. The detection test computes a <em>γ</em> value based on the carrier-phase data, the antenna articulation time history, and the calculations in Equations (3a)-(5). It compares this <em>γ</em> to <em>γ</em><i><sub>th</sub></i>. If <em>γ</em> ≥ <em>γ</em><i><sub>th</sub></i>, then the test indicates that there is no spoofing. If <em>γ</em> &lt; <em>γ</em><i><sub>th</sub></i>, then a spoofing alert is issued.</p>
<p>The exact choice of <em>γ</em><i><sub>th</sub></i> is guided by an analysis of the probability of false alarm. A false alarm occurs if a spoofing attack is declared when there is no spoofing. The false-alarm probability is determined as a function of <em>γ</em><i><sub>th</sub></i> by developing a <em>γ</em> probability density function under the null hypothesis of no spoofing <i>p</i>(<em>γ</em>|<i>H</i><sub>0</sub>). The probability of false alarm equals the integral of <i>p</i>(<em>γ</em>|<i>H</i><sub>0</sub>) from <em>γ</em> = <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-infinity.jpg"><img class="alignnone  wp-image-21649" alt="Eq-infinity" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-infinity.jpg" width="23" height="10" /></a> to <em>γ</em> = <em>γ</em><i><sub>th</sub></i>. This integral relationship can be inverted to determine the <em>γ</em><i><sub>th</sub></i> threshold that yields a given prescribed false-alarm probability</p>
<p>A complication arises because <i>p</i>(<em>γ</em>|<i>H</i><sub>0</sub>) depends on unknown parameters, <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="15" height="21" /></a>  in the case of an unknown UE attitude and 1-D antenna motion. Although sub-optimal, a reasonable way to deal with the dependence of <i>p</i>(<em>γ</em>|<a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="15" height="21" /></a>,<i>H</i><sub>0</sub>) on <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="15" height="21" /></a> is to use the worst-case <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="15" height="21" /></a> for a given <em>γ</em><i><sub>th</sub></i>. The worst-case articulation direction <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rawc.jpg"><img class="alignnone  wp-image-21650" alt="Eq-rawc" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rawc.jpg" width="27" height="20" /></a> maximizes the <i>p</i>(<em>γ</em>|<a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="15" height="21" /></a>,<i>H</i><sub>0</sub>) false-alarm integral. It can be calculated by solving an optimization problem. This analysis can be inverted to pick <em>γ</em><i><sub>th</sub></i> so that the worst-case probability of false alarm equals some prescribed value. For most actual <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="15" height="21" /></a> values, the probability of false alarm will be lower than the prescribed worst case.</p>
<p>Given <em>γ</em><i><sub>th</sub></i>, the final needed analysis is to determine the probability of missed detection. This analysis uses the probability density function of <i>g</i> under the spoofed hypothesis, <i>p</i>(<em>γ</em>|<em>η</em>,<i>H</i><sub>1</sub>). The probability of missed detection is the integral of this function from <em>γ</em> = <em>γ</em><i><sub>th</sub></i> to <em>γ</em> = +<a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-infinity2.jpg"><img class="alignnone  wp-image-21651" alt="Eq-infinity2" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-infinity2.jpg" width="13" height="10" /></a>. The dependence of <i>p</i>(<em>γ</em>|<em>η</em>,<i>H</i><sub>1</sub>) on the unknown dot product <em>η</em> can be handled effectively, though sub-optimally, by determining the worst-case probability of false alarm. This involves an optimization calculation, which finds the worst-case dot product <em>η</em><i><sub>wc</sub></i> that maximizes the missed-detection probability integral. Again, most actual <em>η</em> values will yield lower probabilities of missed detection.</p>
<p>Note that the above-described analyses rely on approximations of the probability density functions <i>p</i>(<em>γ</em>|<a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="15" height="21" /></a>,<i>H</i><sub>0</sub>) and <i>p</i>(<em>γ</em>|<em>η</em>,<i>H</i><sub>1</sub>). The best approximations include dominant Gaussian terms plus small chi-squared or non-central chi-squared terms. It is difficult to analyze the chi-squared terms rigorously. Their smallness, however, makes the use of Gaussian approximations reasonable.</p>
<p>We have developed and evaluated several alternative formulations of this spoofing detection method. One is the case of full 3-D <b><i>b</i></b><i><sub>a</sub></i>(<i>t</i>) antenna motion with unknown UE attitude. The full direction cosines matrix <b><i>A</i></b> is estimated in the modified version of the non-spoofed optimal fit calculations of Equations (3a)-(3c), and the full spoofing direction vector <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-bsp.jpg"><img class="alignnone  wp-image-21652" alt="Eq-bsp" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-bsp.jpg" width="58" height="17" /></a> is estimated in the modified version of Equations (4a)-(4c). A different alternative allows the 1-D motion time history <em>ρ</em><i><sub>a</sub></i>(<i>t</i>) to have an unknown amplitude-scaling factor that must be estimated. This might be appropriate for a UAV drone with a wing-tip-mounted antenna if it induced antenna motions by dithering its ailerons. In fixed-based applications, as might be used by a financial institution, a cell-phone tower, or a power-grid monitor, the attitude would be known, which would eliminate the need to estimate <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="15" height="21" /></a> or <b><i>A</i></b> for the non-spoofed case.</p>
<p><b>Test Results</b></p>
<p>The initial tests of our concept involved generation of simulated truth-model carrier-phase data <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-0jk.jpg"><img class="alignnone  wp-image-21607" alt="Eq-0jk" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-0jk.jpg" width="18" height="25" /></a> using simulated <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg"><img alt="Eq-Bj0" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj0.jpg" width="18" height="21" /></a>, <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg"><img alt="Eq-Bj1" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj1.jpg" width="19" height="22" /></a> , and <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg"><img alt="Eq-Bj2" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-Bj2.jpg" width="22" height="22" /></a> polynomial coefficients, simulated satellite LOS direction vectors <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg"><img class="alignnone size-full wp-image-21611" alt="Eq-rj" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg" width="18" height="20" /></a> for the non-spoofed cases, a simulated true spoofer LOS direction <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rsp.jpg"><img class="alignnone size-full wp-image-21613" alt="Eq-rsp" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rsp.jpg" width="27" height="19" /></a> for the spoofed cases, and simulated antenna motions parameterized by <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img class="alignnone size-full wp-image-21628" alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="21" height="30" /></a> and <em>ρ</em><i><sub>a</sub></i>(<i>t</i>). Monte-Carlo analysis was used to generate many different batches of phase data with different random phase noise realizations in order to produce simulated histograms of the <i>p</i>(<em>γ</em>|<a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="15" height="21" /></a>, <i>H</i><sub>0</sub>) and <i>p</i>(<em>γ</em>|<em>η</em>,<i>H</i><sub>1</sub>) probability density functions  that are used in false-alarm and missed-detection analyses.</p>
<p>The truth-model simulations verified that the system is practical. A representative calculation used one cycle of an 8-Hz 1-D sinusoidal antenna oscillation with a peak-to-peak amplitude of 4.76 centimeters (exactly 1/4 of the L1 wavelength). The accumulation frequency was 1 kHz so that there were <i>M</i><i><sub>j</sub></i> = 125 carrier-phase measurements per satellite per data batch. The number of satellites was <i>L</i> = 6, their <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg"><img alt="Eq-rj" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-rj.jpg" width="18" height="20" /></a> LOS vectors were distributed to yield a geometrical dilution of precision of 3.5, and their carrier-to-noise-density ratios spanned the range 38.2 to 44.0 dB-Hz. The worst-case probability of a spoofing false alarm was set at 10<sup>-5</sup> and the corresponding worst-case probability of missed detection was 1.2 ´ 10<sup>-5</sup>. Representative non-worst-case probabilities of false alarm and missed detection were, respectively, 1.7 ´ 10<sup>-9</sup> and 1.1 ´ 10<sup>-6</sup>. These small numbers indicate that this is a very powerful test. Ten-thousand run Monte-Carlo simulations of the spoofed and non-spoofed cases verified the reasonableness of these probabilities and the reasonableness of the <i>p</i>(<em>γ</em>|<a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="15" height="21" /></a>, <i>H</i><sub>0</sub>) and <i>p</i>(<em>γ</em>|<em>η</em>,<i>H</i><sub>1</sub>) Gaussian approximations that had been used to derive them.</p>
<p>The live-signal tests bore out the truth-model simulation results. The only surprise in the live-signal tests was the presence of significant multipath, which was evidenced by received carrier amplitude oscillations that correlated with the antenna oscillations and whose amplitudes and phases varied among the different received GPS signals. As a verification that these oscillations were caused by multipath, the only live-signal data set without such amplitude oscillations was the one taken in the NASA Wallops anechoic chamber, where one would not expect to find multipath. The multipath, however, seems to have negligible impact on the efficacy of this spoofing detection system.</p>
<p>FIGURES 5 and 6 show the results of typical non-spoofed and spoofed cases from WSMR live-signal tests that took place on the evening of June 19–20, 2012. Each plot shows the spoofing detection statistic <em>γ</em> on the horizontal axis and various related probability density functions on the vertical axis. This statistic has been calculated using a modified test that includes the estimation of two additional unknowns: an antenna articulation scale factor <i>f</i> and a timing bias <i>t</i><sub>0</sub> for the decaying sinusoidal oscillation <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/eq-pa.jpg"><img class="alignnone  wp-image-21654" alt="eq-pa" src="http://www.gpsworld.com/wp-content/uploads/2013/06/eq-pa.jpg" width="224" height="25" /></a>. The damping ratio ζ and the undamped natural frequency <i>w</i><i><sub>n</sub></i> are known from prior system identification tests.</p>
<div id="attachment_21509" class="wp-caption alignnone" style="width: 610px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/In-Fig5.jpg"><img class=" wp-image-21509 " alt="Figure 5. Spoofing detection statistic, threshold, and related probability density functions for a typical non-spoofed case with live data." src="http://www.gpsworld.com/wp-content/uploads/2013/06/In-Fig5.jpg" width="600" height="378" /></a><p class="wp-caption-text">Figure 5. Spoofing detection statistic, threshold, and related probability density functions for a typical non-spoofed case with live data.</p></div>
<p>&nbsp;</p>
<div id="attachment_21510" class="wp-caption alignnone" style="width: 610px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/In-Fig6.jpg"><img class=" wp-image-21510 " alt="Figure 6. Performance of a typical spoofed case with live data: spoofing detection statistic, threshold, and related probability density functions." src="http://www.gpsworld.com/wp-content/uploads/2013/06/In-Fig6.jpg" width="600" height="400" /></a><p class="wp-caption-text">Figure 6. Performance of a typical spoofed case with live data: spoofing detection statistic, threshold, and related probability density functions.</p></div>
<p>The vertical dashed black line in each plot shows the actual value of <em>γ</em> as computed from the GPS data. There are three vertical dash-dotted magenta lines that lie almost on top of each other. They show the worst-case threshold values <em>γ</em><i><sub>th</sub></i> as computed for the optimal and ±2σ estimates of <i>t</i><sub>0</sub>: <i>t</i><sub>0<i>opt</i></sub>, <i>t</i><sub>0<i>opt</i></sub>+2σ<sub><i>t</i><sub>0<i>opt</i></sub></sub>, and <i>t</i><sub>0<i>opt</i></sub><i>-2σ<sub><i>t</i><sub>0<i>opt</i></sub></sub></i>. They have been calculated for a worst-case probability of false alarm equal to 10<sup>-6</sup>. An <i>ad hoc</i> method of compensating for the prototype system’s <i>t</i><sub>0</sub> uncertainty is to use the left-most vertical magenta line as the detection threshold <em>γ</em><i><sub><em></em>th</sub></i>. The vertical dashed black line lies very far to the right of all three vertical dash-dotted magenta lines in Figure 5, which indicates a successful determination that the signals are not being spoofed. In Figure 6, the situation is reversed. The vertical dashed black line lies well to the left of the three vertical dash-dotted magenta lines, and spoofing is correctly and convincingly detected.</p>
<p>These two figures also plot various relevant probability density functions. Consistent with the consideration of three possible values of the <i>t</i><sub>0</sub> motion timing estimate, these are plotted in triplets. The three dotted cyan probability density functions represent the worst-case non-spoofed situation, and the dash-dotted red probability functions represent the corresponding worst-case spoofed situations. Obviously, there is sufficient separation between these sets of probability density functions to yield a powerful detection test, as evidenced by the ability to draw the dash-dotted magenta detection thresholds in a way that clearly separates the red and cyan distributions. Further confirmation of good detection power is provided by the low worst-case probabilities of false alarm and missed detection, the latter metric being 1.6 ´ 10<sup>-6</sup> for the test in Figure 5 and 7 ´ 10<sup>-8</sup> for Figure 6.</p>
<p>The solid-blue distributions on the two plots correspond to the <em>η</em><i><sub>opt</sub></i> estimate and the spoofed assumption, which is somewhat meaningless for Figure 5, but meaningful for Figure 6. The dashed-green distributions are for the <a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg"><img class=" wp-image-21628 alignnone" alt="Eq-tra" src="http://www.gpsworld.com/wp-content/uploads/2013/06/Eq-tra.jpg" width="13" height="18" /></a> estimate under the non-spoofed assumption. The wide separations between the blue distributions and the green distributions in both figures clearly indicate that the worst-case false-alarm and missed-detection probabilities can be very conservative.</p>
<p>The detection test results in Figures 5 and 6 have been generated using the last full oscillation of the respective carrier-phase data, as in Figures 3 and 4, but applied to different data sets. In Figure 3, the last full oscillation starts at <i>t</i> = 3.43 seconds, and it starts at <i>t</i> = 2.11 seconds in Figure 4. The peak-to-peak amplitude of each last full oscillation ranged from 4-6 centimeters, and their periods were shorter than 0.5 seconds. It would have been possible to perform the detections using even shorter data spans had the mechanical oscillation frequency of the cantilevered antenna been higher.</p>
<p><b>Conclusions</b></p>
<p>In this article, we have presented a new method to detect spoofing of GNSS signals. It exploits the effects of intentional high-frequency antenna motion on the measured beat carrier phases of multiple GNSS signals. After detrending using a high-pass filter, the beat carrier-phase variations can be matched to models of the expected effects of the motion. The non-spoofed model predicts differing effects of the antenna motion for the different satellites, but the spoofed case yields identical effects due to a geometry in which all of the false signals originate from a single spoofer transmission antenna. Precise spoofing detection hypothesis tests have been developed by comparing the two models’ ability to fit the measured data.</p>
<p>This new GNSS spoofing detection technique has been evaluated using both Monte-Carlo simulation and live data. Its hypothesis test yields theoretical false-alarm probabilities and missed-detection probabilities on the order of 10<sup>-5</sup> or lower when working with typical numbers and geometries of available GPS signals and typical patch-antenna signal strengths. The required antenna articulation deflections are modest, on the order of 4-6 centimeters peak-to-peak, and detection intervals less than 0.5 seconds can suffice.</p>
<p>A set of live-signal tests at WSMR evaluated the new technique against a sophisticated receiver/spoofer, one that mimics all visible signals in a way that foils standard RAIM techniques. The new system correctly detected all of the attacks. These are the first known practical detections of live-signal attacks mounted against a civilian GNSS receiver by a dangerous new generation of spoofers.</p>
<p><b>Future Directions</b></p>
<p>This work represents one step in an on-going “Blue Team” effort to develop better defenses against new classes of GNSS spoofers. Planned future improvements include 1) the ability to use electronically synthesized antenna motion that eliminates the need for moving parts, 2) the re-acquisition of true signals after detection of spoofing, 3) the implementation of real-time prototypes using software radio techniques, and 4) the consideration of “Red-Team” counter-measures to this defense  and how the “Blue Team” could combat them; counter-measures such as high-frequency phase dithering of the spoofed signals or coordinated spoofing transmissions from multiple locations.</p>
<p><b>Acknowledgments</b></p>
<p>The authors thank the following people and organizations for their contributions to this effort:  The NASA Wallops Flight Facility provided access to their anechoic chamber. Robert Miceli, a Cornell graduate student, helped with data collection at that facility. Dr. John Merrill and the Department of Homeland Security arranged the live-signal spoofing tests. The U.S. Air Force 746th Test Squadron hosted the live-signal spoofing tests at White Sands Missile Range. Prof. Todd Humphreys and members of his University of Texas at Austin Radionavigation Laboratory provided live-signal spoofing broadcasts from their latest receiver/spoofer.</p>
<p><b>Manufacturers</b></p>
<p>The prototype spoofing detection data capture system used an <b>Antcom Corp.</b> (<i>www.antcom.com</i>) 2G1215A L1/L2 GPS antenna. It was connected to an <b>Ettus Research</b> (<i>www.ettus.com</i>) USRP (Universal Software Radio Peripheral) N200 that was equipped with the DBSRX2 daughterboard.</p>
<hr />
<p align="left"><em><b>MARK L. PSIAKI</b> is a professor in the Sibley School of Mechanical and Aerospace Engineering at Cornell University, Ithaca, New York. He received a B.A. in physics and M.A. and Ph.D. degrees in mechanical and aerospace engineering from Princeton University, Princeton, New Jersey. His research interests are in the areas of GNSS technology, applications, and integrity, spacecraft attitude and orbit determination, and general estimation, filtering, and detection.</em></p>
<p align="left"><em><b>STEVEN P. POWELL</b> is a senior engineer with the GPS and Ionospheric Studies Research Group in the Department of Electrical and Computer Engineering at Cornell University. He has M.S. and B.S. degrees in electrical engineering from Cornell University. He has been involved with the design, fabrication, testing, and launch activities of many scientific experiments that have flown on high altitude balloons, sounding rockets, and small satellites. He has designed ground-based and space-based custom GPS receiving systems primarily for scientific applications.</em></p>
<p align="left"><em><b>BRADY W. O’HANLON</b> is a graduate student in the School of Electrical and Computer Engineering at Cornell University. He received a B.S. in electrical and computer engineering from Cornell University. His interests are in the areas of GNSS technology and applications, GNSS security, and GNSS as a tool for space weather research.</em></p>
<h3>VIDEO</h3>
<p>Here is a video (in m4v format) of Cornell University&#8217;s antenna articulation system for the team&#8217;s first prototype spoofing detector tests.</p>
<p><iframe src="http://www.youtube.com/embed/78KGydgRDh0" height="315" width="420" allowfullscreen="" frameborder="0"></iframe></p>
<h3 align="left"><b>FURTHER READING</b></h3>
<p><b>• The Spoofing Threat and RAIM-Resistant Spoofers</b></p>
<p>“Status of Signal Authentication Activities within the GNSS Authentication and User Protection System Simulator (GAUPSS) Project” by O. Pozzobon, C. Sarto, A. Dalla Chiara, A. Pozzobon, G. Gamba, M. Crisci, and R.T. Ioannides, in <i>Proceedings of</i><i> </i><i>ION GNSS 2012</i>, the 25th International Technical Meeting of The Institute of Navigation, Nashville, Tennessee, September 18–21, 2012, pp. 2894-2900.</p>
<p>“<a href="http://www.gpsworld.com/defensesecurity-surveillanceassessing-spoofing-threat-3171/" target="_blank">Assessing the Spoofing Threat</a>” by T.E. Humphreys, P.M. Kintner, Jr., M.L. Psiaki, B.M. Ledvina, and B.W. O’Hanlon in <i>GPS World</i>, Vol. 20, No. 1, January 2009, pp. 28-38.</p>
<p><a href="http://www.navcen.uscg.gov/pdf/vulnerability_assess_2001.pdf" target="_blank"><i>Vulnerability Assessment of the Transportation Infrastructure Relying on the Global Positioning System – Final Report</i></a>. John A. Volpe National Transportation Systems Center, Cambridge, Massachusetts, August 29, 2001.</p>
<p>• <b>Moving-Antenna and Multi-Antenna Spoofing Detection</b></p>
<p><em>“</em>Robust Joint Multi-Antenna Spoofing Detection and Attitude Estimation by Direction Assisted Multiple Hypotheses RAIM<em>” by M. Meurer, A. Konovaltsev, M. Cuntz, and C. Hattich, in </em><i>Proceedings of</i><i> </i><i>ION GNSS 2012</i>, the 25th International Technical Meeting of The Institute of Navigation, Nashville, Tennessee, September 18–21, 2012,<em> pp. 3007-3016.</em></p>
<p>“GNSS Spoofing Detection for Single Antenna Handheld Receivers” by J. Nielsen, A. Broumandan, and G. Lachapelle in <i>Navigation</i>, Vol. 58, No. 4, Winter 2011, pp. 335-344.</p>
<p><em>• </em><b>Alternate Spoofing Detection Strategies</b></p>
<p>“Who’s Afraid of the Spoofer? GPS/GNSS Spoofing Detection via Automatic Gain Control (AGC)” by D.M. Akos, in <i>Navigation</i>, Vol. 59, No. 4, Winter 2012-2013, pp. 281-290.</p>
<p>“Civilian GPS Spoofing Detection based on Dual-Receiver Correlation of Military Signals” by M.L. Psiaki, B.W. O’Hanlon, J.A. Bhatti, D.P. Shepard, and T.E. Humphreys in <i>Proceedings of ION GNSS 2011</i>, the 24th International Technical Meeting of The Institute of Navigation, Portland, Oregon, September 19–23, 2011, pp. 2619-2645.</p>
<p>• <b>Statistical Hypothesis Testing</b></p>
<p><i>Fundamentals of Statistical Signal Processing, Volume II: Detection Theory</i> by S. Kay, published by Prentice Hall, Upper Saddle River, New Jersey,1998.</p>
<p><i>An Introduction to Signal Detection and Estimation</i> by H.V. Poor, 2nd edition, published by Springer-Verlag, New York, 1994.<b></b></p>
<p>Video (in m4v format) of Cornell University&#8217;s antenna articulation system for their first prototype spoofing detector tests.</p>
]]></content:encoded>
			<wfw:commentRss>http://www.gpsworld.com/innovation-gnss-spoofing-detection-correlating-carrier-phase-with-rapid-antenna-motion/feed/</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Following the Team into Danger</title>
		<link>http://www.gpsworld.com/following-the-team-into-danger/?utm_source=rss&#038;utm_medium=rss&#038;utm_campaign=following-the-team-into-danger</link>
		<comments>http://www.gpsworld.com/following-the-team-into-danger/#comments</comments>
		<pubDate>Sat, 01 Jun 2013 12:10:28 +0000</pubDate>
		<dc:creator>GPS World staff</dc:creator>
				<category><![CDATA[Integration with Other Technologies]]></category>
		<category><![CDATA[Navigation]]></category>
		<category><![CDATA[Public Safety]]></category>
		<category><![CDATA[Warfighter]]></category>
		<category><![CDATA[firefighter]]></category>
		<category><![CDATA[first responder]]></category>
		<category><![CDATA[inertial measurement unit]]></category>
		<category><![CDATA[inertial navigation system]]></category>

		<guid isPermaLink="false">http://www.gpsworld.com/?p=21467</guid>
		<description><![CDATA[An Enhanced Personal Inertial Navigation System When a team of firefighters, first responders, or soldiers operates inside a building, in urban canyons, underground, in foliage, or under the forest canopy, the GPS-denied environment presents unique navigation challenges. An enhanced personal inertial navigation system (ePINS), based on a strapdown navigation solution using a mid-grade IMU and [...]]]></description>
				<content:encoded><![CDATA[<h3>An Enhanced Personal Inertial Navigation System</h3>
<p><strong>When a team of firefighters, first responders, or soldiers operates inside a building, in urban canyons, underground, in foliage, or under the forest canopy, the GPS-denied environment presents unique navigation challenges. An enhanced personal inertial navigation system (ePINS), based on a strapdown navigation solution using a mid-grade IMU and wavelet-based motion-classification algorithms, can track positions with errors of less than 2 percent of distance traveled in both indoor and outdoor environments.</strong></p>
<p><em>By Yunqian Ma, Wayne Soehren, Wes Hawkinson, and Justin Syrstad</em></p>
<p>Numerous pedestrian navigation applications are currently available or proposed for development. Some of them include localization for coordinating firefighters, first responders, or soldiers. In these applications, the safety and efficiency of the entire team relies directly on the location and orientation of each team member. Operations in high signal interference areas such as cities, rugged terrain, forest, or indoor spaces deliver intermittent or no GPS signal. An alternative to GPS-based location is required.</p>
<p>In this article, we introduce an enhanced personal inertial navigation system (ePINS) solution specifically designed for environments where GPS is unavailable. ePINS combines an array of state-of-the-art sensors and fusion algorithms into a personal navigation system that provides accurate location information for pedestrian applications.</p>
<div id="attachment_21493" class="wp-caption alignright" style="width: 260px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Ma_concept.jpg"><img class="size-thumbnail wp-image-21493" alt="The ePINS concept." src="http://www.gpsworld.com/wp-content/uploads/2013/06/Ma_concept-250x234.jpg" width="250" height="234" /></a><p class="wp-caption-text">The ePINS concept.</p></div>
<p>The ePINS solution has the following benefits:</p>
<ul>
<li>Accurate positioning in GPS-denied environments;</li>
<li>Small, lightweight unit can be easily carried by first responders, rescue workers, or soldiers;</li>
<li>Ruggedized packaging to withstand difficult first responder and military environments.</li>
</ul>
<p>Features of  the ePINS unit include:</p>
<ul>
<li>State-of-the-art micro-electromechanical systems (MEMS) gyros and accelerometers, barometric altitude sensor, and advanced navigation software;</li>
<li>Advanced motion classification algorithms that accurately identify and measure user activity;</li>
<li>Immunity to magnetic disturbances.</li>
</ul>
<h4>Related Work</h4>
<p>In the field of personal navigation, it is common to find systems that rely on sensors that need infrastructure (for example, Wi-Fi positioning) or sensors that actively emit electro-magnetic radiation (such as Doppler radar). These requirements are major drawbacks for communities such as dismounted soldiers in hostile environments.</p>
<p>Other approaches exploit the so-called Zero-velocity update (ZUPT) mechanism, which resets the inertial measurement unit (IMU) velocity errors during the stationary phase of motion. However, implementation of such schemes relies on sensors embedded in footwear, which is not readily accepted in many user communities.</p>
<p>To address these drawbacks, Honeywell has been developing advanced aiding techniques for personal navigation that do not rely on infrastructure and compute a self-contained, relative-navigation solution based only on passive sensors. One technique that Honeywell has developed uses displacement estimation from human-motion models. This technology has been implemented in the ePINS prototype and shows promising performance.</p>
<p>The human-motion model uses IMU measurements as inputs and was developed to infer distance traveled. It generates a displacement estimate that is used as a measurement in the navigation filtering process. The first version of this model was matured under the DARPA individual Precision Inertial Navigation System (iPINS) program. The iPINS system used an IMU, GPS, barometer, and motion classification to estimate a person’s position in both indoor and outdoor environments. In this system, IMU signal characteristics (e.g., peaks and valleys in the accelerations induced by walking) were exploited to differentiate between walking and running. Honeywell recently expanded the human-motion model to identify more specific motion types using a new wavelet motion classification method.</p>
<h4>System Description</h4>
<p>Figure 1 displays the hardware architecture of the ePINS, a small battery-powered, highly integrated electronic system. The ePINS processing platform is an ARM11-based, i.MX31 system-on-module, paired with support electronics. In addition to the processing platform, the ePINS assembly includes a MEMS IMU, a barometric pressure sensor, a digital magnetometer, and a GPS receiver.</p>
<div id="attachment_21473" class="wp-caption alignnone" style="width: 442px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_figure_1.jpg"><img class=" wp-image-21473 " alt="ePINS hardware architecture." src="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_figure_1.jpg" width="432" height="284" /></a><p class="wp-caption-text">Figure 1. ePINS hardware architecture.</p></div>
<p>The MEMS IMU provides inertial measurements for strapdown navigation. The IMU’s small package size, light weight, low power consumption, and impressive performance make it attractive for use in the ePINS system. The device is less than 5 cubic inches and weighs less than 0.35 pounds. It consumes about 3 watts of power with a typical current draw of 600mA at 5V.</p>
<p>The ePINS software system is shown in Figure 2. The navigation software runs within Honeywell’s Embedded Computing Toolbox and Operating System (ECTOS IIc), which provides a layered, customizable, and reusable software architecture for implementing navigation, guidance, and control software. A Honeywell-developed simulation tool for offline analysis and development of ECTOS-based software was also used in ePINS development and testing.</p>
<div id="attachment_21474" class="wp-caption alignnone" style="width: 442px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_figure_2.jpg"><img class=" wp-image-21474 " alt="Figure 2.  ECTOS IIc hierarchical software structure." src="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_figure_2.jpg" width="432" height="286" /></a><p class="wp-caption-text">Figure 2. ECTOS IIc hierarchical software structure.</p></div>
<p>The ePINS demonstration device can achieve path performance of less 2 percent distance traveled for walking motion after 1 hour of operation, independent of the magnetic environment. Current performance, packaging characteristics, and interfaces are summarized in Table 1.</p>
<div id="attachment_21484" class="wp-caption alignnone" style="width: 442px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_Table_1.jpg"><img class=" wp-image-21484 " alt="table 1  ePINS performance objectives and physical specifications." src="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_Table_1.jpg" width="432" height="223" /></a><p class="wp-caption-text">Table 1. ePINS performance objectives and physical specifications.</p></div>
<h4>Algorithm Description</h4>
<p>Figure 3 depicts the overall sensor integration and data processing scheme used in the ePINS device.</p>
<div id="attachment_21475" class="wp-caption alignnone" style="width: 442px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_figure_3.jpg"><img class=" wp-image-21475 " alt="Figure 3. Sensor integration using the ECTOS extended Kalman filter." src="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_figure_3.jpg" width="432" height="298" /></a><p class="wp-caption-text">Figure 3. Sensor integration using the ECTOS extended Kalman filter.</p></div>
<p><strong>Extended Kalman Filter (EKF). </strong> The EKF estimates the navigation and sensor errors and computes the resets applied to the strapdown navigation solution to increase its accuracy. Error models for the navigation sensors (IMU, barometric altimeter, magnetometer, GPS, and motion classification) are contained in the EKF. For the ePINS device, the virtual measurements from the step-length model and the strapdown navigation solution are fused by the EKF to assist in bounding the time dependent error growth of the strapdown navigator, which in turn helps maintain calibration of the inertial sensors. A key output of the EKF is the navigation confidence, which is an estimate of the accuracy of the navigation solution.</p>
<p>An important aspect of the EKF and step-length modeling is the residual test that the EKF supports. This test provides a reasonableness comparison between the step-length model estimate and the distance predicted by the strapdown navigation system. This capability significantly increases the robustness of the navigation solution, especially when the user is engaged in motions not recognized during motion classification.</p>
<p><strong>Human-Motion Model.</strong> The human-motion model includes two components: wavelet motion classification and step-length model estimation. The wavelet motion classification identifies the type of motion the user is performing, and the step-length model acts as a virtual sensor that quantifies the motion as a distance-traveled estimate.</p>
<p><strong>Wavelet Motion Classification.</strong> Human motions are very diverse and highly irregular. Determining what motion is being performed is a challenging problem of classification. Honeywell’s solution is based on wavelet transformation of IMU data. Predefined, or known, characteristics of a variety of motions (such as walking, running, crawling, etc.) are cataloged and stored to a device’s memory. Estimates of those same characteristics for a user are then computed in real time and compared to the catalog of stored information to find the best match.</p>
<p>Generating the catalog of stored information is an offline task that begins by “segmenting” recorded IMU time domain data into individual steps. An example of the output of the segmentation process is shown in Figure 4.</p>
<div id="attachment_21483" class="wp-caption alignnone" style="width: 442px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Ma_Figure4.jpg"><img class=" wp-image-21483 " alt="Figure 4. Segmentation of the IMU data using the y-axis accelerometer signal." src="http://www.gpsworld.com/wp-content/uploads/2013/06/Ma_Figure4.jpg" width="432" height="302" /></a><p class="wp-caption-text">Figure 4. Segmentation of the IMU data using the y-axis accelerometer signal.</p></div>
<p>Figure 5 displays the segmentation results for two different walking styles (in red and blue) across approximately 15 example steps. As is evident from the graph, walking has characteristics that are common across users, for example, the sharp peaks in the z-axis acceleration caused by foot-ground impacts. Once the data has been segmented, a wavelet transformation on each data channel is performed. Wavelet transformation for many users over many different motion types takes place offline. Subsequently, a wavelet descriptor is built for each motion type based on the transformations into the wavelet domain. With this method, a wide variety of information (that is, descriptors) suitable for input to a classifier is captured about each motion. These descriptors are then cataloged and stored in memory on the ePINS device.</p>
<div id="attachment_21476" class="wp-caption alignnone" style="width: 442px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Ma_Figure_5.jpg"><img class=" wp-image-21476 " alt="Figure 5. Sample steps for two subjects (red) and (blue). " src="http://www.gpsworld.com/wp-content/uploads/2013/06/Ma_Figure_5.jpg" width="432" height="343" /></a><p class="wp-caption-text">Figure 5. Sample steps for two subjects (red) and (blue).</p></div>
<p>Finally, for the online phase, the wavelet descriptor of the incoming IMU data is calculated by performing a wavelet transformation on each data channel. This descriptor is then compared to the pre-computed and stored descriptors to classify the motion. FIGURE 7 shows an example of the motion classifier output, where a running motion was used as an input. The classifier successfully determined the motion type (blue field), frequency and phase of the input motion, depicted by the tallest rectangle in the figure.</p>
<div id="attachment_21478" class="wp-caption alignnone" style="width: 442px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_figure_7.jpg"><img class=" wp-image-21478 " alt="Figure 7. Classification results from a query of running at a certain frequency and phase (depicted by the dark sphere)." src="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_figure_7.jpg" width="432" height="256" /></a><p class="wp-caption-text">Figure 7. Classification results from a query of running at a certain frequency and phase (depicted by the dark sphere).</p></div>
<p><strong>Step-Length Modeling.</strong> Once the current motion is identified, a step-length model specific to that motion is used to aid the navigation algorithms. The model for each motion type is obtained by first collecting data that measures step length and step frequency. From this data, the step-length models can be computed by performing a regression analysis of the step-length vs. step-frequency data. Since the step-length models act as a virtual sensor, the models must be as accurate as possible to achieve better system performance. To attain model accuracy, an accurate data collection method is needed.</p>
<p>For ePINS development, step-length models for multiple users have been identified from step-length and timing information using a precise GPS truth reference system. Step-length regression calculations then determine the step length as a function of step frequency (that is, inverse of the step time period).  An example of GPS truth data and the corresponding regression model are shown in FIGURE 6 for walking motions.</p>
<div id="attachment_21477" class="wp-caption alignnone" style="width: 442px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_figure_6.jpg"><img class=" wp-image-21477 " alt="Figure 6. Step length versus frequency for the walking of subject." src="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_figure_6.jpg" width="432" height="361" /></a><p class="wp-caption-text">Figure 6. Step length versus frequency for the walking of subject.</p></div>
<p>Although basic step-length models are created offline, online calibration of the step-length model can be performed by the EKF if GPS is available during operation. Online calibration tends to increase the overall position accuracy, as variations in the step-length models are likely due to slight variations in biometric differences across humans, terrain features, and even mission plans and duration.</p>
<p><strong>Heading Determination.</strong> Heading initialization is one of the key concerns during system start up. In its current operational use, the ePINS device may perform a dynamic or a static initialization of heading. The static method requires the user to survey the system’s initial heading to an accuracy value that is usually specified by mission performance objectives; the absolute position accuracy is dependent upon the accuracy of the initial heading.</p>
<p>The dynamic method is a general method for heading initialization; it is performed without input from the user, but is possible only when GPS is available. This method of heading initialization does not use any a priori information about heading and requires an EKF implementation with a large-azimuth error model. This method requires an additional period of time in which the heading error uncertainty converges.</p>
<p><strong>User Interface.</strong> During a mission, the user can interact with the navigation system and monitor its output on a display. The current ePINS prototype offers two-way communication via a serial connection. The serial communication is made wireless by the addition of a Bluetooth interface. Users can use this link to monitor the status of the navigation solution and to send commands to the device.</p>
<p>Honeywell has developed an application for the Android platform for this purpose. One of the key features of the interface design is that the navigation system outputs data in a standard NEMA format. Thus, publically available Android applications, not just proprietary applications, can also receive and display the navigation solution output by the ePINS device.</p>
<p>Honeywell’s personal navigation application displays the user’s traveled trajectory in real-time. The application can be adapted to include building floor plans as well as other navigation information.</p>
<h4>Results</h4>
<p>The ePINS prototype has been evaluated both in simulations and indoor/outdoor experiments. The navigation results presented here were obtained in February 2012 at a Honeywell facility (FIGURE 8). First, the user completed the heading calibration, and then online step parameter estimation in the presence of GPS was performed. Once calibration and training was completed, the GPS was disabled to simulate a GPS-denied environment outdoors. The user than transitioned to indoors (with GPS still disabled), and walked a course inside that included walking up and down stairs (FIGURE 9) and ended in a conference room (FIGURE 10).</p>
<div id="attachment_21479" class="wp-caption alignnone" style="width: 442px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Ma_Figure_8.jpg"><img class=" wp-image-21479 " alt="Figure 8. Course for the Honeywell facility demonstration." src="http://www.gpsworld.com/wp-content/uploads/2013/06/Ma_Figure_8.jpg" width="432" height="414" /></a><p class="wp-caption-text">Figure 8. Course for the Honeywell facility demonstration.</p></div>
<div id="attachment_21480" class="wp-caption alignnone" style="width: 442px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Ma_Figure_9.jpg"><img class=" wp-image-21480 " alt="Figure 9. The user walking up stairs." src="http://www.gpsworld.com/wp-content/uploads/2013/06/Ma_Figure_9.jpg" width="432" height="305" /></a><p class="wp-caption-text">Figure 9. The user walking up stairs.</p></div>
<div id="attachment_21481" class="wp-caption alignnone" style="width: 442px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/Ma_Figure_10.jpg"><img class=" wp-image-21481 " alt="Figure 10. The user at the end of the demo." src="http://www.gpsworld.com/wp-content/uploads/2013/06/Ma_Figure_10.jpg" width="432" height="304" /></a><p class="wp-caption-text">Figure 10. The user at the end of the demo.</p></div>
<p>Over these conditions, the ePINS system performed robustly and within performance specifications. Live demonstrations and testing showing similar levels of performance were performed at the 2012 Joint Navigation Conference (JNC) and at military test sites in California and Indiana.</p>
<h4>Summary</h4>
<p>The technical approach of the ePINS solution to the problem of personnel navigation in GPS-denied environments is based on a strapdown navigation solution maintained using a mid-grade IMU and advanced motion-classification algorithms. We integrated an array of sensors and software into a system that provides accurate position information and is suitable for use by first responders, soldiers, and other personnel where GPS is unavailable. ePINS works well for a variety of pedestrian motion types, including walking, running, crawling, walking upstairs, walking downstairs, sidestepping, and walking backwards. The motion classification and modeling method is extensible to other motion types.</p>
<p>We tested the ePINS system in indoor and outdoor environments. FIGURE 11 depicts the future ePINS concept, and TABLE 2 presents its future physical characteristics.</p>
<div id="attachment_21482" class="wp-caption alignnone" style="width: 442px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_figure_11.jpg"><img class=" wp-image-21482 " alt="Figure 11. Future ePINS concept and mounting position." src="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_figure_11.jpg" width="432" height="308" /></a><p class="wp-caption-text">Figure 11. Future ePINS concept and mounting position.</p></div>
<div id="attachment_21485" class="wp-caption alignnone" style="width: 460px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_Table_2.jpg"><img class=" wp-image-21485 " alt="Table 2. Packaging characteristics of the future ePINS." src="http://www.gpsworld.com/wp-content/uploads/2013/06/ma_Table_2.jpg" width="450" height="302" /></a><p class="wp-caption-text">Table 2. Packaging characteristics of the future ePINS.</p></div>
<h4>Acknowledgments</h4>
<p>This article is based on a presentation made at ION GNSS 2012.</p>
<h4>Manufacturers</h4>
<p>The ePINS processing platform uses <a href="http://www.honeywell.com" target="_blank">Honeywell</a> Agile Navigation and Guidance Integrated Electronics support electronics. It includes a Honeywell HG1930 MEMS IMU, a <a href="http://www.bosch-sensortec.com" target="_blank">Bosch</a> Sensortec BMP085 barometric pressure sensor, a Honeywell HMC6343 digital magnetometer, and a <a href="http://www.novatel.com" target="_blank">NovAtel</a> OEMStar GPS receiver.</p>
<hr />
<p><em>Yunqian Ma is a principal scientist at Honeywell Aerospace. He received his Ph.D. degree in electrical engineering from the University of Minnesota, Twin Cities. He is currently the program manager of the GPS-denied navigation program and the next-generation personal navigation program.</em></p>
<p><em>Wayne Soehren is a senior technical manager at Honeywell Aerospace. He was the program manager for the development of Honeywell’s first MEMS-based GPS/INS, which developed the core capability now used in Honeywell’s IGS-2XX family of MEMS-based GPS/INS products. He holds an MSEE from the University of Minnesota.</em></p>
<p><em>Wes Hawkinson is an engineering fellow at Honeywell Aerospace. He holds a BSEE/CE from the University of Wisconsin–Madison.</em><br />
<em> Justin Syrstad is a guidance and navigation scientist. He received a master’s degree in aerospace engineering from the University of Minnesota.</em></p>
]]></content:encoded>
			<wfw:commentRss>http://www.gpsworld.com/following-the-team-into-danger/feed/</wfw:commentRss>
		<slash:comments>1</slash:comments>
		</item>
		<item>
		<title>The System: Galileo Leaves the Building</title>
		<link>http://www.gpsworld.com/the-system-galileo-leaves-the-building/?utm_source=rss&#038;utm_medium=rss&#038;utm_campaign=the-system-galileo-leaves-the-building</link>
		<comments>http://www.gpsworld.com/the-system-galileo-leaves-the-building/#comments</comments>
		<pubDate>Sat, 01 Jun 2013 07:35:13 +0000</pubDate>
		<dc:creator>Alan Cameron</dc:creator>
				<category><![CDATA[Alan Cameron]]></category>
		<category><![CDATA[Galileo]]></category>
		<category><![CDATA[GLONASS]]></category>
		<category><![CDATA[GPS Modernization]]></category>
		<category><![CDATA[The System]]></category>
		<category><![CDATA[IRNSS]]></category>

		<guid isPermaLink="false">http://www.gpsworld.com/?p=21460</guid>
		<description><![CDATA[In the early hours of May 15, Galileo’s first full operational capability (FOC) satellite left manufacturer OHB System AG’s integration hall in Bremen, Germany, after successfully completing integration and system testing. Later that same day, it arrived by road at the European Space Agency’s (ESA’s) technical center at Noordwijk in the Netherlands for a rigorous [...]]]></description>
				<content:encoded><![CDATA[<p>In the early hours of May 15, Galileo’s first full operational capability (FOC) satellite<a href="http://www.gpsworld.com/first-galileo-foc-satellite-heads-to-testing/" target="_blank"> left manufacturer OHB System AG’s integration hall </a>in Bremen, Germany, after successfully completing integration and system testing. Later that same day, it arrived by road at the European Space Agency’s (ESA’s) technical center at Noordwijk in the Netherlands for a rigorous set of tests to check its readiness for launch. The tests will simulate different aspects of launch and space environment. The comprehensive test program will validate the new design and all the FOC satellites to follow.</p>
<p>This first FOC satellite is functionally identical to the first four in-orbit validation (IOV) satellites already in orbit, but has been built by a separate industrial team. Like the other 21 FOC satellites so far procured by ESA, the satellite’s prime contractor is OHB System AG, and the navigation payload was produced by Surrey Satellite Technology Ltd. in Guildford, UK.</p>
<p>Thermal vacuum testing at the European Space Research and Technology Centre (ESTEC) will simulate temperature extremes the satellites must endure in the airlessness of space throughout their 12-year working lifetimes. Without any moderating atmosphere, temperatures can shift hundreds of degrees from sunlight to shadow.</p>
<p>Other activities on the schedule include shaker and acoustic noise testing — simulating the vibration and noise of launch — as well as electromagnetic compatibility and antenna testing, placing the satellite in chambers shielded from all external radio signals to reproduce infinite space and check that its various antennas and electrical systems are interoperable without harmful interference.</p>
<p>“The Galileo FOC satellites provide the same capabilities as the previous IOV satellites, but with improved performance, such as higher transmit power,” explained Giuliano Gatti, the head of the Galileo Space Segment Procurement Office. “They are to all intents a new design that requires a full checkout before getting the green light for launch.”</p>
<p>The second FOC flight model is due to arrive at ESTEC in early June, and the third in the middle of July. The first two satellites are to be placed in orbit on board a Soyuz launcher, with a scheduled lift-off from Kourou in French Guyana this fall, with two more due to follow by the end of the year.</p>
<p>The first four Galileo IOV satellites, launched in 2011 and 2012, were provided by EADS Astrium with Thales Alenia Space Italy responsible for integrating the satellites and Astrium in Portsmouth, UK, providing the navigation payloads. They provided their first navigation fix in March 2013.</p>
<p>The definition, development and in-orbit validation phases of the Galileo programme are being carried out by ESA and co-funded with the European Commission (EC).</p>
<p>The subsequent FOC phase is managed and funded by the EC. The commission has delegated the role of design and procurement agent to ESA for the FOC phase. At the same time as the satellites are being assembled on a production-line basis, ground stations are also being established on European territories around the globe.</p>
<div id="attachment_21128" class="wp-caption alignnone" style="width: 563px"><a href="http://www.gpsworld.com/wp-content/uploads/2013/05/GPSIIF-7.jpg"><img class="wp-image-21128 " alt="Photo credit: Pat Corkery, United Launch Alliance." src="http://www.gpsworld.com/wp-content/uploads/2013/05/GPSIIF-7.jpg" width="553" height="368" /></a><p class="wp-caption-text">Photo credit: Pat Corkery, United Launch Alliance.</p></div>
<h3>GPS Leaves This Earth</h3>
<p>A t 5:38 p.m. Eastern Daylight Time (21:38 UTC) on May 15,  the fourth GPS IIF satellite, Space Vehicle Number (SVN) 66 built by Boeing, <a href="http://www.gpsworld.com/gps-iif-4-successfully-launched-from-cape-canaveral/" target="_blank">ascended towards orbit</a> aboard a United Launch Alliance Atlas V rocket at from Cape Canaveral Air Force Station, Florida.</p>
<p>“The GPS constellation remains healthy and continues to meet and exceed the performance standards to which the satellites were built. Our goal is to deliver sustained, reliable GPS capabilities to America’s warfighters, our allies, and civil users around the world, and this is done by maintaining GPS performance, fielding new capabilities and developing more robust modernized capabilities for the future,” said Colonel Bernie Gruber, director of the U.S. Air Force Space and Missile Systems Center’s GPS Directorate.</p>
<p>The new capabilities of the IIF satellites will provide greater navigational accuracy through improvements in atomic clock technology; a more robust signal for commercial aviation and safety-of-life applications, known as the new third civil signal (L5); and a 12-year design life providing long-term service. These upgrades deliver improved anti-jam capabilities for warfighters and improved security for military and civil users around the world, the Air Force said in a statement.</p>
<p>The IIF-4 satellite is expected to complete testing in August, after which it will be utilized as a reserve or backup satellite. It becomes the fourth satellite in a 12-strong network of GPS IIF spacecraft manufactured by Boeing as lead contractor, the first of which was boosted into orbit in May 2010. The Air Force expects the first of the next-generation GPS IIIA satellites to enter service sometime in 2014.</p>
<h3>System Briefs</h3>
<p><strong>GLONASS.</strong> The GLONASS 747 M-series satellite <a href="http://www.gpsworld.com/glonass-satellite-launched/" target="_blank">launched on April 26</a> has maneuvered into an orbital slot near GLONASS 728, the operational satellite in Plane 1, slot 2. 747 will presumably serve as a reserve until it replaces 728, unless another Plane 1 satellite expires first. The next Russian launch, a GLONASS-M trio, is scheduled for July 1. There are currently 24 operational GLONASS satellites.</p>
<p><strong>IRNSS.</strong> The first Indian Regional Navigation Satellite System satellite is <a href="http://www.gpsworld.com/resources/upcoming-gnss-satellite-launches/" target="_blank">expected to rise</a> at the end of June. The IRNSS plans to orbit of seven: three geostationary and four geosynchronous, providing regional coverage via navigation signals in the L5 and S bands.</p>
]]></content:encoded>
			<wfw:commentRss>http://www.gpsworld.com/the-system-galileo-leaves-the-building/feed/</wfw:commentRss>
		<slash:comments>1</slash:comments>
		</item>
	</channel>
</rss>

<!-- Performance optimized by W3 Total Cache. Learn more: http://www.w3-edge.com/wordpress-plugins/

Page Caching using apc
Object Caching 1655/1726 objects using apc

 Served from: www.gpsworld.com @ 2013-06-11 10:39:44 by W3 Total Cache --