PERFORMANCE OF GLOBAL POSITIONING SYSTEM BLOCK II/IIA/IIR ON-ORBIT NAVSTAR CLOCKS

Similar documents
GLOBAL POSITIONING SYSTEM CONSTELLATION CLOCK PERFORMANCE

ANALYSIS OF ON-ORBIT BEHAVIOR OF GPS BLOCK 11-R TIME KEEPING SYSTEM

On the Line Q Degradation in Hydrogen Masers

KALMAN PLUS WEIGHTS: A TIME SCALE ALGORITHM*

DETECTION OF THE GRAVITATIONAL REDSHIFT OF THE CESIUM FREQUENCY STANDARD AT CRL

An Ensemble of Atomic Fountains

TIME SERIES ANALYSIS OF VLBI ASTROMETRIC SOURCE POSITIONS AT 24-GHZ

The Statistics of GPS

THE LONG-TERM STABILITY OF THE U.S. NAVAL OBSERVATORY S MASERS

VLBA IMAGING OF SOURCES AT 24 AND 43 GHZ

Analysis Comparison between CFD and FEA of an Idealized Concept V- Hull Floor Configuration in Two Dimensions. Dr. Bijan Khatib-Shahidi & Rob E.

Wilson G. Reid SFA, Inc., Largo, MD 20772, USA

Report Documentation Page

Attribution Concepts for Sub-meter Resolution Ground Physics Models

RESPONSE OF A GEOSYNCHRONOUS SPACECRAFT S CRYSTAL OSCILLATOR TO SOLAR FLARES: RESULTS OF A SPACE EXPERIMENT

Report Documentation Page

Satellite Observations of Surface Fronts, Currents and Winds in the Northeast South China Sea

A KALMAN FILTER FOR ATOMIC CLOCKS AND TIMESCALES

Experimental and Theoretical Studies of Ice-Albedo Feedback Processes in the Arctic Basin

Predicting Tropical Cyclone Formation and Structure Change

Swash Zone Dynamics: Modeling and Data Analysis

Real-Time Environmental Information Network and Analysis System (REINAS)

SOLAR FLARES AND PRECISE SATELLITE TIMEKEEPING

System Reliability Simulation and Optimization by Component Reliability Allocation

Assimilation of Synthetic-Aperture Radar Data into Navy Wave Prediction Models

An Examination of 3D Environmental Variability on Broadband Acoustic Propagation Near the Mid-Atlantic Bight

USNO Analysis Center for Source Structure Report

Hurricane Wave Topography and Directional Wave Spectra in Near Real-Time

Diagonal Representation of Certain Matrices

Improvements in Modeling Radiant Emission from the Interaction Between Spacecraft Emanations and the Residual Atmosphere in LEO

Predictive Model for Archaeological Resources. Marine Corps Base Quantico, Virginia John Haynes Jesse Bellavance

AIR FORCE RESEARCH LABORATORY Directed Energy Directorate 3550 Aberdeen Ave SE AIR FORCE MATERIEL COMMAND KIRTLAND AIR FORCE BASE, NM

Analysis of South China Sea Shelf and Basin Acoustic Transmission Data

NAVGEM Platform Support

Effects of Constrictions on Blast-Hazard Areas

1/12 Pacific HYCOM: The End Of A Long Simulation

Report Documentation Page

PIPS 3.0. Pamela G. Posey NRL Code 7322 Stennis Space Center, MS Phone: Fax:

The Effects of Magnetic Storm Phases on F-Layer Irregularities from Auroral to Equatorial Latitudes

Use of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications

SW06 Shallow Water Acoustics Experiment Data Analysis

Award # N J-1716

P. Kestener and A. Arneodo. Laboratoire de Physique Ecole Normale Supérieure de Lyon 46, allée d Italie Lyon cedex 07, FRANCE

The Navy Precision Optical Interferometer for SSA applications: an update

AUTOMATIC FREQUENCY CALIBRATION USING FUZZY LOGIC CONTROLLER

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION REPORT NUMBER

Internal Tide Generation in the Indonesian Seas

Understanding Near-Surface and In-cloud Turbulent Fluxes in the Coastal Stratocumulus-topped Boundary Layers

MODELING SOLAR UV/EUV IRRADIANCE AND IONOSPHERE VARIATIONS WITH MT. WILSON MAGNETIC INDICIES

MTIE AND TDEV ANALYSIS OF UNEVENLY SPACED TIME SERIES DATA AND ITS APPLICATION TO TELECOMMUNICATIONS SYNCHRONIZATION MEASUREMENTS

Experimental and Theoretical Studies of Ice-Albedo Feedback Processes in the Arctic Basin

FRACTAL CONCEPTS AND THE ANALYSIS OF ATMOSPHERIC PROCESSES

REGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT)

Near-Surface Dispersion and Circulation in the Marmara Sea (MARMARA)

Estimating Geometric Aspects of Relative Satellite Motion Using Angles-Only Measurements

The Extratropical Transition of Tropical Cyclones

LAGRANGIAN MEASUREMENTS OF EDDY CHARACTERISTICS IN THE CALIFORNIA CURRENT

Volume 6 Water Surface Profiles

Discovery of Planetary Systems With SIM

Computer Simulation of Sand Ripple Growth and Migration.

Nonlinear Internal Waves: Test of the Inverted Echo Sounder

Development and Application of Acoustic Metamaterials with Locally Resonant Microstructures

Advanced Numerical Methods for NWP Models

INFRARED SPECTRAL MEASUREMENTS OF SHUTTLE ENGINE FIRINGS

ANALYSIS OF LOCALIZED HIGH MAGNETIC SUSCEPTIBILITY ZONES AT JEFFERSON PROVING GROUND, INDIANA

INFRARED SPECTROSCOPY OF HYDROGEN CYANIDE IN SOLID PARAHYDROGEN (BRIEFING CHARTS)

OBSERVABLE RELATIVISTIC FREQUENCY STEPS INDUCED BY GPS ORBIT CHANGES

Scattering of Internal Gravity Waves at Finite Topography

GPS BLOCK IIR RUBIDIUM ATOMIC FREQUENCY STANDARD LIFE TEST

Ocean Acoustics Turbulence Study

HYCOM Caspian Sea Modeling. Part I: An Overview of the Model and Coastal Upwelling. Naval Research Laboratory, Stennis Space Center, USA

Erik L. Swanberg 1 and Steven G. Hoffert 2. Veridian Systems 1, Autometric 2. Sponsored by Defense Threat Reduction Agency

Coastal Ocean Circulation Experiment off Senegal (COCES)

Performance of the Kalman Filter of the Global Positioning System Operational Control Segment During January - March 1997*

CRS Report for Congress

Sensitivity of West Florida Shelf Simulations to Initial and Boundary Conditions Provided by HYCOM Data-Assimilative Ocean Hindcasts

LISA: THE LASER INTERFEROMETER SPACE ANTENNA

TECHNOLOGY TRANSFER OPPORTUNITIES FOR THE CONSTRUCTION ENGINEERING COMMUNITY CONSTRU-CTION CD-IPOIIIOIED BY

A Bound on Mean-Square Estimation Error Accounting for System Model Mismatch

Babylonian resistor networks

Forecasting Tides in Global HYCOM

GPS IIR RUBIDIUM CLOCKS: IN-ORBIT PERFORMANCE ASPECTS

Aviation Weather Routing Tool: A Decision Aid for Manned/Unmanned Aircraft Routing

AN EVALUATION OF VARIOUS SPACE CLOCKS FOR GPS IIF

DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited.

DIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION

Catalytic Oxidation of CW Agents Using H O in 2 2 Ionic Liquids

Crowd Behavior Modeling in COMBAT XXI

Award # N LONG TERM GOALS

Mine Burial Studies with a Large Oscillating Water-Sediment Tunnel (LOWST)

Analysis of Subsurface Velocity Data from the Arctic Ocean

DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited.

Report Documentation Page

On Applying Point-Interval Logic to Criminal Forensics

PREDICTION OF HIGH TEMPERATURE GREASE LIFE USING A DECOMPOSITION KINETIC MODEL

APPLICATION OF THE GSF-1 ALGORITHM TO THE NEAR-OPTIMAL TIMESCALE PREDICTION OF THE HYDROGEN MASER

USMC Enlisted Endstrength Model

Comparative Analysis of Flood Routing Methods

Quantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS)

Grant Number: N IP To compare obtained theoretical results with NPAL experimental data.

Transcription:

33rdAnnual Precise Time and Time nterval (PTT)Meeting PERFORMANCE OF GLOBAL POSTONNG SYSTEM BLOCK /A/R ON-ORBT VSTAR CLOCKS Jay Oaks, Thomas B. McCaskill, Marie M. Largay U.S. Naval Research Laboratory James A. Buisson AE Abstract Analysis of the performance of all on-orbit Navstar space vehicle clocks and Global Positioning System (GPS) monitor station reference clocks is performed by the Naval Research Laboratory (NRL), in cooperation with the GPS Master Control Station, under the sponsorship of the GPS Joint Program Office (JPO). The measurements used in the analysis are collected by multi-channel GPS receivers located at the Air Force and National magery and Mapping Agency (NMA) monitor stations. The offset of each Navstar clock is computed with respect to the Department of Defense Master Clock. The resultant Navstar clock offsets are then used to compute frequency offset, drift offset,frequency stability profiles, and frequency stability histories. The performance of selected cesium and rubidium atomic clocks will be presented. Frequency stability results are calculated using sample times that vary from 5 minutes to several days. The operational Navstar timing signals will be ranked (according to frequency stability) for sample times of 6 hours and day. n addition,the Block R clocks will be compared to results from an ongoing NRL laboratory life-test on two jlight-qualified Block R rubidium atomicfrequency standards. Observed clock behavior and anomalies will be discussed. NTRODUCTON Performance of the on-orbit Navstar atomic clocks in the Global Positioning System (GPS) constellation is summarized using a multi-year database, including data collected to October 200. This ongoing work is sponsored by the GPS Joint Program Office (JPO) and is done in cooperation with the GPS Master Control Station (MCS). The measurements were collected from a network of Air Force (USAF) and National magery and Mapping Agency (NMA) monitor stations. The NMA Washington, D.C. monitor station is located at the United States Naval Observatory (USNO) and is referenced to the Department of Defense (DoD) Master Clock. The results from the NRL analyses of the Navstar spacecraft vehicle (SV) and monitor station (MS) atomic clocks are used by the GPS MCS to aid in the prediction of clock states in the Kalman filter. This aiding results in an improvement in navigation and time transfer performance [l]. The NRL results are reported to the GPS JPO and the MCS, and made available to the GPS working group community via a NRL Web site. Performance and lifetime evaluation tests are being conducted at NFU on Navstar flight-qualified candidate clocks. Currently, two Block W rubidium atomic frequency standards are undergoing life test in the NRL, Precision Clock Evaluation Facility. The results of these analyses are compared with those of the on-orbit clocks and are also available on the NRL, Web site. 469

Report Documentation Page Form Approved OMB No. 0704-088 Public reporting burden for the collection of information is estimated to average hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for nformation Operations and Reports, 25 Jefferson Davis Highway, Suite 204, Arlington VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to a penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number.. REPORT DATE NOV 200 2. REPORT TYPE 3. DATES COVERED 00-00-200 to 00-00-200 4. TTLE AND SUBTTLE Performance of Global Positioning System Block /A/R On-Orbit Navstar Clocks 5a. CONTRACT NUMBER 5b. GRANT NUMBER 5c. PROGRAM ELEMENT NUMBER 6. AUTHOR(S) 5d. PROJECT NUMBER 5e. TASK NUMBER 5f. WORK UNT NUMBER 7. PERFORMNG ORGANZATON ME(S) AND ADDRESS(ES) U.S. Naval Observatory,3450 Massachusetts Avenue NW,Washington,DC,20392-5420 8. PERFORMNG ORGANZATON REPORT NUMBER 9. SPONSORNG/MONTORNG AGENCY ME(S) AND ADDRESS(ES) 0. SPONSOR/MONTOR S ACRONYM(S) 2. DSTRBUTON/AVALABLTY STATEMENT Approved for public release; distribution unlimited. SPONSOR/MONTOR S REPORT NUMBER(S) 3. SUPPLEMENTARY NOTES See also ADM00482. 33rd Annual Precise Time and Time nterval (PTT) Systems and Applications Meeting, 27-29 Nov 200, Long Beach, CA 4. ABSTRACT see report 5. SUBJECT TERMS 6. SECURTY CLASSFCATON OF: 7. LMTATON OF ABSTRACT a. REPORT unclassified b. ABSTRACT unclassified c. THS PAGE unclassified Same as Report (SAR) 8. NUMBER OF PAGES 2 9a. ME OF RESPONSBLE PERSON Standard Form 298 (Rev. 8-98) Prescribed by ANS Std Z39-8

The current tracking network of monitor stations, depicted in Figure, consists of five Air Force MS, and NMA MS. The offset of each MS clock with respect to the DoD Master Clock is computed using Multiple-Path Linked Common-View Time Transfer (LCVTT) [2]. The network links used in the LCVTT computation are shown in Figure 2. The advantage of using multiple paths is that the loss of measurements from one station does not result in the loss of time transfer data for the remaining stations in the network. The offset of each MS clock is combined with the offset of the Navstar SV clocks, observed at each MS, to produce Navstar SV clock offsets with respect to the DoD Master Clock. This process results in multiple measurements of the Navstar SV clocks at each 5-minute epoch. The measurements are averaged at each epoch to get continuous coverage of each Navstar clock with respect to the DoD Master Clock [3]. The continuous coverage measurements are then used to compute the frequency offset, drift offset, frequency stability profiles, and frequency stability histories for the Navstar clocks using the NRL Clock Analysis Software System (CASS). CASS is a complex set of software analysis tools comprised of over 400,000 lines of code that was developed at NRL over the last 8 years. The current active database includes data for all Block va/r SV clocks for 34 Navstar SVs (out of a total of 44 launched). To date, 09 of the 67 Navstar clocks on orbit have been operated and 74 are Block /A/W that are included in the active database. Archived data are maintained for all of the Block Navstar clocks beginning with Navstar, which was launched in 978. CONSTELLATON OVERVEW A summary of the operational Navstar clocks in the GPS constellation as of 30 September 200 is presented in Figure 3. Each Block /ww SV is shown by plane, position in the constellation, and type of clock that was operating on each SV. Fourteen of the 28 clocks operating were rubidium atomic frequency standards (RAFS) and 4 were cesium atomic frequency standards (CAFS). Of the 4 RAFS, six are on the new Block W SVs. The total operating time for each of the Navstar SVs since the SV was inserted into the constellation is shown in Figure 4. Twenty-three SVs have met or have exceeded the required Block VA mean mission duration of 6 years. Eighteen SVs have exceeded 8 years of operation, and nine have exceeded 0 years of operation. Six Block /LA SVs, Navstars 4, 6, 8, 9, 20, and 28, have been decommissioned. Navstar 28 was the only space vehicle that did not meet the required mean mission ddration before being decommissioned. The average lifetime of all SVs in the GPS constellation is 7.5 years. The numbers of clocks that have been placed in operation on each of the active SVs, since the SVs were inserted into the constellation, are shown in Figure 5. Of the active SVs, 0 are on the first clock, 2 are on the second clock, three are on the third clock, and three, Navstars 22,24, and 32, are on the fourth clock. While Navstar 43 shows two clocks placed into operation, one was turned on for test purposes and then switched off. t was a good clock and can be reactivated when the need arises. The Block W SVs are equipped with only three clocks, all RAFS. The Block and A SVs each have four clocks, two CAFS and two RAFS. The operating lifetime, or length of service, of the clocks that were operating as of October 200 is shown in Figure 6. The bars in the front row refer to rubidium clocks, and those in the back row refer to cesium clocks. Four clocks, cesiums on Navstars 3, 5, 23, and 27, have exceeded 8 years of continuous operation, and the cesium clock on Navstar 3 has exceeded 2 years of continuous operation. One of the Block W rubidium clocks, on Navstar 43, has exceeded 4 years of continuous operation. The average age of the currently active CAFS is 6.9 years and the RAFS is 2.0 years, as of October 200. The lifetime of the deactivated clocks ranked by relative order of lifetime is presented in Figure 7. Rubidium clocks are represented in the front row, and cesium clocks are represented in the back row. The notation SVN and SN on the horizontal axis denotes Navstar space vehicle number and clock serial number. The chart shows that the cesium clocks have an average lifetime 470

that is almost three times that of the rubidium clocks, 3.9 years for the 4 cesium clocks compared to.4 years for the 30 rubidium clocks. VSTAR TMNG SGL MEASUREMENTS The phrase Navstar Timing Signal is used rather than Navstar Clock Offset because the AFS output is further modified by the electronics before being broadcast by each Navstar SV to the user. The Block /TA Navstars are equipped with a Frequency Standard Distribution Unit (FSDU), and the Block W Navstars are equipped with a Time Keeping System (TKS) [4]. The TKS provides the capability of tuning the output frequency and drift offset. This capability has been used to adjust the output frequency and drift of the Block R on-orbit RAFS to within a few pp0l2 of the DoD Master Clock. FREQUENCY STABLTY MODELS NRL currently employs two models to measure the frequency stability performance of the Navstar clocks in the time domain. The Allan deviation is normally used in the analysis of cesium clocks that exhibit extremely low aging. The Hadamard deviation adaptively removes the very large drift characteristic of rubidium clocks. These frequency stabilities are computed as a function of sample time to determine the long-term and short-term characteristics of the clocks. Because the Navstar clocks are expected to operate on orbit for a period of years, an analytical method of determining frequency stability history was developed to detect non-stationary behavior and examine frequency stability as a function of time [5]. The frequency stability history is obtained by performing an N-day moving average of the sequence of squared first differences (Allan deviation) or squared second differences (Hadamard deviation) of frequency offset measurements separated by the sample time. The stabilities are computed for a selectable sample time, window width, and data span. The sample times of interest are usually 6 hours and day, and the window width is chosen to be approximately ten times the sample time. The data spans may be over months or years, depending on the type of analysis being performed. MONTHLY FREQUENCY STABLTY RANKNG Six-hour and lday frequency stabilities are calculated for each Navstar clock using month of data. The stability ranking of all operational Navstar SV clocks for the month of September 200 is presented in Figure 8. The chart shows ranking as determined by the -day frequency stability using the Hadamard deviation in the front row. The back row shows the corresponding 6-hour stability for comparison. The table is provided for quick reference to the stability values. The values are listed in order of SV number with the identified clock type and serial number. The Navstar 5 rubidium clock was the most stable clock in the GPS constellation, with a lday stability of.29 pp0i4 for the month of September. The Navstar 35 cesium clock had the least stable -day performance, with a stability of.448 pplo3. All Navstar SV clocks had a drift-corrected -day frequency stability of less than.5 pplo3. The 6-hour frequency stability estimates are more sensitive to measurement noise and systematic effects that may be present in the data. Significant differences can be observed between the 6-hour frequency stability estimates for all Navstar SV clocks as compared to the -day estimates. The 6-hour frequency stability estimates varied from less than 6.5 pp04 to 9 pplo3. 47

FREQUENCY STABLTY HSTORY The frequency stability history of the Navstar 23 cesium clock for a 6-hour sample time is presented in Figure 9. t is computed using the Hadamard deviation with a window width of 96 hours. The 6-hour stability shows a periodic degradation that correlates with the spacecraft vehicle eclipse seasons. A gradual degradation is also seen from a few ~ plo'~ to almost pplo'* over the span of 6 years. The -day sample time frequency stability history in Figure 0 shows a small degradation over the same span and in fact shows good performance of the clock with a stability better than pp03. This implies there is some systematic behavior in the 6-hour case that is not present for a sample time of day. BLOCK R TMNG SGL ALYSS The -day average drift offset of Navstar 4 is presented in Figure with the times of the TKS upload and four unknown frequency breaks labeled on the plot. n the spacecraft timing signal analysis, corrections are not made for unknown anomalies. Hence, no corrections were applied for the unknown frequency breaks in these data. nitially, the corrected -day average drift was -8.4 pp04 per day, with a slope typical of rubidium clocks following activation. On October 200, the drift history approached a value near 2 pplo4 per day. mmediately following a frequency break of 9. ~plo'~, which occurred on 27 October 200, the drift history shows a small dip and then returns to its nominal drift value. The frequency stability history for Navstar 4 is presented in Figure 2, and was computed using the Hadamard deviation for a -day sample time with a loday window width. The -day stability history is well behaved, with all values less than 2 pp04 from activation until the end of October 200. The stability jumped almost an order of magnitude, degrading to a value near.2 pplo3. Following the occurrence of the large frequency break of 9. ~plo'~, the -day frequency stability recovered to a value near.5 pp04. The lday average frequency offset history of Navstar 43, from activation on 9 August 997 to 2 October 200, is presented in Figure 3. Frequency and phase breaks of unknown origin that were detected in the analysis are shown on the graph. To date, a total of six phase breaks and 33 frequency breaks have occurred. The phase offset residuals to a 30" order regression for the Navstar 43 timing signal are presented in Figure 4 using data from 22 January 999 to 20 October 200. A significant change in the phase-offset characteristics began with the first unknown frequency break of 7 pp04 that occurred on 3 July 999. The frequency breaks exhibit a sharp cusp in the phase-offset residuals at the time of the break and a rounded shape at the midpoint between successive breaks. The phase-offset residuals reached a maximum on 2 October 999, which corresponds with a frequency break of -0 pp04. Figure 5 expands the portion of data from 2 May 200 to 2 October 200. Phase offset residuals to a quadratic show that the separation between successive breaks is fairly consistent with a period of 26 days. The period of the frequency breaks and a linear fit to the data is presented in Figure 6. The linear fit exhibits a negative trend, with an overall decrease in the period of 6 days. The cumulative magnitude of the 33 unknown frequency breaks is presented in Figure 7. The breaks began almost two years (688 days) after activation, and all but one has been positive with approximately the same magnitude. As of 2 October 200, the sum of unknown frequency breaks is approximately.65 pplo'*. Frequency breaks of unknown origin have also been observed in one of the two Navstar flight candidate RAFS that has been undergoing a life test in the Precision Clock Evaluation Facility at NRL. The life test for RAFS Serial Number 028 began in April 997. A total of 76 frequency breaks have been detected since the life test began. The cumulative magnitude of the breaks is presented in Figure 8. The characteristics of the data here is similar to that of the Navstar 43 data. Until mid-year 2000, the frequency breaks were predominately positive and have roughly the same magnitude. From that time to present, the 472

RAFS SN028 clock has exhibited a few large negative frequency breaks with significant periods of no breaks at all. This is a dramatic change in the characteristics of the clock. The -day average drift offset for Navstar 44 with respect to the DoD Master Clock is presented in Figure 9. The TKS drift upload that was made by the MCS on March 200 is indicated near the center of the plot. Periodic behavior with a deviation of about f 2 pp0l3 per day is observed in the drift. A Power Spectral Density (PSD) of the -day average drift offset was computed monthly and is presented for periods of maximum and minimum drift deviation. Figure 20 is the PSD where the drift deviation is a maximum and shows a peak at 2.5 days. Figure 2 shows the PSD with a less predominant peak at 3 days for the span of time where the drift deviation is a minimum. The period of the drift deviation taken at monthly intervals is plotted in Figure 22 and is a good fit to a linear regression with an increasing trend. The data show an increase of almost day in the period of drift deviation in the time span of about a year. The 4-day average drift offset of the Navstar 46 timing signal with respect to the DoD Master Clock is presented in Figure 23 for the time span 25 December 999 to 8 October 200. A smaller drift deviation than was exhibited in Navstar 44 is seen with a period of about 2 days. The large increase in drift at turn-on is the activation transient that is typical in rubidium clocks. Figure 24 shows the drift deviation on an expanded scale, beginning with data after the activation transient. The drift deviation is f 2 pp04, an order of magnitude smaller than that of Navstar 44. The Navstar 54 timing signal -day average drift offset with respect to the DoD Master Clock is presented in Figure 25. A TKS drift correction upload was made by the MCS as shown on 24 July 200. A total of four frequency breaks of unknown origin have been detected in the Navstar 54 timing signal since activation. The first break occurred 28 days after activation on 7 September 200. At that same time the drift exhibited changes from previous normal behavior. A drift variation can be seen in September and again in October of almost pplo3 per day. This variation correlates with a degradation in the -day stability during that same time span, as shown in Figure 26. nitially, the -day frequency stability was near.9 pp0l4 during the first month of operation and approaches pp04 until the onset of the unknown frequency breaks and drift anomalies. The frequency stability history performance degrades, with a peak that exceeds 4 pp04 at the end of the data span. The frequency-stability profiles using the Allan deviation are shown in Figure 27 for all 28 clocks in the GPS constellation, as of October 200. The effect of the large drift characteristic of the Block MA rubidium clocks,can be seen in the rapid degradation in stability after a -day sample time. Frequencystability profiles are shown in Figure 28 for the same 28 Navstar clocks, but corrected for drift using the Hadamard deviation. The profiles vary by an order of magnitude for the different technology clocks. The best clocks are the Block R rubidium with a lday stability of about 2 pp04. The profiles for Navstar 22 and Navstar 30 at the top of the graph are poor because the clocks have been activated recently and they are probably still undergoing initial turn-on transitions. The stability profiles for the Block R Navstars are shown in Figure 29 for closer examination. Five of the six clocks (Navstars 4,43,46,5, and 54) have frequency stability estimates at day that are well below the GPS Block W specification of 6 pp04. The Navstar 44 rubidium clock, with a frequency stability of 6.4 pp04 at day, is close to the specification. The Navstars 44 and 46 rubidium clocks show a similar degradation in stability sample times of 3-6 days. The stability profile for Navstar 5 stands out as the best by far, with a 6-day stability approaching 5 pplo5. A summary of the systematic behavior found in the Navstar Block R timing signals is presented in Figure 30. 473

CONCLUSONS The -day frequency stability of all of the Navstar clocks in the GPS constellation, as of October 200, was better than.5 pp03. The best clock was the Navstar 5 Block R rubidium, which had a -day stability of.3 pp04. The average lifetime for all currently activated Navstar atomic clocks is 4.4 years. The oldest active clock in the constellation is the Navstar 3 Block cesium, which has been continuously operated for more than 2 years. The oldest Block R rubidium clock is the Navstar 43 rubidium (S/N 6), which has been operated for more than 4 years. REFERENCES [l] S. Hutsell, W. Reid, J. Crum, H. Mobbs, and J. Buisson, 997, Operational Use of the Hadamard Variance in GPS, in Proceedings of the 28th Annual Precise Time and Time nterval (PTT) Applications and Planning Meeting, 3-5 December 996, Reston, Virginia, USA (U.S. Naval Observatory, Washington, DC), pp. 20-24. [2] W. Reid, 2000, Multiple-PathLinked Common-View Time Transfer, in Proceedings of the 3lSt Annual Precise Time and Time nterval (PTT) Systems and Applications Meeting, 7-9 December 999, Dana Point, California, USA (U.S. Naval Observatory, Washington, DC), pp. 43-53. [3] W. Reid, 997, Continuous Observation of Navstar Clock OfSsetfrom the DoD Master Clock Using Linked Common View-Time Transfer, in Proceedings of the 28th Annual Precise Time and Time nterval (PTT) Applications and Planning Meeting, 3-5 December 996, Reston, Virginia, USA (U.S. Washington, DC), pp. 397-408. [4] M. Epstein and T. Dass, 2002, Management of Phase and Frequency f o r GPS R Satellites, in these Proceedings. [5] T. McCaskill, 997, Analysis of the Frequency Stability History of GPS Navstar Clocks, in Proceedings of the 997 EEE nternational Frequency Control Symposium, 28-30 May 997, Orlando, Florida, USA (EEE Publication 97CH3606), pp. 295-303. 474

VSTM 43 TlMlNG SGL ONE-DAY AVERAGE FREOUENCY OFFSEJ from od)meaer clock Using RAFS Sene No 6 S-AUG-97 to 2-OCT-0 VSTAR 43 TlMlNG SGL 3O-MlNLl7E PHASE OFFSET RESDUALS from OOD Master clock Uang RAFS senel No 6 22-JAN-SS to 20-OCT-0 Figure 4. Figure 3. VST4R 43 TlMlNG SGL W E OmET from h h t W clock h, l g RAFS Send No 6 2-MAY-0 to 2-OCT-0-24 6 2 W UNEAR FT AND PEROD OF VSTAR 43 REpErmVE FREOUENCY BREAKS from Do0 Master C/o& Uang RAFSsenelNo 6? 52070 62090 520 6230 6260 6270 MODHE0 JUUAN DATE 6290 O 9 5220 Figure 6. CUMULATVE CORRECTlON N THE FREOUENCY OFFSET OF RAFS S/N 02s Naval Research Laboratory Recismn CMck EvahJabonFmhry Weeks -236 CUMULATVE CORRECllON N THE FREWENCY ORSET of VSlAR 43 TMNG sk;hw RAFsserlalNo E 60800 6rdoo Si00 6 b 6E00 9 MoDlFlED &LMN DATE Figure 5. G O O, 6300 6400 65JO 6800 6700 6800 6900 52WO 6200 62200 6800 62doo 60200 604oO 508oO 50800 6xx)O 6200 6400 6 W 6800 62000 62200 62iOO MODFED JJUAN DATE MoDlFlED JUUAN DATE Figure 8. Figure 7. 477

SPECTRUM OF VSTAR 44 77MlNG ONE-DAY AVERAGE DRFT OFFSET from Do0 Master Clock Using RAFS Serial No. 9 6-AUG-00 to 0-SEP-00 VSTM 44 TMNG SGL CORRECTED ONE-DAY AVERAGE DRFT OrrSET from 000 Master Clock U m g RAFS Serial No. 9 l6-aug-00 0 0-OCT-0 f p & ~ y Yx x ) q + 000, 8 li -. -6 00-8 6760 6800 6860 6900 6960 62000 62060 6200 6260 52200 MODFED JUUAN DATE 03 QA 06 08 0. 0 8 09 6600 6TKx) 0.7 00. 09 Figure 22. ~ ' ' ' ' ' ~ " ' t ~, ', 6800 06 MODFED JULAN DATE VSMR 46 TMNG SGL FOUR-DAY AVERAGE DRFT OFFSET from Dd) Master Clock U s r n ~ RAFS Sarial No 2 RAFSSenalNo 27 26DEC-99 to 8-0'3-07 6530 05. 0 6l760 6800 6850 6900 6960 6'2000 6'2060 6200 6'260 67.200 0 VSJAR 46 TMNG sk;l FOUR-DAY AVERAGE DRFT OFFSn from Dd) Master COck Uang ' 64 PEROD OF MOMHLY SPECTRAL PEAK USNG ONE-DAY AVERAGE DRFT OmET Figure 2. J 03 FOR VSAR 44 RUBUDUM CLOCK (Sh'9) FREQUENCY (CYCLES PER DAY).. ;. 02 Figure 20. SPECTRUM OF VSTM44 77MlNG ONE-DAY AVERAGE DRFT OFFSET from DoD Master Clock Usrng RAFS Serial No 9 0-MAY-0 to 0-JUN-0 02 0 FREOUENCY CYCLES PER DAY) Figure 9. 0.0 0 6900 26-FEE-00 to 78-00-07, 62000 6200 i 67.200 MODflED JULAN DATE Figure 23. Figure 24. 478

VSTAR 64 TMNG SGL ONE-DAY AVERAGE DRFT ORStT fmm Do0 Mestar Clock Uwng RAFS Senel No 44 -FEE-0 to 8-OCT-0 FREOUENCY STABLTY HSTORY OF VflM 64 TMNG SGL OFFSET f r m DOD Master Clock Using RAFS Serial No 44 S e W Tune day Window Width 0 days 200,.. 6QM) 62060 62000 -.. 6200 moo 6250 6960 62060 62iQO MODRED JUUAN DATE sxxx) MODFED JUUAN DATE Figure 26. Figure 25. FREQUENCY STABLTY OF VSTAR CLOCKS (HADAMARD) April 200 to October 200 FREOUENCY STABlLTY OF VSTAR CLOCKS (ALLAN) April 200 to October 200 4 8 - -v 0 8 r r - T F 4 m 8 6 7? - unarm, eam,. 0 SAMPLE TlME (DAYS) / Figure 28. FREQUENCY STABlLlTy OF VSAR Block R Timing Signal Systematics BLOCK R RUBlDlUM CLOCKS (HADAMARD) April 200 to October 200 [ Satellite SVN 4 SVN44 SVN44ATP SVN46 Anomalv Drift SVN 5 0-'6L, 0 SAMPLE TME (DAYS) 0 00 Figure 27. 0 6250 0 00 SVN 54 SN 28 /Lab) SN30(Lab) SAMPLE TlME (DAYS) Figure 29. Drift Drift Drift Periodic Total Breaks 2.07days 2.52davs None None None 47 davs Figure 30. 479.62days Drift Freauencv 2.32days. 7 6 ~

QUESTONS AND ANSWERS JM CAMPARO (The Aerospace Corporation): Do you have any idea what those anomalies might be due to that you are seeing? JAY OAKS: No. We present all the data and try to correlate it with things that we know about, but we have not come up with anything, and we presented it more or less to the working community to try to come up with reasons. Again, most of these things that we are seeing are below what the operational people care about. But they are interesting to us because there are things that if we want to improve clocks, if we can figure out what s causing them, maybe we can make clocks better in the future. VCTOR RENHARDT (Boeing Satellite Systems): Do you have access to the telemetry data from the spacecraft to see if it is correlating with environmental things? OAKS: We have been off and on about trying to get the telemetry data. One of the problems with the current telemetry data is granularity or resolution - we cannot really correlate very well. With the laboratory clocks, we get that down pretty well, and we can correlate, the anomalies that we see in the clock to all of the different things that are monitored. But for the most part, when we ve tried to look at the telemetry data, it is too coarse to do a lot of good. REWHARDT: Yes, most telemetry data are okay if you don t have to use them. 480