Seeing the Light

These days, amateur astronomers are far better equipped — at least electronically — than they were just a decade ago. Many of them have acquired CCD cameras at modest expense and have learned well how to use them to estimate the brightness of variable stars, a process called CCD photometry. This makes for more precise measurements than the “old way,” which was to look at the star, possibly assisted by a telescope or binoculars, but with no other instrument to estimate brightness than the human eye. The target star is compared to nearby stars of known brightness (comparison stars) and a “judgement call” is made. This process is called visual photometry.


Perhaps you’re thinking “At last! With this innovative new technology, amateur astronomers can finally produce data which is of value comparable to that obtained by professionals!”

You’d be wrong. Not about amateurs making observations of comparable value to professionals, they can. But you’d be mistaken to think that they couldn’t do so before.

How can that be? How can “eyeballing” a star and making a “guesstimate” of its brightness compare to the precision of professional-grade high-tech photometry by professional astronomers at professional observatories? After all, visual observers only report their estimates to the nearest tenth of a magnitude, while professionals report to the nearest thousandth of a magnitude (a “millimagnitude”, abbreviated “mmag”). The fundamental resolution of professional photometry is a hundred times finer than visual photometry. How can visual observers compete with that?

Let’s take a close look at an interesting star which has been observed by amateurs and professionals alike, a Cepheid-type variable known as X Cygni, or just “X Cyg” for short.

I managed to find 385 brightness estimates of this star by professional astronomers, published in the scientific literature, which are available online through the McMaster Cepheid Photometry and Radial Velocity Data Archive. Here’s a graph of the brightness over time:

It doesn’t look very informative. That’s because Cepheid-type variable change brightness so quickly that the pattern of fluctuations isn’t easily visible in this graph. This particular star shows a periodic oscillation with period 16.386 days, and since the graph covers over 7,000 days the cycles are squeezed too thin to be visible. Besides, professionals can’t monitor individual stars all the time, so for most of the 7,000 days shown on the graph, they just weren’t looking at it.

When patterns aren’t readily visible, it’s time for analysis! A natural choice for things that oscillate periodically is Fourier analysis. It produces what’s called a spectrum for the data, and here’s the particular version called the amplitude spectrum:

For each frequency of oscillation, it computes the size (amplitude) of oscillation which best matches the data. The strongest response is at frequency f=0.06103 cycles per day, which corresponds ot a period 1/f=16.386 days.

There’s also response at most other frequencies, but that turns out to be due to the main period. Its influence leaks to other frequencies, a phenomenon known as aliasing which we examined in another post. There are some fancy-shmancy ways to eliminate most of the aliasing in order to “clean up” the spectrum, one of which produces this result:

That’s quite a bit clearer! Now there are only 8 frequencies at which we note strong response, and they’re all multiples of the fundamental oscillation frequency of this star. They’re the “harmonics” of its oscillations, and they combine to make the particular form of its cycles.

It’s also worth noting that in between those genuine oscillation frequencies, the level of the spectrum is extremely low. The average value is only about 5 mmag, so if there were another oscillation even a little stronger than that we’d be able to detect it. That’s how good the professional data is!

And what is the form of its cycles? The easiest way to see it is to fold the light curve. Essentially, we wrap time around in circles, with one circuit exactly equal to the period of the oscillation. Doing so should make each moment of maximum brightness align with all the other maxima, and each minimum align with the other minima, and give us a good idea of the shape, or waveform, of the cycles of this star.

That enables us to plot the folded light curve. It’s customary to plot two folded light curves laid end-to-end. When we do so we get this:

Not only can we clearly and easily see the shape of the oscillation (the waveform), we can even note some interesting fine details, like the “bump” in the light curve at about phase -0.2 cycles. Such bumps are seen in the light curves of many Cepheid-type variables. Note that the bump itself is considerably smaller than a tenth of a magnitude, how could visual photometry could reveal that? Professional photometry has shown us the shape of this star’s variations in fine detail, and nailed down the period and amplitude of this star with precision.

Take that, you old-style amateurs with your visual photometry!

What about that old-style amateur visual photometry? It’s certainly not as precise as professional work. But when I scan the database of the AAVSO for X Cyg, I find a lot more observations by amateurs doing visual photometry than the professionals had. In fact, there are 22,745 visual observations of this one star.

I took the data from a single one of those amateur observers. His observer ID code is “LX” and he was Wayne Lowder (may he rest in peace), formerly president of the AAVSO, a very skilled observer, and a genuinely wonderful person. He alone did 4,852 visual observations of X Cyg, over 12 times as many as all the professional data I was able to find combined. Here’s his raw data:

It doesn’t look like much either. But again, the cycles are squeezed together, and there’s so much data that different points are plotted on top of each other. Time for analysis!

Let’s do Fourier analysis again. It gives this:

This amplitude spectrum is better than what we had with the published photometry. That’s because the professionals observed so intermittently, there were so many large gaps between their observing runs that aliasing was a worse problem for their data.

But this amplitude spectrum is not as good as the “cleaned” version from professionals. So, let’s apply the same cleaning technique to the amateur data:

Well well! The cleaned spectrum is almost as good as that from published photometry after all. There are 8 frequencies at which we have significant response, and they’re all multiples of the fundamental oscillation frequency of this star. We can also expand the y-axis and superimpose both spectra, leaving out the signal frequencies so we can focus on the background “noise” level, with that from the published photometry in blue and from visual photometry in red, to see how they compare:

The background (“noise”) level is only about 5 mmag for published photometry. But it’s a bit less than 10 mmag for visual photometry. That means we could detect oscillations with amplitudes only a little bigger than that. Which is impressive, since it’s smaller than the resolution of visual photometry. It turns out that we can detect features smaller than the data resolution because of the statistical power of large numbers. Amateurs contribute so much data, that its statistical power is magnified greatly.

What about the shape of the cycles? We can make a folded light curve with the visual data, and I’ll superimpose (in blue) the folded light curve we have previously from published photometry:

The visual data shows a lot more variation — a lot more “scatter” — than the published photometry. But there’s a lot more of it. We can take advantage of its plenitude by averaging the data over intervals 0.05 cycles wide. When we do that with both data sets they compare thusly:

In spite of the greater variation in visual data, their greater numbers make the averaged light curve nearly as precise as that from professionals. We can even see fine details which are smaller than a tenth of a magnitude, like that “bump” on the light curve at phase around -0.2 cycles. There’s so much visual photometry we can even average it over smaller intervals, only 1/50th of a cycle:

There is indeed enough information in the visual data (because there’s so much of it) to see the detailed structure of the waveform.

We do notice that there’s an offset between the two light curves. Wayne Lowder’s data is consistently about 0.17 magnitudes fainter than the published photometry. If we offset his data by this amount we can see just how well the two different sources of data agree:

Yes, Virginia, visual data really can compete with professional photometry. But to do so requires large quantities of data. Experience and skill on the part of the observer doesn’t hurt either! It’s a lot of effort to accumulate nearly 5,000 visual observations of just one star (and Wayne Lowder observed a lot more than just this one star!), but its scientific value is tremendous. The value of the data set taken as a whole is comparable to that from professionals.

Amateur astronomy, whether by high-tech or just by eye, is a rewarding task, both for the observer and the analyst. If I have anything to say about the matter, it will not be a thankless task.

About these ads

17 responses to “Seeing the Light

  1. Is there some sort of central repository for this kind of data? If one person’s observations are nearly as good as a combination of professional datasets, I imagine combining everyone’s data would be even better.

    [Response: AAVSO]

  2. Yes, thank you Wayne (and others). It is simply fascinating that this level of accuracy can be obtained by one observer. I’m also intrigued by the consistency of his error as well. Almost makes me wonder if his comparison star’s magnitude was mis-labeled!

  3. Gavin's Pussycat

    Ah, this brings back great memories… I never contributed variable star observations, but I was along in Moonwatch! Endless nights under the stars…

    KenM, I suspect the magnitude systematic offset is simply due to Wayne Lowder’s eyes not being a properly calibrated ‘V’ channel of a UVB photometer, and the colours of X Cygni and his comparison stars being different.

  4. My wife bought me some binoculars for my birthday back in March, and a month or so later my mother-in-law sent me a small telescope that was going unused at her place, and as a result I’ve been beginning to get into the amateur astronomy thing. I joined a local club and have been to a couple of talks and observation events. So I’m quite new to it all. But it astounds me already how much I notice that I wouldn’t have looked for before, what I can see that I wouldn’t have seen before. Getting some training in observation really opens up a world of detail. I looked at the tutorial on the AAVSO site, and that looks like a great way to start building up one’s observation chops.

    From this post, and the concept of trained observers generally, it feels like all sorts of comparisons and analogies to expertise in other fields (like climate science) are begging to be made. Which could apply to further discussions about the ease of misleading people who haven’t developed a feel for the field, if one were so inclined (and I’ve heard some people actually do this–what’s up with that?).

    But I’ve got to go to work, so I don’t have the time to chase it all down. Thanks for this, Tamino!

  5. KenM, I’m guessing but it’s possible that the consistent offset is a result of different spectra from the target star and the reference stars causing different responses in the professional telescope and the visual observer’s eye.

    • Gavin's pussycat

      Ed, that was my idea too. And looking at the comparison curve: is it just me, or does the colour of this star change at brightness maximum and minimum?

  6. Thanks Tamino, that was a fascinating post. Very eye-opening, too, in terms of how it’s still possible to do useful science with relatively cheap & relatively primitive instruments.

    Then again, it’s how *all* science was done before the invention of the modern, expensive instruments.

    • I forgot to add – I observed the Venus Transit yesterday with a pinhole camera arrangement – a cardboard tube ~1200mm long, foil with a pinhole on one end, a piece of white paper on the other, and a viewing hole cut in the side. Worked brilliantly, and while it would have been nice to have some magnification, I’d say that some useful science could have been done w.r.t. timing of first / last contact. The photos I took with the point&shoot camera I had with me looked pretty ordinary, especially as everything was hand-held, but then again – holy crap, I took photos of a transit of Venus with an instrument I built myself in 10 minutes! :-)

  7. Tamino, do you have any recommendations for tools for advanced amateurs for doing scientific work? I learned how to use IRAF to do crowded field photometry of globular clusters around 1996 or so in undergrad, but haven’t touched it since then. Is that still the gold standard for doing scientific analysis?

    [Response: Photometry isn’t my bag (just data analysis). But I do know that AAVSO has a lot of photometric observers, and the current director (Arne Henden) is one of the most skilled around. I’d suggest searching the AAVSO website (http://www.aavso.org/), and if you don’t find what you’re looking for you could email them. Their reason for existing is to help folks who want to do astronomy.]

  8. But what if some of the more dedicated amateur astronomers are aware of the general behavior of cepheid type stars. This could aid in their visual photometry but possibly lead to systematic biases.

    I’ve always found that the best way to eliminate systematic biases in observations is to to have 2-3 stiff drinks before making any observations :)

  9. Tamino, on this subject, you may be interested in my most recent post on observations of the Venus transit. (Some broader climate discussion included)

    http://climatephys.org/2012/06/07/the-venus-transit/

  10. Two things stand out in my educational career that changed my world view. Oh, sure, I read the classics, I studied music and art. I write about design, and I’m published, and with a new book on the way. But you know what? Calculus and Statistics shed more light on the real world than any other subject could ever do. These two subjects reveal a fundamental wondrous fact about the universe that we find ourselves in: it’s mathematically describable. This realization was a huge leap for me that continues to inform my thinking, now 40 years later.

    My statistics teacher had the same ability to logically present data: complex data presented in a way that is easy to understand.

    So thanks for these posts. They bring back great memories! I do wish that some of the underlying analyses were included (formulas). But it’s a nit. These essays are terrific examples of how we bring understanding to our real world.

  11. Whether intentional or not, this is a good analogy to how one can estimate temperature trends to a tenth of a degree per century from source data having one-degree resolution and even worse accuracy.

  12. Does the X Cyg light curve change any over time?

    in particular, did it change over the period in which Wayne Lowder observed it?

    [Response: Not that I can tell. Cepheid-type variables tend to have remarkably stable oscillations. But some of them do show period changes, and a *very* few show amplitude changes (the most famous case being Polaris).]

  13. Indeed I did, and thanks.