EAON circular # 5

The timing of visual observations :

1/ Personal equation and accuracy

2/ Timing with stopwatch

3/ Timing with tape recorder

4/ Other methods

5/ Notes on time-keeping sources.

If you want to make an visual observation of any scientific interest, it is necessary to subtract your reaction time (personal equation = P.E.) from your timings and to give the estimated accuracy of your measures. Very often, observers underestimate their P.E. and the uncertainty of their results. Here are some cues to estimate them, at least crudely.

PART ONE

Personal equation, and accuracy

Everyone has his own P.E. : at each timing, you must subtract it to find the "true" time. But the P.E. varies a little at every timing : this fluctuation generates an uncertainty in the "true" time and decreases the accuracy of the timing.

Here is a simple method to measure the P.E. with a digital stopwatch. With a piece of paper, hide the digits of seconds : SS.ss, and start the stopwatch. As soon as you see the minute digit changing, stop your stopwatch and look at the hidden digits : you have a first measurement of your P.E. Repeat the measure a number of times to find the most probable value and its usual variation.

When doing your experiments, be careful to isolate yourself from noise, and especially rhythmic noises, since one can unconsciously take up and follow the beat.

As an example, I made about 30 tests of this method (on 2001 May 04 from UT 09 35 to 10 15), P.E. measured in 1/100th s :

40 37 46 43 43 43 37 56 40 37

46 40 37 40 40 37 144 40 68 40

34 59 40 43 40 53 34 43 34 31 31

At the 17th measurement, I was inattentive and forgot to scrutinize the display! Such a misfortune sometimes happens when watching an occultation, too ! Could 1.5s as P.E. be a reasonable value in such circumstances ???

So, I discarded this value and did one more test, which one can't do when monitoring a real event !

Well, from these values, it is obvious that this stopwatch is running by 1/30th of second ; we find :

31(x2) 34(x3) 37(x5) 40(x8) 43(x5) 46(x2) ...All these values are spaced by about 0.03s . Test your own stopwatch, one never knows ! You have to take it in account for any timing made with such a peculiar stopwatch.

Without any computation, you will observe that all these "correct" values span from 0.31s to 0.68s, so one could crudely estimate that the P.E. is : 0.5s ± 0.2s .

On the other hand, if you observe that the most frequent value is 0.4s, it could be better to estimate the P.E. to 0.4s, with an accuracy of +0.3/-0.1s . Then when you get a timing t , the "true" time is t - 0.4s, -0.3/ +0.1s (= between t - 0.7s and t - 0.3s).

If you prefer to make some statistic computation, the mean of the 30 "correct" measurements is 0.42s (rounded from 0.417s), the estimated standard deviation is 0.08s (rounded from 0.081s).

[more advanced : it is better to take the estimated standard deviation, denoted s or s(n-1) in hand calculators. With uncertainty in the range of 0.1/0.2s, it makes no sense to keep all the digits at the end of computation]

Thus the error bar at 95% confidence is ± (0.08x2) = ± 0.16s = ± 0.2s (rounded). One can estimate that the P.E. is 0.4s, with an accuracy of ± 0.2s .

[more advanced : however it is not valid to use a gaussian law (or rather here a Student's law) for a random variable which is not at all gaussian. The reaction time law is bounded and strongly skewed, because nerve impulse travels at a finite speed of a few m/s between brain and finger. Furthermore, the distribution of reaction time could be bi-modal, with values clustered around two main ones, a fast one and a slower other one. And it depends if you are tense or inattentive, or tired ... Hence a symetric error bar is almost certainly wrong]

How many tests will you have to do ? Well, more is better : 100 is better than 30, and 200 better than 100. You can achieve this in several sessions. Don't try to take a value after every 10 seconds : this lap may be too short, and you could unconsciously anticipate the change in the digit of tens (I tried it, and found too small values !)...

[more advanced : draw a graph of the values. If there's a hint of multiple clustering : i/ make more tests, ii/ you could try to fit a log-gaussian law to each cluster]

Remember that there are causes of unaccuracy other than the intrinsic fluctuation of P.E. . For example, if the stopwatch rounds the values to multiples of .03s it introduces a further unaccuracy by at least .015s

[more advanced : according the way it rounds, to the nearest value, or to the smaller, or greater one. It still is NOT a gaussian variable, but a uniform one with constant probability].

You can estimate crudely the total uncertainty as the sum of the spans :

± 0.16 ± 0.015 = ± 0.175s, rounded to 0.2s (three digits would be delusive).

[more advanced : IF the random variables are independent (it is true here), the whole standard deviation (sd) is the quadratic sum of the various sds = the square root of the sum of the squared sds. Note that here the sd of the rounding uniform variable is 0.015 s / SQRT(3). Thus the total sd is SQRT( 0.006593 + 0.000075) = 0.081656s , hence a gaussian error bar ± 0.1633s... which is not valid since the random variables are NOT gaussian...]

Are there other methods to measure your P.E. ?

Yes, but it requires more hardware. The best one could use an artificial star seen through an eyepiece, lighting on and off at random times, while it is properly recorded (for example with a chart recorder) with your pressing on a button.

According to your time-keeping source(*) and your method of recording, you can have more uncertainties to take into account.

All these uncertainties increase the error bar !

You might also be surprised by an unexpected phenomenon (blinks, occultation of an unknown double star, etc..), and have a supplementary time lag : do not underestimate the resulting uncertainty!

[more advanced : For the variables that are gaussian, it is right to compute the quadratic sum of their standard deviations, and then the 95% confidence error bar. For the others variables, it is conservative to add the half-span to the previous 95% (or 0.975 single-sided probability) error bar].

-------------------------------------------------------------------------------------------------------------------------------------------------------------

(*) 2 short stories, both outside Europe : 15 years ago, I discovered that a local phone time-signal was exactly 1 minute fast. A decade ago in another country, two observers synchronized their stopwatch from local radio braodcasting station. The observations were positive, but timed 15 seconds late ...

PART TWO

A positive asteroidal occultation observation requires TWO absolute timings :

- the disappearance (or immersion)

- the reappearance ( or emersion ).

It is necessary to link what you see to a reliable time-keeping source. You generally need a medium (a temporary time-keeper, very often a stopwatch or a tape recorder) to make this link, so there are several sources of possible errors :

- your reaction time (personal equation),

- the unaccuracy of the medium,

- the comparison of the medium with the time source,

- the uncertainty of the time-keeping source itself.

Every point needs attention and testing.

Timing with stopwatches

First of all, verify how reliable are the stopwatches you use !

How to start or stop a stopwatch at a precise time ?

Listen to the bips of a time source, take up the rhythm (3 or 4 beats are enough) and press the button synchronously with the following beat (just like a musician in a band).

How to measure the drift of a stopwatch ?

Start the stopwatch at a precise time and let it run. Hours later, stop the running stopwatch in the same way, and compare its measured time with the precise time elapsed. This will give you a correction factor if you need it.

Example : suppose you start the stopwatch at UT 19 36 04.00 and stop it one day later at UT 19 36 04.00 , and then the stopwatch says it measured 24h00m01.55s. Then your correction factor is 86400/86401.55 = 0.999982 . It means almost nothing if you run such a stopwatch 10 minutes, but it must be taken into account if it runs 1 hour or more.

Warning ! Recent experiments (march and april 2002) show that some digital stopwatches can be fast as much as 0.58 second per an hour, and even more ! Repeat your experiments to be sure that such a rate is constant and not varying at random.

Also, there could be a time delay when you press the "start button" ( as much as 0.16s in one experiment).

[More advanced : You still can use the stopwatch if it varies at constant rate but you have to take it into account]

How precise is your synchronous starting or stopping of the stopwatch ?

Test it in the same way, by measuring (many times) precise short rounds (of 1 minute for instance). You should find that a measured minute is between 59.95s and 1m00.05s .

[More advanced : Your mean error is the standard deviation (sd) of the differences between precise and measured durations. Note that the mean of these differences should be zero : it is approximately OK if the absolute value of this mean is less than 2 sd divided by the square root of the number of measurements (the so-called standard error, for it is here a gaussian variate).]

Is there a time delay in the display ?

Experiments show that digital LCD very often have a time delay (increasing with cold). To test it, hide the digits of tenths and hundredths of seconds with a piece of paper. Start the stopwatch in rhythm with a reliable time-keeping source. Then (time source turned off not to be disturbed) look at the display, take up rhythm of the seconds digit, and stop the watch. The interval measured should be an integer number of seconds : generally it is not. Take note of the difference (only fractions of a second, generally positive, maybe negative sometimes). Repeat this test many times (for example 50 times). Then compute the mean and standard deviation of these numbers. The mean is an estimate of the time delay of your display. The standard error of this estimation is equal to the standard deviation divided by the square root of the number of experiments (~7 if you made 50 tests). If not negligible, take it into account for the estimation of your P.E.

EXAMPLE : with the stopwatch used in PART ONE, I found a mean time delay of 0.05s , the standard deviation is 0.035s, the standard error is (for 50 tests) 0.035/7 = 0.005s . This stopwatch has a time delay of 0.05s ± 0.01 at 2 sd level. Thus my P.E. is 0.40 - 0.05 = 0.35s

There are TWO kinds of stopwatches :

- simple run-and-stop stopwatch

- stopwatch with memories (intermediate or lap timing).

One timing with one stopwatch (without memory)

NOTE : (one timing is enough only for regular -no grazing - lunar occultation)

There are TWO methods to proceed :

- start the stopwatch before and stop it when you see the phenomenon.

- start the stopwatch when you see the phenomenon and stop it afterwards, when comparing with the time-keeping source.

The risk with the first method is to stop accidentally the running clock BEFORE the observation. Then reverse to second method ! Make a subtraction instead of an addition ...

The error comes from the variability of your personal equation (maybe greater if the phenomenon takes you by surprise) and the uncertainty of your comparison with precise time-source which is (in 95% of occurences) 2 times your mean error (see above).

If the stopwatch has memories, use them to repeat several times the comparison with precise time one or a few minutes later. Then compute the mean of timings (you can't do that without memories : here is a drawback of simple stopwatches).

[More advanced : you can compound quadratically these uncertainties since they are stochastically independent. Take the square root of the sum of uncertainties squared : if P.E. is +0.25 / -0.15s and rhythmic timing is ± 0.07s, you'll get :

SQRT[ 0.25² + 0.07² ] =0.2596 and SQRT[ 0.15² + 0.07² ] = 0.1655 hence an error bar +0.26 / -0.17s ]

Two timings :

You need TWO simple stopwatches (and it can be confusing if you're timing a very short occultation : here is another drawback) or use a spared memory to take the reappearance in intermediate timing. Reduce each timing separately, as indicated above.

Don't forget that you may estimate differently your P.E. at disappearnce (you might be surprised) and at reappearance (you might be tense and attentive, since reappearance comes soon after disappearance). Some observers prefer to time the disappearance with one stopwatch, and the duration of occultation with the other. The drawback is that both timings could suffer from the "surprise lag" effect.

stopwatch with memories

Only one stopwatch is enough to carry all the work to be done. Start the stopwatch BEFORE the beginning of observation according to your time-keeping source (radio signal or reliable phone time service) and do all the timings with the "intermediate time" button. Spare some memories to compare again your stopwatch with the time-keeping source AFTER the observation is done.

NOTE : If you have only one memory, time the disappearance with it and the reappearance by stopping the stopwatch. Or start the stopwatch at immersion, take intermediate time at emersion, and compare with time source when stopping.

[More advanced : If you took some intermediate times in comparison with the time-keeping source before and after the observation, you can compute a regression line and have a better accuracy for the event timings.The best accuracy is reached when the observation time is at midtime between the start and the end of the stopwatch run. The timing of events are isolated estimates (not averaged ones) from the regression line : the uncertainty is more important than for the estimation of a mean ]

¤ --- ¤ --- ¤

PART THREE

Timing with tape recorder

In this standard method, one marks vocally the event while simultaneously recording time signals. Radio signals like DCF or WWV give a bip every second.

For the vocal mark, it is the consonnant "T" or "P" and "I" or "E" vowel sound that are the shortest : so "TIP" or "PIP" give the most accurate vocal mark. Of course the timing to consider is the beginning of the mark.

If you have no adequate radio receiver, you can use a stopwatch (in the "countdown" mode) which sounds a bip every 5 or 10 seconds and start it according to a reliable time signal (like the hourly signals of some broadcasting stations, or wired phone time service. DO NOT trust mobil phone or FM band stations before testing them : sometimes the signals could come from satellite relay with unknown time delay).

Remember that a tape recorder is NOT a clock. Recording ONE time signal on the tape before the observation is a bad way of timing an event. An observer struck by a sudden stopwatch failure tried that "technique" to save his observation of the 1998.03.21 occultation of HIP 28954 by 39 Laetitia. While listening to the tape many times he found that the timing was varying from "start + 20 minutes" to "start + 22 minutes" : no absolute timing was possible.

Here are 24 successive timings (in april 2001, temperature +24°C) of the same 60 seconds interval (recorded in march 2001, temperature +6°C) :

58.59 58.72 58.03 58.87 58.81 58.81 59.21 59.47

59.31 59.52 59.21 57.95 59.03 57.92 59.25 59.04

59.00 59.11 59.03 59.07 58.39 58.91 59.00 58.65

Never trust a tape recorder as a clock !

After the recording is made, one must read it up to time the events. Listen to the tape and count the recorded bips to know which one is just before and which one is just after every mark of events. Thus you know (in the favourable case) the second of the event (don't forget to subtract your P.E.). Now you must refine :

Use a stopwatch to time some of the recorded bips AND the mark of the event. Start the stopwatch as usual by taking up the rhythm of the recorded bips (not possible if one bip every 10 seconds), and stop it at the mark of event. THEN, don't forget that the stopwatch stops NOT at the time of the mark BUT AFTER YOUR TIME DELAY ! To retrieve the correct time of event, you have to subtract TWO times the P.E. from the stopwatch indication : a first time to find the time of the mark, and a second time to find the time of event.

EXAMPLE : I start the stopwatch from the recorded bip corresponding to UT 19h01m04s and stop it when I hear the mark of disappearance. The stopwatch says 20.55s, and I must subtract 0.35s of reaction time (± 0.1s, since I have no surprise effect : I have already listened to this tape, I know that the mark is coming now) : the mark time is thus UT 19h01m04s + 20.55s - 0.35s = 19h01m24.20s ± 0.1s . The event time is 19h01m24.20s - 0.35s = 19h01m23.85s ± 0.14s [0.14=SQRT(0.1²+0.1²)]

If I had been surprised by the disappearance when observing, then my reaction time when recording could have been rather 0.5s (± 0.2s). In this case the disappearance time is UT 19h01m24.20s - 0.50s = 19h01m23.70s . The total uncertainty from P.E. would be SQRT(0.1² + 0.2²) = 0.224s (better rounded to 0.23s) .

It is much better to use a stopwatch with memories to read your tape, so you can time a bip BEFORE the mark, the mark itself, and a bip AFTER the mark. The tape speed is varying a little at every reading (see above !), due to changing motor temperature, rewinding of the tape, etc ... measuring the time between bips of known interval allows to correct the timing of the event mark from this factor.

In the example above, the stopwatch gives 58.52s between the bips of UT19h01m04s and 19h02m04s . So the intervals measured on the tape have to be multiplicated by 60/58.52 = 1.0253 and the mark time is not 20.55s after the 19h01m04s bip, but 20.55*1.0253 = 21.07 s : the mark time is UT19h01m04s + 21.07s - 0.35s (this 0.35 P.E. is live, not recorded) = UT19h01m24.72s , and the disappearance time would be UT19h01m24.72s - 0.5s = UT19h01m24.2s (rounded, ± 0.45s)

(Note : this example comes from a real experiment made today -20020520- with an old cassette record : the 1998.03.21 HIP 28954 occultation by 39 Laetitia when I was NOT surprised at disappearance, thanks to a Martin Federspield update.)

[More advanced : The uncertainty of the P.E. in listening is divided by the square root of the number of readings, BUT NOT the uncertainty of P.E. at the time of event, which is recorded and occured only one time ]

If your tape recorder has two speeds (or more), you can record at high speed and read at low speed. Be sure of monitoring reading times as indicated above : the variation in low speed can be much more important ! Remember that "low" speed is not necessarily "half" speed. The vantage is that the P.E. and its uncertainty when reading is always the same : so it is divided by ~2 with respect to the recording. If you can reach ± 0.1s when reading, it gives you ± 0.05s for the recording (to compound with the ± 0.1s when timing the event).

If you have no stopwatch (try to have one) you might test the "ear method" : if you recorded a bip every second, you know the time of marks at better than 1s . You then can try to estimate if the mark is near the preceding bip, or near the following bip, or in the middle (it it is not coincident with a bip). So you can have a refined time with an additional uncertainty of 0.25s (here, there is no listening P.E. , and no averaging is possible).You also could try to use simultaneously stopwatch and tape timing.

All this can take a lot of time ! It is the price to pay for scientific accuracy ...

Acccuracy

Sources of error are :

1/ Persistency at disappearance (the retina -?- property which allows us to see movies),

2/ Variability of P.E.,

3/ Uncertainty in the signals of the time-keeping source,

4/ Uncertainty in the start (or stop) in rhythm of the stopwatch,

5/ Uncertainty in the rounding of digits of the stopwatch,

6/ Variation of the running stopwatch or of the tape speed.

1/ is unknown, 2/ is the main source of error, 3/ should be negligible, 4/ & 5/ are small, 6/ is to take into account. You can compound all them quadratically and take an error half-bar of at least 2 standard deviations ...

¤ --- ¤ --- ¤

PART FOUR

Other methods

Chart recorder : ticks are made on a paper tape every second (from an astronomical clock) and every time the observer marks an event by pressing on a button. Formerly standard technique in observatories with meridian refractor.

Programmable hand-held calculator : some of them (such as HP 71 and HP 75) can be programmed so that they give a reliable timing and record the time when pressing a key. One must calibrate carefully the device before using it.

Timing from a computer clock : computer internal clocks were never designed to reach the time accuracy of a good wristwatch (see further). DO NOT trust this method.

Radio signal receivers with memories : I know a lucky observer who owns a digital radio controlled stopwatch with a thousand memories ... (don't forget the P.E.)

Eye-ear method : The observer is listening to a time signal (one bip every second) while watching in the scope, and counts every second to know the time of events. Furthermore he must estimate the fraction of a second between the bips and the occurrence of the event. To count the seconds from the bips while watching an unexpected event through a scope may be too demanding from the observer : according to David Dunham (The Astronomical Journal MAY 1990), there is NO known example of a correct timing by this "method". Dr. Eberhard Bredner (personal communication, 2002) replies that he had been aware of a few exceptions. Theoretically, with some training, one should reach an accuracy better than 0.2s (there is NO P.E.) if one could solve the problem of counting the seconds, for example by coupling this method with the tape recording one. The eye-ear estimation must be made "live". If reappearance occurs a few seconds after disappearance, it could be somewhat confusing...

PART FIVE

Taking time from a short wave radio receiver (tuned on WWV from Boulder, Colorado, or other transmitters in Europe) : if the signal is reflected from the ionosphere, the waves time of propagation is not exactly known (generally less than 0.05s).

Taking time from a long wave radio receiver (transmitters in United Kingdom, Germany, Switzerland ...). These receivers can sound a bip every second : they generally have a time delay, tested from 0.022s to 0.062s, which must be taken into account.

Taking time from phone time service : it could be tested against another reliable time source such as a radio time signal receiver. DO NOT use a mobil phone : a simple test shows that it is at least 0.2s late. If it comes from a local phone network, the time signals could come down from satellite relays ...

Taking time from radio controlled clock : BE SURE that the clock is continuously receiving the radio signal (generally it is NOT) and update correctly the display. Check that the internal antenna (a ferrite bar) stands in the right orientation (if not, it could not receive enough strong a signal to update). Most clocks updates only one time a day (a few minutes around 3 a.m.) and drift can accumulate : a badly updated radio "controlled" clock was catched with a display several minutes (!) in error. Remember that metallic and concrete buildings can shield the radio receiver.

Check also how the display is updated : some clocks display the nearest second (uncertainty : ± 0.5s) , some other clocks display a default rounded time (uncertainty from delay 0 to 0.999... s ! ). A time delay in the display can be temperature dependent.

[More advanced : Take notice that such uncertainties are not gaussian ]

Taking time from radio broadcast station : signals are sparse and could be unreliable (see the "short story" at the end of PART ONE).

Taking time from a GPS display : signals received by the machine are excellent, but the display is not reliable. Some models have a time delay in the display, others can sometimes round the display to the next second in advance.

Taking time from GPS satellite signals : one must use a GPS receiver with the 1pps (pulse per second) output. Japanese astronomers (Hiroyuki Geshiro at Oishikohgen Obs. & Shinji Toyomasu at Misato Obs.) developped a "Geshiro clock for GPS". A more simple circuitry can provide short "ticks" recorded on a tape. This could be one of the best methods.

Taking time from a computer : as said above, the internal clock is unreliable. There are methods supposed to keep up the computer clock with a very precise timing through internet. Tests have been made which show an uncertainty sometimes up to 1 second. But the Geneva observatory have consistently excellent results with NTP (and an atomic clock).

Taking time from an atomic clock : of course ... (I saw ONE report of this kind).

For EAON, Raymond DUSSER, 2003.