zotero/storage/ZYFCE23H/.zotero-ft-cache

2547 lines
48 KiB
Plaintext
Raw Permalink Normal View History

2024-08-27 21:48:20 -05:00
I I I
I
I I
I
I
I
I
I
I
I I I I I I
'I
PERFORMANCE AND RESULTS OF PORTABLE CLOCKS IN AIRCRAFT
by J.C. Hafele
Professor Hafele is an Assistant Professor of Physics, Washington University, St. Louis, Mi.ssouri..
1. 0 INTRODUCTION During the first two weeks of October, 1971, R. E. Keating of the U.S.
Naval Observatory and I flew twice around the world on regularly scheduled commercial flights, once eastward and once westward, with four HewlettPackard 5061A cesium beam clocks. For about a week before the first trip, between the trips, and for about a week after the second trip we recorded a continuous phase comparison between each clock and the MEAN (USNO) time scale. Time differences between the four clocks were also recorded with a time interval counter at regular intervals before, during, and after each trip, thereby permitting evaluation of the flying mean time scale. The experiment was conducted for two reasons: (1) to compare the known performance of similar clocks under fairly well controlled laboratory known performance of similar clocks under fairly well controlled laboratory conditions, and (2) to try to detect relativistic effects on the time recorded by clocks during terrestrial circumnavigations. It was performed with the interest and complete cooperation of Dr. G. M. R. Winkler and others of the Naval Observatory Time Service Division, and with financial support from the Office of Naval Research. A major portion of the credit for the experiment goes to Richard Keating, whose steadfast desire to perform such an experiment generated the necessary motivation. Although the
-261-
I
I
performance of these clocks was found to be noticeably degraded during the
trips, we also found that the ensemble of four flying clocks produces con-
I
sistent time differences which are difficult to explain without invoking
relativistic effects.
I
2.0 A PRECISION LONG-TERM TIME REFERENCE - MEAN (USNO)
I
When a reliable means of counting the accumulated number of
periods of an oscillator is combined with a cesium beam frequency standard,
I
the assembly approaches the realization of an ideal standard clock. In
fact, a time interval of one second is now, by definition, exactly 9 , 19 2, 631, 77 O accumulated periods of an "ideal" cesium beam frequency standard. 1 Experi-
I
ence shows, however, that the times recorded (periods accumulated from a
I
common starting instant) by two "real" cesium beam clocks having no inten-
12 tional frequency offset differ by as much as several parts in 10 . Though
I
this performance is truly remarkable, certainly beyond any expectations of
a decade or so ago, for many purposes it is desirable to establish a stable,
I
continuous atomic time scale which is independent of the characteristics
of any particular cesium beam clock. It is not difficult to accept the premise
I
that the long-term average or mean of the times indicated by two or more
properly functioning cesium beam clocks is more reliable (or stable) than the
I
individual times indicated. Ideally, the mean of a very large (infinite)
number of clocks would provide a perfectly reliable time scale; however,
I
budget limitations prevent realization of this goal. So we compromise by
taking the mean of an ensemble of as many clocks as possible, or as is
I
necessary to produce the desired reliability. The mean time scale at the
U.S. Naval Observatory (MEAN (USNO)) represents a suitable average of the
I
1Winkler, G.M.R., R.G. Hall, and D.B. Percival (1970), "The U.S. Naval
I
Observatory Clock Time Reference and the Performance of a Sample of Atomic
Clocks, " Metrologia Vol. 6, p. 12 6.
I
I
-262-
I
I
c:
I
of the time recorded by an ensemble of more than 15 cesium beam clocks. 2
I
As some of these clocks deteriorate or fail, they are replaced by standby
clocks, thereby permitting an indefinite extension of the time scale into
I
the future. The time indicated by MEAN(USNO) is believed by many to be
the closest approximation we have today to an ideal long-term atomic time
I
scale.
The basis for confidence in this time scale is the belief that both
I
long- and short-term differences (or fluctuations) between the mean and
each clock of the ensemble are random. If one clock of the ensemble
I
increases its rate (relative to the mean), it will not be long before some
other member on the average decreases. In this way, ultimately, the mean
I
is correctly maintained. In fact, the stability of the mean is actually
improved by taking into account rate changes for each member of the
I
3 ensemble.
The MEAN (USNO) time comes from a "paper" clock; that is, the
I
mean is calculated after the fact and is not instantaneously available for
time comparisons. This presents no real problem, however, because one of
I
the members of the ensemble is chosen as the primary reference or "master"
for the ensemble and all clocks are compared with this master. A secondary
I
master is kept in reserve in case the primary master fails. This MASTER (USNO)
time is tracked to follow the mean as closely as possible and is physically
I
available for instantaneous time comparisons. Since the difference between
MASTER(USNO) and MEAN(USNO) is calculated at regular intervals, it is an
I
easy matter to convert to an equivalent comparison with MEAN (USNO) after
the intercomparison measurements are completed.
I
I
I
I
-263-
I
I
I
3. 0
LABORATORY PERFORMANCE OF CESIUM BEAM CLOCKS
Now that we have a long-term reference time scale, it is straight-
I
forward to evaluate the performance of any particular clock simply by comparing the time it indicates with the time of the reference. There is
I
extensive literature on evaluation of the performance of precision atomic oscillators, particularly on the subject of short-term time and frequency
I
stability (see Barnes, J.A., @t al and Allan, D. W., et al and bibliographies therein). 4 There seems to be a consensus in this literature that
I
time domain stability is best described in terms of fractional frequency 2
fluctuations, with the quantitative measure being the Allan variance cr (r).
I
Cesium beam clocks, like all clocks to a greater or lesser extent, are subject to short-term FM noise and this noise limits the precision with which the time
I
or frequency can be determined during short averaging times ( several c.ays or less). In most cases, however, this type of noise is not the limiting factor
I
for long-term stability. Long-term performance appears to be governed more
,1
or less by instantaneous and randomly occurring rate (or frequency) changes.
The characteristic long-term performance of cesium beam clocks can be described in terms of u (T) , the square root of the Allan variance. As the
I
averaging time T increases from zero, u ( T) characteristically decreases
on a log-log plot with a slope of -½. However, beyond a certain T , depend-
I
ing on the particular clock, the slope changes rather abruptly to zero, and for
I
longer averaging times u(T) remains constant and in some cases eventually
4Ibid.; Barnes, J.A., et al. (1971) "Characterization of Frequency Stability, 11
I
IEEE Trans. Instr. Meas. IM-20 (May 1971), p. 105; Allan, D.W., B.E. Blair, D. D. Davis, and H. E. Machlan (1971) 11 Precision and Accuracy of Remote
I
Synchronization via Portable Clocks, LORAN C, and Network Tele vis ion
Broadcasts," Proc. 25th Ann. Symp. Freg. Control, 26-28 April 1971 (Electronic Industries Assoc.) p. 195 .; Martin, D. (1971) "Frequency Stability Measure-
I
ments by Computing Counter System, 11 Hewlett-Packard Journal Vol. 23
(Novemberl971) p. 9.
I
I
-264-
I
I c:::
I
5
increases.
The averaging time for which
u ( T) decreases with slope -½
I
gives the time for which the frequency fluctuations are random about a constant mean. If the frequency fluctuations remained indefinitely random
I
about a constant mean; that is, if u ( r) continued to decrease indefinitely
with slope -½, one could achieve indefinite long-term stability with a single
I
clock simply by increasing the averaging time. In other words, if this were the case, once the average rate relative to MEAN(USNO) was established
I
the average rate could be projected indefinitely into the future. This is not
the case, however.
I
This change in the slope for CT ( r) indicates that "quasipermanent"
changes in the average rate for the clock occur randomly with the average
I
time interval between changes comparable to the averaging time for the slope
change. By quasipermanent I want to suggest,or imply, that the clock changes
I
(unpredictably) to a new average rate every so often. Uncertainty in the
onset and magnitude for each new quasipermanent rate is related to the
I
minimum value of the Allan variance.
Table I illustrates this long-term performance. The table lists
I
quasipermanent rate changes for five cesium beam clocks over a 40-day
period at the U.S. Naval Observatory. The data were taken from continuous
I
phase comparison records with MEAN(USNO) as reference. Approximate rate
changes are relative to MEAN (USNO). N is the number of changes during this
I
period, the '1 R's are the signed magnitudes of the changes, 6 R is the
average rate change, and T is the approximate average duration between rate
I
changes.
I
5 Winkler, G. M. R. , R. G. Hall,and D. B. Percival; Winkler, G. M. R. (1970) "Requirements and Performance for Today's Atomic Standards, " PTT!
I
Strategic Planning Conference, Vol. 1, December 19 70, p. 12 9.
I
I
-265-
I
I
I
TABLE I. QUASIPERMANENT RATE CHANGES FOR FIVE CESIUM BEAM CLOCKS
CLOCK (s. no.)
N (in 40 d)
6R (nsec/day)
-
-
6R
T
I
(nsec/day) (days)
140
3
-25; +28; -54
3911
I
-17
1471
2
+22; -19
-17
10
I
-17
20
+l
13
I
60
3
+17; -12; +17
+7
10
283
6
-33 I +21; -50; +57; -41; +22
-4
7
I
The cause for these quasipermanent rate changes is not understood in detail. 6 Both the time of their occurrence and their direction appear to occur at random.
I
Though they are quite unpredictable in nature and, of course, must average to zero over sufficiently long periods with properly functioning clocks (very many
I
T periods) they limit the precision with which the rate of any particular clock can be projected into the future.
I
As can be seen in the table, the average time between changes and their magnitudes depends on the particular clock. T for clocks at USNO
I
has been found to range between 1 and 40 days. 7 "Good" clocks are those which go longer than average between rate changes. One of the purposes of
I
the experiment was to study deterioration, if any, in the performance of four good clocks under traveling conditions.
I
4. 0 RELATIVISTIC EFFECTS DURING TERRESTRIAL CIRCUMNAVIGATION
I
In 1905 Einstein laid a radical new basis for the concepts of space and
time. Though Newton's absolute time had proved adequate for most practical
I
6 Ibid.
I
7 Winkler, G.M.R., R. G. Hall, and D.B. Percival (1970) "The U.S. Naval Observatory Clock Time Reference and the Performance of a Sample of Atomic
I
Clocks, " Metrologia, Vol. 6, p. 12 6.
I
I
-266-
I
I
I
I
purposes, Einstein produced convincing arguments against it. Absolute time contains an element of mystery which is incompatible with precisely
I
defined scientific quantities. Consequently, Einstein defined a new empirical basis for time by accepting a definition which states, in effect,
I
that "time is that which is indicated by a clock," and then proceeded to develop his relativity theories on that basis. Einstein's relativity has
I
proved to be completely compatible with all relevant observations; in fact, no definitive test ever performed has disproved it. The results of our flying
I
clock experiments, at least at the present state of analysis, offer no exceptions.
I
The special theory of relativity predicts that a moving clock will run slow compared with similar clocks distributed at rest and suitably synchron-
I
ized in an inertial reference system. In addition, the general theory predicts that a clock in a stronger gravitational field will run slow compared
I
with a similar clock in a weaker field. Most physicists believe these pre-
dictions have been verified through studies of lifetimes of elementary parti-
I
cles and studies employing the Mossbauer effect, but a small vociferous
minority and indeed, I would say, most of the general public, are reluctant to
I
accept the prediction that these effects also apply to ordinary time indicating
clocks, particularly to biological clocks. For an excellent book devoted
I
almost entirely to the "clock paradox, 11 see L. Marder' s "Time and the
Space Traveler, 11 printed by Allen and Univin, London, 1971.
I
Relativistic effects on times recorded by clocks are vanishingly
small at ordinary speeds, certainly too small to detect with ordinary
I
chronometers. However, the development of portable cesium beam clocks, with stabilities of better than 1 part in 1012 over a period of days, makes a
I
direct test with time recording clocks feasible. I have shown that the pre-
dicted kinematic effects on the time recorded by clocks during equatorial
I
circumnavigations at ordinary jet aircraft speeds and altitudes are enhanced
by the Earth's rotation to a level which is probably above the threshold for
I
-267-
I
I
I
detection with cesium beam clocks. 8 Although the assumption of an
equatorial circumnavigation at constant ground speed and altitude is not
I
essential, it does simplify somewhat the calculations for estimating the magnitude of expected relativistic effects. For an equatorial circumna V'.""
I
igation with constant ground speed v (m/sec) and altitude h (m), the predicted relativistic time gain for the flying clock over a similar reference
I
clock kept at "rest" on the Earth's surface is given by9
I
.ill! 2R0v + v2
c2
2c 2
(1)
I
where r and r O are the respective times recorded by the flying and
I
ground clocks; R (m' is the Earth's radius and f! (rad/sec) its angular speed; g (rr,/sec 2) is the surface value of the acceleration of gravity; and c (m/sec)
I
is the speed of light. In Equation 1, the ground speed is positive for eastward and negative for westward circumnavigations. At latitudes higher than
I
the equatorial plane, the term in Equation 1, which depends linearly on the ground speed and gives an east-west directional asymmetry, becomes pro-
I
portional to the cosine of the latitude, lO so for latitudes as high as 60°, the directional dependence is reduced by less than one-half. Although
I
Equation 1 is correct only to lowest order in the speed and altitude, it represents an extremely accurate approximation for jet speeds and altitudes
I
because only these lowest order terms are detectable in this case.
I
8 Hafele, J.C. (1970), "Relativistic Behavior of Moving Terrestrial Clocks," Nature, Vol. 227, p. 270; Hafele, J.C. (1971), "Reply to Schlegel,"
I
Nature Phy. Sci. Vol. 229, p. 238; Hafele, J.C. (1971),"An Empirical
Resolution of the Relativistic Clock Paradox," Bull. Am. Phys. Soc., Vol. 16, p. 611; Hafele, J.C., "Relativistic Time for Terrestrial Circumnavigations,"
I
Am. J. Phys. (in press, to appear January 1972).
9Ibid.
I
10cutler, L.S. (1970), "Correction to 'Flying Clock' Comparisons Extended to
East Europe, Africa and Australia," Hewlett-Packard Journal, Vol. 21 (March 1970), p. 10.
I
-268-
I
I
c:::
I The actual time gain ~ r is a bit more instructive than the time
I
.......
ratio of Equation 1 and it follows from multiplication of Equation 1 by r0 •
rr:
i:::::
Because standard clocks keep the same time while sitting on the ground
I
anywhere on Earth (at average sea level and to this order of approximation), 11
only the actual time in flight during a trip contributes to relativistic effects.
I
However, ground time does contribute to increasing the random, unpredictable
time offset and therefore to the threshold for detection of relativistic effects.
I
(Relativistic effects were not detected during previous flying clock trips
because they accumulate only while the clocks are in flight, and for those
I
trips most of the time was spent on the ground.) Suppose for the moment
that ground time, for example, for refueling stops, is negligibly small compared
I
with the time it takes to fly around the world. Then the time recorded by the
ground clock during the circumnavigation is given by
I
I
I
Solving for AT and inserting this value for T O in Equation 1 gives
I
21rRgh _L 2 lvl
C
..::L IVI
-1r.21L JVJ •
C
I
Figure 1 is a graph of this equation showing AT versus v for altitudes of O, 10, and 20 kilometers. The dots in Figure 1 correspond to the cruising
I
speed and altitude for a Boeing 707. 12
I
11Cocke, W. J. (1966), "Relativistic Corrections for Terrestrial Clock
Synchronization," Phys. Rev. Letters, Vol. 16 (19 66) p. 662. ·
I
12 Schaefer, H. J. ·(1971), "Radiation Exposure in Air Travel," Science,
Vol. 173 (1971), p. 780.
I
I -269-
I
0 0 0
0
QJ (I)
..E..
~
PREDICTED RELATIVISTIC TIME GAIN .1 T
vs GROUND SPEED v
ALTITUDE LATiTUDE
= h = o0
EASTWARD v >0
WESTWARD v < 0
11111mm .1 T THRESHOLD
I
N
-..J
v(m/sec)
0
I
0 0
LO
I
Figure 1. PREDICTED RELATIVISTIC TIME GAIN 6:r vs. GROUND SPEED v
FOR EQUATORIAL CIRCUMNAVIGATIONS
I
c-·
I
Personnel of the Naval Observatory's Time Service Di vision have
I
a considerable background of practical experience with cesium beam clocks,
both in the laboratory and under traveling conditions. Their experience is
I
that at the end of trips the time difference between the flying clocks and
MEAN{USNO) is random with a zero center gaussian distribution having a
I
spread
of
about
60
nanoseconds
per
day
of
13 trip.
Since this spread was for
all portable clocks, "good" clocks may be expected to provide a somewhat
I
lower threshold. A conservative estimate for the threshold of detection of
relativistic effects with "good" clocks is given by
I
I
D-T
threshold
60 nsec
day
x T O (days) •
I
The hatched lines in Figure 1 represent this threshold. Figure 1 shows there
is a variety of available conditions with jet aircraft which predict relativistic
I
time differences exceeding this threshold.
Commercial around-the-world flights do not, of course, maintain
I
constant altitude, latitude, or ground speed. In this case, it is necessary
to perform a numerical integration of the relativistic equations. The necessary
I
calculation is given by
I
- ~light
1 2c 2
(2R!lcosHr) cos6(r) v(T) +v2(r))]
dr,
I
path
(2)
where, for each interval of the summation, A is the latitude, (J is the
I
azimuth or bearing of the plane's velocity relative to east, and the rest of
the symbols have the same meaning as for Equation 1 (v is the unsigned
I
magnitude of the ground speed in Equation 2; the azimuth 8 accounts for the
direction) •
I
13 Winkler, G.M.R., R.G. Hall, and D.B. Percival; Winkler, G.M.R. (1970)
I
"Requirements and Performance for Today's Atomic Standards, " PTTI Strategic Planning Conference, Vol. 1, December 1970, p. 129.
-271-
I
I
I
5. 0 PREDICTED RELATIVISTIC TIME DIFFERENCES
Necessary data for each around-the-world flight were provided by
I
the various flight crew captains. (We are greatly indebted to Pan Am, American Airlines, and TWA, and their flight captains, for this extraordinary
I
service, and to Pan Am and TWA for sending a company official as escort.)
I
In most cases the flight paths were traced on appropriate flight maps with
the time (GMT), altitude, and ground speed recorded regularly at navigation
I
check points. The least accurate quantities are probably the times over the
check points, which appear to have been recorded to only the nearest three or
I
four minutes in some cases. Our own records of liftoff and touchdown times,
recorded to the nearest minute, are in good agreement with the corresponding
I
times provided by the flight captains.
The information on the flight maps divided the eastward circumnaviga-
I
tion into 125 intervals and the westward circumnavigation into 108 intervals.
The latitud8 and longitude with the corresponding time for each check point
I
permit calculation of the average latitude, azimuth relative to east, and ground
speed for each interval. The average altitude for each interval was estimated
I
from the altitudes recorded at the end points. These calculations and the
summation indicated in Equation 2 were carried out with the Washington
I
University IBM 3 60/50 computer. Table II summarizes the data for each
direction, and gives the predicted relativistic time gains based on Equation 2.
I
6 .Q AMBIENT CONDITIONS
In addition to recordin~ time differences between the four flying
I
clocks at regular intervals during each trip, we recorded at irregular intervals the temperature, pressure, and relative humidity with a small
I
aneroid desk-top weather station; the compass direction with a small pocket
I
compass; and the ambient magnetic field with a small pocket magnetometer.
Before and after each trip the clocks were kept in an unused darkroom at
I
I
-272-
I
I
c::
I
TABLE II. SUMMARY OF AROUND-THE-WORLD FLIGHTS AND PREDICTED
RELATIVISTIC TIME GAINS FOR FLYING CLOCKS
I
r··r
EASTWARD
WESTWARD
I
Day GMT Location
04 1 1930
2 USNO D
Day GMT Location 13 1940 USNO D
I
05 0012 Dulles D (Pan Am 747)
232,2 Dulles D (TWA 707)
0656 London A3 (*Pan Am 707) 14 0400 Los Angeles (*TWA 707)
0814 D
0503 D
I
0909 Frankfurt A l 036 D
1248 Istanbul A
1014 1313 2015
Honolulu A D Guam A
I
1357 1513 1619
D Beirut A D
2113 D 15 0006 Okinawa A
0107 D
1813 Tehran A
0209 Taipei A
I
1940 D 2241 New Delhi A
0303 D 0413 Hong Kong A
06 0000 D
1248 D
I
0333 Bangkok A 0513 D
1514 Bangkok A 1632 D
0745 Hong Kong A
2006 Bombay A
I
0855 D 1216 Tokyo A (*Pan Am 747)
2115 D
16 0403 Tel Aviv A (*T'-.M.ll 707) I
1432 D
0509 D
I
2110 Honolulu A 2314 D 07 0350 Los Angeles A (*AA 707)
0645 0733 0903
Athens A D Rome A
0447 D
1001 D
I
0713 Dallas A 0753 D
1138 Paris A 1425 D
0959 Dulles A
1557 Shannon A **
I
1255 USNO return
1706 D 2338 Boston A
I
17 0118 D 0226 Dulles A
0400 USNO return
Trip time
65 .42 hours
80. 33 hours
I
Avg ground speed 243. meters/sec
Avg altitude
8. 90 kilometers
Avg latitude
34. degrees N
218. meters/sec 9. 3 6 kilometers
31. degrees N
I
Rel time gain 1 October, 1971
-40. nsec (loss)
+275. nsec
I
2 3
D
-
Depart
A - Arrive
* Indicates clocks transferred to a different aircraft; 707 and 747 indicate
Boeing aircraft type.
I
** Indicates an unscheduled fuel stop.
-273-
I
I
I
the Time Service Building where the temperature was held between 68 and
73°F. Temperatures during the trips averaged about 75°F with extremes
I
between 70 and 85°F. During each flight between landing points, the cabin
pressure dropped from 30 to 25 in Hg, except in one or two cases where it
I
was not lowered that much. For each trip the relative humidity dropped
from 60% to about 30o/0. Noticeably rough flying conditions occured only
I
during the flight between Athens and Rome on the westward trip. At no time
during either trip did any of the clocks lose regulation control; that is, the
I
lights indicating momentary loss of feedback control did not light up.
The four clocks were carried as two assemblies, each assembly
I
consisting of a bottom clock, a middle H-P K02-5060A battery pack and
charging unit, and a top clock. To suppress any possible magnetic coup-
I
ling between the clocks, additional magnetic shielding was placed between
each clock and the battery pack of each assembly. Each assembly was
I
rigidly bolted together.
For the eastward trip, the clocks were oriented in the airplane with
I
the front of the clocks towards the front of the plane between Washington
and London and between Tokyo and Washington, and with the front of the
I
clocks towards the back of the plane from London to Tokyo. The clocks
were oriented with the front of the clocks towards the front of the airplane
I
during the entire westward trip. The clocks were strapped in seats in the passenger compartment.
·I
When we discovered that one of the four batteries in one of the
battery packs was dead, there was much concern on several occasions
I
when electric power for charging the batteries was not available, because this abnormal situation caused the batteries to discharge rapidly. Fortu-
I
nately, two of the clocks were never in jeopardy on this account, and our intercomparison data indicate no noticeable effect from low battery voltage
I
for the clocks in question.
I
-274-
I
I
I
I
Near the end of the data period after the clocks were returned to
I
the Naval Observatory, two environmental tests were performed. First,
rr:
r.::::
the clocks were reoriented to see if changes in the direction of the
I
Earth's magnetic field caused noticeable changes in the clock rates; none
were found. Secondly, line power charging the battery pack with the dead
I
battery was interrupted for one hour to see if discharge of the remaining
I
three batteries caused any noticeable change in the rates; again, none was observed.
I
Although environmental studies have shown that cesium beam clocks are susceptible to such effects as temperature changes and AC magnetic
I
fields, it is important to note that no known effect consistently increases or decreases the rates for these clocks. That is, all known environmental
I
effects cause rate changes that are random in direction. This observation gives additional confidence in the mean of an ensemble. (Strong DC mag-
I
netic fields cause predictable changes, but the clocks are triply shielded against the Earth's weak field.)
I
7. 0 MEASURED TIME DIFFERENCES
A strip chart recorder with associated phase comparison electronics
I
was used to record a continuous phase comparison between each of the four
flying clocks and MASTER(USNO) before and after each of the trips. In
I
addition, at periodic intervals before and after the trips, absolute time
differences between the one second pulses or "ticks" were measured to
I
the nearest nsec with a Hewlett-Packard computing counter. 14 The time
differences between the ticks of each clock and MEAN(USNO) via MASTER(USNO)
I
provide a calibration of the continuous phase graphs. Thus the rate of each
clock relative to MEAN(USNO) before and after each trip is established.
I
14 Martin, D. (1971), "Frequency Stability Measurements by Computing System,"
I
Hewlett-Packard Journal, Vol. 23 (November 1971) p. 9.
I -275-
I
I
I
Identically the same electronic arrangement was used for all time inter-
comparisons.
I
Six intercomparison times between the four clocks (three of them are
redundant but serve as a check on the measurements) were recorded with
I
the same computing counter at hourly intervals for a short period before,
during, and for a short period after each trip, with additional readings taken
I
shortly before and after each touchdown and liftoff. These interc;:..nparison
data permit evaluation of the mean of the flying ensemble (MEAN(FLYING))
I
in much the same way as MEAN (USNO) is determined. Although evaluation
of MEAN (FLYING) is not yet completed, we expect this analysis of the data
I
to produce the highest level of confidence in our results, and we hope to
be able to report them soon.
I
The results reported here depend only on the intercomparison data
with MEAN (USNO) before and after the trips. The final analysis is not
I
expected to change these results significantly.
An early and very preliminary estimate of the results was based on
I
extrapolation of the rate for each clock immediately before the trips. The
difference between this extrapolated time (to the end of the trip) and the
I
observed time was averaged among the four clocks for each trip. Even with
this very simple approach, consistent time differences were found. But
I
. even a casual observer would be quick to criticize this approach by pointing
out that it is necessary that the rate after each trip be the same as before
I
the trip, because it assumes no quasipermanent changes occurred during the
trip. This approach is difficult to defend because quasipermanent changes
I
did occur for each clock during each trip; it represents the lowest level of
analysis and produced the least confidence. I wish to present here a
I
somewhat higher level of analysis.
I
I
I
-276-
I
I
c:::
I
Let rmean and 1c be the respective times indicated by MEAN(USNO)
I
and one of the clocks of the flying ensemble, and define the rate (difference) Ras
I
R
Tc - 1'mean =
I
C
rmean
'Y mean
I
Figure 2 shows measured values of 6:rc versus 'l'mean for each clock taken
before, between, and after the trips . (The curves in Figure 2 are displaced
I
vertically from the original data to improve clarity; this displacement of course
had no effect on the relative rates.) The vertical lines correspond to the
I
times the clocks were removed from and returned to the Naval Observatory.
Of course, no comparisons with MEAN (USNO) were possible during the trips.
I
The slope of each trace in Figure 2 gives R • A slope of 30° for this graph
C
-12
corresponds to a rate of approximately 10 nsec/hour (2. 8 x 10 sec/sec).
I
It can be seen that all the clocks agree with MEAN (USNO) to within 3 x 10-l 2 sec/sec. Notice also that numerous quasipermanent rate changes
I
occurred while the clocks were sitting in the laboratory. Moreover, rate
changes that are noticeably larger than those typical in the laboratory
I
occurred for each clock during at least one of the trips, except for clock
447. This result suggests that these clocks cannot be expected to perform
I
under traveling conditions as well as they do in the laboratory; this is not
particularly new information.
I
A careful study of the data shown in Figure 2 indicates that the rate
for each clock immediately before and after each trip can be evaluated to
I
perhaps better than 10% (except for the case of clock 3 61 after the eastward
trip). The following is an attempt to justify the approach used to derive the
I
results presented in this report.
During both trips each clock suffered one or more quasipermanent
I
rate changes. On the basis of "the least number of assumptions, assume
only one occurred. This change is then given by the difference in the rates
I
-277-
I
30001 7000
6000
361 IC - (MEAN)
5000
I
I
EASTBOUND
I
I
~
WESTBOUND
], 4000
u
I r-..:>
1-
<3
-..J
co
I
3000
120 IC - A. I (MEAN)
2000
IOOO
o! , 1 , , , I, , I , , , , l , , I :=::::::-::, ~
C
40
80 120 150 200 240 280 320 360 400 440 480 520 560 600 640
'f' MEAN (HOURS)
Figure 2. MEASURED VALUES OF f..7c VS T MEAN FOR EACH CLOCK
TAKEN BEFORE I BETWEEN I AND AFTER TRIPS
I
c::
I
I
immediately before and after the trips. Without the inflight intercomparison
data, we have no way of knowing at what time during the trip this rate
I
change occurred. The least error accrues with the assumption that it occurred
at the midpoint of the trip. Thus the average inflight rate (assuming no rela-
I
tivistic effects) is half the sum of the initial and final rates. Extrapolation
with this average rate then gives the final times the clocks would have read
I
if there were no relativistic effects. Figure 3 illustrates this approach and
defines the symbols used.
I
Table III lists for each clock and each trip the initial and final rates,
the average rate, the rate change, the observed initial and final time dif-
I
ference between each clock and MEAN (USNO), and the difference between
the extrapolated time and the observed final time difference. Notice that the
I
the average rate change for each trip is less than 1 nsec/hr (3 x 10- 13 sec/
sec). Hence MEAN (FLYING) can be expected to produce considerably more
I
reliable results . Most people {myself included) would be reluctant to agree
that the time gained by any one of these clocks is indica.tive of anything, but
I
the rather striking consistency between all four clocks must be taken seri-
ously. The averages of the four final time differences are consistently nega-
I
tive or near zero for the eastward trip and consistently positive for the west-
ward trip. Corresponding standard deviations are also listed in Table III.
I
Figure 4 summarizes our results at this stage of the analysis. The
numbers in the blocks indicate the serial numbers for the corresponding
I
clocks. The average of the observed time gains, with corresponding standard
deviations, and the predicted time gains are shown at the bottom of Figure 4.
I
It is amusing to notice that the values for the eastward trip are in excellent
agreement despite our expectation that we would not be able to detect a
I
definite effect in this direction. On the other hand, the effect for the
I
westward trip was predicted to be considerably larger and detectable, but the observed value is more than one standard deviation below the predicted
I
value. Hopefully our final analysis will improve the situation.
-279-
I
I
I
TABLE III. FLYING CLOCK DATA - RATES REIATIVE TO U.S. NAVAL
OBSERVATORY MEAN
-
I
CLOCK
Ri
Rf
R
ti R ti·q
ti T f
ti T
(ser. no.) (ns/hr) tns/hr) (ns/hr) (ns/hr) (ns)
(ns)
(ns)
I
Eastward flight
r = 65. 42 hr
I
361
+2.66 +4 .38 +3.52 +1. 72 1790 191U -110
408
-1. 78 +3.22 +0.72 +5.00 -20
30
+3
I
120
-4.50 -8.89 -6.70 -4.39 -290 -780
-52
447
-7.16 -8.41 -7.78 -1.25 -1140 -1705
-56
--
I
ti R = +0. 27 ns/hr eri1R = 3.8 ns/hr
~T = -54 ns (T = 46 ns
i1T
I
Westward flight
I
r=80.33hr
3 61
+6.89 +3.97 +5.43 -2.93 2880 3390
+74
I
408
+4.84 +2.16 +3.50 -2.68 490 980 +209
120
-8.88 -4.56 -6.72 +4.31 -2100 -2400 +240
I
447
-7.17 -9.42 -8.30 -2.25 -2840 -3390
-+11-6
I
KR = -0. 89 ns/hr
--x;- = +160 ns
(T = 2.6 ns/hr i1R
(TtiT = 78 ns
I
-13
1 ns/hr = 2. 78 x 10 sec/sec
I
(T ~ [ l: {t.r3-t.r2_2t
(standard deviation)
I
I
I
-280-
I
I
-------------------
I
I
N
I
CX) I-'
t--0
I
II
t-,.0
<]
I
I
I
I
I
I
2
~--~---n II
I
I
h.i T L---R~f-----
I~ Ti
I
II- T lI FLIGHT TIME
I
I
I
I
I
T MEAN---a► Figure 3. AVERAGE RATE METHOD FOR t6.'t'
I
I
WESTWARD
I
~~ @~ I
I
-200 -100 0 100 200 300 400
I
L\ 'T(nsec)
I
I
EASTWARD
@~~
I
-200 -100 0 100 200 300 400
I
I
L\1'(nsec)
I
I
• I •I I EAST I
I WEST
I I
-200 -100
0 100 200 300 400
I L\4J'(nsec) KEY:
THEORY
I
EXP
I
I
Figure 4. TIME GAINED BY FLYING CLOCKS
I
-282-
I
I
I
I
8. 0 SUMMARY AND CONCLUSIONS
I
Portable cesium beam clocks (model 5061A) cannot be expected to
perform as well under traveling conditions as they do in the laboratory.
I
Our results show that quasipermanent rate changes as large as 5 nsec/hr
(120 nsec/day) may occur during trips with clocks that have shown consid-
I
erably better performance in the laboratory. Of course, such changes
reduce the utility of these clocks. For example, if a flying clock changes
I
rate by 5 nsec/hr shortly after the beginning of a two-week trip, and no
other significant changes occur, synchronizations with this clock shortly
I
before the end of the trip would be off by 1. 6µ sec. However, our results
also suggest that the average of four flying clocks permits synchronization
I
with an uncertainty of less than 1 nsec/hr (24 nsec/day), assuming no
intercomparison data are recorded. With intercomparison data, it should be
I
possible to reduce the uncertainty even further.
Although the final analysis of our data is not yet completed, we
I
have established, with an intermediate level of analysis, that portable
cesium beam clocks are capable of showing relativistic effects with rela-
I
tively inexpensive commercial jet flights. The results of this analysis are
in reasonable agreement with theoretical predictions. However, those who
I
doubt the validity of conventional relativity theory, and there are many
people in this category, probably will not be converted by the results shown
I
in Figure 4. Indeed, the difference between theory and measurement in
Figure 4 is disturbing, and if our final analysis does not improve agreement,
I
an improved version of this experiment should be given serious considera-
tion. The standard deviation on the measurement could be reduced consid-
I
erably, probably by a factor of ten, with such improvements as the use of
dual beam clocks and circumnavigations with less ground time. In any
I
event, this experiment verifies unequivocally the existence of the predicted
east-west directional asymmetry; only more precise magnitudes remain to be
I
established.
I -283-
I
I
DISCUSSION
I
I
MR. CHI: Is your predicted time drift, or time change, nanoseconds per
day, the initial and final? Are they about the same length of time to determine those values?
I
DR. HAFELE: What I did was to go from the time, either starting or stopping of the flight, to the first rate change that was obvious in the laboratory.
I
Now all of these included at least twenty-four hours. Many of the clocks
changed either a day before the flight (there was a rate change in the laboratory} or a day after the flight (there was a rate change in the labora-
I
tory} and that's why, particularly in the case of 361 after the eastbound
flight, it is quite uncertain what the rate is after the flight. Hopefully,
I
the intercomparison data will clear this question up a lot. What we'll have
is a mean flying rate before the flight and a mean flying rate after the flight, which should be equal; if these four clocks are working randomly, if these
I
effects are random, they will scatter about a constant mean. However,
there will be an offset in the mean due to the relativistic effect. The theory of relativity cannot induce any permanent rate changes; it's only while the
I
clock is moving that these changes occur.
DR. ALLEY: I think that Professor Hafele and Mr. Kea ting and Dr. Winkler
I
and others are to be commended for carrying out, for the first time, an experi-
ment in which a recording clock has actually been carried in a closed path in space-time and returned. I would, however, caution members of the
I
audience, or remind them, that specia 1 relativity has been verified in many,
many experiments, as has the frequency change of transmitters, like gamma ray emitting atoms and so on. I think the primary significance of this lies
I
in the ,first time of returning a clock and comparing actual elapsed time.
DR. WINKLER: I would also like to make a comment here and attempt a
I
clarification of the additional data processing which is still to be done, and
why it hasn't been done yet. I feel very guilty about the fact that we haven't been able to provide better support to this case. What has happened is that
I
the data set is available in digital form, but there are a couple of errors
which are quite evident when you look at the chart. Charts are automatically produced from the data set, and Professor Hafele's assistant, Mr. Keating,
I
had to be sent on temporary analysis duty immediately after their return.
Therefore we felt it was wise to let the same people complete the analyses who have collected the data.
I
I
-284-
I
I
I
I
The analyses will be completed when the remaining 9 0 percent of
I
the information is processed. I think what has been represented is not
"'"I'"'.
more than 10 percent of what is actually available by way of measurements.
The basic idea is to keep account, by using the internal measurements,
I
of the measured frequency variations of each standard and to be able to connect the average rate across a flight duration to what you find after
you return. If the rate that you predict is altered by all corrections follow-
I
ing from the hourly measurements made during the flight, and that result is very close to measurements after return, I think you will have a very
great confidence in the result. So I am confident that data processing
I
will produce a much more objective and much greater confidence in the final result, because we use an automatic procedure in which judgments
are left to the computer.
I
DR. HAFELE: I wonder if I could respond to Professor Alley's comment.
He said that the special theory had been thoroughly proved by all ki.nds
I
of experiments. Well, I think that in the same respect there's never been an experiment done by anybody on either the special or the general theory
of relativity which disproves either one. The general theory just makes
I
some interesting predictions that you can't test. Does a clock on the ceiling run slower than a clock on the floor? We don't know for sure, but
it looks as though when you send gamma rays up from a radioactive nucleus,
I
they are absorbed only if you doppler shift the upper nucleus. Does that prove that a clock on the ceiling runs slower than a clock on the floor?
Many people will say "yes, it has to, and there's no point in doing the
I
experiment." But then there are a lot of people who don't buy that argument. So the special theory has. been tested in the same way that the general
theory has been tested so far.
I
DR. ALLEY: I just want to make clear that I do regard such experiments
as worthwhile. In fact, there are people who have questioned the·
I
assumption that clocks actually run at different rates in different gravitational potentials. But that is an integral part of the structure of general relativity
which has been tested to a certain extent. However, convincing demonstra-
I
tions of this, I feel, are very much in order.
MR. GATTERER: For what it's worth department, I would not accuse Dr. Bonanomy
I
of saying that this observation is worth very much. He may not even wish to be quoted in this gathering, but he's had recent experience in transporting
portable clocks around in cars, and it is his opinion that, or at least it is his
I
suspicion that, the environment in an aircraft is a little more devastating to a cesium standard than it is in a car. His experience indicates that his
closing error is smaller than, for example, the Naval Observatory's portable
I
clock trips. It is entirely possible that the environment that may be to blame is a relativistic change, although I think he thinks that the loading and the
I -285-
I
I
unloading process is damaging; therefore your data are varied, perhaps,
I
in knocks and bumps. DR. WINKLER: If you look at the differential phase plots available on the
I
little paper and follow these differential phase plots which have been using the data throughout, the laboratory stay, the westward trip, the eastward
I
trip and so on, you cannot discern any substantial change in the appearance
of these phase plots. There are clock pairs--incidentally not in the same package--which have an extremely uniform performance, and there are
I
instances in which you can identify very well the moment and the magnitude
of the frequency change. I don't know if we have labeled them, but cesium 12 0 has suffered on the westward trip, a very definite and very accurately
I
determined frequency shift. You question the effect of the environment.
Well, that has been the main purpose of conducting the experiment; that is why we got the eight thousand dollars from ONR, to study that effect.
I
I think we have plenty of information on it, some of it can be inspected
visually, and I meant you to do that. I think the aircraft environment is somewhat better than in a car, it depends whose car. I think that is a very
I
wise comment because Dr. Bonanomy's car may have been particularly
bumpy.
I
DR. RUEGER: Were you instrumented in the three dimensions of the accelera-
tional field during the flights?
I
DR. HAFELE: No, we don't know what the accelerational fields were, except
for those that we felt. The plane jiggled a little bit between Athens and Rome, and there was some rough weather. When we were sitting in the seat
I
we felt about as heavy as we always feel.
DR. RUEGER: But the circling of an aircraft over an airport would be an
I
asymmetrical field you'd get on this instrument.
DR. HAFELE: Very slight. You mean that the force down would not be
I
balanced by upward force; you mean that there would be a sidewise force on
the clock. Very slight, I suppose that, yes, if you turned the clocks over it changes them. Well, I'm not the expert on how these clocks perform
I
under various environmental conditions.
DR. RUEGER: Well, we knew that this was a very important parameter in
I
the carrying of crystal clocks, and we've always been very careful to block
them up square and do other things to keep them from being the least bit out of square.
I
I
I
-286-
I
I
c::
I
DR. HAq;LE: I see what you mean. Yes, the clocks sat in the seat and
.....,
I
they weren't square in the seat, they were tilting back a bit in the seat, so they weren't sitting level in the plane. It was almost impossible to
do that; we weren't rigged to do that. We just plopped them down on the
I
seat and strapped them in.
MR. COCHRAN: Has anyone checked the gravitational effect of a high-
I
altitude ground level against sea level? Or would that be within the accuracy of the clock?
I
DR. HAFELE: You can see from the plot (Figure 1) the difference between altitude zero and altitude 10 kilometers. With a high altitude balloon you
I
can easily go to 100 kilometers. Nobody's ever done such an experiment. DR. RUEGER: There was, I think, an experiment done with crystal clocks
I
in eccentric orbits. If you massage the data rather exhaustively, they could show there was indeed an influence on what we call the gravitational
redshift which at the surface of the earth is in the order of two parts in
I
ten of the thirteen per kilometer. This is an effect of the difference of gravitational potential.
I
DR. HAFELE: You"re talking about frequency now, though; I'm talking about time. That experiment was done with frequency on a crystal oscillator
in a satellite.
I
DR. RUEGER: Frequency is synonymous with time for most people. Anyway,
the experiment did show the periodic effect with the period of the orbit and
I
I think there was a confirmation in that of, perhaps, 20 percent. Is that reasonable?
I
DR. HAFELE: I looked atthose data and you believe what you wanted in there. If you don't want to believe it you don't need to. If you want to, you see it
there.
I
DR. WINKLER: I would like to add an entirely different comment here. The
experiments may have been more practical and may have more practical
I
relevance than appears at this stage. Carrying clocks around the earth in a non-inertial frame of reference is essentially the same thing as sending
an electromagnetic wave around the earth one way and the other way. If
I
you remember the discussions about utilizing the Defense Communications Satellite System for high precision synchronization, if we have a network
which depends on an internal high precision synchronization to within a
I
I -287-
I
I
few tenths of nanoseconds, then we will have to expect that the network
I
will have some border difficulties. This depends on which way around
we have propagated the precision timing, and this is a test which we have considered doing. We have discussed it with some people in DCA;
I
there are some difficulties of an entirely practical nature. On the otrer
hand, it has practical implications for the design of the systems and the operational routines to be employed.
I
DR. HAFELE: The magnitude of that effect is, I believe, about 250 nanoseconds difference in time. It's a big effect--a quarter of a microsecond.
I
MR. BARTKO: Why don't you take the predicted offset and get the clock out of synchronization and then try to bring it back in synchronization by
I
going around the world, reversing the experiment so to speak?
DR. HAFELE: I'm not sure I understand. You want to predict where the
I
plane's going to go so you can make a prediction of what the relativistic
effect will be, and then you change the setting of the clock before you 1a ke off? And then come around and see if it reads zero? What is gained
I
by that? MR. BARTKO: Well, that's the reverse of what you did originally.
I
DR. HAFELE: I must admit I don't see the gain. I suppose that changing something in one of these clocks is one of the worst things you· can do;
I
just reaching in there and turning a knob, it takes it two weeks to settle down.
I
DR. KLEPCZYNSKI: I don't believe you could make a prediction in advance as to what the retardation or advance would be, because you wouldn't
I
know what altitude or velocity you 're flying at.
DR. HAFELE: That would be a problem.
I
I
I
I
I
-288-
I
I