That was both educational and enthralling, Mark!
Crystal oscillators are something we treat as a "black box" and just hope they work. Sometimes that's good enough, sometimes not.
I've worked on a design for a clock that needed to keep time within 1 second per day at 20C and within 4 seconds per day between -30C and 65C. That's 11 ppm and 44 ppm respectively.
Crystals are typically specified with a manufacturing tolerance of 30-50ppm and a temperature sensitivity depending on the specific angle of cut. The circuit board, microprocessor and capacitors are going to have manufacturing tolerances as well, adding up to the same order of variation. The crystals also change their specs a little as they age.
So on the production line we had an Agilent Function Generator calibrated to 1 ppm and periodically sent off for certification. It generated a 1Hz reference signal that was fed into a timer pin on the micro. It measured the period of that signal to better than 0.5ppm, calculated a correction factor and stored that in its EEPROM. About once per second it would add that correction factor to the calculated time, as well as adding in a cubic correction factor derived from an on-board temperature sensor.
If you find you really need to know how these things work, start here:
http://www.ieee-uffc.org/frequency_control/teaching.asp
Especially, starting at page 664:
http://www.ieee-uffc.org/frequency_control/teaching/pdf/fcdevices.pdf
Tom