PHASE and AMPLITUDE VARIATION of CHANDLER WOBBLE by John
PHASE AND AMPLITUDE VARIATION
OF CHANDLER WOBBLE
by
John Alexander Linton
B.A. Sc., University of British Columbia, 1970
A Thesis Submitted in Partial Fulfilment of the Requirements for the Degree of
Master of Science
in the Department
of
Geophysics
We accept this thesis as conforming to the
required standard
The University of British Columbia
May, 1973 In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of British Columbia, I agree that the Library shall make it freely available for reference and study.
I further agree that permission for extensive copying of this thesis for scholarly purposes may be granted by the Head of my Department or by his representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission.
Department of G>e
The University of British Columbia Vancouver 8, Canada
Date Z7o^7 ABSTRACT Normally, the wobble of the earth has been dealt with in a manner that assumes that the two main periodic components have constant phase and amplitude. An initial assumption of this thesis is that both these parameters can vary with time. A technique of predictive filtering is used to determine the Chandler component of the wobble from basic latitude measurements at the five I. L. S. observatories. A simple analytic proce• dure is employed to obtain the phase and amplitude variation of the periodic Chandler motion. The results indicate that major changes in both phase and amplitude occur in the period 1922. 7 to 1949. These changes are possibly associated with earthquake activity, although there is nothing to indicate that there is a correlation between individual earthquakes and events in the Chandler motion. The calculated period of the Chandler wobble is 437 days and the damping time is so uncertain that a value approaching infinity is not unlikely. ii TABLE OF CONTENTS Chapter Page I INTRODUCTION 1 EARLY HISTORY OF WOBBLE RESEARCH 2 The Annual Wobble 3 The Chandler Wobble 7 Chandler Damping 10 Chandler Excitation 12 SUMMARY 16 II THEORY OF EARTH WOBBLE 17 Choice of Axes 19 Perturbation Equations 19 Rotational Deformation 23 Special Solutions to Equation (2.9) 24 III THE LATITUDE SERVICES 30 Method of Observation 32 Reduction of Data by the I. L. S. . . . 33 Instantaneous Pole Path 34 The Problem , . 35 IV PREPARATION OF EQUALLY SPACED LATITUDE MEANS AT EACH OBSERVATORY . . 37 Time Variable Noise Removal 45 V FREQUENCY DOMAIN FILTERING TO OBTAIN CHANDLER LATITUDE VARIATION 56 Motivation 56 Predictive Filtering 58 The Effect of Phase and Amplitude Shifts on the Fourier Spectrum 62 Band Reject Filtering of Annual and Noise 70 iii TABLE OF CONTENTS (Continued) Chapter Page VI THE PHASE AND AMPLITUDE RESULTS FOR THE 5 I. L. S. OBSERVATORIES 78 Phase Results 87 Amplitude Results 87 Interpretation with Respect to Established Theories 90 Conclusions 95 BIBLIOGRAPHY 96 APPENDIX I 99 iv LIST OF TABLES Table Page I Annual Components of the Rotation Pole in Units of ". 01 5 II Estimates of Chandler Period 5a III Earthquakes of magnitude greater than 7. 9 (after Richter, 1958) 9Qa V LIST OF FIGURES Figure Page 1 Power spectrum of Polar Motion and Latitude Data (taken from Munk and McDonald, 1960) 4 2 Power spectrum of Polar Motion (taken from Rudnick, 1956) 4 3 Chandler polar motion amplitude characteristics (taken from Yumi, 1970) 9 4 Chandler period variation (taken from Iijimi, 1971) 9 5 Power spectrum showing split Chandler peak (taken from Rochester and Pedersen, 1972) 11 6 Smylie and Mansinha's fitting of smooth arcs to observed polar motion (taken from Scientific American, Dec. 1971) 15 7 Orientation of reference axes . 20 8 Polar motion corresponding to damped free wobble of the earth 26 9 Response of damped earth as a function of frequency 27 10 Longitude of I. L. S. observatories 31 11 Configuration of star pair with respect to zenith and celestial pole 32 12 Filter to Remove Random Noise 46 13 Output from variable noise filter (solid line) and constant noise filter (solid line with circles). The dotted line is the original series .... 49 14-18 Plots of basic smoothed data from each of the observatories 51-55 vi LIST OF FIGURES (Continued) Figure Page 19 Fourier transform of 10 years of a perfect sinusoid (period = 435 days). The dotted vertical line represents transform of an infinite sinusoid of the same frequency 57 20 Fourier transform of synthetic wobble data of Figure 22 59 21 Method for Predictive Extension of a Data Series 61 22 Synthetic latitude variation composed of 2 main frequency components (annual and Chandler) plus random noise 63 23 Prediction of time series in Figure 22 to five times its original length 64 24 The Fourier amplitude spectrum of the extended series (Figure 23) 65 25 Chandler variation used to create the series of Figure 22. Composed of sine wavelets all of the same frequency (1/435 days"*) but each varying in phase and amplitude 66 26 Fourier amplitude spectrum of Figure 25. The dashed line is the spectrum of a sinusoid with constant phase and amplitude 67 27 Synthesis of the Fourier transform of a truncated sinusoid 69 28 Effect of band-reject filtering of annual on the spectrum of the extended synthetic latitude variation (Figure 23 and Figure 24) 72 29 Chandler variation resulting from band- reject filtering of the synthetic data. The dashed line represents the known synthetic Chandler variation 73 vii LIST OF FIGURES (Continued) Figure Page 30 Phase variation of the synthetic Chandler. The solid line is the phase of the synthetic Chandler motion recovered by band-reject filtering and the dashed line is the phase of the known original input Chandler 76 31 Amplitude variation of the synthetic Chandler. Solid and dashed lines represent the filtered and actual versions respectively . . 77 32-36 Chandler latitude variation for each of the observatories 79-83 37 Chandler amplitude variation for all observatories 85 38 Chandler phase variation for all observatories 86 39 Amplitude variation averaged over number of observatories 88 40 Phase variation averaged over number of observatories 89 41 Chandler amplitude variation after Guinot (1972) 91 42 Chandler phase variation after Guinot (1972) 91 43 Fourier amplitude spectrum of simplified undamped Chandler motion 94 viii ACKNOWLEDGEMENT I wish to thank Doug Smylie and Tad Ulrych for their generous support and encouragement. 1 INTRODUCTION At first glance, the description of the earth's rotation about its axis would seem a simple task. That is. if the earth is treated as a rigid, rotating sphere, one could simply apply the relevant classical laws of physics as developed by Euler. However, the real earth is neither rigid nor spherical. These deviations from the 'ideal' body result in complex perturbations in the rota• tion of the earth. The bulging of the earth at the equator produces the varia• tion in the direction of the axis of angular momentum in space(precession and nutation). Another effect is the variation in the length of the day. The most important perturbation to geophysicists, and the subject of this thesis, is the periodic variation in the orientation of the instantaneous rotation axis with respect to the geographic frame of the earth - or alternately, the movement of the rotational pole (point near North pole where instantaneous rotation axis pierces surface) over a horizontal plane parallel to the equator. This motion is called wobble. Most measurements of wobble are in terms of the motion of the rotational pole in units of distance - the complete planar motion being described by two direction components. Succeeding chapters of this thesis will deal with the history of earth wobble research, the classical theory, and the pole path data, all with a view of justifying a new approach to the analysis of wobble data. 2 EARLY HISTORY OF WOBBLE RESEARCH Since the history of research involving the earth's wobble is thoroughly documented by Munk and McDonald (I960) in their book on the rotation of the earth, only the highlights are discussed here. It was Euler (1765) who first predicted that a rigid sphere would, in addition to its normal rotation, wobble freely about its axis of greatest rotational inertia. He pre• dicted that the period of the wobble should be 305 days. Pe'ters, Bessel, Kelvin and Newcomb all searched without success (or with false success) for a motion of the predicted period. In 1888, Kdstner discovered a variation in the constant of aberration that he attributed to a . 2 second annual variation in latitude. To show that this was in fact a variation due to wobble, measure• ments were made at stations in Berlin and Waikiki - roughly 180 degrees apart - and the measurements showed the expected opposite change in latitude. S. C. Chandler in 1891 announced that the wobble, in addition to an annual component, also contains a motion of period 428 days. Newcomb in the same year showed that the 428 day period could, in fact, be the free wobble period lengthened from 305 to 428 days by the non-rigidity of the earth. These dis• coveries and the lack of confidence in the observational evidence to that time led to the establishment of the International Latitude Service (I. L. S. ) which functioned from 1899. 0 to 1962. 0 when it was replaced by the International Polar Motion Service. Since the time the I. L. S. was established, the major thrust 3 of research has been centered on two aspects of the wobble problem: 1. Analysis of the data and determination in detail of the characterist• ics of the two main frequency components (Chandler and annual). 2. Explanation of these characteristics by a plausible model for the real earth. The annual wobble The annual wobble has, for convenience's sake, been universally treated as a spectral delta function which, in more familiar terms, is just simple circular or elliptical motion of the rotational pole. On first glance at the power spectrums depicted in Fig. 1 & 2 (both are taken from Munk and McDonald), this assumption would seem to be justified by the narrowness of the annual peak. The power spectra for perfect circular motion of period one year would have a peak width of approximately . 04 cycles/year compared to the peak width in these plots of . 03 to . 06 cycles/year. Unfortunately, the simplicity this implies is not quite achieved. Walker and Young (1957) point out. in their exhaustive analysis of I. L. S. data that, "the components of the annual motion vary according to type of series and interval, there being noticeable differences in both the amplitudes and phases of the several compo• nents. " Table I gives various estimates of the amplitude for each of the two directional components of annual polar motion - rrij towards Greenwich, rr^ towards ninety degrees east of Greenwich, The amplitude varies by almost a factor of two between minimum and maximum. This variation means that if the annual wobble is to be modelled as simple circular motion, the model I- -E&gure 1 i Power spectrum of Polar Motion and Latitude Data (taken from Munk and McDonald, I960) Figure 2 : Power spectrum of Polar Motion (taken from Rudnlck,1956) 5 TABLE I ANNUAL COMPONENTS OF THE ROTATION POLE IN UNITS OF ". 01 # Source Interval mi m2 Jeffreys 1952 1892-1938 - 3.6cos® -8. 5 sin© 7. Ocos ® -2. 9s Pollak 1927 1890-1924 -3. 7 -8.9 7. 0 -3.9 Rudnick 1957 1891-1945 -3. 2 -8. 2 6. 7 -2. 8 Walker & Young 1899-1954 -6.4 -7. 1 7. 0 -4.6 1957 1900-1934 -5. 5 -7. 0 7. 5 -4.6 1900-1920 -4.8 -6. 0 6.6 -3. 7 Jeffreys 1940 1912-1935) and ) -3. 2 -7. 8 5. 6 -1.6 Markowitz 1916-1940) 1942 * ® is the longitude of the mean sun measured from the beginning of the year. # after Munk & McDonald, I960. 5a TABLE II ESTIMATES OF CHANDLER PERIOD^ Time Interval Values of the Period, Years Source 1890-1915 1.13, 1.20, 1.27 Witting (1915) 1890-1924 1.20 Pollak (1927) 1890-1924 1.19 Stumpff (1927) 1890-1922 1.08, 1.14, 1.19, 1.27 Wahl (1938) 1922-1938 1.13 II 1892-1933 1. 223 Jeffreys (1940) 1908-1921 1. 202 n 1900-1940 1.108, 1.170, 1.208, 1.250 Labrouste (1946) 1891-1945 1. 196 Rudnick (1956) 1900-1920 1. 193 Walker & Young (1957) 1890-1924 1. 191 Danjon & Guinot (1954) 1929-1953 1. 172 11 1897-1957 1. 186 Pachenko (I960) 1897-1922 1.190 it 1930-1957 1.180 n 1891 -1952 1. 193 Arato (1962) * after Wells, 1972. 6 must be fitted to as short a data length as possible - short enough that the wobble does not depart significantly from circular motion over the interval. Yashkov (1965) and Iijimi (1971), using harmonic analysis on ten and twelve year data intervals respectively, claim to detect variations in period as high as , 05 years from the mean period of one year. Considering the vagaries of harmonic analysis, these results are somewhat questionable. Wells (1972) discovered that the annual component of the latitude variation at an observa• tory (two or more latitude observatories are required to produce a pole path - see Chapter 2) is very inconsistent among observing stations in both phase and amplitude. He goes on to suggest that "the annual polar motion as deter• mined from the usual least squares solution from the latitude series, does not accurately reflect the true annual wobble of the earth. " If the annual latitude variation is markedly different from station to station, one must assume that a portion of the annual variation at any station is due to local effects. This implies that the annual latitude variation at a station is not entirely due to wobble of the earth. What the exact cause of these local variations might be is a matter of conjecture - possibly variations in local vertical, temperature, and 'seeing' . Whatever the cause, it is extremely desirable that, if one is going to remove the annual component in order to examine the Chandler component, the removal should be done on the latitude series before they are used to calculate a pole path. Also, if one wishes to 'see' a real annual wobble, the local effects must somehow be considered. Despite this, the best practical model for the annual component of wobble remains a spectral delta function adapted to vary in amplitude with time. However, care must 7 be taken in equating the annual wobble that best fits this model to real wobble of the earth. The annual wobble has, ever since its discovery, been attri• buted principally to a seasonal shift in air mass over the surface of the earth - along with less significant variations in snow, vegetation and ocean mobility. Spitaler (1901), Rosenhead (1929). and Munk and Hassan (1960) all conclude that the seasonal air mass shift has the phase and magnitude to account for the observed wobble. However, on a careful reading of the Munk and Hassan paper, we find little evidence to support their conclusion. The excitation function required to produce the observed wobble and the excitation function calculated from meteorological records differ by more than an order of magnitude - even if looking at just the annual component of each of these functions. Their subsequent conclusion that "the annual is principally due to seasonal air shifts", is, at best, premature. If air shifts are the cause of the annual wobble, then variations in climate from year to year would reflect themselves as variations in the frequency, amplitude and phase of this compo• nent. This emphasizes our earlier conclusion that the annual wobble should be treated as time-variable in phase and amplitude. The Chandler wobble The Chandler wobble represents the free motion (Jeffreys, 1957) of the earth in response to some perturbation of the instantaneous rotation axis away from its equilibrium position. Since free motion is sensitive to the physical characteristics of the system involved, detailed knowledge of the 8 Chandler motion enables one to draw conclusions about the rigidity, shape and non-elasticity of the real earth. The essential feature of the Chandler wobble is that its spectral representation resembles that of a damped oscillator - that is, the spectral peak is spread over frequencies near the peak value (see Fig. 1 & 2). Table II (taken from Wells, 1972) shows the Chandler period (period at peak) during various epochs as calculated by the authors indicated. This table illustrates that the Chandler wobble might have time-variable amplitude and frequency. Plots by Yumi (1970) and Iijimi (1971) (Fig, 3 & 4), indicate how the period and amplitude of the two Chandler components vary with time. The amplitude variation is especially remarkable - between ten and thirty feet; It should be pointed out that these authors have used various forms of harmonic analysis - sometimes-on data lengths as short as ten years (Yumi) - and therefore results that indicate a shift in frequency are subject to large errors due to resolution problems of the technique. Estimates of the damping time (time for signal to decay to 1/e of its original value) vary from two years (Walker and Young, 1957) to ninety years (Pachenko). These estimates of damping time lead to a Q factor (dimensionless measure of energy dissipation) of thirty to fifty, which is several times smaller than Q values calculated from seismic data. Fellgett (I960) calculated a damping time of 12.4 years but concluded that the value is "uncertain by a factor of at least ten. and damping times as short as 2. 5 years or as long as several hundred years are not excluded. " The 9 « g -io\ a. E o —i aoL l i L j i L _j 1 i 1993 ISOS 1317 1929 19*/ 19S3 I96ST Figure 3 : Chandler polar motion amplitude characteristics (taken from Yuml, 1970) 1M\ i—i—i i r i—i—r T— — — l.20\ a> 1.16 T3 o 1.14- ** L rn. J JL i i i J_ 190072 2436 Figure 4 : Chandler period variation (taken from Iijiml, 1971) 10 apparent disparity between Q factors still remains. One solution (Wells. 1972) would be a frequency dependent Q. That is, the earth would have low damping for seismic frequencies and high damping for wobble frequencies. There is no evidence to support this hypothesis. Some (Melchior, Niccolini, 1957) explain the broadened spectral peak of the Chandler component in terms of a variable frequency - of which we have already some evidence. This explanation requires a variation of four percent in the Chandler period. This is not apparent in the observations. Yashkov (1965) carried out harmonic analysis on polar motion data and produced a split Chandler peak. Fig. 5 shows a similar split peak obtained by Rochester and Pedersen (1972) from a seventy year record of monthly mean pole positions. Colombo and Shapiro (1968) assume the split peak represents two real neighbouring frequencies and on that basis resolve the Q factor disparity. Unfortunately for them, Fedorov and Yatskiv (1965) proved that this apparent splitting could in fact be caused by a simple phase shift in a pure harmonic of the Chandler frequency. Thus, the general characteristics of the Chandler wobble can be summarized as: a time-variable wide-band periodic function of central period 435 days and a spectral width typical of a damped system with a damping time estimated to be ten to thirty years. Chandler damping It is generally accepted that the apparent damping of the Chandler wobble, as indicated by its spectral width, is a real damping due to inelastic response of the earth to the free wobble. The major problem 11 100 90 80 ct70 o •"60 O x 50 o cc 40 UJ * ? 30 BAND WIDTH 20 10 LA • /V. VY. .5 .7 .8 .9 1.0 I.I FREQUENCY (CPY) Figure 5 s Power spectrum showing split Chandler peak (taken from Rochester and Pedersen, 1972) TOWARD 90 DEGREES Figure 6 : Smylle and Manslnha's fitting of smooth arcs to observed polar motion (taken from Scientific American, Dec. 1971) 12 has been to explain the discrepancy between the damping required by this wobble characteristic and that observed in seismic work. The conclusions reached by Munk and McDonald in this regard are: 1. Solid friction in the mantle could produce a Q factor of 100 to 200 - as opposed to a Q factor from wobble data of 30 to 50. 2. The additional damping could be due to: (a) the oceans. (b) the lower mantle having a viscosity that is especially important on a time scale of the order of the Chandler period. (c) core-mantle interaction. They conclude by stating that "the situation is appallingly uncertain". There the matter rests (or perhaps wobbles). Chandler excitation Even without observational evidence that the free wobble is damped, it can be assumed that damping does occur (as in all real systems). Thus, left undisturbed, the Chandler wobble would eventually dissipate com• pletely. Therefore, an excitation mechanism is required. The form of this excitation has been the cause of much conjecture in recent years. It was thought that irregularities in the annual atmospheric loading were the cause of the excitation until Munk and Hassan (I960) proved that this mechanism was several times too small. Today, two possible causes of the excitation remain: 1. Core-mantle interaction of some type. 13 2. Redistribution of mass due to earthquakes and subsequent perturbation of the moment of inertia. In regard to the first of these causes, the main interaction occurring at the core-mantle boundary is thought to be electro-magnetic coupling. Munk and Hassan (I960) and Rochester and Smylie (1965) found that this mechanism fails by several orders of magnitude. However, Rochester (1970), Stacey (1970) and Runcorn (1970) argue that effects of secular variation and preces• sion in concert with coupling could possibly account for observed wobble. Until the physical parameters at the core-mantle boundary are better defined, electromagnetic coupling will remain a somewhat confusing issue. The form of the excitation mechanism is likely to be either random catastrophic events that cause sharp discontinuities in the free (circular) motion of the rotational pole or a periodic forcing function of frequency equal to the observed wobble frequency. Earthquake excitation of the wobble has been considered for many years (Cecchini, 1928), but calculations based on simple models yield a negligible effect, (Munk and McDonald, I960). Press (1965) applied dis• location theory to earth models and produced much larger displacement fields than had been anticipated. Smylie and Mansinha (1967) used Press's tech• niques to calculate theoretical displacement fields of earthquakes and the associated change in the moment of inertia. They concluded that the pertur• bations were large enough to maintain the Chandler wobble. If the excitation is due to random earthquakes, then the rotational pole path will be smooth, circular arcs (spirals if damping is considered) separated by sharp breaks. 14 and therefore, Smylie and Mansinha (1968) fitted circular arcs in a least squares fit to the observed Chandler pole path. As can be seen in Figure 6, they allow random shifts in the origin, radius and phase of the motion - a shift occurring when the data and the fitted curve no longer agree within some preset circle of validity. They found excellent correlation between these 'breaks' and major earthquakes. In fact, the 'breaks' tended to occur just before the earthquake so that if the effect was real, earthquakes could be predicted. Haubrich (1970) found that a re-analysis of the same data produced quite different events and so he concluded that the strong correlation found by Mansinha and Smylie was merely a coincidence. Later work by Smylie et al (1970) using a deconvolution method tends to support the original correlation - though not so strongly. Smylie and Mansinha (1971) calculated a pole path displacement caused by the 1964 Alaskan earthquake on a real earth model. Their result of one foot compares well with similar work by Dahlen (1971). Displacements of this order could quite easily sustain the Chandler wobble. However, so far, observational evidence of such shifts is not conclusive. It should be kept in mind that both the damping and excitation mechanisms of the Chandler wobble are based on the assumption of a damping time of about twenty years. An order of magnitude shift in this value, which is perfectly possible according to Fellgett (I960), would have an enormous effect on all the models that are now used to describe the Chandler wobble. For example, if the damping time is actually 200 years, the conflict between seismic and wobble damping factors would no longer exist and the excitation required to produce the wobble would be far smaller - allowing a variety of mechanisms. now thought to be insignificant, to affect the Chandler wobble. 16 SUMMARY The most important facts that can be gleaned from this brief history are that: 1. both the annual and Chandler components are time variable in amplitude and perhaps in frequency. 2. no assumption should be made about the damping time of the Chandler motion. 3. the Chandler may have random sudden changes in phase and amplitude. 17 CHAPTER II THEORY OF EARTH WOBBLE The rotation of bodies about their centre of mass is governed by Euler's equation of motion which states that the time rate of change of angular momentum L is equal to the applied torque P if all measure• ments are made in an inertial system. ~F~ _ T~ 1 inertial Since it is convenient to use a coordinate system attached to the rotating earth (angular velocity *w in space), the transform « • — i i ~~~ = ~~~ 4- :w x L L L inertial rotating is used. The equation of motion is now « 7^ = w x L~ 4- L~ (2.1) where the time derivative is taken in the rotating frame. Total angular momentum is given by = (7)] L /{{tYi [ r x fvv"x r 4- v (r)) p dV earth where v(r) is velocity of material at point "r with respect to the rotating frame and p(r) is the density of material at that point. The term r x (w x T) can be expanded by the vector identity ~ 18 rx(wxr) = w(r.r)- r(w.r) which becomes r x ( w x r ) = e. [r2 &.. - x. x. "1 w. 1 L ij 1 J J J when the standard summation convention (repeated subscripts indicate summa• tion over that subscript) is used. The e"- are unit vectors in the coordinate directions and c$ij is the Kronecker delta. The angular momentum is now p L = w. fr) dV 1 L 1 J eartarth ij J J 4- fff r x v ( r ) p ( r ) dV earth ' e. I., w. 4- J 1 iJ J where I.. = fff \ r2 cf'.. - x. x. 1 p ( r ) dV (Z. 2) XJ earth L *J 1 J J ~ is the inertia tensor of the earth, and = JyT 7xv(r) p (r) dV earth ' is the angular momentum relative to the rotating frame. Substituting equation 2.2 into equation 2. 1 gives the Liouville equation: where t£ is the normal alternating tensor. 19 Choice of Axes The most convenient choice of axes is one in which the coordin• ate system is attached to the observatories - relative motion of the observa• tories due to nonrigidity of the earth being the only complication. The axes are oriented so that the origin is at the centre of mass, the x^ axis is reasonably close to the axis of rotation and the x^ axis passes through the Greenwich meridian. See Fig. 7. If steady, equilibrium rotation about the x^ axis is assumed, the above coordinate system will be aligned with the principal axes of inertia - which are defined as the axes that make I.. = 0 if i 4 j. Thus, if the principal moments of inertia for the earth are denoted A, A and C, then l C I11 = I22 = A AND 33 = where two of the principal moments are equal because of cylindrical symmetry. Perturbation Equations The development of the equations governing the case where the axis of rotation is perturbed from its equilibrium position follows the method of Munk and McDonald (1960). The departures from uniform, diurnal rotation are almost undetectable in practice, so that first order perturbation theory should work extremely well. Therefore, let 20 Z (geographic North pole) instantaneous rotation axis — fixed in space except for effects of Iml precession earth s centre of mass Y (90° E of Greenwich) (to Greenwich) Figure 7 : Orientation of reference axes 21 m m ) w = f\ 1 4- JX 4- e^f\. (1 4- 3 J\ - angular velocity of earth; rri|, m^. << 1 Ill5sA+Cll+bll hi- A+C22 + B22 X33 * C + C33 + B33 X ! =B + C l = I = B + C 13 = 31 13 13 Z3 32 23 23 I12 = I21 = b12 + c12 b is the contribution to I., from rotational deformation above that due to ij XJ steady, diurnal rotation, and c-j is deformation due to other causes. Then, the angular momentum components are L l = hi W1 + I12W2 + I13W3+ A = A-ami+ j\ (b13 + c13) +• x1 wl + L2 = hi *22 W2 + J23 W3 + ^ 2 = AAm2 + A(b23+c23)+ Xz wj L3 = I31 + I32 w2 4- I33 w3 4- Xz A A = C (1 + m3) 4- (c33 4- b33 ) 4- /3 where terms like (b.j 4- c.j)m k J\ are neglected in comparison to those like (b.j 4- C.J)JTL . Substituting into the Liouville equation (2. 3) yields 22 P ~ £ + A J\ m 4- _TL (b , 4- c ) 4- C yL2m_ 11 1 13 13 2 _i2 -jT_ Am2- Jl?(b23 4- c23) - yi (2.4) p = A jT.m 4- JT (b 4- c ) 4- i +AAm 2 ^ 23 2 3 ^ * + SL /I ii 2 (b13 4- c13) 4- - C Jl m (2.5) = CJl m3 4- yi (c33 4- b33) 4- /3 (2.6) where, in addition to the previous approximation, we now neglect products of lj and mj. Now, let m = rrij 4- i r = p + i p — 1 2 i - i ± a x z By multiplying equation 2. 5 by i = (V-l ), and adding the result to equation 2. 4. we obtain m 71 A - i m -A.2 ( C - A.) = _Q - -ij/l - ( cl3 4- i c23)_fL - ( bu 4- i b'23 )_/! 2 2 - i ( c13 4- i c23)/l - i ( b^ + i b23)yi_ (2.7) 23 This is the equation of motion for the wobble - m gives the complex angular displacement of the axis of rotation away from "e^. or. if one wishes, the motion of the instantaneous pole of rotation over the surface of the earth, with . 1 second of angular motion equivalent to ten feet of polar motion. The re• maining equation m = P 3 77T [ 3 " ^3 - < '33+ b33 >^-] <2'8) governs the change in length of day. Rotational Deformation Since the terms b.. are in fact functions of m, it is useful at !J — this point to obtain the exact functional relationship. It can be shown that k_ d m m, _TL e. c. b I 12 = l* 3 G k d m 7 J\ 2 b.23 " 3 G B = k2 d5ml^2 13 3~G where d equals the radius of the earth, G is the gravitational constant and k is the Love number-defined such that if the symmetric earth is subjected to a disturbing force field expressible as the gradient of a second order spherical harmonic potential U^. the additional potential V at the deformed surface is V = k2 U-,. These results are based on the assumption of linear str ess-strain relations, and are derived on page 25 of Munk and McDonald 24 (1960). Substituting for b.^ in equation 2.7 gives m/L(A+ — ) - i rn j\ ( C - A - _2 ) - 3G - 3G = H - ( + 1 1 (c + 1 c ) 1 " " ^13 ' 13 23 ^ If we neglect k^d _/l /3G compared to A, and allow (C-A)/A to be replaced c 2 5 2 by (C-A)/C and k2 -i [JZ -1 -ii^ 2 2J - i (C13 4- i C23) J1 (2.9) ^ A5 Z k?d y, C - A H = (3 The quantity H determines the rate of precession and from direct observa• tion has the value H = 1 /305. 5. Special Solution to Equation 2. 9 A. Free Wobble For this case, the right hand side of (2.9) is zero and so we have m - i G~Q m = 0 i This represents a circular motion of the instantaneous pole of rotation with (T angular velocity 0 and radius A. For a rigid earth with the same preces- sional constant H as the real earth, the angular velocity is given by i CTQ = J\ H and this corresponds to a period of 305. 5 days. No periodicity of this value has been detected in the wobble. For a non-rigid earth, the angular velocity is 5 2 k?d il cr = AH- — o 3GC The only unknown in this equation is k^. and Takeuchi (1950) has computed k^ = . 281 for one of Bullen's earth models. Using this value, the period of the free wobble is 2 TT 435 days 5 2 H - (k2d JT. /3GC) which is close to the observed value for the Chandler wobble. However, the good agreement is probably more due to good luck than anything else, since no consideration is made for the effect of oceans or for core-mantle coupling. According to Munk and McDonald (I960), these two effects seem to just balance each other with regard to the period of the wobble. B. Damped Free Wobble Once again, the right hand side of equation 2.9 is zero. To allow for damping of the motion, 0s is allowed to be complex: cr = w. + i/t o ° 26 w is the angular frequency of the free wobble as calculated previously, and o f is the damping time. Equation 2.9 then yields the solution i w t . - t/Z Q m A e e and the motion of the pole is a west to east spiral as shown in Fig. 8. Fig. 8. Polar motion corresponding to damped free wobble of the earth. C. Harmonic Excitation with Damping Allow to be complex as above, and assume the right i