Got a PlayStation Move and an old vinyl turntable ? You can try to detect Earth's rotation, find geographic north, and determine your latitude.

1. Overview In this project we follow in the footsteps of Léon Foucault who in 1851 produced the first convincing experimental evidence of Earth's diurnal rotation, thus confirming a crucial element of the heliocentric model of the solar system. For lack of a 28 kg pendulum and a 70-meter-high ceiling to suspend it from, we use a table-top apparatus whose main components are a PlayStation Move motion-sensing videogame controller and a 33/45 RPM record player. With an acquisition time of several hours, our device has little practical use, but this is an enlightening physics experiment that anyone can easily replicate. Technical details and source code are provided. Readers who prefer to skip the hardware tinkering phase will find raw datasets as well. According to our review (see Appendix A, Historical context ), all earlier experimental confirmations of Earth's rotation required expensive or dedicated, carefully crafted instruments - including the pendulum, despite its apparent simplicity. With a few improvements (see Section 8, “ Perspectives ”), our approach will allow millions of individuals to confirm Earth's rotation for themselves at negligible cost and effort, without advanced technical skills, using only household items - a democratization of the Copernican revolution made possible by the consumer electronic industry.

2. Disclaimer This project was developed independently and is not endorsed by any third party. Claims of performance contained in this article are derived from our own tests and measurements on specific units only. The PS Move is marketed as a peripheral with Bluetooth connectivity (Bluetooth is an open standard which promotes interoperability between computers and wireless devices), and we believe this experiment uses it exactly as intended (as a motion-sensing peripheral). There is no need to modify the device or circumvent any protection. Still, following these instructions might void your warranty, damage your hardware, or worse. Use at your own risk !

3. Gyroscopes Gyroscopes used to be synonymous with sophisticated navigation systems for aircrafts, nuclear submarines, rockets and ballistic missiles. Figure 1. Mechanical gyroscope for gravitation research, 2004

In recent years, thanks to the industrial-scale production of microelectromechanical systems, they have become ubiquitous. MEMS gyroscopes can now be found in millions of video cameras, GPS receivers, mobile phones and videogame controllers. Due to their principle of operation, gyroscopes measure their rate of rotation, not with respect to the ground, but with respect to an inertial frame of reference aligned with the distant stars. Since the Earth itself rotates with respect to the celestial sphere, a stationary MEMS gyroscope will sense that rotation. For mainstream applications whose purpose is to track the orientation of a vehicle or hand-held device with respect to the ground, Earth's rotation is merely a source of small measurement errors. However, there are industrial applications for devices sensitive enough to detect it. By measuring not just the rate of rotation but also the direction of its axis, one can determine true north (also known as geographic north) and latitude. The gyrocompass has become a crucial instrument in environments where neither the night sky nor Earth's magnetic field can be observed reliably, such as underseas (submarine navigation) and underground (drilling, pipe mapping). Unlike the directional gyroscope, the gyrocompass does not require manual initialization and does not drift. Sophisticated devices are used for such applications. It was unclear, when we started this project, whether inexpensive consumer-grade MEMS gyroscopes could detect Earth's rotation, be it only qualitatively.

4. Back-of-the-envelope calculations



PlayStation Move



PlayStation Move



PlayStation Move



PlayStation Move

Affiliate links According to [IFIXIT], the gyroscope in the PS MOVE is a chip labeled Y5250H. We found this suspicious because the marking is reminiscent of a single-axis yaw gyroscope from ST (presumably a custom variant of their LY5150ALH) whereas the PS Move is known to measure rotations on all three axes. We disassembled our unit and discovered a slightly different board layout, with an additional dual-axis ST PR425A gyroscope where IFIXIT had photographed an unidentifiable shielded device. Figure 2. Yaw and Pitch/Roll gyroscopes in a PS Move (CECH-ZCM1E)



From earlier experiments (see [LINMCTOOL]), we know that the PS Move is a regular Bluetooth-compliant peripheral, and that its designers were kind enough not to add any obfuscation or encryption on top of the standard. As a result, sensor data can be retrieved by any Bluetooth-enabled computer. Data is transmitted wirelessly as Bluetooth HID reports. Gyroscope data transmitted over Bluetooth suggests a sampling rate of about 176 Hz with a full scale range of about ±3000 °/s and a dynamic range of 16 bits, providing a resolution of 0.09 °/s. This is most likely achieved by oversampling since internal hardware details imply that six analog sensors are multiplexed toward two ADCs capable of 12 bits at up to 1 MHz each. Measurements from a stationary gyroscope has a standard deviation of about 0.45 °/s, which translates to a noise density of 0.05 °/s/√Hz, consistent with the state of the art in consumer-grade MEMS gyroscopes. Figure 3. Data from the X-axis gyroscope of a stationary PS Move

The Earth rotates 360 ° in one sidereal day of 23 h 56 min 4 s, i.e. at a rate of 0.0042 °/s. This is 20 times smaller than the gyroscope resolution and well below the noise floor. Obviously, some post-processing of the gyroscope data will be required.

5. Principle of operation 5.1. Naive approach: averaging and common-mode rejection Under the usual assumptions about measurement noise, we should be able to measure Earth's rotation with acceptable uncertainty simply by averaging a few million samples from a stationary gyroscope aligned in the north-south direction. The resulting value would reflect Earth's rotation plus an unwanted offset called the zero-rate bias. We would subtract that offset after estimating it by performing a second measurement with the same gyroscope aligned east-west (where it is not affected by Earth's rotation). Unfortunately the zero-rate bias of a gyroscope cannot be calibrated out so easily because it tends to drift slowly and unpredictably. Bias stability is a crucial parameter of the performance of a gyroscope. For obvious reasons, devices with excellent bias stability are even classified as weapons components under ITAR. With inexpensive consumer-grade MEMS gyroscopes, the bias drifts within minutes. Therefore we cannot expect consistent results across successive measurements lasting several hours each. Figure 4. Bias drift over 11 hours.



5.2. Reduction of bias and low-frequency noise by homodyning In [MST2005], a stepper motor mechanically toggles the orientation of a gyroscope back and forth 180 ° every 1.8 s, and the gyroscope is sampled at rest for a fraction of a second at each extremal orientation. The bias value remains stable over such short time scales. Thus, it is effectively removed by subtracting one set of measurements from the other. [BMR] analyzes the benefits of slewing gyroscopes and accelerometers in sophisticated patterns with two degrees of freedom. Our approach is similar. Instead of building a custom motorized mount and control system, we spin a gyroscope continuously on an off-the-shelf vinyl record player. Figure 5. Homodyning

From a signal processing perspective, we modulate the signal of interest mechanically so as to shift it out of a noisy part of the spectrum, and then we demodulate at the same frequency with the software equivalent of a lock-in amplifier. When applicable, this is a popular technique for dealing with non-uniform noise spectra, and with low frequency noise in particular. Figure 6. Noise spectrum

The output of the demodulator is a complex number which provides both the rate of Earth's rotation (multiplied by the cosine of latitude) and the direction of its axis with respect to the phase origin of the modulation.

6. Implementation 6.1. Turntable



Direct-Drive Turntable



Direct-Drive Turntable



Direct-Drive Turntable



Direct-Drive Turntable

Affiliate links Record players are useful as amateur scientific instruments because audiophiles tend to have high expectations regarding mechanical construction, speed accuracy and stability, vibrations, magnetic shielding of motors, etc. For our first attempt we used the least expensive unit we could find, a belt-driven model. After measuring a rotation rate of 46 RPM instead of 45, we concluded that it did not even have a crystal-controlled synchronous motor, and we upgraded to an entry-level direct-drive model (130 EUR retail). The direct-drive unit proved satisfactory. Since we had to implement robust synchronization mechanisms in software eventually, a precise turntable does not appear to be a crucial requirement anymore, and we encourage readers to try and replicate the experiment with any unit they may have access to. The PS Move is mounted on the platter in its natural orientation. Either the X-axis (pitch) or the Y-axis (roll) gyroscope can be used to sense Earth's rotation. Figure 7. Axes



6.2. Effect of magnetic fields Early versions of our apparatus failed to produce conclusive results. It took some time to discover that gyroscope measurements were being affected by Earth's weak magnetic field. This was confirmed by measuring the response of a stationary gyroscope to a permanent magnet spinning nearby. We found no mention of such an effect in publications about MEMS gyroscopes. Actually, magnetic fields could be affecting not the MEMS sensors themselves, but the associated electronic circuitry.



Regulated Linear Power Supply 0-15V 0-3A



Regulated Linear Power Supply 0-15V 0-3A



Regulated Linear Power Supply 0-15V 0-3A



Regulated Linear Power Supply 0-15V 0-3A

Affiliate links Magnetic shielding might help address this issue, but mu-metal is hard to procure and would interfere with Bluetooth communications. Instead we use a Helmholtz coil to cancel Earth's magnetic field in the vicinity of the gyroscope. We spent more time than necessary building this coil; a crude version could be made by wrapping electrical wire around a cardboard box. Figure 8. Coil: 40 cm diameter, 2x10 turns, 400-600 mA



6.3. Clocking The PS Move, the host computer and the turntable each have their own clock, and we have no direct way of synchronizing them. The host computer is synchronized to a global time reference via NTP. But since the operating system and Bluetooth communication stack introduce unpredictable latency, high resolution time-stamping of incoming packets has little value. We only take advantage of NTP synchronization to accurately measure sampling rates and rotation rates over long periods (several hours).

The reporting rate of the PS Move seems to be slightly affected by temperature (at least 20 ppm/K in one of our tests). Samples are transmitted with what looks like a time-stamp field, but we could not find a use for it and we simply assumed that the sensors are sampled at regular intervals. The sequence number field allows detection of packet loss on the wireless link.

The turntable appears to have a precise and stable rotation rate, but this is not good enough to keep it in-phase with the computer over several hours. Thus, a crucial component of our homodyne signal processing is a phase-locked loop which synchronizes the software-based demodulator with the mechanical modulation performed by the turntable. In our current implementation, the PLL is simply derived from another sensor of the PS Move, with narrow-band filtering around 0.56 or 0.75 Hz. If a magnetometer is used as reference, the PLL will be aligned with the local residual magnetic field. (Despite the action of the Helmholtz coil, the magnetometers detect Earth's field because they spin 80 mm off-axis.) If an accelerometer is used as reference, the PLL can be aligned in a specific direction by deliberately tilting the turntable. Figure 9. Phase-locked loop

6.4. Effect of linear acceleration Being mechanical devices with vibrating structures, MEMS gyroscopes tend to respond not only to rotation rates, but also to linear accelerations. Manufacturers know how to mitigate this with pairs of counterbalanced proof masses, and any residual effect is small enough to not even be mentioned in datasheets of gyroscopes marketed for mainstream applications. However, even aircraft-grade gyroscopes claim sensitivities as high as 0.05 °/s/g. This implies that a linear acceleration of only 0.08 g may cause an unwanted response as large as Earth's rotation's. In our apparatus, the gyroscope experiences linear accelerations from gravity and from the rotational motion of the turntable. We carefully align the rotation axis of the turntable vertically. Hence, in the rotating frame of reference, gravity and centrifugal forces are stationary, and we hope that their effect is cancelled by the homodyning process. 6.5. Effect of fluctuations in the rotation rate of the turntable The rotation of the turntable will be sensed by any gyroscope whose axis is not perfectly orthogonal to that of the turntable. The effects of its DC component and random fluctuations will be cancelled by the homodyning process, but any synchronous fluctuation will directly affect the final result. Synchronous fluctuations may be caused by: friction in the turntable bearings

an unbalanced payload spinning around a non-vertical rotation axis

inhomogeneous stiffness of rubber belts or similar components

the switched nature of magnetic forces in the motor. In Figure 10, “ Fluctuation of turntable rotation rate ” the 16th-harmonic fluctuation strongly suggests that the motor has 16 magnetic poles. The fundamental component of about 1000 ppm is a more serious concern. We have not determined its cause. Figure 10. Fluctuation of turntable rotation rate

It should be possible to damp fluctuations mechanically with a flywheel attached to the platter via flexible mounts, but our attempts with household materials proved unsuccessful. We also tried to estimate the relative amplitude of fluctuations and subtract their predicted impact on the demodulated output, with limited success. Eventually we chose to reduce the effect of fluctuations as much as possible by carefully aligning the gyroscope axis with the rotation plane. Eventually we chose to reduce the effect of fluctuations as much as possible by carefully aligning the gyroscope axis with the rotation plane. 6.6. Coil alignment The coil is aligned with the horizontal component of Earth's magnetic field. The vertical component needs not be cancelled because it is stationary in the vicinity of the gyroscope, at the center of the rotating frame. Figure 11. Aligning the inactive coil with Earth's magnetic field

The sensing gyroscope is located 22 mm behind the "PS" button, in the upper half of the device. It must spin at the center of the coil. 6.7. Horizontal alignment of the turntable Any inclination of the turntable with respect to the horizontal plane causes gravity in the rotating frame of reference to wobble in a way that can be measured by the accelerometers of the PS Move, using the homodyning technique again. Repeatability is better than 0.05 °; a rigorous estimation of accuracy has not been carried out. Any inclination is corrected manually by inserting strips of paper (about 100 µm thick) under the feet of the turntable. 6.8. Gyroscope axis alignment To ensure that the sensing gyroscope is orthogonal to the axis of the turntable, we simply measure its response to the rotation of the turntable (ideally, it should be zero). This is done by averaging the signal and subtracting the zero-rate bias, all within one minute. Repeatability is better than 0.05 °, but we do not make any claim about accuracy here either. 6.9. Software 6.9.1. Data acquisition linmctool (see [LINMCTOOL]) receives data from the PS Move: # wget http://www.pabr.org/linmctool/linmctool-20110206.c # gcc --std=gnu99 -Wall linmctool-20110206.c -lusb -o linmctool # ./linmctool --timestamp | bzip2 > data.dat.bz2 6.9.2. Processing copernitron-20110210.m is a script for GNU Octave / MATLAB with the following minimal functionality: Parsing data files recorded by linmctool .

Homodyne demodulation synchronized with the X- and Y-axis magnetometers.

Estimation of alignment errors between the turntable and the horizontal plane.

Simplified estimation of alignment errors between the PS Move and the rotating platter.

Estimation of the residual magnetic field. This actually reflects the radial and tangential components of the field on a circumference about 80 mm from the axis, not the field that the gyroscopes experience.

Measurement of Earth's rotation rate, in rotations per sidereal day, with a confidence interval accounting for sensor noise and gyroscope alignment errors.

Estimation of the direction of true north with respect to the residual magnetic field (which might differ from magnetic north).

7. Datasets and results The figures below show measurements of the rate and direction of Earth's rotation obtained from each dataset and each horizontal gyroscope (gX and gY). The demodulator outputs a complex number which we normalize to the expected value (blue circle: one rotation per sidereal day) and display with the PLL phase origin up. Note that error margins (green and red circles) account for known effects only: sensor noise and gyroscope alignment errors. Even with 17 hour long recordings, the process is not entirely reliable. Unfortunately we cannot perform a control experiment with Earth's rotation switched off and all other things equal; but several facts suggest that the signal actually reflects Earth's rotation rather than experimental artifacts or other directional phenomena: The two horizontal gyroscopes, gX and gY, provide two independent results from each dataset. We can check whether the two results are in agreement (up to error margins and a 90 ° phase difference).

Repeatability is tested by processing independent fragments of each dataset (green circles).

We can slightly alter the alignment of the turntable with the horizontal plane and with the gyroscope and verify that the demodulated signal still points toward the same azimuth.

Similarly, we can rotate the chassis of the turntable by 90 ° to ensure that the signal does not reflect a preferential direction of mechanical friction.

Across all datasets, even with poor alignment and large residual magnetic fields, the measured rotation rate is within ±50 % of the expected value. Each dataset is about 100 MB. To download, please change the URL extension to ".dat.bz2". stationary-magnet-off.dat. This is 11 hours of noise from a stationary device.

Figure 12. dd45eeehelm.dat. Residual magnetic north is up.

dd45eeehelm.dat. The X-axis gyroscope is poorly aligned; the Y-axis gyroscope yields a reasonable margin of error. dataset = dd45eeehelm latitude = 48 rpm_nominal = 45 count = 9894401 duration = 5.5933e+04 Sample rate: 176.898 Hz (ignoring packet loss). mX: Residual mag. north 2.3 µT at -0 ° CW from PLL origin. mY: Residual mag. north 2.3 µT at 1 ° CW from PLL origin. aX: Turntable tilt 0.011 ° up at 78 ° CW from PLL origin. aY: Turntable tilt 0.013 ° up at 72 ° CW from PLL origin. Rate fluctuation: 1057 ppm. Using gyro gX. Range: 12.175 °/s mean: -0.637 °/s stddev: 0.474 °/s. Tilt: 0.135 ± 0.021 °. Earth rotation rate: 0.668 rotations per sidereal day ±0.309 with 95% confidence. True north: -10.4 ° clockwise from PLL origin. Using gyro gY. Range: 15.219 °/s mean: 0.005 °/s stddev: 0.406 °/s. Tilt: -0.001 ± 0.021 °. Earth rotation rate: 0.859 rotations per sidereal day ±0.105 with 95% confidence. True north: 31.5 ° clockwise from PLL origin.

Figure 13. dd33alignedB.dat. Residual magnetic north is up.

dd33alignedB.dat. The coil current is deliberately set slightly higher than needed (530 mA) in order to confirm that the experiment detects the direction of true north, rather than that of the local magnetic field via unintended effects. dataset = dd33alignedB latitude = 48 rpm_nominal = 33.333 count = 9894401 duration = 5.6350e+04 Sample rate: 175.589 Hz (ignoring packet loss). mX: Residual mag. north 2.0 µT at -1 ° CW from PLL origin. mY: Residual mag. north 2.0 µT at 1 ° CW from PLL origin. aX: Turntable tilt 0.006 ° up at 53 ° CW from PLL origin. aY: Turntable tilt 0.006 ° up at 65 ° CW from PLL origin. Rate fluctuation: 1008 ppm. Using gyro gX. Range: 4.834 °/s mean: -0.535 °/s stddev: 0.474 °/s. Tilt: 0.153 ± 0.029 °. Earth rotation rate: 0.750 rotations per sidereal day ±0.269 with 95% confidence. True north: -144.7 ° clockwise from PLL origin. Using gyro gY. Range: 4.566 °/s mean: 0.258 °/s stddev: 0.386 °/s. Tilt: -0.074 ± 0.029 °. Earth rotation rate: 0.547 rotations per sidereal day ±0.173 with 95% confidence. True north: -137.2 ° clockwise from PLL origin.

Figure 14. dd45aligned.dat. Residual magnetic north is up.

dd45aligned.dat. The turntable and sensors are carefully aligned. We suspect that natural fluctuations of Earth's magnetic field make the direction of the residual field unstable, yielding an unreliable PLL. dataset = dd45aligned latitude = 48 rpm_nominal = 45 count = 9894401 duration = 5.6344e+04 Sample rate: 175.607 Hz (ignoring packet loss). mX: Residual mag. north 0.9 µT at -1 ° CW from PLL origin. mY: Residual mag. north 0.9 µT at -1 ° CW from PLL origin. aX: Turntable tilt 0.002 ° up at 93 ° CW from PLL origin. aY: Turntable tilt 0.003 ° up at 147 ° CW from PLL origin. Rate fluctuation: 1049 ppm. Using gyro gX. Range: 4.745 °/s mean: -0.629 °/s stddev: 0.460 °/s. Tilt: 0.133 ± 0.021 °. Earth rotation rate: 0.771 rotations per sidereal day ±0.302 with 95% confidence. True north: -74.4 ° clockwise from PLL origin. Using gyro gY. Range: 4.029 °/s mean: 0.333 °/s stddev: 0.381 °/s. Tilt: -0.071 ± 0.021 °. Earth rotation rate: 0.521 rotations per sidereal day ±0.199 with 95% confidence. True north: -68.9 ° clockwise from PLL origin.

Figure 15. dd45tiltup115.dat. Azimuth 115 °E is up.

dd45tiltup115.dat. The turntable is deliberately tilted up toward azimuth 115 °E in order to locate true north with respect to a known direction. This tilt seems to affect phase measurements in a way that is not understood yet. dataset = dd45tiltup115 latitude = 48 rpm_nominal = 45 count = 8199671 duration = 4.6695e+04 Sample rate: 175.602 Hz (ignoring packet loss). *** PLL TILT *** mX: Residual mag. north 1.1 µT at -39 ° CW from PLL origin. mY: Residual mag. north 1.0 µT at -40 ° CW from PLL origin. aX: Turntable tilt 0.065 ° up at -1 ° CW from PLL origin. aY: Turntable tilt 0.062 ° up at 1 ° CW from PLL origin. Rate fluctuation: 1013 ppm. Using gyro gX. Range: 4.834 °/s mean: -0.310 °/s stddev: 0.476 °/s. Tilt: 0.066 ± 0.021 °. Earth rotation rate: 1.225 rotations per sidereal day ±0.214 with 95% confidence. True north: -167.7 ° clockwise from PLL origin. Using gyro gY. Range: 4.029 °/s mean: 0.190 °/s stddev: 0.374 °/s. Tilt: -0.040 ± 0.021 °. Earth rotation rate: 0.884 rotations per sidereal day ±0.158 with 95% confidence. True north: -60.8 ° clockwise from PLL origin.

8. Perspectives Hardware improvements . Our experiments are typically started in the evening and left to run for up to 17 hours. The PS Move is attached to a cardboard column with rubber bands, and the rubber feet of the turntable stand on a wooden table. All of these materials may expand or contract in response to changes in ambient temperature or humidity. Mechanical alignments could be improved with more suitable materials. Mechanical modelling of vibrating MEMS gyroscopes . We implicitly assumed that our gyroscopes respond instantaneously to variations of the rate and direction of their rotation, ignoring the action of the Euler force on the vibrating structure. We also ignored the impact of fluctuating linear accelerations. On the one hand, this is not unreasonable because MEMS gyroscopes typically vibrate at frequencies three or four orders of magnitude higher than that of the turntable. On the other hand, homodyne demodulation might be sensitive enough to detect periodic dynamic effects. A better model could improve results and explain the anomalies in Figure 15, “ dd45tiltup115.dat. Azimuth 115 °E is up. ”. Binary modulation instead of continuous rotation . For this project, a major challenge was to require only household items and minimize construction work. This led to using a game controller, a turntable, and harmonic modulation of the signal of interest. But the approach of [MST2005], where the gyroscope is sampled at rest in two opposite orientations, avoids many of the complications of continuous rotation. Readers are encouraged to devise mechanical contraption made from common items which would achieve suitable motions (crank and connecting rod ? windshield wiper motor ? LEGO ?). Note that the motor needs not be controlled by the computer: we can use the magnetometers to trigger sampling whenever the gyroscope reaches a proper orientation. Modelling the effect of magnetic fields . The Helmholtz coil is the only component in our device which requires some hardware tinkering, and this detracts from the simplicity of the project. Magnetic shielding should be investigated. As an alternative, it might be possible to understand and model how magnetic fields affect measurements and use that knowledge to cancel their effect in software. Unknown preprocessing ? . What we considered as raw sensor data throughout this article might actually be preprocessed in undocumented ways inside the PS Move. For example, we noticed a peak at 88 Hz in the noise spectrum. This suggests that the two samples in each Bluetooth packet are processed differently. Measuring latitude . Latitude can be derived from the known value of Earth's rotation rate and a local measurement of its horizontal component. Alternatively, we could spin a gyroscope in a vertical meridian plane. Multiple models of the PS Move ? . As mentioned in Section 4, “ Back-of-the-envelope calculations ”, there are at least two versions of PS Move hardware (redesigns are routine throughout the lifetime of industrial products). On the one hand, we cannot guarantee that the experiment will work with versions other than CECH-ZCM1E; On the other hand, it might actually work better with those which have built-in magnetic shielding around the gyroscope. Other devices with gyroscopes . As far as we know the PS Move has the gyroscopes with the highest dynamic range (16 bits) in its market segment. But it also has an unusually wide full scale range (±3000 °/s), which is excellent for motion tracking, but suboptimal for our project. Many MEMS gyroscopes have a secondary "amplified" analog output with a reduced range. According to [WIIBREW], the Wii MotionPlus has gyroscopes with only 14 bit dynamic range, but it does support a "slow" mode with a narrow range. Therefore it should be able to detect Earth's rotation with comparable performance. If it were discovered that the gyroscopes of the PS Move can be configured to operate with a narrower range too, acquisition time might be reduced. Mobile app . We are investigating whether the gyroscopes found in mobile phones are suitable for this experiment. Detect tidal forces ? . Both the Moon and the Sun cause gravity at Earth's surface to fluctuate by about 0.1 ppm. One could attempt to detect this effect with inertial MEMS sensors, using similar techniques.