Physical Geology Research Insights
Physical Geology Research Insights
Author(s): W. Hopkins
Source: Philosophical Transactions of the Royal Society of London , 1839, Vol. 129 (1839),
pp. 381-423
Published by: Royal Society
JSTOR is a not-for-profit service that helps scholars, researchers, and students discover, use, and build upon a wide
range of content in a trusted digital archive. We use information technology and tools to increase productivity and
facilitate new forms of scholarship. For more information about JSTOR, please contact [email protected].
Your use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
https://about.jstor.org/terms
Royal Society is collaborating with JSTOR to digitize, preserve and extend access to
Philosophical Transactions of the Royal Society of London
XX. Researches in Physical Geology. .By W. HOPKINS, Esq. M1A. F.R.S., Fellow
of the Royal Astronomical Society, of the Geological Society, and of the Cambridge
Philosophical Society.-First Series.
On the Phenomena of Precession and Nutation, assuming the Fluidity of the Interior
of the Earth.
BEFORE I proceed to the discussion of the question which forms the principal sub-
ject of the present communication, I shall offer some general remarks on the refri-
geration of the globe, as introductory not only to this memoir, but to others which I
hope hereafter to bring under the notice of the Society.
In the first place, we may observe that there are two distinct processes of cooling,
of which one belongs to bodies which are either solid or imperfectly fluid, and is
termed cooling by conduction, and the other to masses in that state of more perfect
fluidity which admits of a free motion of the component particles among thernselves.
In this case the cooling is said to take place by circulation or convection. The na-
ture of the former process has been ascertained with considerable accuracy by ex-
periment, and the laws of the phenomena have been made the subject of mathematical
investigation, but of the exact laws of cooling by the latter process we are compara-
tively ignorant. It is manifest, however, that since time must be necessary for the
transmission of the hotter and lighter particles from the central to the superficial
parts of the mass, as well as for that of the colder and heavier particles in the oppo-
site direction, the temperature must increase with the depth beneath the surface;
and, moreover, that this increase will be the more rapid, the more nearly the fluidity
of the mass approaches that limit at which this process of cooling would cease, and
that by conduction begin, since the rapility of circulation would constantly diminish
as the fluidity should approximtate to that limit. But still, even in this limniting case,
it seems probable that the ten(lency to produce an equality of temperature throughout
the mass will be much greater, and consequently the rate of increase of temperature
in approaching the centre much less, than if the cooling of the mass had proceeded
by conduction during the same time, the conductive power being very small.
If ttle matter composing thle globe was originally in a high state of fluidity from
heat, the process of cooling would undoubtedly, in the first instance, be by circula-
tion. The manner in which t;he transition will take place fiom this mode of refri-
which would then begin to cool by conduction, rapidly at the surface expo
low temnperature of the planetary space, and extremely slowly in the central
account of the small conductive power of the matter composing the ear
quently the globe would consist, after a certain time, of an exterior solid
interior fluid matter, of which the fluidity would increase in approaching th
where it might still approach to that more perfect fluidity which admits of
convection. With reference, however, to the mechanical action of any f
ducing either mnotion or hydrostatic pressure in the interior mass, the w
might, as an approximation, be considered perfectly fluid. No attempt ha
made to determine the present probable thickness of the earth's crust, assum
have been originally in a state of fluidity, on account of the difficulty alread
tioned, arising from our ignorance of the influence of high temperature i
solidification, compared with that of great pressure in promoting it. All that h
hitherto been determined on the subject is, that the present state of the earth's surf
may be consistent with the existence of a solid crust, of which the thickness is sm
compared with the earth's radius.
Let us now recur to the other case above mentioned, that in which the increase
pressure in descending towards the centre of the mass is supposed to have a grea
effect in promoting solidification than the inciease of tempe ature in preventing
Supposing the Inass to have been first in a state in wich every pa was cooling b
convection, this process would first cease, and that of cooling by conduction be
at the centre, while the superior portion would still continue to cool by convecti
so that these two processes would for a time be going on simultaneously in differ
parts of the mass. It is manifest, however, that the central portion, coolirng by c
duction, would constantly increase, while the exterior portion, cooling by convect
would constantly diminish, so that at length no part of the imass would be cooling b
the latter process. Before it should reach this stage of the refrigeration the cent
portion of a mnass so large as the earth mighlt become perfectly solid, so that at
instant when the circulation should entirely cease, the whole might consist of a s
central nucleus, surrounded by the external portion still in a state of fusion, an
which the fluidity would vary continuously from the solidity of the nucleus to
fluidity of the sutface, where, at the instant we are speaking of, it would be just slc
as not to admit of circulation,
When the mass should have arrived at this stage of the cooling, a change would take
place in the process of solidification, whichl it is important to remnark. Trhe superficial
parts of the mass must in all cases cool the most rapidly, and now (in consequence of
the imperfect fluidity) being no longer able to descend, a crust will be formed on the
surface, from whlich the process of solidification will proceed far more rapidly down-
wards, than upwards on the solid nucleus. Consequently, then, our globe would
arrive at that state, according to the mode of cooling we are now considering, in
which it would be composed of a solid shell, and a solid central nucleus, with matter
in a state of fusion between thlem, the fluidity of which, however, would necessarily
be less than that which might exist in the fused mnass very near the centre in the
case previously considered.
With respect to the thickness of the shell which may be consistent with the pre-
sent appearances of the earth's surface, the same conclusion will hold as in the former
case, i. e. it may be smnall compared with the eartth's radius. What would be the radius
of the solid nucleus at the instant of the first incrutstation of thle surface, or that which
would correspond to any assigned thickness of the exterior shell, it is quite impossible
to determnine from the want of all experimental evidence respecting the tendency of
great pressure to promote solidification at very high temperatures, and our ignorance
of the temperature at which the superficial incrustation of a large mass would begin,
when exposed to the temperature of the planetary space. It is, therefore, manifestly
MDCCCXXXIX. 3 D
impossible to decide by any such reasoning as the above, whether the exterior shell
solid nucleus are now united, or are separated by matter still in a state of fusion
Upon the whole, reasoning suctl as the above can lead us to nothing more defi
than the following conclusions respecting the actual state of the earth, assumin
to have once been in a state of perfect fluidity.
(1.) It may consist of a solid exterior shell and an internal mass in a state of fusio
of which the fluidity is greatest at the centre; and it is possible that the thick
of the shell lmay be small compared with its radius, and the fluidity at the centre
approximate to that which would admit of cooling by convection.
(2.) It may consist of an exterior shell, and a central solid nucleus, with matter in
state of fusion between them. The thickness of the shell, as well as the radius of
solid nucleus, may possibly be small compared with the radius of the earth. T
fluidity of the intervening mass must necessarily be considerably more imperfect t
that which would just admit of cooling by circulation.
(3.) Tile earth may be solid from the surface to the centre.
It appears then that the direct investigation of the manner of the earth's refriger
tion, assuming its original fluidity from heat, still leaves us in a state of perfect un
tainty as to the actual condition of its central parts, not from any imperfection in
mathemnatical part of the investigation, but from the want of the experimental dete
nation of values which it must ever be found extremely difficult, if not impossible
obtain with accuracy. Under these circumstances, we are naturally led to consi
fluidity may consist withl- the observed phenomena of the precessional niotion of th
pole. These phenomena have been shown to be perfectly in accordance with the
ternal solidity of the earth under certain lhypotlheses, which may be deemed perfec
reasonable, respecting the lalw of density; but so far firom any attempts lhaving bee
hitherto made to determine what would be the precessional motion on the supposition
of interior fluidity, I am not aware that the problem has been before suggested. I sh
now proceed to its solution, which forms the principal object of the presenit memoir
On Precession and Nutataion; assuming the Fluidity of the terior of the Earth
In the peset memoir shall investigate the pamount of the luni-solar precessio
and nutation, assuming the earth to consist of a solid spheroidal shell filled with flui
To present the problemn under its most simple form I shall first suppose the solid sh
to be bounded by a determinate inner spheroidal surface, of which the ellipticity
equal to that of the outer surface, the change from the solidity of the shell to t
fluidity of the included mass not being gradual but abrupt. I shall also here sup
pose both the shell and fluild homogeneous and of equal density. From ttlis I p
pose in a future memoir to pass to the case in which the earth is considered
heterogeneous.
?. Statement of the Problem.
1. If S denote the position of the sun, A the centre of the earth, A P its axis of in
stantaneous rotation, the sun's attraction tends to produce an angular velocity of
earth about an axis througlh A, and perpendicular to the plane S A P. The movi
force producing this rotation (supposing the earth a homogeneous spheroid),
3 1J, 4 r 2 n2 - C2) sin 2 A*
- 2 r3 *15 a12 c, (a2 C2) sin 2
whelre
= absolute force of the sun's attraction.
a, = equatorial radius.
c,- polar radius.
Also the moment of inertia of the spheroid about the axis of this rotatory mot
I. The Attraction of the Sun on the Solid Shell.-This will produce effects of pre-
cisely the same kind as if the spheroid were solid, but with different coefficients
(Art. 2.), and therefore, if the motion of I' depended on this cause alone, we should have
dx
dt - A- Bl cos 2 (n t + X)
dy
dt= D sin2 (nt +),
(where n t + x is the longitude of the sun at the time t), these bei
expressions which give the precessional motion of the pole, and
as far as they depend on the sun's action.
II. The Attraction of the JiMoon on the Shell.-This alone woul
dx
dt= A - B. cos 2 (n' t + X')
would be equal at all points similarly situated with respect to the line just ment
through the centres of the sun and earth. If the interior surface of the shell be sp
roidal, but of small eccentricity, very nearly the same effect will be produced
pressure in this case will be exactly equal at points similarly situated with respe
a plane through the sun and the axis of rotation (A P) of the shell, and will c
quently tend to communicate a rotatory motion to the shell about an axis per
dicular to this plane and through the earth's centre; i. e. about the same axis as
about which the attraction mentioned in (I.) tends to communicate a rotatory motio
Also the effects of this pressure must recur with recurring positions of the sun ex
in the same manner as the effects of the sun's attraction just alluded to. Hence
terms depending on this cause will be of the same form as those in (I.), as wil
fact, be proved to be the case when we come to investigate their exact value.
will, therefore, give us
dx
dt = A3- B3cos 2 (nt + t )
dy __
dt D3 sin 2 (n t ).
IV. The Interior Pressure on the Shellfrorn the Attraction of the M
Mass.-This will give us terms similar to those arising fiom the s
this cause alone, therefore, we should have
dx
d-t A4 - B4 cos 2 (n' t + 1')
dy
-y- D4 sin 2 (n t + X).
V. The Interior Pressure on the Shell from the Centrifugal Force of the Fluid Mass.
-When P and P' do not coincide, the interior fluid mass will tend, from the effect of
centrifugal force, to assume a form different from that of the interior surface of the
solid shell. Thus normal pressures will be produced on the interior surface of the
shell; and they will manifestly act symmetrically with respect to a plane through P,
P' and A the centre of the earth, i. e. through the axes of rotation of the solid shell
and of the fluid mass. Consequently the tendency of these pressures will be to com-
municate a motion of rotation to the shell about an axis through A, and perpendi-
cular to this plane; and the consequent motion of P, if this force alone acted on the
shell, would be perpendicular to P' P, the axis of rotation of the shell having, fiom
this cause, an angular velocity in space = - (Art. 1.) c", being the quantity analogous
to vc in the article referred to; or since P and P' are supposed to be on the surface of
a sphere whose radius is unity, -X will be the linear velocity of P perpendicular to P' P.
Now when we come to the calculation of the quantities involved in these investiga-
tions, we shall I
find w
that - = I
y7 sin 2 3, where yr is a constant quantity depending on
~~~
the diurnal angular velocity (w), and on the magnitudes and ellipticity of the f
spheroid and solid shell; and where 3 the angle P A P', or = the line P' P. C
sequently,
t7 == Y sin 2. P' P,
or, since P' P will always be extremely small, the linear velocity of P perpendic
to P P,
= -2 Y1 .P' P;
and resolving this in directions parallel and perpendicular to M' M, we have (; being
the angle which P' P makes with the axis of x)
dy -2 y. P P cos = 2 , (x - x')
6. If we now take the sumnn of the different terms which express the effects of th
d- '
several causes affecting the motion of I, we obtain for the complete values of and
dtt
d = (A1 + A2, + A3 + A4) -- (B1 + B3) cos 2 (n t + x) - (B2 + B4) cos 2 (n' t + ')
-2 7Y (y -y');
d _
^dyt~ =(D + D3) sin 2 (n t + X) + (D2 + D4) sin 2 (n't
+ 2 (- x') ;
or putting
A 4- A2 + A3+ A = A
Bl + B3 = B
B2+ B4 = B'
D1 + D3= D
D2 + D4= D'
dx
d t + 2 y, (y -- y') = A - B cos 2 (n t + ;) - B' cos 2 (n' t + /')
dy -* (A.)
dt -_ 2 y, (x - x')
?. Motion of t
7. When any accel
which the whole sur
viz. that X d x + Y d
III. Centrifugal Force.-In investigating the equations of motion for the solid shel
it has been assumed (Arts. 4. 5.) that the spheroidal axis of the shell will not generally
be coincident with the axis A B' of rotation, which is now proved to be true, since the
disturbing forces of the sun and moon, while they produce a motion of the shell, cause
no motion by their immediate action in the fluid. Let B' A b' (fig. 2.) be the axis of ro-
tation of the interior fluid, and suppose the spheroidal axis first to coincide with it, th
clotted ellipse then representing the section of the interior surface of the solid shell.
'The shell, its form being supposed coincident with that of equilibrium of the fluid
will, in this case, produce no constraint on the fluid motion; but conceive the shel
to be afterwards brought into the position represented by the continuous ellipse, A B
being its spheroidal axis, while B' b' shall still represent the instantaneous axis of ro
tation of the fluid. It is manifest that the planes of rotatory motion of the fluid par-
tides near B' and b' can no longer, as in the former case, be perpendicular to A B', bu
must be constrained to move in planes very nearly parallel to the tangent planes
B' and b'; and it is also sufficiently obvious that whatever effect is produced on
T3'
_ B '._\~' B;___ __ZLT' Fig. 2.
X 13 _ ---_ z
1/i,-: ______ ....
= --- IIN ==
1Is \ 74
r b
This is so small that we shall still consider the angular velocity uniform, which
be proved in the sequel to be correct to the degree of approximation to which it
requisite to carry our investigations.
MDCCCXXXIX. 3 E:
=2 w2 G r cos t sin 0
= ? r sin 2 8, (k = 2 62 g 3);
3 r = Z cos . r
= r2 cos2 0
dR
d-r =k rcos2 ,
d.Or
dr
-= kr ( +- cos 2 ).
Substituting these values in the above equation,
d. f r
k r cos 2 = kr (l + cos 2)+ -) -
d. Olr
d *d r r --krn
' r - - r2 + , () ;
3O' -- - r;
Now suppose forces L- r and - - r equal and in opposite directions to act on each
particle perpendicular to r, together with Z. The motion produced by Z will not be
k
affected by this superposition. But the forces Z (= k x) and - - r are in equilibrium,
and therefore the motion produced by Z must be the same as that which would be
produced by -- r, acting perpendicular to r.
13, Since the motion we are considering is in space of two dimensions, the surface
of the fluid must be defined by some plane curve, if the particular form of which the
result at which we have just arrived is quite independent, being subject only to the
condition that no part of the fluid surface shall be free. Let us suppose the curve to
be a circle, of which the centre is the origin of coordinates. The angular accelerating
force on each particle := -, and is, therefore, the same for every particle. Also the
3 E 2
reaction of the surface would, in this case, have no effect on the angular mo
the fluid. Consequently the angular velocity generated by Z in a unit of time
k
14. If the boundary of the fluid be an ellipse of which the centlre is the origin and
the eccentricity very small, the same result will manifestly be very approximately
true.
This last is the case, in which it was necessary to determine the angular motio
(Art. 11.). It follows that the angular velocity of the internal fluid mass round th
k
axis of y, which would be generated by the force Z in a unit of time = --, or (substi-
tuting the proper value of k) it = W2 s 3, neglecting quantities of the order s2 3.
15. This angular velocity will be coimpounded with that about the axis of z (w).
Now if we again suppose the fluid mass to be spherical, it would manifestly move
precisely as if it were solid, since the angular velocities w and w2 s S about the axes of
z and y respectively are common to all the particles of the mass, and the axis of in-
stantaneous rotation would consequently hlave an angular motion in space perpen-
dicular to the plane B A B' (fig. 2.) and = -= w 3. If the fluid mass be sphe-
roidal, as in our actual case, the ellipticity being small the same result will be very
approximately true.
We may now proceed to the formation of the differential equations for the motion
of the instantaneous axis of rotation of the interior fluid, or of the point P' (fig. 1.).
M-B
Ml =+ B+^-
yD Ml-A=
Blt +'B'
2Dt
?
N=D+ + N-D'-. +
Then
dx
d- + 2 yy L - M cos 2 (n t + X) -M' cos 2 (n' t + X'),
Eachl term in C1 and C corresponding to the several terms in <> (t) and T (t) may be
determined separately.
Let P (t) = L, t (t) -L'; then
d1
dt =2 Ly. cos 2 y + L sin 2y t,
L LI
.'. C= ...sin 2
.. C= -- 2 K (sint cos - t t)
C
C2= 2-2K (cos +
K ..K\ 2 4'
y+ sin rt t).
-2 K -
e4Ksin 2 +
- 2 tJ t C3;
cos 2 yt) +
LL'L
C2 = cos 2 y t + 2 sin 2 y
M+N M -N
-4(y+n) cos2 (y + nt + )-4(y -) cos 2 (-n t-X)
MI + Ml'- N'
4 (y + ') COS 2 ( +
--Kcose siney )
-2 K ( 4y- + ~yj ) + C4,
where c3 and c4 are arbitrary constants.
Substituting these values in equations (C.), we obtain after reduction,
L' nM-yN i n'M' -yN' N,
*x -= -- e- + , (f-_ n~) sin 2 (n t r+ k) + 2 (
n M-y Nn B -- .2 D- y D --y r B
MB y - B B /y B , D)?T
_1__ _ B- Je B *l_ .
' (y2-n2) - n D n
Similarly,
n_MI-YNf
M'-yN' f. 1I Bt _y
B /y
1 B'
B \ )D
D
"2 ( _ - n12) 2 '
ryM-nN yB + ynr D
-?9- n . -Yn-,D D
7n
Similarly,
y M' -_ N _ fi ri (y B'\"1f
y_ -y- yU
.YXD' r-(-n -: Di-, IS- if2'
,B , 'Y 7~ Bt D~
MDCCCXXXIX. 3F
D I- Fy B\D
r-) eCos2X - (i - ) Cos2
l-+
yrsint + ,2 (, ,R
i --B-'_l isi sIin
' D'c s 22 y
t t
+ (v _ B{ D At
L' osa2
+ A t sin 2 Y , (F._
y 'D sin 2 eA9(r
susiin7 s B1D
BDsin 2 x
7D )y B\ D
Y__--cos2 CS 2 (2 t + X .cost )2
_4 scns 2 y t (F)
+-"II)^
2 r6
rrz Cosi2-os
osi2'
( y~-n'" Yi2 lX'
_~--) ( --
1e r .D-
sin .. A
D
.Y _ -(cos2 t c4 2 +, cos 2 n +
+ DI I'V B DI A cosS 2 t ~ _- \~n ~/
y2=-
/ -2yBn
- _\D) cos 2 (n t + x)
Cos 2 B - ACO 2 C
cos 2 t) .
--
We have now to determine the values of A, B, D, B', D', r^ and 72, for which pur-
pose (Art. 6.) we must find the values of the following quantities:
AD1 B B D1
AC B2 D-
A3 B3 D3
A4 134 D4
Y1 72'
3p,iA 4 D4\
r=3
----- ' 5cl(a
- r 15 a-)sin2A,
a (a1 - 2) sin 2 A.,
a1
where q =- the ratio of the outer to the inner radius of the shell.
The moment of inertia of the shell
=: ,- a5 (q5- I),
3F2
and therefore c- 3 p
s= 1o r A
--asin2
w 2r3
6 ~
= Ti- sin 2 A (T = one year)
3 r
= ^sin 2 A
B2 = i( s + ~ (cos2 I sin 2 I,
D.2=--sc
Dt Igcsls2
( s I sin 2 ipo T,
21. Numerical value of A, B3 and D3.-Let the interior surface of the shell be re-
ferred to three rectangular co-oldinates xy z, the spheroidal axis being now that of S ;
and let p denote the normial pressure at the point x y z; j and % the angles which the
normal makes with lines parallel to the axes of x and z respectively. Then if
X = p . S cos ,, Z = p. S cos ,
the moment of the normal pressures with respect to the axis ofy
= : (Zx - Xz)
= 2 (x cos ~ - z cos S) p. S
= Z (- z -PCos ,
P S cos ~.
But 1
cos
COSC? = .. ' ^y r ,I2+ d z\
^ #-
dz
dx
= 1 ( - a) .)xp Scos
2 e I xp. S cos ,
where
a - c
X -= 3 2 x,
Y =--r3 y',
13
and
. p - - 2 "?.p
,,in2A+ ~
xzsin2A +z2cos2A
~~~~~~ S~ - --
Since 2 7 cos 2+ . d p = 0, each of these integrals except the last will vanish be-
tween the limits p- = 0 and p = 2 T. Consequently the moment about the axis of y
=3 r
- s a5sin23
3 A cos2
0 to cos2 d4
d 0 to
= _ 2 }
4-r 5 siN 2 o A54
5 rZ
1 a
q--- . .. (Art. 19.)
Ience we have A3 --~5-
-A 1'
3 - 1'
D DIs~
22. falues of A4, B4, and D4.-In the same manner we obtain
AA
B4-- _ B
D4 -- __
D2
23. Value of yl.-We have seen (Art. 5. V.) that the angullar motion of the solid
shell produced by the centrifugal force of the fluid will be about an axis through the
centre of the spheroid, and perpendicular to the plane passing through the spheroidal
axis of the shell and the axis of rotation of the fluid. Let this axis be now taken for
that of y', and the axis of irotation of.the fluid for that of z', and let x' y' z' be the co-
ordinate of any particle of the fluid; then p denoting the fluid pressure there,
-= o2 sin 2 P a5 fsin3 S2 co s 2 c p d 6 d
4,r
= c 2 ssin 23 a5.
Dividing this by the moment of inertia (= a5 (qS 1) we have (if t be the an-
gular velocity generated in a unit of time by the force we are now considering)
c (i _ o) dsin 2 3,
and
-- = (q 1) sin 2;
and therefore (Art. 5, V.)
Y1 = 2 (qs-1)
7g2r
since if tl = one day, o = 2-r
25. Substituting the values of A, B, B', D, and D' (Art. 6.) we obtain
These give us
- cos I,
B I) cos I
7'Y1 _1 _-- l
YI , 36626- 186t
7 = 36626, =q 8'q
-'= 27'32, =- 70.
Also, taking the ellipticity which the earth would have ha
fluid and homogeneous, we have s = '004 nearly; which al
- = -68.
Hence it appears
7 that- can
18. never exceed
7 or '047, a
may be gre
thickness of the earth's crust,
M DcccxxxIB x D
C - ,2 8 '+ 72 -( )
'-4CCCXX I Y .
MDCCCXXXIX. 3 G
Its value depends on that of -, which may either be much less than unity, equal to
unity, or greater than unity, according to the value of the ratio (q) which the outet
bears to the inner radius of the shell (Art. 25.).
1. Let the shell be thin, or q nearly = 1; then will- be small (Art. 25.). Conse
quently the coefficient becomes
( (l _7) B
Cy en
_ _ 72B
fy 2n
q5-1 B
q5 * OBn
3a
=-- -- sin I cos I.
2. Let the thickness of the shell be such that - = I nearly. Then (Art. 25.
q5- 1 1
2q q5 I neally,
2q - -- 1 nearly,
q = \e7-- ~/ nearly,
and
(3-71)-= 1'3 nearly,
which determines the corresponding value of the ratio of the inner and
the shell. I shall reserve this case for a distinct consideration in the
= - sin I cos I
4v
(q being considerably greater than unity). This value is identical with that found in
the former case.
Y2 B_ q5--1 Bt
-"7 q fn1-- q- ,a t nt
=- 4(+) T Tc 21 Isin 2 i
The coefficient of cos 2 y t consists of two parts, of which the latter (since - is sm
reduces itself to
D- sin 2 ';
and in like manner, when - is small, (or the thickness of the shell small) the first
part becomes
|. , sin2 X,
and the whole coefficient becomes
{% B/ ry A
{ C(os 2 COS 2 X'- }
q A
3q v e , ,.. y f 1 , 1 v
= - sin [ cos I { + 1 * os2 T
v
y
is small,
7y2 D
7 in
3e r
- -- T4Y sin I.
This is also true when q becomes considerable.
The coefficient of cos 2 (n' + '-) becomes
% Dr
7y 2 r
The numerical values of the coefficients of sin 2 y t and cos 2 y t are respecti
the same as those of cos 2 7 t and sin 2 y t in the expression for x.
The constant term becomes
.Sffs.y ,f i , 1 Q 1 1)
+ y sin cos -T + 1-' cos 2r + C
$Q s a ~-
y = -- v sin
1 3 f sin I sinIsin i . t
'--sin 2 2 - 2 sM2 sin 2 i
+^~ qs * 2si"n2A- y^as-+-i) sin2A
sin I cos I cos Q I sin i
cos2 ?, +- (. i ) cos2X '
t4'c- + qcos.2 r -+ C;
+ sin I os ssin I cos si I cos i c 5 1 t
~-.'.. v2 .'. /2( r+ 1) J
when C and C' are small constant terms whose values are
qS-- i tl
q e
x = - 271
G cos 2 y t +- 7H . sin 2 y t,
/ ' 7'
y =71
-72
71
G sin 2 yt- 72 H cos
or i = - G2 H2
7'
cos 2 (y + K),
Referring to equation (E) (Art. 18.) we find that when y - n is very small, w
have (taking what then become the most important periodical terms)
x = - (i - ) - I sin 2 (n t + )
7i /e B D D
- n (D- l) sin 2 X cos 2 y t
Now -- 7__
7y --n tI 72+n
Y+ n 7-n (since
ry-n y n nearly)
also
B
-_ = cos I.
D-vc
~/' D
Hence putting - (1 -cos I) = h, we have
x -- - -- sin 2 (n t + X)
=- n { sin 2 (y t -+ ) - sin 2 (n t + X)
sin 2 (y - n) t
-2 h. ) cos 2 (n t + X) very nearly,
n
0
an expression which assumes the fornm - when y = n accurately. To put it under a
more convenient form, assume t' a particular value of t, such that
2 (y -n) t' = some multiple of x
2 m r ;
let y -- n = s
For the clearer interpretation of th
shell, and therefore y and s, to re
x -2 h 2ss1
sin l(t-
(tt
1Si
t cos
t") 2 (n t'
Since s, is supposed very small compared with n and the product s tI" may be co
sidered as nearly constant for any one year, in which time sin 2 (n t" + L) will pa
through its period, and the solar nutation for that year will depend on
s sin 2 sl t"
Consequently from the time when t = t' or t" = O, this nutation will increase eve
year till 2 s t" = , after which it will again decrease. We should thus have a secula
inequality in the solar nutation, of which the whole period would be -, and of w[ic
In the actual case in whicb y constantly decreases, suppose that at the time t" firo
the time t',
s =; y - n == s - r e',
r denoting the rate of decrease of y, and being taken constant during a small aug-
mentation in the thickness of the shell. Then shall we have
dt
3~BK~~~~ut~* d 1
But r= sQ t (Art. 25.)
a1 w
-as5 -a a1 T
dLy a1ba4 ira da
d /_ aa4 d a
dt-5 (a5 - a5) 2 at
qb re 1 da
= 5 (_ 1)tt * a at
Let a = increase of thickness in time T (one year); then
da
a=t T,
and q5 = ? a
r = 5 (' i 2T)T
aa I
=3 '35-
a T t,
by the substitution of numerical values (Art. 26. 2.). Hence r may be known when
a is determined.
l (-cos)D.
h--_ (1 c- os I) D.
= (qii ]v-2)T
= (q5 1'464- 2) .
h
Hence the greatest value of the secular inequality =
q0 '?4 \"
\q_ I5 1464- -
5 1-464- 2
Qs- 1
This gives 1
Ts = '77 nearly,
or ea = -77 al.
Also when y = n accurately, we ha
qt6
q5l 1'464 - 2 = 0,
which gives
-q =-'768472,
or, if a' be the corresponding value of a,
a' = '768472 . a1.
MDCCCXXXIX. 3 H
Therefore
a- a'-= 6 miles approximately.
The period would be about 125 years.
In order that these numerical results tnay be approximately true, the variation of
a - a' during the period of 125 years must be small compared with a - a'. If we
suppose the thickness of our shell to ilncrease at the same rate as that of the earth's
crust in the process of its solidification, this will probably be true.
Again, if we suppose the inequality to amount to about 1000", we obtain
a - a' = 130 feet nearly;
and the period, supposing a constant, would be about 25,000 years, i. e. in one fourth
of that period the part of solar nutation dependent on the term we are discussing would
pass from zero to about 1000 seconds. If, however, we suppose the solid shell of our
spheroid to increase in thickness at the samre rate as the crust of the earth, the dif-
ference between a and a' would possibly not remain within the value just mentioned
for nearly so long as 6000 years; in which case, supposing the inequality to be zero
when a - a' should equal about 130 feet, it could never afterwards amolunt to nearly
1000 seconds; nor could it have been previously so great, because its previous values
must have corresponded to values of a - a' less than the above value. Our investi-
gation, however, does not tell us whether 120 or 130 feet would be near the value of
a - a' the last time the secular inequality should vanish before a became = a', and
consequently we cannot say with certainty that 1000 seconds would be the extreme
limit to which the inequality would attain. In fact, the exact determination of this
limit woutld require the very accurate determination of a as a function of t, which
cannot be known in the case of the earth's crust without an accurate knowledge of
the conductive power of the matter which constitutes it. From the small value, how.
ever, of a - a' and great length of the period corresponding to the maximum of 1000"
for the inequality we have been considering, it may perhaps be deemed extremnely
improbable that it should ever exceed that value in the case of the earth. The dura-
tion of time for which the effect of the cause we ate discussing on solar nutation
would be sensible to observation would be, that necessary for the thickening of the
earth's crust so to increase that a - a' should pass from + (6 or 8 miles) to - (6 or
8 miles), and might therefore be approximately determined if the quantity denoted
by r in this article were known.
about the axis of y, and that which, combined with the angular mot
causes the angular motion in space of this latter axis. The value
found on the hypothesis of w being constant, or of the rotatory mo
being the samne as if the sections of the inner surface of the shell were
of being ellipses of small eccentricity (Art. 9.); and the pressure on th
of the shell depending on the centrifugal force has been calculated on
pothesis. It will be necessary therefore to examine the errors thus co
33. Let us conceive a closed cylinder entirely filled with fluid, which
fornrmly about the axis of the cylinder with a velocity w, and is not act
external force. If the form of the cylinder be then changed withou
volume, so that each section perpendicular to its axis shall become an
eccentricity instead of being circular, it is manifest from the conditions
that the angular motion of the fluid, though no longer uniform, will
about the same axis, as in the circular cylinder. Consequently if p' b
at any point on the surface of the elliptical cylinder, u' the velocity of t
point, we shall have
p' C - 0'2;
and if p and v be any corresponding values of p' and v' (wlich may be taken
mean values)
P=C-- V2,
and
ar2 - V2
and -1
,lr = 1 + ? :
whence
v'2 - . = v2 s' ( - 2 cos2 0');
and
We have here taken r' and 4' as the polar coordinates of the elliptical section of the
surface, but it is evident that this expression for the fluid pressure will be equally tlr
for any point of the fluid of which r' and 4' are the co-ordinates, p and r being taken
with reference to an ellipse passing through the point (r' d') similar to the elliptic sec-
tion of the cylinder, and similarly situated.
The case which presents itself in our actual problem is analogous to the one ju
considered, so far as regards the elliptical form of the sections made by the planes of
rotatory motion, these planes being parallel to the tangent plane at B' (fig. 2.). Th
common ellipticity of these sections is s (3 (Art. 9.), and, therefore, in finding the effe
on the shell, of the pressure arising from the centrifugal force on the fluid in article
23., we ought to have used for p the pressure as found in the last article, p', or
p + -2 r2 3 cos 2 4'.
This, however, would only introduce into p a term of the order sP , and which, there
fore, may be omitted, as shown in the investigation just referred to (Art. 23.). Al
taking this expression for p' as applicable to any point of the fluid (as explained i
the preceding paragraph), it is easily seen that the force Z (= w2 g P . x) will only
altered, in consequence of the ellipticity of the sections of the shell, by a quantit
small compared with itself, and which may, therefore, be neglected. Our result
then, will be quite accurate to the degree of approximation to which we have pr
ceeded, assuming the parallelism of the planes of rotatory motion to the tangent plan
at B'. I shall now proceed to this point.
34. It has been shown (Art. 15.) that the angular velocity of A B' in space = . s 3
and also (Art. 28.) that the angular velocity of A B is of the same order as that A
i. e. of the order (3. The angular motion of A B' (Arts. 12.... 15.) is due entire
to the obliquity of the planes of rotatory motion of the fluid particles, and the above
value of it is calculated on the hypothesis of these planes of rotation being parallel
the tangent plane at B'. It is easy to see, however, that if this hypothesis be not
curate, the value of y2 (Art. 16.) and therefore of the angular velocity of A
(Art. 15.) will still be quantities of the same order respectively as the calculated values
which have been given, so long as the planes of rotation shall make angles with planes
perpendicular to A B' which, instead of being = to i (Art. 8.), shall be merely of t
same order of magnitude. Without assuming, therefore, the accuracy of the ab
hypothesis, we mnay still assert that the angular velocities of A B and A B' will vneces-
sarily be of the order (3.
The positions of the planes of rotatory motion will be affected by the change of
position of the shell, or of A B, and also by that of A B'. It will be convenient to
consider these cases separately, first, supposing A B' fixed while A B moves, and
then taking A B fixed and A B' in motion.
In the assumption, that A B' shall be at rest, it is meant that it shall here be con-
sidered as unaffected by the angular motion which has been investigated, and which
is due to the obliquity of the planes of rotatory motion. Our first object will be to
examine whether any motion will be communicated to A B' as the direct and imine-
diate consequence of the motion of the shell, and independently of centrifugal force
in the fluid, to which the previously calculated velocity of A B' is entirely due; also
A B' ought strictly to be considered as the line of quiescent Jfuid particles, in which
sense it will not necessarily be a straight line, as we have hitherto considered it in
calculating the effects of centrifugal force on the fluid. It will, therefore, be neces-
sary to examine the degree of its deviation from rectilinearity.
Suppose the shell to be at first in the position represented by the dotted line
(fig. 2.) and then to be brought into that represented by the continuots line, A B,
coinciding at first with A B'. Then while A B moves through the angle B' A B (3),
the normal motion (N N") at any point (N) cannot exceed a quantity of the order s 3,
as is easily shown*. Also it is evident that (considering only the velocity due to the
displacement of the shell) the ratio
vel. of fluid particle at N
vel. of point B of the shell
N N"
must be a quantity of the same order as -Br B i. e. of the order s; and it is easily seen
that for any particle in the interior of the fluid the motion cannot exceed a quantity of
that order. Also the conditions of symmetry will evidently require that the particle
at A should remain at rest.
If the spheroidal axis, instead of moving from A B' to A B, move frorn A B to A 38,
the same conclusion respecting the ratio of the velocity of any fluid particle to the
velocity of B will still be true, as is easily seen.
Let vl be the velocity of B, u that of a particle Q, from the cause we are consideringg,
the distance of Q from A being r, and its distance from the axis A B' = g. Since v
will be of the order s v1 let v = k - vl, where k is a numerical quantity, the value of
which may depend on the position of the particle. Also the velocity of Q from the
motion of rotation round A B' = g. Consequently, if Q be so situated that these
* N v must manifestly vanish with e as well as with nd the expression for it must, therefore, involve
some power of e as a factor.
velocities are impressed upon it along the same line and in opposite directions, th
whole velocity of Q will
= k a v1 - r,
and if this = 0, Q will be a point in the line of quiescent paLticles. This gives us
A _ke uI'
r co a
Now if A B should move from A B' through an angle 3, it is easy to show that
tangent plane at B' must move through an angle of thle order s (; and it is easi
seen likewise that if A B move in any other direction, as from A B to A 3 throu
an angle (', the angle moved through by the tangent plane at B' will be of the or
(3'. In every case, therefore,
ang. vel. of the tang. plane = k' s . ang. vel. of A B
=k S2 (
where k is some finite numerical quantity. Consequently, since the angle describe
by thle mean plane of rotation in a unit of time cannot be greater than this, th
angle, and therefore the instantaneous constr aining force, must be of the order 12 (3
36. Let us now consider the relation between the constraining force and the ang
which the instantaneous mean plane of rotation niakes with the instantaneous tangent
plane at B'. Let i denote, as heretofore, the angle between the tangent plane and
plane perpendicular to A B', ' that between the tangent plane and the mean plan
Now instead of the shell mloving on continuously, let us conceive its motion to ce
at any instant, and consider its action, when thus at rest, on the fluid mass. If
take a fluid particle near to B' and in contact with the surface, B' may be consider
as the centre of its rotatory motion, provided its distance from that point be not less
than a quantity of the order s 3 (since the angular displacement of B' cannot exce
a quantity of the order S2 3). Consequently, if the motion of the plane should cea
at any instant, as above supposed, it is manifest that the plane of motion of this par-
ticle must be immediately constrained to coincide with the tangent plane at B', i.
the constraining force upon it must have been such as to change its plane of motio
through an angle of the order ' in a very short space of time. If we take a particl
in contact with the surface, rather more remote from B', the same conclusion mu
be approximately true, though a somewhat greater time may be necessary to produ
an equal change in the position of its plane of rotation. And similarly if we take
particle in contact with the shell at any point, for instance, between B' and N', th
reaction of the surface must produce a similar effect on its plane of rotation; an
moreover, it is easily seen that if the shell be supposed to retnain thus at rest for
whole revolution, for example, the effect produced in that time must be of the sam
order of magnitude as that for particles near to B'. Precisely the same effects mu
take place about b' and between b' and L', whence it will necessarily follow th
similar effects and of the same order of magnitude must be produced on the plane
of motion of the particles constituting the interior part of the fluid intermediate to
the portions N' n' and L' ' of the surface. Similar effects must also be produced on
the portion of the fluid exterior to that just specified, though these effects may de-
crease in magnitude as the particles are situated nearer to C and c.
Hence then it follows that (taking, for- the greater distinctness, one day for the uni
of time) if B, and therefore the tangent plane at B', were to remain at rest for a unit
of time, the constraining force, estimated as above described, arising from the reaction
of the shell on the fluid, would be such as to cause the mean plane of rotation to mo
through an angle of the same order of magnitude as the instantaneous angle betw
that plane and the tangent plane at B', i. e. of the order /. But it is evident tha
the tangent plane at B', instead of retnaining at rest, as we have liere conceived i
do, have its actual motion from the instantaneous mean plane, the whole effect
one day on the plane of rotation of a particle near to B' orI b' must be greater than
the surface had remained at rest; and the same conclusion must also be true for
ticles more remote fiom B' or b'. Consequently the angle through which the m
plane of rotation moves in one unit of time, must, a fortiori, in the actual motion b
of the same order as /', i. e. the constraining force, estimated by this angle, must b
the order /. But it has been already shown (Art. 35.) that this force must be of
order s2 P. Consequently / must be of the order 2 3 ; or, since =- 2 s (, the an
between the mean plane of rotation and the tangent plane at B' is a small quantit
a higher order than k, which proves the truth of our assumption, in the previous in
vestigations, of the coincidence of these planes to the required degree of approx
tion.
37. We have hitherto considered B to nmove while B' remains at rest; let us n
consider B' to move while B remains at rest. Suppose A B' to move through
angle (' in its motion in space which has been previously investigated, B' then. comin
to 33'. If the shell were spherical, the angle between the tangent planes at B' and
respectively would = 3', and in the spheroid the angle between these planes can d
from P' only by a quantity of the order s P'. Consequently, in order that the m
plane of rotation should be always parallel to the tangent plane at the extremit
the axis of rotation of the fluid, it must move through an angle of the same order
that (f') described by that axis; whereas when A B moves througli an angle (',
corresponding angular motion of the mean plane of rotation is (as we have sho
only of the order s P'. We must examine how this angular motion of the mean p
is produced when A B' is in motion.
While the axis of instantaneous rotation in a rigid body changes its position in
body, the instantaneous planes of rotatory motion necessarily retain their perpendi
larity to it, and therefore the angular motion of those planes is equal to that of
axis. Now we have shown (Art. 15,) that the change in the position of A B' is
duced in a manner exactly similar to that in a solid body, so that the same ca
produces simultaneously the angular motion of A B', and an equal angular motio
the planes of rotation ; whence it is easily seen that the Ymean plane of rotatory mot
when the axis lhas moved to 3', cannot, on this account alone, deviate from paralleli
with the tangent plane at 33' by a quantity greater than of the order s 3'. Co
quently the additional angular velocity of the mean plane of motion necessary to pr
serve it in parallelism with the tangent plane cannot exceed a quantity of the or
g . ang vel. of A B',
=K K 2 .
I. The precession will be the same, whatever be the thickness of the shell, as if the
whole earth were homogeneous and solid.
II. The lunar nutation will be the same as for the homogeneous spheroid to such
a degree of approximation that the difference is inappreciable to observation.
III. The solar nutation will be sensibly the same as for the homogeneous spheroid,
unless the thickness of the shell be very nearly of a certain value, something less than
one-fourth of the earth's radius, in which case this nutation might become much
greater than for the solid spheroid.
IV. In addition to the aove motons of pecession and nutation, the pole of the
earth would have a small circular motion, depending entirely on the internal fluidity.
'l'The radius of the circle thus described would be greatest when the thickness of the
shell should be least; but the inequality thus produced would not for the smallest
thickness of the shell exceed a quantity of the same order as the solar nutation; and
for any but the most inconsiderable thickness of the shell would be entirely inappre-
ciable to observation.
In niy next communication I propose to consider the case in which both the sol
shell and the inclosed fluid mass are of variable density.
W. HOPKINS.
Cambridge,
November 19, 1838.
MDCCCXXXIX. 3 i